/[MITgcm]/MITgcm/lsopt/lsopt_top.F
ViewVC logotype

Diff of /MITgcm/lsopt/lsopt_top.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.4 by heimbach, Thu Sep 9 15:51:26 2004 UTC revision 1.7 by heimbach, Sun Mar 12 08:21:47 2006 UTC
# Line 125  c--------------------------------------- Line 125  c---------------------------------------
125        mm    = nn        mm    = nn
126        mupd  = nupdate        mupd  = nupdate
127    
128        REAL_BYTE = 4        REAL_BYTE = 4
129        isize = REAL_BYTE        isize = REAL_BYTE
130    
131  c-----------------------------------------  c-----------------------------------------
# Line 279  cph) Line 279  cph)
279  cph(  cph(
280        print *, 'pathei: vor lswri '        print *, 'pathei: vor lswri '
281  cph)  cph)
282           call lswri( isize, iiter, nn, xx, gg, lphprint )  cph         call lswri( isize, iiter, nn, xx, gg, lphprint )
283    
284  cph(  cph(
285        print *, 'pathei: vor gnorm0 '        print *, 'pathei: vor gnorm0 '
# Line 377  c--------------------------------------- Line 377  c---------------------------------------
377              print iform, iiter, tact, ifunc, ff, gnorm0              print iform, iiter, tact, ifunc, ff, gnorm0
378       $           , DNRM2( nn, xx, 1 ), 0.       $           , DNRM2( nn, xx, 1 ), 0.
379    
380              write(94,'(i5,2x,1pe11.4,4(2x,1pe8.1))')  Cml            write(94,'(i5,2x,1pe11.4,4(2x,1pe8.1))')
381       &           iiter, ff, gnorm0, tact,  Cml     &           iiter, ff, gnorm0, tact,
382       &           DNRM2( nn, xx, 1 ), 0.  Cml     &           DNRM2( nn, xx, 1 ), 0.
383    
384              if ( itmax .EQ. 0 ) then              if ( itmax .EQ. 0 ) then
385                 ifail = 10                 ifail = 10
# Line 450  c Line 450  c
450        call dostore( nn, xx, .true., 1 )        call dostore( nn, xx, .true., 1 )
451        call dostore( nn, gg, .true., 2 )        call dostore( nn, gg, .true., 2 )
452  cph(  cph(
453        call lswri( isize, iiter, nn, xx, gg, lphprint )  cph      call lswri( isize, iiter, nn, xx, gg, lphprint )
454  cph)  cph)
455    
456        gnorm = DNRM2( nn, gg, 1 )        gnorm = DNRM2( nn, gg, 1 )
# Line 461  c--------------------------------------- Line 461  c---------------------------------------
461        print iform, iiter, tact, ifunc, ff, gnorm        print iform, iiter, tact, ifunc, ff, gnorm
462       $     , DNRM2( nn, xx, 1 ), tact*DNRM2( nn, dd, 1 )       $     , DNRM2( nn, xx, 1 ), tact*DNRM2( nn, dd, 1 )
463    
464        write(94,'(i5,2x,1pe11.4,4(2x,1pe8.1))')  Cml      write(94,'(i5,2x,1pe11.4,4(2x,1pe8.1))')
465       &     iiter, ff, gnorm, tact,  Cml     &     iiter, ff, gnorm, tact,
466       &     DNRM2( nn, xx, 1 ), tact*DNRM2( nn, dd, 1 )  Cml     &     DNRM2( nn, xx, 1 ), tact*DNRM2( nn, dd, 1 )
467    
468  c-----------------------------------------  c-----------------------------------------
469  c check output mode of ifail  c check output mode of ifail
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22