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

Diff of /MITgcm/pkg/ctrl/ctrl_set_unpack_xz.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 58  c     == local variables == Line 58  c     == local variables ==
58        integer bi,bj        integer bi,bj
59        integer ip,jp        integer ip,jp
60        integer i,j,k        integer i,j,k
61        integer ii        integer ii,jj,kk
62        integer il        integer il
63        integer irec,iobcs,nrec_nl        integer irec,iobcs,nrec_nl
64        integer itlo,ithi        integer itlo,ithi
# Line 229  cgg       And now back-calculate what io Line 229  cgg       And now back-calculate what io
229                     do i = imin,imax                     do i = imin,imax
230                      if ( globmskxz(i,bi,ip,bj,jp,k,iobcs) .ne. 0. ) then                      if ( globmskxz(i,bi,ip,bj,jp,k,iobcs) .ne. 0. ) then
231                         cbuffindex = cbuffindex + 1                         cbuffindex = cbuffindex + 1
232                         globfld3d(i,bi,ip,j,bj,jp,k) = cbuff(cbuffindex)                         jj=mod((j-1)*nr+k-1,sny)+1
233                           kk=int((j-1)*nr+k-1)/sny+1
234                           globfld3d(i,bi,ip,jj,bj,jp,kk) =
235         &                      cbuff(cbuffindex)
236  #ifdef ALLOW_NONDIMENSIONAL_CONTROL_IO  #ifdef ALLOW_NONDIMENSIONAL_CONTROL_IO
237                         globfld3d(i,bi,ip,j,bj,jp,k) =                         globfld3d(i,bi,ip,jj,bj,jp,kk) =
238       &                      globfld3d(i,bi,ip,j,bj,jp,k)/       &                      globfld3d(i,bi,ip,jj,bj,jp,kk)/
239  # ifdef CTRL_UNPACK_PRECISE  # ifdef CTRL_UNPACK_PRECISE
240       &                      sqrt(weightfldxz(i,bi,ip,bj,jp,k,iobcs))       &                      sqrt(weightfldxz(i,bi,ip,bj,jp,k,iobcs))
241  # else  # else
# Line 240  cgg       And now back-calculate what io Line 243  cgg       And now back-calculate what io
243  # endif  # endif
244  #endif /* ALLOW_NONDIMENSIONAL_CONTROL_IO */  #endif /* ALLOW_NONDIMENSIONAL_CONTROL_IO */
245                      else                      else
246                         globfld3d(i,bi,ip,j,bj,jp,k) = 0. _d 0                         globfld3d(i,bi,ip,jj,bj,jp,kk) = 0. _d 0
247                      endif                      endif
248                     enddo                     enddo
249                    enddo                    enddo
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22