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

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

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

revision 1.10 by jmc, Wed Jun 22 00:27:47 2005 UTC revision 1.11 by jmc, Sat Oct 22 20:17:44 2005 UTC
# Line 223  C- a hack to prevent Water-Vapor vert.tr Line 223  C- a hack to prevent Water-Vapor vert.tr
223            CALL GAD_C2_IMPL_R( bi,bj,k, iMin,iMax,jMin,jMax,            CALL GAD_C2_IMPL_R( bi,bj,k, iMin,iMax,jMin,jMax,
224       I                        dTtracerLev, rTrans,       I                        dTtracerLev, rTrans,
225       U                        b5d, c5d, d5d,       U                        b5d, c5d, d5d,
226       I                        myThid)       I                        myThid )
227           ELSEIF (advectionScheme.EQ.ENUM_FLUX_LIMIT) THEN           ELSEIF ( advectionScheme.EQ.ENUM_UPWIND_1RST
228         &       .OR. advectionScheme.EQ.ENUM_DST2 ) THEN
229              diagonalNumber = 3
230              CALL GAD_DST2U1_IMPL_R( bi,bj,k, iMin,iMax,jMin,jMax,
231         I                        advectionScheme, dTtracerLev, rTrans,
232         U                        b5d, c5d, d5d,
233         I                        myThid )
234             ELSEIF ( advectionScheme.EQ.ENUM_FLUX_LIMIT ) THEN
235            diagonalNumber = 3            diagonalNumber = 3
236            CALL GAD_FLUXLIMIT_IMPL_R( bi,bj,k, iMin,iMax,jMin,jMax,            CALL GAD_FLUXLIMIT_IMPL_R( bi,bj,k, iMin,iMax,jMin,jMax,
237       I                        dTtracerLev, rTrans, localTijk,       I                        dTtracerLev, rTrans, localTijk,
238       U                        b5d, c5d, d5d,       U                        b5d, c5d, d5d,
239       I                        myThid)       I                        myThid )
240           ELSEIF (advectionScheme.EQ.ENUM_UPWIND_3RD .OR.           ELSEIF ( advectionScheme.EQ.ENUM_UPWIND_3RD
241       &           advectionScheme.EQ.ENUM_CENTERED_4TH) THEN       &       .OR. advectionScheme.EQ.ENUM_CENTERED_4TH
242         &       .OR. advectionScheme.EQ.ENUM_DST3 ) THEN
243            diagonalNumber = 5            diagonalNumber = 5
244            CALL GAD_U3C4_IMPL_R( bi,bj,k, iMin,iMax,jMin,jMax,            CALL GAD_U3C4_IMPL_R( bi,bj,k, iMin,iMax,jMin,jMax,
245       I                        advectionScheme, dTtracerLev, rTrans,       I                        advectionScheme, dTtracerLev, rTrans,
246       U                        a5d, b5d, c5d, d5d, e5d,       U                        a5d, b5d, c5d, d5d, e5d,
247       I                        myThid)       I                        myThid )
248             ELSEIF ( advectionScheme.EQ.ENUM_DST3_FLUX_LIMIT ) THEN
249              diagonalNumber = 5
250              CALL GAD_DST3FL_IMPL_R( bi,bj,k, iMin,iMax,jMin,jMax,
251         I                        dTtracerLev, rTrans, localTijk,
252         U                        a5d, b5d, c5d, d5d, e5d,
253         I                        myThid )
254           ELSE           ELSE
255            STOP 'GAD_IMPLICIT_R: Adv.Scheme in Impl form not yet coded'            STOP 'GAD_IMPLICIT_R: Adv.Scheme in Impl form not yet coded'
256           ENDIF           ENDIF
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22