/[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.1.2.2 by heimbach, Thu Apr 4 10:20:16 2002 UTC revision 1.1.2.3 by heimbach, Fri Dec 6 01:47:35 2002 UTC
# Line 151  cph( Line 151  cph(
151        print *,'ph-opt 1 ', nvartype, nvarlength, filensx, filensy        print *,'ph-opt 1 ', nvartype, nvarlength, filensx, filensy
152  cph)  cph)
153    
154        read( funit ) (((nWetcTile(i,j,k), i=1,nsx), j=1,nsy),        read( funit ) (nWetcGlobal(k), k=1,nr)
155       &     k=1,nr)        read( funit ) (nWetsGlobal(k), k=1,nr)
156        read( funit ) (((nWetsTile(i,j,k), i=1,nsx), j=1,nsy),        read( funit ) (nWetwGlobal(k), k=1,nr)
157       &     k=1,nr)        read( funit ) (nWetvGlobal(k), 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.1.2.2  
changed lines
  Added in v.1.1.2.3
  ViewVC Help
Powered by ViewVC 1.1.22