/[MITgcm]/MITgcm/pkg/generic_advdiff/gad_dst3_adv_y.F

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

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

-revision 1.1 by adcroft,
Tue Sep  4 14:53:11 2001 UTC
+revision 1.10 by jmc,
Sun Oct 22 01:08:04 2006 UTC
 Line 3 
 C $Name$
  #include "GAD_OPTIONS.h"
-       SUBROUTINE GAD_DST3_ADV_Y(
+ CBOP
-      I           bi,bj,k,deltaT,
+ C !ROUTINE: GAD_DST3_ADV_Y
-      I           vTrans, vVel,
-      I           tracer,
+ C !INTERFACE: ==========================================================
+       SUBROUTINE GAD_DST3_ADV_Y(
+      I           bi,bj,k,deltaTloc,
+      I           vTrans, vFld,
+      I           maskLocS, tracer,
       O           vT,
       I           myThid )
- C     /==========================================================\
+ C !DESCRIPTION:
- C     | SUBROUTINE GAD_DST3_ADV_Y                                |
+ C  Calculates the area integrated Meridional flux due to advection of a
- C     | o Compute Meridional advective Flux of Tracer using      |
+ C  tracer using 3rd-order Direct Space and Time (DST-3) Advection Scheme
- C     |   3rd Order DST Sceheme                                  |
- C     |==========================================================|
+ C !USES: ===============================================================
        IMPLICIT NONE
  C     == GLobal variables ==
  #include "SIZE.h"
  #include "GRID.h"
+ #include "EEPARAMS.h"
+ #include "PARAMS.h"
  #include "GAD.h"
  C     == Routine arguments ==
+ C !INPUT PARAMETERS: ===================================================
+ C  bi,bj             :: tile indices
+ C  k                 :: vertical level
+ C  deltaTloc         :: local time-step (s)
+ C  vTrans            :: meridional volume transport
+ C  vFld              :: meridional flow
+ C  tracer            :: tracer field
+ C  myThid            :: thread number
        INTEGER bi,bj,k
-       _RL deltaT
+       _RL deltaTloc
        _RL vTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
-       _RL vVel(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
+       _RL vFld  (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
+       _RS maskLocS(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
        _RL tracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
-       _RL vT    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
        INTEGER myThid
+ C !OUTPUT PARAMETERS: ==================================================
+ C  vT                :: meridional advective flux
+       _RL vT    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
  C     == Local variables ==
+ C !LOCAL VARIABLES: ====================================================
+ C  i,j               :: loop indices
+ C  vLoc              :: velocity [m/s], meridional component
+ C  cfl               :: Courant-Friedrich-Levy number
        INTEGER i,j
        _RL Rjm,Rj,Rjp,cfl,d0,d1
+       _RL vLoc
+ #ifdef OLD_DST3_FORMULATION
+       _RL psiP,psiM,thetaP,thetaM
+       _RL smallNo
+ c     _RL Rjjm,Rjjp
+       IF (inAdMode) THEN
+        smallNo = 1.0D-20
+       ELSE
+        smallNo = 1.0D-20
+       ENDIF
+ #endif
        DO i=1-Olx,sNx+Olx
         vT(i,1-Oly)=0.
-Line 41 
 C     == Local variables ==
+Line 75 
 C     == Local variables ==
        ENDDO
        DO j=1-Oly+2,sNy+Oly-1
         DO i=1-Olx,sNx+Olx
-         Rjp=(tracer(i,j+1)-tracer(i,j))*maskS(i,j+1,k,bi,bj)
+         Rjp=(tracer(i,j+1)-tracer(i, j ))*maskLocS(i,j+1)
-         Rj =(tracer(i,j)-tracer(i,j-1))*maskS(i,j,k,bi,bj)
+         Rj =(tracer(i, j )-tracer(i,j-1))*maskLocS(i, j )
-         Rjm=(tracer(i,j-1)-tracer(i,j-2))*maskS(i,j-1,k,bi,bj)
+         Rjm=(tracer(i,j-1)-tracer(i,j-2))*maskLocS(i,j-1)
-         cfl=vVel(i,j,k,bi,bj)*deltaT*recip_dyc(i,j,bi,bj)
+         vLoc = vFld(i,j)
-         d0=(2.-abs(cfl))*(1.-abs(cfl))*oneSixth
+ c       vLoc = vTrans(i,j)*recip_dxG(i,j,bi,bj)
-         d1=(1.-cfl)*(1.+cfl)*oneSixth
+ c    &       *recip_drF(k)*_recip_hFacS(i,j,k,bi,bj)
+         cfl=ABS(vLoc*deltaTloc*recip_dyC(i,j,bi,bj))
+         d0=(2.-cfl)*(1.-cfl)*oneSixth
+         d1=(1.-cfl*cfl)*oneSixth
+ #ifdef OLD_DST3_FORMULATION
+         IF ( ABS(Rj).LT.smallNo .OR.
+      &       ABS(Rjm).LT.smallNo ) THEN
+          thetaP=0.
+          psiP=0.
+         ELSE
+          thetaP=(Rjm+smallNo)/(smallNo+Rj)
+          psiP=d0+d1*thetaP
+         ENDIF
+         IF ( ABS(Rj).LT.smallNo .OR.
+      &       ABS(Rjp).LT.smallNo ) THEN
+          thetaM=0.
+          psiM=0.
+         ELSE
+          thetaM=(Rjp+smallNo)/(smallNo+Rj)
+          psiM=d0+d1*thetaM
+         ENDIF
+         vT(i,j)=
+      &   0.5*(vTrans(i,j)+ABS(vTrans(i,j)))
+      &      *( Tracer(i,j-1) + psiP*Rj )
+      &  +0.5*(vTrans(i,j)-ABS(vTrans(i,j)))
+      &      *( Tracer(i, j ) - psiM*Rj )
+ #else /* OLD_DST3_FORMULATION */
          vT(i,j)=
-      &   0.5*(vTrans(i,j)+abs(vTrans(i,j)))
+      &   0.5*(vTrans(i,j)+ABS(vTrans(i,j)))
-      &      *( Tracer(i,j-1) + d0*Rj + d1*Rjm )
+      &      *( Tracer(i,j-1) + (d0*Rj+d1*Rjm) )
-      &  +0.5*(vTrans(i,j)-abs(vTrans(i,j)))
+      &  +0.5*(vTrans(i,j)-ABS(vTrans(i,j)))
-      &      *( Tracer(i, j ) - d0*Rj + d1*Rjp )
+      &      *( Tracer(i, j ) - (d0*Rj+d1*Rjp) )
+ #endif /* OLD_DST3_FORMULATION */
         ENDDO
        ENDDO

 Legend:



Removed from v.1.1
 


changed lines


 
Added in v.1.10
 Legend:



Removed from v.1.1
 


changed lines


 
Added in v.1.10
-Removed from v.1.1
+Added in v.1.10

	ViewVC Help
Powered by ViewVC 1.1.22