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

Diff of /MITgcm/optim/optim_writedata.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.3 by heimbach, Fri Dec 6 01:42:25 2002 UTC
# Line 114  cph( Line 114  cph(
114           print *, 'pathei: nsx ', nsx           print *, 'pathei: nsx ', nsx
115           print *, 'pathei: nsy ', nsy           print *, 'pathei: nsy ', nsy
116                    
117           print *, 'pathei: nWetcTile ',           print *, 'pathei: nWetcGlobal ',
118       &        (((nWetcTile(i,j,k), i=1,nsx), j=1,nsy), k=1,nr)       &        (nWetcGlobal(k), k=1,nr)
119           print *, 'pathei: nWetsTile ',           print *, 'pathei: nWetsGlobal ',
120       &        (((nWetsTile(i,j,k), i=1,nsx), j=1,nsy), k=1,nr)       &        (nWetsGlobal(k), k=1,nr)
121           print *, 'pathei: nWetwTile ',           print *, 'pathei: nWetwGlobal ',
122       &        (((nWetwTile(i,j,k), i=1,nsx), j=1,nsy), k=1,nr)       &        (nWetwGlobal(k), k=1,nr)
123             print *, 'pathei: nWetvGlobal ',
124         &        (nWetvGlobal(k), k=1,nr)
125           print *, 'pathei: ncvarindex ',           print *, 'pathei: ncvarindex ',
126       &        (ncvarindex(i), i=1,maxcvars)       &        (ncvarindex(i), i=1,maxcvars)
127           print *, 'pathei: ncvarrecs ',           print *, 'pathei: ncvarrecs ',
# Line 144  c--   Write the header. Line 146  c--   Write the header.
146        write( funit ) bjg        write( funit ) bjg
147        write( funit ) nsx        write( funit ) nsx
148        write( funit ) nsy        write( funit ) nsy
149        write( funit ) (((nWetcTile(i,j,k), i=1,nsx), j=1,nsy),        write( funit ) (nWetcGlobal(k), k=1,nr)
150       &     k=1,nr)        write( funit ) (nWetsGlobal(k), k=1,nr)
151        write( funit ) (((nWetsTile(i,j,k), i=1,nsx), j=1,nsy),        write( funit ) (nWetwGlobal(k), k=1,nr)
152       &     k=1,nr)        write( funit ) (nWetvGlobal(k), k=1,nr)
       write( funit ) (((nWetwTile(i,j,k), i=1,nsx), j=1,nsy),  
      &     k=1,nr)  
153    
154  cgg(    Add OBCS Mask information into the header section for optimization.  cgg(    Add OBCS Mask information into the header section for optimization.
155  #ifdef ALLOW_OBCSN_CONTROL  #ifdef ALLOW_OBCSN_CONTROL
156            write(funit) ((((nWetobcsn(i,j,k,iobcs), k=1,nr),            write(funit) ((nWetobcsnGlo(k,iobcs), k=1,nr),iobcs= 1,nobcs)
      &          iobcs= 1,nobcs), i=1,nsx) , j=1,nsy)  
157  #endif  #endif
158  #ifdef ALLOW_OBCSS_CONTROL  #ifdef ALLOW_OBCSS_CONTROL
159            write(funit) ((((nWetobcss(i,j,k,iobcs), k=1,nr),            write(funit) ((nWetobcssGlo(k,iobcs), k=1,nr),iobcs= 1,nobcs)
      &          iobcs= 1,nobcs), i=1,nsx) , j=1,nsy)  
160  #endif  #endif
161  #ifdef ALLOW_OBCSW_CONTROL  #ifdef ALLOW_OBCSW_CONTROL
162            write(funit) ((((nWetobcsw(i,j,k,iobcs), k=1,nr),            write(funit) ((nWetobcswGlo(k,iobcs), k=1,nr),iobcs= 1,nobcs)
      &          iobcs= 1,nobcs), i=1,nsx) , j=1,nsy)  
163  #endif  #endif
164  #ifdef ALLOW_OBCSE_CONTROL  #ifdef ALLOW_OBCSE_CONTROL
165            write(funit) ((((nWetobcse(i,j,k,iobcs), k=1,nr),            write(funit) ((nWetobcseGlo(k,iobcs), k=1,nr),iobcs= 1,nobcs)
      &          iobcs= 1,nobcs), i=1,nsx) , j=1,nsy)  
166  #endif  #endif
167  cgg)  cgg)
168    
# Line 195  cph) Line 191  cph)
191                       do k = 1,ncvarnrmax(icvar)                       do k = 1,ncvarnrmax(icvar)
192                          cbuffindex = 0                          cbuffindex = 0
193                          if (ncvargrd(icvar) .eq. 'c') then                          if (ncvargrd(icvar) .eq. 'c') then
194                             cbuffindex = nwetctile(bi,bj,k)                             cbuffindex = nWetcGlobal(k)
195                          else if (ncvargrd(icvar) .eq. 's') then                          else if (ncvargrd(icvar) .eq. 's') then
196                             cbuffindex = nwetstile(bi,bj,k)                             cbuffindex = nWetsGlobal(k)
197                          else if (ncvargrd(icvar) .eq. 'w') then                          else if (ncvargrd(icvar) .eq. 'w') then
198                             cbuffindex = nwetwtile(bi,bj,k)                             cbuffindex = nWetwGlobal(k)
199                            else if (ncvargrd(icvar) .eq. 'v') then
200                               cbuffindex = nWetvGlobal(k)
201    
202  cgg(   O.B. points have the grid mask "m".  cgg(   O.B. points have the grid mask "m".
203                          else if (ncvargrd(icvar) .eq. 'm') then                          else if (ncvargrd(icvar) .eq. 'm') then
# Line 209  cgg    From "icvrec", calculate what iob Line 207  cgg    From "icvrec", calculate what iob
207                            iobcs= icvrec - igg*nobcs                            iobcs= icvrec - igg*nobcs
208  #ifdef ALLOW_OBCSN_CONTROL  #ifdef ALLOW_OBCSN_CONTROL
209                            if (icvar .eq. 11) then                                                if (icvar .eq. 11) then                    
210                               cbuffindex = nwetobcsn(bi,bj,k,iobcs)                               cbuffindex = nWetobcsnGlo(k,iobcs)
211                            endif                            endif
212  #endif  #endif
213  #ifdef ALLOW_OBCSS_CONTROL  #ifdef ALLOW_OBCSS_CONTROL
214                            if (icvar .eq. 12) then                            if (icvar .eq. 12) then
215                               cbuffindex = nwetobcss(bi,bj,k,iobcs)                               cbuffindex = nWetobcssGlo(k,iobcs)
216                            endif                            endif
217  #endif  #endif
218  #ifdef ALLOW_OBCSW_CONTROL  #ifdef ALLOW_OBCSW_CONTROL
219                            if (icvar .eq. 13) then                            if (icvar .eq. 13) then
220                               cbuffindex = nwetobcsw(bi,bj,k,iobcs)                               cbuffindex = nWetobcswGlo(k,iobcs)
221                            endif                            endif
222  #endif  #endif
223  #ifdef ALLOW_OBCSE_CONTROL  #ifdef ALLOW_OBCSE_CONTROL
224                            if (icvar .eq. 14) then                            if (icvar .eq. 14) then
225                               cbuffindex = nwetobcse(bi,bj,k,iobcs)                               cbuffindex = nWetobcseGlo(k,iobcs)
226                            endif                            endif
227  #endif  #endif
228                          endif                          endif
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22