/[MITgcm]/MITgcm/pkg/ctrl/ctrl_set_pack_xz.F
ViewVC logotype

Diff of /MITgcm/pkg/ctrl/ctrl_set_pack_xz.F

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

revision 1.6 by heimbach, Thu Oct 30 19:09:05 2003 UTC revision 1.7 by heimbach, Thu Nov 6 22:05:08 2003 UTC
# Line 1  Line 1 
 C  
 C $Header$  
 C $Name$  
1    
2  #include "CTRL_CPPOPTIONS.h"  #include "CTRL_CPPOPTIONS.h"
3    
   
4        subroutine ctrl_set_pack_xz(        subroutine ctrl_set_pack_xz(
5       &     cunit, ivartype, fname, masktype,weighttype,       &     cunit, ivartype, fname, masktype,weighttype,
6       &     weightfld, lxxadxx, mythid)       &     weightfld, lxxadxx, mythid)
# Line 36  c     == global variables == Line 32  c     == global variables ==
32  #include "GRID.h"  #include "GRID.h"
33    
34  #include "ctrl.h"  #include "ctrl.h"
 #include "cost.h"  
35    
36  #ifdef ALLOW_ECCO_OPTIMIZATION  #ifdef ALLOW_ECCO_OPTIMIZATION
37  #include "optim.h"  #include "optim.h"
# Line 77  cgg( Line 72  cgg(
72        character*(80) weightname        character*(80) weightname
73  cgg)  cgg)
74    
75        _RL     cbuff      ( snx*nsx*npx*nsy*npy )        real*4     cbuff      ( snx*nsx*npx*nsy*npy )
76        _RL     globfldxz  ( snx,nsx,npx,nsy,npy,nr )        _RL     globfldxz  ( snx,nsx,npx,nsy,npy,nr )
77        _RL     globfld3d  ( snx,nsx,npx,sny,nsy,npy,nr )        _RL     globfld3d  ( snx,nsx,npx,sny,nsy,npy,nr )
78        _RL     globmskxz  ( snx,nsx,npx,nsy,npy,nr,nobcs )        _RL     globmskxz  ( snx,nsx,npx,nsy,npy,nr,nobcs )
# Line 211  cgg    Get rid of any sensitivity to bar Line 206  cgg    Get rid of any sensitivity to bar
206                 do ip = 1,nPx                 do ip = 1,nPx
207                  do bi = itlo,ithi                  do bi = itlo,ithi
208                   do i = imin,imax                   do i = imin,imax
209                      jj=mod((j-1)*nr+k-1,sny)+1
210                      kk=int((j-1)*nr+K-1)/sny+1
211                    if (globmskxz(i,bi,ip,bj,jp,k,iobcs)  .ne. 0. ) then                    if (globmskxz(i,bi,ip,bj,jp,k,iobcs)  .ne. 0. ) then
212                       cbuffindex = cbuffindex + 1                       cbuffindex = cbuffindex + 1
                      jj=mod((j-1)*nr+k-1,sny)+1  
                      kk=int((j-1)*nr+K-1)/sny+1  
213  #ifdef ALLOW_NONDIMENSIONAL_CONTROL_IO  #ifdef ALLOW_NONDIMENSIONAL_CONTROL_IO
214                       if (lxxadxx) then                       if (lxxadxx) then
215                          cbuff(cbuffindex) =                          cbuff(cbuffindex) =
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22