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

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

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

revision 1.23 by jmc, Wed Jan 7 21:35:00 2004 UTC revision 1.24 by edhill, Mon Mar 29 03:33:51 2004 UTC
# Line 51  C !USES: =============================== Line 51  C !USES: ===============================
51  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
52    
53  C !INPUT PARAMETERS: ===================================================  C !INPUT PARAMETERS: ===================================================
54  C  bi,bj                :: tile indices  C bi,bj            :: tile indices
55  C  iMin,iMax,jMin,jMax  :: loop range for called routines  C iMin,iMax        :: loop range for called routines
56  C  kup                  :: index into 2 1/2D array, toggles between 1 and 2  C jMin,jMax        :: loop range for called routines
57  C  kdown                :: index into 2 1/2D array, toggles between 2 and 1  C kup              :: index into 2 1/2D array, toggles between 1|2
58  C  kp1                  :: =k+1 for k<Nr, =Nr for k=Nr  C kdown            :: index into 2 1/2D array, toggles between 2|1
59  C  xA,yA                :: areas of X and Y face of tracer cells  C kp1              :: =k+1 for k<Nr, =Nr for k=Nr
60  C  uTrans,vTrans,rTrans :: 2-D arrays of volume transports at U,V and W points  C xA,yA            :: areas of X and Y face of tracer cells
61  C  rTransKp1            :: 2-D array of volume transport at W pt, interface k+1  C uTrans,vTrans    :: 2-D arrays of volume transports at U,V points
62  C  maskUp               :: 2-D array for mask at W points  C rTrans           :: 2-D arrays of volume transports at W points
63  C  uVel, vVel, wVel     :: 3 components of the velcity field (3-D array)  C rTransKp1        :: 2-D array of volume trans at W pts, interf k+1
64  C  diffKh               :: horizontal diffusion coefficient  C maskUp           :: 2-D array for mask at W points
65  C  diffK4               :: bi-harmonic diffusion coefficient  C uVel,vVel,wVel   :: 3 components of the velcity field (3-D array)
66  C  KappaRT              :: 3-D array for vertical diffusion coefficient  C diffKh           :: horizontal diffusion coefficient
67  C  Tracer               :: tracer field  C diffK4           :: bi-harmonic diffusion coefficient
68  C  tracerIdentity       :: identifier for the tracer (required for KPP and GM)  C KappaRT          :: 3-D array for vertical diffusion coefficient
69  C  advectionScheme      :: advection scheme to use  C Tracer           :: tracer field
70  C  calcAdvection        :: =False if Advec terms computed with multiDim scheme  C tracerIdentity   :: tracer identifier (required for KPP,GM)
71  C  implicitAdvection    :: =True if vertical Advec term is computed implicitly  C advectionScheme  :: advection scheme to use
72  C  myThid               :: thread number  C calcAdvection    :: =False if Advec computed with multiDim scheme
73    C implicitAdvection:: =True if vertical Advec computed implicitly
74    C myThid           :: thread number
75        INTEGER bi,bj,iMin,iMax,jMin,jMax        INTEGER bi,bj,iMin,iMax,jMin,jMax
76        INTEGER k,kUp,kDown,kM1        INTEGER k,kUp,kDown,kM1
77        _RS xA    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RS xA    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
# Line 92  C  myThid               :: thread number Line 94  C  myThid               :: thread number
94        INTEGER myThid        INTEGER myThid
95    
96  C !OUTPUT PARAMETERS: ==================================================  C !OUTPUT PARAMETERS: ==================================================
97  C  gTracer              :: tendancy array  C gTracer          :: tendancy array
98  C  fVerT                :: 2 1/2D arrays for vertical advective flux  C fVerT            :: 2 1/2D arrays for vertical advective flux
99        _RL gTracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL gTracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
100        _RL fVerT (1-OLx:sNx+OLx,1-OLy:sNy+OLy,2)        _RL fVerT (1-OLx:sNx+OLx,1-OLy:sNy+OLy,2)
101    
102  C !LOCAL VARIABLES: ====================================================  C !LOCAL VARIABLES: ====================================================
103  C  i,j                  :: loop indices  C i,j              :: loop indices
104  C  df4                  :: used for storing del^2 T for bi-harmonic term  C df4              :: used for storing del^2 T for bi-harmonic term
105  C  fZon                 :: zonal flux  C fZon             :: zonal flux
106  C  fmer                 :: meridional flux  C fmer             :: meridional flux
107  C  af                   :: advective flux  C af               :: advective flux
108  C  df                   :: diffusive flux  C df               :: diffusive flux
109  C  localT               :: local copy of tracer field  C localT           :: local copy of tracer field
110        INTEGER i,j        INTEGER i,j
111        _RL df4   (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL df4   (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
112        _RL fZon  (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL fZon  (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
Legend:
Removed from v.1.23  
changed lines
  Added in v.1.24
  ViewVC Help
Powered by ViewVC 1.1.22