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

Diff of /MITgcm/pkg/ctrl/ctrl_set_unpack_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 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,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 iob Line 229  cgg      And now back-calculate what iob
229                     do bi = itlo,ithi                     do bi = itlo,ithi
230                      if ( globmskyz(bi,ip,j,bj,jp,k,iobcs) .ne. 0. ) then                      if ( globmskyz(bi,ip,j,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)                         ii=mod((i-1)*nr+k-1,snx)+1
233                           kk=int((i-1)*nr+k-1)/snx+1
234                           globfld3d(ii,bi,ip,j,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(ii,bi,ip,j,bj,jp,kk) =
238       &                      globfld3d(i,bi,ip,j,bj,jp,k)/       &                      globfld3d(ii,bi,ip,j,bj,jp,kk)/
239  # ifdef CTRL_UNPACK_PRECISE  # ifdef CTRL_UNPACK_PRECISE
240       &                      sqrt(weightfldyz(bi,ip,j,bj,jp,k,iobcs))       &                      sqrt(weightfldyz(bi,ip,j,bj,jp,k,iobcs))
241  # else  # else
# Line 240  cgg      And now back-calculate what iob Line 243  cgg      And now back-calculate what iob
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(ii,bi,ip,j,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