/[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.6 by mlosch, Fri Apr 27 13:58:14 2012 UTC
# Line 41  c     == local variables == Line 41  c     == local variables ==
41        _RL            ff        _RL            ff
42        _RL            dfminFrac        _RL            dfminFrac
43    
44  #if defined (DYNAMIC)  #ifdef DYNAMIC
45        _RL            vv(nn)        _RL            vv(nn)
 #elif defined (USE_POINTER) || (MAX_INDEPEND == 0)  
       _RL            vv  
       pointer (pvv,vv(1))  
46  #else  #else
47        integer nmax        integer nmax
48        parameter( nmax = MAX_INDEPEND )        parameter( nmax = MAX_INDEPEND )
# Line 95  c--   Read control parameters from file. Line 92  c--   Read control parameters from file.
92        open(unit = modeldataunit,file = 'data.ctrl',        open(unit = modeldataunit,file = 'data.ctrl',
93       &     status = 'old', iostat = errio)       &     status = 'old', iostat = errio)
94        if ( errio .lt. 0 ) then        if ( errio .lt. 0 ) then
95         stop ' stopped in optim_numbmod while opening data.ctrl'         stop ' stopped in optim_readparms while opening data.ctrl'
96        endif        endif
97    
98        do while ( .true. )        do while ( .true. )
# Line 118  c--   Read optimization parameters from Line 115  c--   Read optimization parameters from
115        open(unit = modeldataunit,file = 'data.optim',        open(unit = modeldataunit,file = 'data.optim',
116       &     status = 'old', iostat = errio)       &     status = 'old', iostat = errio)
117        if ( errio .lt. 0 ) then        if ( errio .lt. 0 ) then
118         stop ' stopped in optim_numbmod while opening data.optim'         stop ' stopped in optim_readparms while opening data.optim'
119        endif        endif
120    
121        do while ( .true. )        do while ( .true. )
# Line 158  c     otherwise we (ab-)use the file ctr Line 155  c     otherwise we (ab-)use the file ctr
155          fmin = ff          fmin = ff
156         endif         endif
157        endif        endif
158        if ( fmin.le.0.0 .and. optimcycle.eq.0 ) then        if ( optimcycle.eq.0 ) then
159         print*, ' OPTIM_READPARMS: fmin = ', fmin, ' should be > 0'         if ( fmin.le.0.0 ) then
160         stop 'S/R OPTIM_READPARMS: ABNORMAL END'          print '(A,E12.6,A)', ' OPTIM_READPARMS: fmin = ', fmin,
161         &       ' should be > 0'
162            stop 'S/R OPTIM_READPARMS: ABNORMAL END'
163           endif
164           if ( fmin.gt.ff ) then
165            print '(A,E12.6,A,E12.6)', ' OPTIM_READPARMS: fmin = ', fmin,
166         &       ' should be < ff = ', ff
167            stop 'S/R OPTIM_READPARMS: ABNORMAL END'
168           endif
169        endif        endif
170    
171  c--   Do some final printout.  c--   Do some final printout.
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22