/[MITgcm]/MITgcm/optim/optim_readdata.F
ViewVC logotype

Diff of /MITgcm/optim/optim_readdata.F

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

revision 1.4 by heimbach, Fri Mar 7 05:20:50 2003 UTC revision 1.9 by mlosch, Mon Jul 25 08:24:38 2011 UTC
# Line 37  cgg   Include ECCO_CPPOPTIONS because th Line 37  cgg   Include ECCO_CPPOPTIONS because th
37  cgg   have headers with options for OBCS masks.  cgg   have headers with options for OBCS masks.
38  #include "ECCO_CPPOPTIONS.h"  #include "ECCO_CPPOPTIONS.h"
39    
 #include "ecco.h"  
40  #include "ctrl.h"  #include "ctrl.h"
41  #include "optim.h"  #include "optim.h"
42  #include "minimization.h"  #include "minimization.h"
# Line 75  c     == local variables == Line 74  c     == local variables ==
74        integer funit        integer funit
75    
76        integer cbuffindex        integer cbuffindex
77        _RL     cbuff( sNx*nSx*nPx*sNy*nSy*nPy )        real*4 cbuff( sNx*nSx*nPx*sNy*nSy*nPy )
78    
79        character*(128) fname        character*(128) fname
80    
81        integer         filei  c      integer         filei
82        integer         filej  c      integer         filej
83        integer         filek  c      integer         filek
84    c      integer         fileiG
85    c      integer         filejG
86    c      integer         filensx
87    c      integer         filensy
88        integer         filenopt        integer         filenopt
       integer         fileig  
       integer         filejg  
       integer         filensx  
       integer         filensy  
89        _RL             fileff        _RL             fileff
90    
91  cgg(  cgg(
# Line 97  cgg) Line 96  cgg)
96    
97  c     == end of interface ==  c     == end of interface ==
98    
99          print *, 'pathei-lsopt in optim_readdata'
100    
101  c--   The reference i/o unit.  c--   The reference i/o unit.
102        funit = 20        funit = 20
103    
# Line 129  c--   Read the data. Line 130  c--   Read the data.
130    
131  c--   Generate file name and open the file.  c--   Generate file name and open the file.
132        write(fname(1:128),'(4a,i4.4)')        write(fname(1:128),'(4a,i4.4)')
133       &     dfile,'_',expId(1:10),'.opt', nopt       &     dfile,'_',yctrlid(1:10),'.opt', nopt
134        open( funit, file   = fname,        open( funit, file   = fname,
135       &     status = 'old',       &     status = 'old',
136       &     form   = 'unformatted',       &     form   = 'unformatted',
# Line 139  c--   Generate file name and open the fi Line 140  c--   Generate file name and open the fi
140  c--   Read the header.  c--   Read the header.
141        read( funit ) nvartype        read( funit ) nvartype
142        read( funit ) nvarlength        read( funit ) nvarlength
143        read( funit ) expId        read( funit ) yctrlid
144        read( funit ) filenopt        read( funit ) filenopt
145        read( funit ) fileff        read( funit ) fileff
146        read( funit ) fileiG        read( funit ) fileiG
# Line 147  c--   Read the header. Line 148  c--   Read the header.
148        read( funit ) filensx        read( funit ) filensx
149        read( funit ) filensy        read( funit ) filensy
150    
 cph(  
       print *,'ph-opt 1 ', nvartype, nvarlength, filensx, filensy  
 cph)  
   
151        read( funit ) (nWetcGlobal(k), k=1,nr)        read( funit ) (nWetcGlobal(k), k=1,nr)
152        read( funit ) (nWetsGlobal(k), k=1,nr)        read( funit ) (nWetsGlobal(k), k=1,nr)
153        read( funit ) (nWetwGlobal(k), k=1,nr)        read( funit ) (nWetwGlobal(k), k=1,nr)
154  #ifdef ALLOW_CTRL_WETV  #ifdef ALLOW_CTRL_WETV
155        read( funit ) (nWetvGlobal(k), k=1,nr)        read( funit ) (nWetvGlobal(k), k=1,nr)
156  #endif  #endif
157    #ifdef ALLOW_SHIFWFLX_CONTROL
158          read(funit) (nWetiGlobal(k), k=1,nr)
159    c     read(funit) nWetiGlobal(1)
160    #endif
161    
162  cgg(    Add OBCS Mask information into the header section for optimization.  cgg(    Add OBCS Mask information into the header section for optimization.
163  #ifdef ALLOW_OBCSN_CONTROL  #ifdef ALLOW_OBCSN_CONTROL
# Line 184  cph( Line 185  cph(
185  cph      if (lheaderonly) then  cph      if (lheaderonly) then
186           print *, 'pathei: nvartype ', nvartype           print *, 'pathei: nvartype ', nvartype
187           print *, 'pathei: nvarlength ', nvarlength           print *, 'pathei: nvarlength ', nvarlength
188           print *, 'pathei: expId ', expId           print *, 'pathei: yctrlid ', yctrlid
189           print *, 'pathei: filenopt ', filenopt           print *, 'pathei: filenopt ', filenopt
190           print *, 'pathei: fileff ', fileff           print *, 'pathei: fileff ', fileff
191           print *, 'pathei: fileiG ', fileiG           print *, 'pathei: fileiG ', fileiG
# Line 200  cph      if (lheaderonly) then Line 201  cph      if (lheaderonly) then
201       &        (nWetwGlobal(k),  k=1,nr)       &        (nWetwGlobal(k),  k=1,nr)
202           print *, 'pathei: nWetvGlobal ',           print *, 'pathei: nWetvGlobal ',
203       &        (nWetvGlobal(k),  k=1,nr)       &        (nWetvGlobal(k),  k=1,nr)
204    #ifdef ALLOW_SHIFWFLX_CONTROL
205             print *, 'pathei: nWetiGlobal ',
206         &        (nWetiGlobal(k), k=1,nr)
207    #endif
208    #ifdef ALLOW_OBCSN_CONTROL
209             do iobcs=1,nobcs
210              print *, 'pathei: nWetobcsnGlo (iobcs=', iobcs,')',
211         &         (nWetobcsnGlo(k,iobcs), k=1,nr)
212             enddo
213    #endif
214    #ifdef ALLOW_OBCSS_CONTROL
215             do iobcs=1,nobcs
216              print *, 'pathei: nWetobcssGlo (iobcs=', iobcs,')',
217         &         (nWetobcssGlo(k,iobcs), k=1,nr)
218             enddo
219    #endif
220    #ifdef ALLOW_OBCSW_CONTROL
221             do iobcs=1,nobcs
222              print *, 'pathei: nWetobcswGlo (iobcs=', iobcs,')',
223         &         (nWetobcswGlo(k,iobcs), k=1,nr)
224             enddo
225    #endif
226    #ifdef ALLOW_OBCSE_CONTROL
227             do iobcs=1,nobcs
228              print *, 'pathei: nWetobcseGlo (iobcs=', iobcs,')',
229         &         (nWetobcseGlo(k,iobcs), k=1,nr)
230             enddo
231    #endif
232           print *, 'pathei: ncvarindex ',           print *, 'pathei: ncvarindex ',
233       &        (ncvarindex(i), i=1,maxcvars)       &        (ncvarindex(i), i=1,maxcvars)
234           print *, 'pathei: ncvarrecs ',           print *, 'pathei: ncvarrecs ',
# Line 272  c--   Read the data. Line 301  c--   Read the data.
301                                cbuffindex = nWetwGlobal(k)                                cbuffindex = nWetwGlobal(k)
302                             else if (ncvargrd(icvar) .eq. 'v') then                             else if (ncvargrd(icvar) .eq. 'v') then
303                                cbuffindex = nWetvGlobal(k)                                cbuffindex = nWetvGlobal(k)
304    #ifdef ALLOW_SHIFWFLX_CONTROL
305                               else if (ncvargrd(icvar) .eq. 'i') then
306                                  cbuffindex = nWetiGlobal(k)
307    #endif
308  cgg(   O.B. points have the grid mask "m".  cgg(   O.B. points have the grid mask "m".
309                             else if (ncvargrd(icvar) .eq. 'm') then                             else if (ncvargrd(icvar) .eq. 'm') then
310  cgg    From "icvrec", calculate what iobcs must be.  cgg    From "icvrec", calculate what iobcs must be.
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22