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

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

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

revision 1.3 by heimbach, Tue Jun 24 16:07:07 2003 UTC revision 1.4 by heimbach, Thu Jul 24 22:00:18 2003 UTC
# Line 59  c     == local variables == Line 59  c     == local variables ==
59        integer bi,bj        integer bi,bj
60        integer ip,jp        integer ip,jp
61        integer i,j,k        integer i,j,k
62        integer ii        integer ii,kk
63        integer il        integer il
64        integer irec,iobcs,nrec_nl        integer irec,iobcs,nrec_nl
65        integer itlo,ithi        integer itlo,ithi
# Line 213  cgg    Get rid of any sensitivity to bar Line 213  cgg    Get rid of any sensitivity to bar
213                   do j = jmin,jmax                   do j = jmin,jmax
214                    if (globmskyz(bi,ip,j,bj,jp,k,iobcs)  .ne. 0. ) then                    if (globmskyz(bi,ip,j,bj,jp,k,iobcs)  .ne. 0. ) then
215                       cbuffindex = cbuffindex + 1                       cbuffindex = cbuffindex + 1
216                         ii=mod((i-1)*nr+k-1,snx)+1
217                         kk=int((i-1)*nr+k-1)/snx+1
218  #ifdef ALLOW_NONDIMENSIONAL_CONTROL_IO  #ifdef ALLOW_NONDIMENSIONAL_CONTROL_IO
219                       if (lxxadxx) then                       if (lxxadxx) then
220                          cbuff(cbuffindex) =                          cbuff(cbuffindex) =
221       &                       globfld3d(i,bi,ip,j,bj,jp,k) *       &                       globfld3d(ii,bi,ip,j,bj,jp,kk) *
222  #ifdef CTRL_PACK_PRECISE  #ifdef CTRL_PACK_PRECISE
223       &                       sqrt(weightfldyz(bi,ip,j,bj,jp,k,iobcs))       &                       sqrt(weightfldyz(bi,ip,j,bj,jp,k,iobcs))
224  #else  #else
# Line 224  cgg    Get rid of any sensitivity to bar Line 226  cgg    Get rid of any sensitivity to bar
226  #endif  #endif
227                       else                       else
228                          cbuff(cbuffindex) =                          cbuff(cbuffindex) =
229       &                       globfld3d(i,bi,ip,j,bj,jp,k) /       &                       globfld3d(ii,bi,ip,j,bj,jp,kk) /
230  #ifdef CTRL_PACK_PRECISE  #ifdef CTRL_PACK_PRECISE
231       &                       sqrt(weightfldyz(bi,ip,j,bj,jp,k,iobcs))       &                       sqrt(weightfldyz(bi,ip,j,bj,jp,k,iobcs))
232  #else  #else
# Line 232  cgg    Get rid of any sensitivity to bar Line 234  cgg    Get rid of any sensitivity to bar
234  #endif  #endif
235                       endif                       endif
236  #else /* ALLOW_NONDIMENSIONAL_CONTROL_IO undef */  #else /* ALLOW_NONDIMENSIONAL_CONTROL_IO undef */
237                       cbuff(cbuffindex) = globfld3d(i,bi,ip,j,bj,jp,k)                       cbuff(cbuffindex) = globfld3d(ii,bi,ip,j,bj,jp,kk)
238  #endif /* ALLOW_NONDIMENSIONAL_CONTROL_IO */  #endif /* ALLOW_NONDIMENSIONAL_CONTROL_IO */
239                    endif                    endif
240                   enddo                   enddo
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22