/[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.1.2.2 by heimbach, Fri Dec 6 01:47:35 2002 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 = 8        REAL_BYTE = 4
129        isize = REAL_BYTE        isize = REAL_BYTE
130    
131  c-----------------------------------------  c-----------------------------------------
# Line 255  c--------------------------------------- Line 255  c---------------------------------------
255  c--- start if cold start ---  c--- start if cold start ---
256           if (lphprint) then           if (lphprint) then
257              print '(a)', 'pathei-lsopt: cold start'              print '(a)', 'pathei-lsopt: cold start'
             print *  
258           end if           end if
259    
260             print *, 'pathei-lsopt vor simul', nn
261             print *, 'pathei-lsopt xx(1), gg(1) ', xx(1), gg(1)
262    
263           call simul( indic, nn, xx, ff, gg )           call simul( indic, nn, xx, ff, gg )
264  cph(  
265           print *, 'pathei: nach simul:   nn, ff = ', nn, ff           print *, 'pathei: nach simul:   nn, ff = ', nn, ff
266           print *, 'pathei: nach simul:   xx, gg = ', xx(1), gg(1)           print *, 'pathei: nach simul:   xx(1), gg(1) = ', xx(1), gg(1)
 cph)  
267    
268           do i = 1, nn           do i = 1, nn
269              xdiff(i) = 1.              xdiff(i) = 1.
# Line 278  cph) Line 279  cph)
279  cph(  cph(
280        print *, 'pathei: vor lswri '        print *, 'pathei: vor lswri '
281  cph)  cph)
282           call lswri( 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 376  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 449  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( 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 460  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.1.2.2  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22