/[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.3 by heimbach, Fri Dec 6 01:42:25 2002 UTC revision 1.4 by heimbach, Thu Sep 9 15:51:26 2004 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 )           call lswri( isize, iiter, nn, xx, gg, lphprint )
283    
284  cph(  cph(
285        print *, 'pathei: vor gnorm0 '        print *, 'pathei: vor gnorm0 '
# 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 )        call lswri( isize, iiter, nn, xx, gg, lphprint )
454  cph)  cph)
455    
456        gnorm = DNRM2( nn, gg, 1 )        gnorm = DNRM2( nn, gg, 1 )
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22