/[MITgcm]/MITgcm_contrib/mlosch/optim_m1qn3/optim_readparms.F
ViewVC logotype

Diff of /MITgcm_contrib/mlosch/optim_m1qn3/optim_readparms.F

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

revision 1.1 by mlosch, Thu Apr 26 11:10:06 2012 UTC revision 1.3 by mlosch, Fri Apr 27 09:09:39 2012 UTC
# Line 95  c--   Read control parameters from file. Line 95  c--   Read control parameters from file.
95        open(unit = modeldataunit,file = 'data.ctrl',        open(unit = modeldataunit,file = 'data.ctrl',
96       &     status = 'old', iostat = errio)       &     status = 'old', iostat = errio)
97        if ( errio .lt. 0 ) then        if ( errio .lt. 0 ) then
98         stop ' stopped in optim_numbmod while opening data.ctrl'         stop ' stopped in optim_readparms while opening data.ctrl'
99        endif        endif
100    
101        do while ( .true. )        do while ( .true. )
# Line 118  c--   Read optimization parameters from Line 118  c--   Read optimization parameters from
118        open(unit = modeldataunit,file = 'data.optim',        open(unit = modeldataunit,file = 'data.optim',
119       &     status = 'old', iostat = errio)       &     status = 'old', iostat = errio)
120        if ( errio .lt. 0 ) then        if ( errio .lt. 0 ) then
121         stop ' stopped in optim_numbmod while opening data.optim'         stop ' stopped in optim_readparms while opening data.optim'
122        endif        endif
123    
124        do while ( .true. )        do while ( .true. )
# Line 158  c     otherwise we (ab-)use the file ctr Line 158  c     otherwise we (ab-)use the file ctr
158          fmin = ff          fmin = ff
159         endif         endif
160        endif        endif
161        if ( fmin.le.0.0 .and. optimcycle.eq.0 ) then        if ( optimcycle.eq.0 ) then
162          if ( fmin.le.0.0 ) then
163         print*, ' OPTIM_READPARMS: fmin = ', fmin, ' should be > 0'         print*, ' OPTIM_READPARMS: fmin = ', fmin, ' should be > 0'
164         stop 'S/R OPTIM_READPARMS: ABNORMAL END'         stop 'S/R OPTIM_READPARMS: ABNORMAL END'
165        endif        endif
166          if ( fmin.lt.ff ) then
167           print*, ' OPTIM_READPARMS: fmin = ', fmin,
168         &      ' should be < ff = ', ff
169           stop 'S/R OPTIM_READPARMS: ABNORMAL END'
170          endif
171    
172  c--   Do some final printout.  c--   Do some final printout.
173        print*        print*
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22