/[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.5 by jmc, Sun Feb 9 02:01:49 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    c     IF (.FALSE.) THEN
66        IF (staggerTimeStep) THEN        IF (staggerTimeStep) THEN
67          phxFac = pfFacMom          phxFac = pfFacMom
68          phyFac = pfFacMom          phyFac = pfFacMom
# 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 (.FALSE.) THEN
83    c     IF (staggerTimeStep) THEN
84           DO j=jMin,jMax
85            DO i=iMin,iMax
86             pf(i,j) = pf(i,j)+phi_hyd(i,j,k)
87            ENDDO
88           ENDDO
89          ENDIF
90    
91  C--   Zonal velocity coriolis term  C--   Zonal velocity coriolis term
92  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.5
  ViewVC Help
Powered by ViewVC 1.1.22