/[MITgcm]/MITgcm/pkg/mdsio/mdsio_readfield.F
ViewVC logotype

Diff of /MITgcm/pkg/mdsio/mdsio_readfield.F

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

revision 1.20 by dimitri, Wed Dec 22 07:27:07 2004 UTC revision 1.21 by jmc, Fri Feb 11 03:04:43 2005 UTC
# Line 82  C Local variables Line 82  C Local variables
82        Real*8 global   (Nx,Ny)        Real*8 global   (Nx,Ny)
83        _RL    local(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL    local(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
84  #ifdef ALLOW_EXCH2  #ifdef ALLOW_EXCH2
85        integer domainHeight,domainLength,tby,tgx,tny,tnx,tn        integer domainHeight,domainLength,tgy,tgx,tny,tnx,tn
86  #endif /* ALLOW_EXCH2 */  #endif /* ALLOW_EXCH2 */
87        COMMON /GlobalLo/ mpi_myXGlobalLo, mpi_myYGlobalLo        COMMON /GlobalLo/ mpi_myXGlobalLo, mpi_myYGlobalLo
88        INTEGER mpi_myXGlobalLo(nPx*nPy)        INTEGER mpi_myXGlobalLo(nPx*nPy)
# Line 180  cafe             call print_error( msgbu Line 180  cafe             call print_error( msgbu
180  cafe             stop 'ABNORMAL END: S/R MDSREADFIELD'  cafe             stop 'ABNORMAL END: S/R MDSREADFIELD'
181  cafe        endif  cafe        endif
182          tn = W2_myTileList(bi)          tn = W2_myTileList(bi)
183          tby = exch2_tbasey(tn)          tgy = exch2_tyglobalo(tn)
184          tgx = exch2_txglobalo(tn)          tgx = exch2_txglobalo(tn)
185          tny = exch2_tny(tn)          tny = exch2_tny(tn)
186          tnx = exch2_tnx(tn)          tnx = exch2_tnx(tn)
187          if (fileIsOpen) then          if (fileIsOpen) then
188           do k=1,nNz           do k=1,nNz
189            do j=1,tNy            do j=1,tNy
190  cafe          write(*,fmt='(1X,A,I3,A,I3,A,I3,A,I3,A,I3,A,I3)') 'tby=', tby,  cafe          write(*,fmt='(1X,A,I3,A,I3,A,I3,A,I3,A,I3,A,I3)') 'tgy=', tgy,
191  cafe     &        ', tgx=', tgx,  cafe     &        ', tgx=', tgx,
192  cafe     &            ', tnx=',tnx, ', tny=', tny, ', j=',j,', tn=',tn  cafe     &            ', tnx=',tnx, ', tny=', tny, ', j=',j,', tn=',tn
193                            
194              irec = domainLength*tby + (tgx-1)/tnx + 1 +              irec = domainLength*(tgy-1) + (tgx-1)/tnx + 1 +
195       &             domainLength*(j-1) +       &             domainLength*(j-1) +
196       &             domainLength*domainHeight*tny*(k-1) +       &             domainLength*domainHeight*tny*(k-1) +
197       &             domainLength*domainHeight*tny*nNz*(irecord-1)       &             domainLength*domainHeight*tny*nNz*(irecord-1)
Legend:
Removed from v.1.20  
changed lines
  Added in v.1.21
  ViewVC Help
Powered by ViewVC 1.1.22