/[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.4 by heimbach, Fri Mar 7 05:20:50 2003 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)  #ifdef ALLOW_CTRL_WETV
158        read( funit ) (((nWetwTile(i,j,k), i=1,nsx), j=1,nsy),        read( funit ) (nWetvGlobal(k), k=1,nr)
159       &     k=1,nr)  #endif
160    
161  cgg(    Add OBCS Mask information into the header section for optimization.  cgg(    Add OBCS Mask information into the header section for optimization.
162  #ifdef ALLOW_OBCSN_CONTROL  #ifdef ALLOW_OBCSN_CONTROL
163            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)  
164  #endif  #endif
165  #ifdef ALLOW_OBCSS_CONTROL  #ifdef ALLOW_OBCSS_CONTROL
166            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)  
167  #endif  #endif
168  #ifdef ALLOW_OBCSW_CONTROL  #ifdef ALLOW_OBCSW_CONTROL
169            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)  
170  #endif  #endif
171  #ifdef ALLOW_OBCSE_CONTROL  #ifdef ALLOW_OBCSE_CONTROL
172            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)  
173  #endif  #endif
174  cgg)  cgg)
175        read( funit ) (ncvarindex(i), i=1,maxcvars)        read( funit ) (ncvarindex(i), i=1,maxcvars)
# Line 196  cph      if (lheaderonly) then Line 192  cph      if (lheaderonly) then
192           print *, 'pathei: filensx ', filensx           print *, 'pathei: filensx ', filensx
193           print *, 'pathei: filensy ', filensy           print *, 'pathei: filensy ', filensy
194                    
195           print *, 'pathei: nWetcTile ',           print *, 'pathei: nWetcGlobal ',
196       &        (((nWetcTile(i,j,k), i=1,nsx), j=1,nsy), k=1,nr)       &        (nWetcGlobal(k),  k=1,nr)
197           print *, 'pathei: nWetsTile ',           print *, 'pathei: nWetsGlobal ',
198       &        (((nWetsTile(i,j,k), i=1,nsx), j=1,nsy), k=1,nr)       &        (nWetsGlobal(k),  k=1,nr)
199           print *, 'pathei: nWetwTile ',           print *, 'pathei: nWetwGlobal ',
200       &        (((nWetwTile(i,j,k), i=1,nsx), j=1,nsy), k=1,nr)       &        (nWetwGlobal(k),  k=1,nr)
201             print *, 'pathei: nWetvGlobal ',
202         &        (nWetvGlobal(k),  k=1,nr)
203           print *, 'pathei: ncvarindex ',           print *, 'pathei: ncvarindex ',
204       &        (ncvarindex(i), i=1,maxcvars)       &        (ncvarindex(i), i=1,maxcvars)
205           print *, 'pathei: ncvarrecs ',           print *, 'pathei: ncvarrecs ',
# Line 267  c--   Read the data. Line 265  c--   Read the data.
265                          do k = 1,ncvarnrmax(icvar)                          do k = 1,ncvarnrmax(icvar)
266                             cbuffindex = 0                             cbuffindex = 0
267                             if (ncvargrd(icvar) .eq. 'c') then                             if (ncvargrd(icvar) .eq. 'c') then
268                                cbuffindex = nwetctile(bi,bj,k)                                cbuffindex = nWetcGlobal(k)
269                             else if (ncvargrd(icvar) .eq. 's') then                             else if (ncvargrd(icvar) .eq. 's') then
270                                cbuffindex = nwetstile(bi,bj,k)                                cbuffindex = nWetsGlobal(k)
271                             else if (ncvargrd(icvar) .eq. 'w') then                             else if (ncvargrd(icvar) .eq. 'w') then
272                                cbuffindex = nwetwtile(bi,bj,k)                                cbuffindex = nWetwGlobal(k)
273                               else if (ncvargrd(icvar) .eq. 'v') then
274                                  cbuffindex = nWetvGlobal(k)
275  cgg(   O.B. points have the grid mask "m".  cgg(   O.B. points have the grid mask "m".
276                             else if (ncvargrd(icvar) .eq. 'm') then                             else if (ncvargrd(icvar) .eq. 'm') then
277  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 280  cgg    From "icvrec", calculate what iob
280                               iobcs= icvrec - igg*nobcs                               iobcs= icvrec - igg*nobcs
281  #ifdef ALLOW_OBCSN_CONTROL  #ifdef ALLOW_OBCSN_CONTROL
282                               if (icvar .eq. 11) then                                                   if (icvar .eq. 11) then                    
283                                 cbuffindex = nwetobcsn(bi,bj,k,iobcs)                                 cbuffindex = nWetobcsnGlo(k,iobcs)
284                               endif                               endif
285  #endif  #endif
286  #ifdef ALLOW_OBCSS_CONTROL  #ifdef ALLOW_OBCSS_CONTROL
287                               if (icvar .eq. 12) then                               if (icvar .eq. 12) then
288                                 cbuffindex = nwetobcss(bi,bj,k,iobcs)                                 cbuffindex = nWetobcssGlo(k,iobcs)
289                               endif                               endif
290  #endif  #endif
291  #ifdef ALLOW_OBCSW_CONTROL  #ifdef ALLOW_OBCSW_CONTROL
292                               if (icvar .eq. 13) then                               if (icvar .eq. 13) then
293                                 cbuffindex = nwetobcsw(bi,bj,k,iobcs)                                 cbuffindex = nWetobcswGlo(k,iobcs)
294                               endif                               endif
295  #endif  #endif
296  #ifdef ALLOW_OBCSE_CONTROL  #ifdef ALLOW_OBCSE_CONTROL
297                               if (icvar .eq. 14) then                               if (icvar .eq. 14) then
298                                 cbuffindex = nwetobcse(bi,bj,k,iobcs)                                 cbuffindex = nWetobcseGlo(k,iobcs)
299                               endif                               endif
300  #endif  #endif
301  cgg)  cgg)
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22