/[MITgcm]/MITgcm/pkg/ptracers/ptracers_init_varia.F
ViewVC logotype

Diff of /MITgcm/pkg/ptracers/ptracers_init_varia.F

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

revision 1.2 by jmc, Wed Oct 25 20:41:21 2006 UTC revision 1.3 by jmc, Thu Oct 26 00:29:33 2006 UTC
# Line 70  C     Now read initial conditions and al Line 70  C     Now read initial conditions and al
70        IF (nIter0.EQ.PTRACERS_Iter0) THEN        IF (nIter0.EQ.PTRACERS_Iter0) THEN
71         DO iTracer = 1, PTRACERS_numInUse         DO iTracer = 1, PTRACERS_numInUse
72          IF ( PTRACERS_initialFile(iTracer) .NE. ' ' ) THEN          IF ( PTRACERS_initialFile(iTracer) .NE. ' ' ) THEN
          _BEGIN_MASTER( myThid )  
73           CALL READ_FLD_XYZ_RL(PTRACERS_initialFile(iTracer),' ',           CALL READ_FLD_XYZ_RL(PTRACERS_initialFile(iTracer),' ',
74       &        pTracer(1-Olx,1-Oly,1,1,1,iTracer),0,myThid)       &        pTracer(1-Olx,1-Oly,1,1,1,iTracer),0,myThid)
          _END_MASTER(myThid)  
75           _EXCH_XYZ_R8(pTracer(1-Olx,1-Oly,1,1,1,iTracer),myThid)           _EXCH_XYZ_R8(pTracer(1-Olx,1-Oly,1,1,1,iTracer),myThid)
76          ENDIF          ENDIF
77         ENDDO         ENDDO
# Line 84  C     Apply mask Line 82  C     Apply mask
82         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
83          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
84           DO k=1,Nr           DO k=1,Nr
85            DO j=1,sNy            DO j=1-Oly,sNy+OLy
86             DO i=1,sNx             DO i=1-Olx,sNx+Olx
87              IF (maskC(i,j,k,bi,bj).EQ.0.)              IF (maskC(i,j,k,bi,bj).EQ.0.)
88       &           pTracer(i,j,k,bi,bj,iTracer)=0. _d 0       &           pTracer(i,j,k,bi,bj,iTracer)=0. _d 0
89             ENDDO             ENDDO
# Line 95  C     Apply mask Line 93  C     Apply mask
93         ENDDO         ENDDO
94        ENDDO        ENDDO
95    
       _BARRIER  
   
96  C Read from a pickup file if needed  C Read from a pickup file if needed
97        IF (nIter0.GT.PTRACERS_Iter0) THEN        IF (nIter0.GT.PTRACERS_Iter0) THEN
98  C     Suffix for pickup files  C     Suffix for pickup files
# Line 105  C     Suffix for pickup files Line 101  C     Suffix for pickup files
101         ELSE         ELSE
102          WRITE(suff,'(A10)') pickupSuff          WRITE(suff,'(A10)') pickupSuff
103         ENDIF         ENDIF
104         CALL PTRACERS_READ_PICKUP( nIter0,myThid )         CALL PTRACERS_READ_PICKUP( nIter0, myThid )
105        ENDIF        ENDIF
106    
107  #endif /* ALLOW_PTRACERS */  #endif /* ALLOW_PTRACERS */
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22