/[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.2 by heimbach, Fri Nov 15 04:03:25 2002 UTC revision 1.6 by heimbach, Thu Sep 9 15:58:42 2004 UTC
# Line 75  c     == local variables == Line 75  c     == local variables ==
75        integer funit        integer funit
76    
77        integer cbuffindex        integer cbuffindex
78        _RL     cbuff( sNx*nSx*nPx*sNy*nSy*nPy )        real*4 cbuff( sNx*nSx*nPx*sNy*nSy*nPy )
79    
80        character*(128) fname        character*(128) fname
81    
82        integer         filei  c      integer         filei
83        integer         filej  c      integer         filej
84        integer         filek  c      integer         filek
85    c      integer         fileig
86    c      integer         filejg
87    c      integer         filensx
88    c      integer         filensy
89        integer         filenopt        integer         filenopt
       integer         fileig  
       integer         filejg  
       integer         filensx  
       integer         filensy  
90        _RL             fileff        _RL             fileff
91    
92  cgg(  cgg(
# Line 97  cgg) Line 97  cgg)
97    
98  c     == end of interface ==  c     == end of interface ==
99    
100          print *, 'pathei-lsopt in optim_readdata'
101    
102  c--   The reference i/o unit.  c--   The reference i/o unit.
103        funit = 20        funit = 20
104    
# Line 147  c--   Read the header. Line 149  c--   Read the header.
149        read( funit ) filensx        read( funit ) filensx
150        read( funit ) filensy        read( funit ) filensy
151    
152  cph(        read( funit ) (nWetcGlobal(k), k=1,nr)
153        print *,'ph-opt 1 ', nvartype, nvarlength, filensx, filensy        read( funit ) (nWetsGlobal(k), k=1,nr)
154  cph)        read( funit ) (nWetwGlobal(k), k=1,nr)
155    #ifdef ALLOW_CTRL_WETV
156        read( funit ) (((nWetcTile(i,j,k), i=1,nsx), j=1,nsy),        read( funit ) (nWetvGlobal(k), k=1,nr)
157       &     k=1,nr)  #endif
       read( funit ) (((nWetsTile(i,j,k), i=1,nsx), j=1,nsy),  
      &     k=1,nr)  
       read( funit ) (((nWetwTile(i,j,k), i=1,nsx), j=1,nsy),  
      &     k=1,nr)  
158    
159  cgg(    Add OBCS Mask information into the header section for optimization.  cgg(    Add OBCS Mask information into the header section for optimization.
160  #ifdef ALLOW_OBCSN_CONTROL  #ifdef ALLOW_OBCSN_CONTROL
161            read(funit) ((((nWetobcsn(i,j,k,iobcs), k=1,nr),        read( funit ) ((nWetobcsnGlo(k,iobcs), k=1,nr),iobcs= 1,nobcs)
      &          iobcs= 1,nobcs), i=1,nsx) , j=1,nsy)  
162  #endif  #endif
163  #ifdef ALLOW_OBCSS_CONTROL  #ifdef ALLOW_OBCSS_CONTROL
164            read(funit) ((((nWetobcss(i,j,k,iobcs), k=1,nr),        read( funit ) ((nWetobcssGlo(k,iobcs), k=1,nr),iobcs= 1,nobcs)
      &          iobcs= 1,nobcs), i=1,nsx) , j=1,nsy)  
165  #endif  #endif
166  #ifdef ALLOW_OBCSW_CONTROL  #ifdef ALLOW_OBCSW_CONTROL
167            read(funit) ((((nWetobcsw(i,j,k,iobcs), k=1,nr),        read( funit ) ((nWetobcswGlo(k,iobcs), k=1,nr),iobcs= 1,nobcs)
      &          iobcs= 1,nobcs), i=1,nsx) , j=1,nsy)  
168  #endif  #endif
169  #ifdef ALLOW_OBCSE_CONTROL  #ifdef ALLOW_OBCSE_CONTROL
170            read(funit) ((((nWetobcse(i,j,k,iobcs), k=1,nr),        read( funit ) ((nWetobcseGlo(k,iobcs), k=1,nr),iobcs= 1,nobcs)
      &          iobcs= 1,nobcs), i=1,nsx) , j=1,nsy)  
171  #endif  #endif
172  cgg)  cgg)
173        read( funit ) (ncvarindex(i), i=1,maxcvars)        read( funit ) (ncvarindex(i), i=1,maxcvars)
# Line 196  cph      if (lheaderonly) then Line 190  cph      if (lheaderonly) then
190           print *, 'pathei: filensx ', filensx           print *, 'pathei: filensx ', filensx
191           print *, 'pathei: filensy ', filensy           print *, 'pathei: filensy ', filensy
192                    
193           print *, 'pathei: nWetcTile ',           print *, 'pathei: nWetcGlobal ',
194       &        (((nWetcTile(i,j,k), i=1,nsx), j=1,nsy), k=1,nr)       &        (nWetcGlobal(k),  k=1,nr)
195           print *, 'pathei: nWetsTile ',           print *, 'pathei: nWetsGlobal ',
196       &        (((nWetsTile(i,j,k), i=1,nsx), j=1,nsy), k=1,nr)       &        (nWetsGlobal(k),  k=1,nr)
197           print *, 'pathei: nWetwTile ',           print *, 'pathei: nWetwGlobal ',
198       &        (((nWetwTile(i,j,k), i=1,nsx), j=1,nsy), k=1,nr)       &        (nWetwGlobal(k),  k=1,nr)
199             print *, 'pathei: nWetvGlobal ',
200         &        (nWetvGlobal(k),  k=1,nr)
201           print *, 'pathei: ncvarindex ',           print *, 'pathei: ncvarindex ',
202       &        (ncvarindex(i), i=1,maxcvars)       &        (ncvarindex(i), i=1,maxcvars)
203           print *, 'pathei: ncvarrecs ',           print *, 'pathei: ncvarrecs ',
# Line 267  c--   Read the data. Line 263  c--   Read the data.
263                          do k = 1,ncvarnrmax(icvar)                          do k = 1,ncvarnrmax(icvar)
264                             cbuffindex = 0                             cbuffindex = 0
265                             if (ncvargrd(icvar) .eq. 'c') then                             if (ncvargrd(icvar) .eq. 'c') then
266                                cbuffindex = nwetctile(bi,bj,k)                                cbuffindex = nWetcGlobal(k)
267                             else if (ncvargrd(icvar) .eq. 's') then                             else if (ncvargrd(icvar) .eq. 's') then
268                                cbuffindex = nwetstile(bi,bj,k)                                cbuffindex = nWetsGlobal(k)
269                             else if (ncvargrd(icvar) .eq. 'w') then                             else if (ncvargrd(icvar) .eq. 'w') then
270                                cbuffindex = nwetwtile(bi,bj,k)                                cbuffindex = nWetwGlobal(k)
271                               else if (ncvargrd(icvar) .eq. 'v') then
272                                  cbuffindex = nWetvGlobal(k)
273  cgg(   O.B. points have the grid mask "m".  cgg(   O.B. points have the grid mask "m".
274                             else if (ncvargrd(icvar) .eq. 'm') then                             else if (ncvargrd(icvar) .eq. 'm') then
275  cgg    From "icvrec", calculate what iobcs must be.  cgg    From "icvrec", calculate what iobcs must be.
# Line 280  cgg    From "icvrec", calculate what iob Line 278  cgg    From "icvrec", calculate what iob
278                               iobcs= icvrec - igg*nobcs                               iobcs= icvrec - igg*nobcs
279  #ifdef ALLOW_OBCSN_CONTROL  #ifdef ALLOW_OBCSN_CONTROL
280                               if (icvar .eq. 11) then                                                   if (icvar .eq. 11) then                    
281                                 cbuffindex = nwetobcsn(bi,bj,k,iobcs)                                 cbuffindex = nWetobcsnGlo(k,iobcs)
282                               endif                               endif
283  #endif  #endif
284  #ifdef ALLOW_OBCSS_CONTROL  #ifdef ALLOW_OBCSS_CONTROL
285                               if (icvar .eq. 12) then                               if (icvar .eq. 12) then
286                                 cbuffindex = nwetobcss(bi,bj,k,iobcs)                                 cbuffindex = nWetobcssGlo(k,iobcs)
287                               endif                               endif
288  #endif  #endif
289  #ifdef ALLOW_OBCSW_CONTROL  #ifdef ALLOW_OBCSW_CONTROL
290                               if (icvar .eq. 13) then                               if (icvar .eq. 13) then
291                                 cbuffindex = nwetobcsw(bi,bj,k,iobcs)                                 cbuffindex = nWetobcswGlo(k,iobcs)
292                               endif                               endif
293  #endif  #endif
294  #ifdef ALLOW_OBCSE_CONTROL  #ifdef ALLOW_OBCSE_CONTROL
295                               if (icvar .eq. 14) then                               if (icvar .eq. 14) then
296                                 cbuffindex = nwetobcse(bi,bj,k,iobcs)                                 cbuffindex = nWetobcseGlo(k,iobcs)
297                               endif                               endif
298  #endif  #endif
299  cgg)  cgg)
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22