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

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

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

revision 1.6 by jmc, Tue Dec 5 22:21:50 2006 UTC revision 1.7 by jmc, Wed Apr 4 01:39:06 2007 UTC
# Line 8  C !ROUTINE: GAD_DST2U1_ADV_X Line 8  C !ROUTINE: GAD_DST2U1_ADV_X
8    
9  C !INTERFACE: ==========================================================  C !INTERFACE: ==========================================================
10        SUBROUTINE GAD_DST2U1_ADV_X(        SUBROUTINE GAD_DST2U1_ADV_X(
11       I           bi,bj,k, advectionScheme,       I           bi,bj,k, advectionScheme, calcCFL,
12       I           deltaTloc, uTrans, uFld,       I           deltaTloc, uTrans, uFld,
13       I           tracer,       I           tracer,
14       O           uT,       O           uT,
# Line 30  C  bi,bj             :: tile indices Line 30  C  bi,bj             :: tile indices
30  C  k                 :: vertical level  C  k                 :: vertical level
31  C  advectionScheme   :: advection scheme to use: either 2nd Order DST  C  advectionScheme   :: advection scheme to use: either 2nd Order DST
32  C                                                or 1rst Order Upwind  C                                                or 1rst Order Upwind
33    C  calcCFL           :: =T: calculate CFL number ; =F: take uFld as CFL
34    C  deltaTloc         :: local time-step (s)
35  C  uTrans            :: zonal volume transport  C  uTrans            :: zonal volume transport
36  C  uFld              :: zonal flow  C  uFld              :: zonal flow / CFL number
37  C  tracer            :: tracer field  C  tracer            :: tracer field
38  C  myThid            :: thread number  C  myThid            :: thread number
39        INTEGER bi,bj,k        INTEGER bi,bj,k
40        INTEGER advectionScheme        INTEGER advectionScheme
41          LOGICAL calcCFL
42        _RL deltaTloc        _RL deltaTloc
43        _RL uTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL uTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
44        _RL uFld  (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL uFld  (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
# Line 49  C  uT                :: zonal advective Line 52  C  uT                :: zonal advective
52  C !LOCAL VARIABLES: ====================================================  C !LOCAL VARIABLES: ====================================================
53  C  i,j               :: loop indices  C  i,j               :: loop indices
54  C  rLimit            :: centered (vs upwind) fraction  C  rLimit            :: centered (vs upwind) fraction
 C  uLoc              :: velocity [m/s], zonal component  
55  C  uCFL              :: Courant-Friedrich-Levy number  C  uCFL              :: Courant-Friedrich-Levy number
56        INTEGER i,j        INTEGER i,j
57        _RL uLoc, uCFL, xLimit, uAbs        _RL uCFL, xLimit, uAbs
58  CEOP  CEOP
59    
60        xLimit = 0. _d 0        xLimit = 0. _d 0
# Line 62  CEOP Line 64  CEOP
64         uT(1-Olx,j)=0.         uT(1-Olx,j)=0.
65         DO i=1-Olx+1,sNx+Olx         DO i=1-Olx+1,sNx+Olx
66    
67          uLoc = uFld(i,j)          uCFL = uFld(i,j)
68          uCFL = ABS( uLoc*deltaTloc          IF ( calcCFL ) uCFL = ABS( uFld(i,j)*deltaTloc
69       &                  *recip_dxC(i,j,bi,bj)*recip_deepFacC(k) )       &                  *recip_dxC(i,j,bi,bj)*recip_deepFacC(k) )
70    
71  c       uT(i,j) =  c       uT(i,j) =
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22