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

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

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

revision 1.11 by jmc, Tue Dec 5 22:21:50 2006 UTC revision 1.12 by jmc, Wed Apr 4 01:39:06 2007 UTC
# Line 8  C !ROUTINE: GAD_DST3_ADV_Y Line 8  C !ROUTINE: GAD_DST3_ADV_Y
8    
9  C !INTERFACE: ==========================================================  C !INTERFACE: ==========================================================
10        SUBROUTINE GAD_DST3_ADV_Y(        SUBROUTINE GAD_DST3_ADV_Y(
11       I           bi,bj,k,deltaTloc,       I           bi,bj,k, calcCFL, deltaTloc,
12       I           vTrans, vFld,       I           vTrans, vFld,
13       I           maskLocS, tracer,       I           maskLocS, tracer,
14       O           vT,       O           vT,
# Line 31  C     == Routine arguments == Line 31  C     == Routine arguments ==
31  C !INPUT PARAMETERS: ===================================================  C !INPUT PARAMETERS: ===================================================
32  C  bi,bj             :: tile indices  C  bi,bj             :: tile indices
33  C  k                 :: vertical level  C  k                 :: vertical level
34    C  calcCFL           :: =T: calculate CFL number ; =F: take vFld as CFL
35  C  deltaTloc         :: local time-step (s)  C  deltaTloc         :: local time-step (s)
36  C  vTrans            :: meridional volume transport  C  vTrans            :: meridional volume transport
37  C  vFld              :: meridional flow  C  vFld              :: meridional flow / CFL number
38  C  tracer            :: tracer field  C  tracer            :: tracer field
39  C  myThid            :: thread number  C  myThid            :: thread number
40        INTEGER bi,bj,k        INTEGER bi,bj,k
41          LOGICAL calcCFL
42        _RL deltaTloc        _RL deltaTloc
43        _RL vTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL vTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
44        _RL vFld  (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL vFld  (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
# Line 51  C  vT                :: meridional advec Line 53  C  vT                :: meridional advec
53  C     == Local variables ==  C     == Local variables ==
54  C !LOCAL VARIABLES: ====================================================  C !LOCAL VARIABLES: ====================================================
55  C  i,j               :: loop indices  C  i,j               :: loop indices
 C  vLoc              :: velocity [m/s], meridional component  
56  C  vCFL              :: Courant-Friedrich-Levy number  C  vCFL              :: Courant-Friedrich-Levy number
57        INTEGER i,j        INTEGER i,j
58        _RL Rjm,Rj,Rjp,vCFL,d0,d1        _RL Rjm,Rj,Rjp,vCFL,d0,d1
       _RL vLoc  
59  #ifdef OLD_DST3_FORMULATION  #ifdef OLD_DST3_FORMULATION
60        _RL psiP,psiM,thetaP,thetaM        _RL psiP,psiM,thetaP,thetaM
61        _RL smallNo        _RL smallNo
# Line 79  c     _RL Rjjm,Rjjp Line 79  c     _RL Rjjm,Rjjp
79          Rj =(tracer(i, j )-tracer(i,j-1))*maskLocS(i, j )          Rj =(tracer(i, j )-tracer(i,j-1))*maskLocS(i, j )
80          Rjm=(tracer(i,j-1)-tracer(i,j-2))*maskLocS(i,j-1)          Rjm=(tracer(i,j-1)-tracer(i,j-2))*maskLocS(i,j-1)
81    
82          vLoc = vFld(i,j)          vCFL = vFld(i,j)
83          vCFL = ABS( vLoc*deltaTloc          IF ( calcCFL ) vCFL = ABS( vFld(i,j)*deltaTloc
84       &                  *recip_dyC(i,j,bi,bj)*recip_deepFacC(k) )       &                  *recip_dyC(i,j,bi,bj)*recip_deepFacC(k) )
85          d0=(2.-vCFL)*(1.-vCFL)*oneSixth          d0=(2.-vCFL)*(1.-vCFL)*oneSixth
86          d1=(1.-vCFL*vCFL)*oneSixth          d1=(1.-vCFL*vCFL)*oneSixth
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22