/[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.3 by adcroft, Thu Aug 30 00:40:37 2001 UTC revision 1.4 by adcroft, Tue Sep 4 14:53:11 2001 UTC
# Line 103  C-    Advective flux in X Line 103  C-    Advective flux in X
103         CALL GAD_U3_ADV_X(bi,bj,k,uTrans,localT,af,myThid)         CALL GAD_U3_ADV_X(bi,bj,k,uTrans,localT,af,myThid)
104        ELSEIF (advectionScheme.EQ.ENUM_CENTERED_4TH) THEN        ELSEIF (advectionScheme.EQ.ENUM_CENTERED_4TH) THEN
105         CALL GAD_C4_ADV_X(bi,bj,k,uTrans,localT,af,myThid)         CALL GAD_C4_ADV_X(bi,bj,k,uTrans,localT,af,myThid)
106          ELSEIF (advectionScheme.EQ.ENUM_DST3 ) THEN
107           CALL GAD_DST3_ADV_X(
108         &       bi,bj,k,deltaTtracer,uTrans,uVel,localT,af,myThid)
109          ELSEIF (advectionScheme.EQ.ENUM_DST3_FLUX_LIMIT ) THEN
110           CALL GAD_DST3FL_ADV_X(
111         &       bi,bj,k,deltaTtracer,uTrans,uVel,localT,af,myThid)
112        ELSE        ELSE
113         STOP 'GAD_CALC_RHS: Bad advectionScheme (X)'         STOP 'GAD_CALC_RHS: Bad advectionScheme (X)'
114        ENDIF        ENDIF
# Line 167  C-    Advective flux in Y Line 173  C-    Advective flux in Y
173         CALL GAD_U3_ADV_Y(bi,bj,k,vTrans,localT,af,myThid)         CALL GAD_U3_ADV_Y(bi,bj,k,vTrans,localT,af,myThid)
174        ELSEIF (advectionScheme.EQ.ENUM_CENTERED_4TH) THEN        ELSEIF (advectionScheme.EQ.ENUM_CENTERED_4TH) THEN
175         CALL GAD_C4_ADV_Y(bi,bj,k,vTrans,localT,af,myThid)         CALL GAD_C4_ADV_Y(bi,bj,k,vTrans,localT,af,myThid)
176          ELSEIF (advectionScheme.EQ.ENUM_DST3 ) THEN
177           CALL GAD_DST3_ADV_Y(
178         &       bi,bj,k,deltaTtracer,vTrans,vVel,localT,af,myThid)
179          ELSEIF (advectionScheme.EQ.ENUM_DST3_FLUX_LIMIT ) THEN
180           CALL GAD_DST3FL_ADV_Y(
181         &       bi,bj,k,deltaTtracer,vTrans,vVel,localT,af,myThid)
182        ELSE        ELSE
183         STOP 'GAD_CALC_RHS: Bad advectionScheme (Y)'         STOP 'GAD_CALC_RHS: Bad advectionScheme (Y)'
184        ENDIF        ENDIF
# Line 234  C-    Compute vertical advective flux in Line 246  C-    Compute vertical advective flux in
246          CALL GAD_U3_ADV_R(bi,bj,k,rTrans,tracer,af,myThid)          CALL GAD_U3_ADV_R(bi,bj,k,rTrans,tracer,af,myThid)
247         ELSEIF (advectionScheme.EQ.ENUM_CENTERED_4TH) THEN         ELSEIF (advectionScheme.EQ.ENUM_CENTERED_4TH) THEN
248          CALL GAD_C4_ADV_R(bi,bj,k,rTrans,tracer,af,myThid)          CALL GAD_C4_ADV_R(bi,bj,k,rTrans,tracer,af,myThid)
249  c       CALL GAD_C2_ADV_R(bi,bj,k,rTrans,tracer,af,myThid)         ELSEIF (advectionScheme.EQ.ENUM_DST3 ) THEN
250    c       CALL GAD_DST3_ADV_R(
251    c    &       bi,bj,k,deltaTtracer,rTrans,wVel,tracer,af,myThid)
252            STOP 'GAD_CALC_RHS: GAD_DST3_ADV_R not coded yet'
253           ELSEIF (advectionScheme.EQ.ENUM_DST3_FLUX_LIMIT ) THEN
254    c       CALL GAD_DST3FL_ADV_R(
255    c    &       bi,bj,k,deltaTtracer,rTrans,wVel,tracer,af,myThid)
256            STOP 'GAD_CALC_RHS: GAD_DST3FL_ADV_R not coded yet'
257         ELSE         ELSE
258          STOP 'GAD_CALC_RHS: Bad advectionScheme (R)'          STOP 'GAD_CALC_RHS: Bad advectionScheme (R)'
259         ENDIF         ENDIF
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22