/[MITgcm]/MITgcm/pkg/mom_fluxform/mom_cdscheme.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_fluxform/mom_cdscheme.F

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

revision 1.3.6.1 by heimbach, Fri Mar 7 04:46:40 2003 UTC revision 1.4 by jmc, Sat Feb 8 02:13:02 2003 UTC
# Line 8  C !ROUTINE: MOM_CDSCHEME Line 8  C !ROUTINE: MOM_CDSCHEME
8    
9  C !INTERFACE: ==========================================================  C !INTERFACE: ==========================================================
10        SUBROUTINE MOM_CDSCHEME(        SUBROUTINE MOM_CDSCHEME(
11       I        bi,bj,k,dPhiHydX,dPhiHydY,       I        bi,bj,k,phi_hyd,dPhiHydX,dPhiHydY,
12       I        myThid)       I        myThid)
13    
14  C !DESCRIPTION:  C !DESCRIPTION:
# Line 27  C     == Global variables == Line 27  C     == Global variables ==
27  C !INPUT PARAMETERS: ===================================================  C !INPUT PARAMETERS: ===================================================
28  C  bi,bj                :: tile indices  C  bi,bj                :: tile indices
29  C  k                    :: vertical level  C  k                    :: vertical level
30    C  phi_hyd              :: hydrostatic pressure
31  C     dPhiHydX,Y        :: Gradient (X & Y dir.) of Hydrostatic Potential  C     dPhiHydX,Y        :: Gradient (X & Y dir.) of Hydrostatic Potential
32  C  myThid               :: thread number  C  myThid               :: thread number
33        INTEGER bi,bj,k        INTEGER bi,bj,k
34          _RL phi_hyd(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
35        _RL dPhiHydX(1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RL dPhiHydX(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
36        _RL dPhiHydY(1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RL dPhiHydY(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
37        INTEGER myThid        INTEGER myThid
# Line 60  C     Adams-Bashforth weighting factors Line 62  C     Adams-Bashforth weighting factors
62        ab05   = -0.5 - abEps        ab05   = -0.5 - abEps
63    
64  C-- stagger time stepping: grad Phi_Hyp is not in gU,gV and needs to be added:  C-- stagger time stepping: grad Phi_Hyp is not in gU,gV and needs to be added:
65        IF (staggerTimeStep) THEN        IF (.FALSE.) THEN
66    c     IF (staggerTimeStep) THEN
67          phxFac = pfFacMom          phxFac = pfFacMom
68          phyFac = pfFacMom          phyFac = pfFacMom
69        ELSE        ELSE
# Line 76  C     Pressure extrapolated forward in t Line 79  C     Pressure extrapolated forward in t
79       &  +ab05*(etaNm1(i,j,bi,bj)*Bo_surf(i,j,bi,bj) )       &  +ab05*(etaNm1(i,j,bi,bj)*Bo_surf(i,j,bi,bj) )
80         ENDDO         ENDDO
81        ENDDO        ENDDO
82          IF (staggerTimeStep) THEN
83           DO j=jMin,jMax
84            DO i=iMin,iMax
85             pf(i,j) = pf(i,j)+phi_hyd(i,j,k)
86            ENDDO
87           ENDDO
88          ENDIF
89    
90  C--   Zonal velocity coriolis term  C--   Zonal velocity coriolis term
91  C     Note. As coded here, coriolis will not work with "thin walls"  C     Note. As coded here, coriolis will not work with "thin walls"
Legend:
Removed from v.1.3.6.1  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22