/[MITgcm]/MITgcm/pkg/generic_advdiff/gad_advection.F
ViewVC logotype

Diff of /MITgcm/pkg/generic_advdiff/gad_advection.F

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

revision 1.8 by adcroft, Fri Sep 28 16:49:54 2001 UTC revision 1.10 by jmc, Tue Mar 5 15:17:45 2002 UTC
# Line 81  C  myThid               :: thread number Line 81  C  myThid               :: thread number
81        INTEGER bi,bj        INTEGER bi,bj
82        INTEGER advectionScheme        INTEGER advectionScheme
83        INTEGER tracerIdentity        INTEGER tracerIdentity
84        _RL Gtracer(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr,nSx,nSy)        _RL Tracer(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr,nSx,nSy)
85        _RL myTime        _RL myTime
86        INTEGER myIter        INTEGER myIter
87        INTEGER myThid        INTEGER myThid
88    
89  C !OUTPUT PARAMETERS: ==================================================  C !OUTPUT PARAMETERS: ==================================================
90  C  gTracer              :: tendancy array  C  gTracer              :: tendancy array
91        _RL Tracer(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr,nSx,nSy)        _RL gTracer(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr,nSx,nSy)
92    
93  C !LOCAL VARIABLES: ====================================================  C !LOCAL VARIABLES: ====================================================
94  C  maskUp               :: 2-D array for mask at W points  C  maskUp               :: 2-D array for mask at W points
# Line 256  CADJ STORE localTij(:,:)  = comlev1_bibj Line 256  CADJ STORE localTij(:,:)  = comlev1_bibj
256         CALL GAD_DST3FL_ADV_X(         CALL GAD_DST3FL_ADV_X(
257       &       bi,bj,k,deltaTtracer,uTrans,uVel,localTij,af,myThid)       &       bi,bj,k,deltaTtracer,uTrans,uVel,localTij,af,myThid)
258        ELSE        ELSE
259         STOP 'GAD_ADVECTION: adv. scheme incompatibale with mutli-dim'         write(0,*) advectionScheme
260           STOP 'GAD_ADVECTION: adv. scheme incompatibale with multi-dim'
261        ENDIF        ENDIF
262    
263        DO j=1-Oly,sNy+Oly        DO j=1-Oly,sNy+Oly
# Line 411  C-    Surface "correction" term at k>1 : Line 412  C-    Surface "correction" term at k>1 :
412          DO i=1-Olx,sNx+Olx          DO i=1-Olx,sNx+Olx
413           af(i,j) = af(i,j)           af(i,j) = af(i,j)
414       &           + (maskC(i,j,k,bi,bj)-maskC(i,j,k-1,bi,bj))*       &           + (maskC(i,j,k,bi,bj)-maskC(i,j,k-1,bi,bj))*
415       &             rTrans(i,j)*localTijk(i,j,k)       &             rTrans(i,j)*tracer(i,j,k,bi,bj)
416    c    &             rTrans(i,j)*localTijk(i,j,k)
417          ENDDO          ENDDO
418         ENDDO         ENDDO
419        ELSE        ELSE
420  C-    Surface "correction" term at k=1 :  C-    Surface "correction" term at k=1 :
421         DO j=1-Oly,sNy+Oly         DO j=1-Oly,sNy+Oly
422          DO i=1-Olx,sNx+Olx          DO i=1-Olx,sNx+Olx
423           af(i,j) = rTrans(i,j)*localTijk(i,j,k)           af(i,j) = rTrans(i,j)*tracer(i,j,k,bi,bj)
424    c        af(i,j) = rTrans(i,j)*localTijk(i,j,k)
425          ENDDO          ENDDO
426         ENDDO         ENDDO
427        ENDIF        ENDIF
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22