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

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

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

revision 1.36 by mlosch, Mon Oct 10 05:53:49 2005 UTC revision 1.37 by jmc, Sat Oct 22 19:59:45 2005 UTC
# Line 343  CADJ &     comlev1_bibj_k_gad_pass, key= Line 343  CADJ &     comlev1_bibj_k_gad_pass, key=
343  #endif  #endif
344  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
345    
346          IF (advectionScheme.EQ.ENUM_FLUX_LIMIT) THEN          IF ( advectionScheme.EQ.ENUM_UPWIND_1RST
347         &     .OR. advectionScheme.EQ.ENUM_DST2 ) THEN
348              CALL GAD_DST2U1_ADV_X( bi,bj,k, advectionScheme,
349         I                           dTtracerLev(k),uTrans,uVel,localTij,
350         O                           af, myThid )
351            ELSEIF (advectionScheme.EQ.ENUM_FLUX_LIMIT) THEN
352            CALL GAD_FLUXLIMIT_ADV_X( bi,bj,k, dTtracerLev(k),            CALL GAD_FLUXLIMIT_ADV_X( bi,bj,k, dTtracerLev(k),
353       I                              uTrans, uVel, maskLocW, localTij,       I                              uTrans, uVel, maskLocW, localTij,
354       O                              af, myThid )       O                              af, myThid )
# Line 485  CADJ &     comlev1_bibj_k_gad_pass, key= Line 490  CADJ &     comlev1_bibj_k_gad_pass, key=
490  #endif  #endif
491  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
492    
493          IF (advectionScheme.EQ.ENUM_FLUX_LIMIT) THEN          IF ( advectionScheme.EQ.ENUM_UPWIND_1RST
494         &     .OR. advectionScheme.EQ.ENUM_DST2 ) THEN
495              CALL GAD_DST2U1_ADV_Y( bi,bj,k, advectionScheme,
496         I                           dTtracerLev(k),vTrans,vVel,localTij,
497         O                           af, myThid )
498            ELSEIF (advectionScheme.EQ.ENUM_FLUX_LIMIT) THEN
499            CALL GAD_FLUXLIMIT_ADV_Y( bi,bj,k, dTtracerLev(k),            CALL GAD_FLUXLIMIT_ADV_Y( bi,bj,k, dTtracerLev(k),
500       I                              vTrans, vVel, maskLocS, localTij,       I                              vTrans, vVel, maskLocS, localTij,
501       O                              af, myThid )       O                              af, myThid )
# Line 697  CADJ &     = comlev1_bibj_k_gad, key=kke Line 707  CADJ &     = comlev1_bibj_k_gad, key=kke
707  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
708    
709  C-    Compute vertical advective flux in the interior:  C-    Compute vertical advective flux in the interior:
710           IF (vertAdvecScheme.EQ.ENUM_FLUX_LIMIT) THEN           IF ( advectionScheme.EQ.ENUM_UPWIND_1RST
711  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|       &      .OR. advectionScheme.EQ.ENUM_DST2 ) THEN
712               CALL GAD_DST2U1_ADV_R( bi,bj,k, advectionScheme,
713         I                            dTtracerLev(k),rTrans,wVel,localTijk,
714         O                            fVerT(1-Olx,1-Oly,kUp), myThid )
715             ELSEIF (advectionScheme.EQ.ENUM_FLUX_LIMIT) THEN
716             CALL GAD_FLUXLIMIT_ADV_R( bi,bj,k, dTtracerLev(k),             CALL GAD_FLUXLIMIT_ADV_R( bi,bj,k, dTtracerLev(k),
717       I                               rTrans, wVel, localTijk,       I                               rTrans, wVel, localTijk,
718       O                               fVerT(1-Olx,1-Oly,kUp), myThid )       O                               fVerT(1-Olx,1-Oly,kUp), myThid )
Legend:
Removed from v.1.36  
changed lines
  Added in v.1.37
  ViewVC Help
Powered by ViewVC 1.1.22