/[MITgcm]/MITgcm_contrib/dgoldberg/streamice/streamice_advect_2dtracer.F
ViewVC logotype

Diff of /MITgcm_contrib/dgoldberg/streamice/streamice_advect_2dtracer.F

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

revision 1.1 by dgoldberg, Sat Jun 8 22:15:33 2013 UTC revision 1.2 by dgoldberg, Wed Jun 12 20:48:08 2013 UTC
# Line 120  CADJ STORE trac2d  = comlev1, key=ikey_d Line 120  CADJ STORE trac2d  = comlev1, key=ikey_d
120        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
121         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
122          DO j=1-3,sNy+3          DO j=1-3,sNy+3
123           DO i=1-1,sNx+1           DO i=1,sNx
124            Gi = (myXGlobalLo-1)+(bi-1)*sNx+i            Gi = (myXGlobalLo-1)+(bi-1)*sNx+i
125            Gj = (myYGlobalLo-1)+(bj-1)*sNy+j            Gj = (myYGlobalLo-1)+(bj-1)*sNy+j
126              IF (((Gj .ge. 1) .and. (Gj .le. Ny))
127         &       .or.STREAMICE_NS_PERIODIC) THEN
128    
129    
130            IF (STREAMICE_hmask(i,j,bi,bj).eq.1.0 .or.            IF (STREAMICE_hmask(i,j,bi,bj).eq.1.0 .or.
131       &        STREAMICE_hmask(i,j,bi,bj).eq.2.0) THEN       &        STREAMICE_hmask(i,j,bi,bj).eq.2.0) THEN
# Line 145  CADJ STORE trac2d  = comlev1, key=ikey_d Line 148  CADJ STORE trac2d  = comlev1, key=ikey_d
148       &      * time_step_loc       &      * time_step_loc
149  #endif  #endif
150            ENDIF            ENDIF
151              ENDIF
152           ENDDO           ENDDO
153          ENDDO          ENDDO
154         ENDDO         ENDDO
# Line 169  CADJ STORE trac2d  = comlev1, key=ikey_d Line 173  CADJ STORE trac2d  = comlev1, key=ikey_d
173    
174        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
175         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
176          DO j=1-1,sNy+1          DO j=1,sNy
177           DO i=1-1,sNx+1           DO i=1,sNx
178            Gi = (myXGlobalLo-1)+(bi-1)*sNx+i            Gi = (myXGlobalLo-1)+(bi-1)*sNx+i
179            Gj = (myYGlobalLo-1)+(bj-1)*sNy+j            Gj = (myYGlobalLo-1)+(bj-1)*sNy+j
180    
181            IF (STREAMICE_hmask(i,j,bi,bj).eq.1.0 .or.            IF (STREAMICE_hmask(i,j,bi,bj).eq.1.0 .or.
182       &        STREAMICE_hmask(i,j,bi,bj).eq.2.0) THEN       &        STREAMICE_hmask(i,j,bi,bj).eq.2.0) THEN
183             IF (Gi.eq.34.and.Gj.eq.94) THEN  !           IF (Gi.eq.34.and.Gj.eq.94) THEN
184              print *, "GOT HERE YFLUX", ytracflux(i,j,bi,bj),  !            print *, "GOT HERE YFLUX", ytracflux(i,j,bi,bj),
185       &        ytracflux(i,j+1,bi,bj),trac2d(i,j,bi,bj),  !     &        ytracflux(i,j+1,bi,bj),trac2d(i,j,bi,bj),
186       &        vtrans(i,j,bi,bj), vtrans(i,j+1,bi,bj),  !     &        vtrans(i,j,bi,bj), vtrans(i,j+1,bi,bj),
187       &        bcmasky(i,j,bi,bj)  !     &        bcmasky(i,j,bi,bj)
188             ENDIF  !           ENDIF
189    
190  #ifdef STREAMICE_TRACER_AB  #ifdef STREAMICE_TRACER_AB
191             GAD_trac_2d(i,j,bi,bj) = GAD_trac_2d(i,j,bi,bj) -             GAD_trac_2d(i,j,bi,bj) = GAD_trac_2d(i,j,bi,bj) -
# Line 207  CADJ STORE trac2d  = comlev1, key=ikey_d Line 212  CADJ STORE trac2d  = comlev1, key=ikey_d
212    
213        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
214         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
215          DO j=1-OLy,sNy+OLy          DO j=1,sNy
216           DO i=1-OLx,sNx+OLx           DO i=1,sNx
217    
218            trac2d(i,j,bi,bj) = trac2d(i,j,bi,bj) + time_step_loc *            trac2d(i,j,bi,bj) = trac2d(i,j,bi,bj) + time_step_loc *
219       &      GAD_trac_2d(i,j,bi,bj)       &      GAD_trac_2d(i,j,bi,bj)
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22