/[MITgcm]/MITgcm/model/src/diags_phi_rlow.F
ViewVC logotype

Diff of /MITgcm/model/src/diags_phi_rlow.F

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

revision 1.4 by gforget, Mon Oct 26 21:48:43 2009 UTC revision 1.5 by jmc, Fri Jan 15 19:25:21 2010 UTC
# Line 6  C $Name$ Line 6  C $Name$
6  CBOP  CBOP
7  C     !ROUTINE: DIAGS_PHI_RLOW  C     !ROUTINE: DIAGS_PHI_RLOW
8  C     !INTERFACE:  C     !INTERFACE:
9        SUBROUTINE DIAGS_PHI_RLOW(        SUBROUTINE DIAGS_PHI_RLOW(
10       I                       k, bi, bj, iMin,iMax, jMin,jMax,       I                       k, bi, bj, iMin,iMax, jMin,jMax,
11       I                       phiHydF, phiHydC, alphRho, tFld, sFld,       I                       phiHydF, phiHydC, alphRho, tFld, sFld,
12       I                       myTime, myIter, myThid)       I                       myTime, myIter, myThid)
13  C     !DESCRIPTION: \bv  C     !DESCRIPTION: \bv
14  C     *==========================================================*  C     *==========================================================*
15  C     | S/R DIAGS_PHI_RLOW                                      C     | S/R DIAGS_PHI_RLOW
16  C     | o Diagnose Phi-Hydrostatic at r-lower boundary  C     | o Diagnose Phi-Hydrostatic at r-lower boundary
17  C     |   = bottom pressure (ocean in z-coord) ;  C     |   = bottom pressure (ocean in z-coord) ;
18  C     |   = sea surface elevation (ocean in p-coord) ;  C     |   = sea surface elevation (ocean in p-coord) ;
19  C     |   = height at the top of atmosphere (in p-coord) ;  C     |   = height at the top of atmosphere (in p-coord) ;
20  C     *==========================================================*  C     *==========================================================*
# Line 32  C     == Global variables == Line 32  C     == Global variables ==
32    
33  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
34  C     == Routine Arguments ==  C     == Routine Arguments ==
35  C     bi,bj      :: tile index              C     bi,bj      :: tile index
36  C     iMin,iMax,jMin,jMax :: Loop counters  C     iMin,iMax,jMin,jMax :: Loop counters
37  C     phiHydF    :: hydrostatic potential anomaly at middle between  C     phiHydF    :: hydrostatic potential anomaly at middle between
38  C                   2 centers k & k+1 (interface k+1)  C                   2 centers k & k+1 (interface k+1)
39  C     phiHydC    :: hydrostatic potential anomaly at cell center  C     phiHydC    :: hydrostatic potential anomaly at cell center
40  C                  (atmos: =Geopotential ; ocean-z: =Pressure/rho)  C                  (atmos: =Geopotential ; ocean-z: =Pressure/rho)
41  C     alphRho    :: Density (z-coord) or specific volume (p-coord)  C     alphRho    :: Density (z-coord) or specific volume (p-coord)
42  C     tFld       :: Potential temp.  C     tFld       :: Potential temp.
43  C     sFld       :: Salinity  C     sFld       :: Salinity
44  C     myTime :: Current time  C     myTime     :: Current time
45  C     myIter :: Current iteration number  C     myIter     :: Current iteration number
46  C     myThid :: Instance number for this call of the routine.  C     myThid     :: my Thread Id number
47        INTEGER k, bi,bj, iMin,iMax, jMin,jMax        INTEGER k, bi,bj, iMin,iMax, jMin,jMax
48        _RL phiHydF(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL phiHydF(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
49        _RL phiHydC(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL phiHydC(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
# Line 67  CEOP Line 67  CEOP
67        IF ( buoyancyRelation .EQ. 'OCEANIC' ) THEN        IF ( buoyancyRelation .EQ. 'OCEANIC' ) THEN
68    
69  C----- Compute bottom pressure deviation from gravity*rho0*H  C----- Compute bottom pressure deviation from gravity*rho0*H
70  C      Start from phiHyd at the (bottom) tracer point and add Del_h*g*rho'  C      Start from phiHyd at the (bottom) tracer point and add Del_h*g*rho'
71  C      with Del_h = distance from the bottom up to tracer point  C      with Del_h = distance from the bottom up to tracer point
72    
73         IF (integr_GeoPot.EQ.1) THEN         IF (integr_GeoPot.EQ.1) THEN
# Line 105  C  --  Finite Difference Form Line 105  C  --  Finite Difference Form
105    
106  C  --  end if integr_GeoPot = ...  C  --  end if integr_GeoPot = ...
107         ENDIF         ENDIF
108            
109  C  -- end buoyancyR = Oceanic (z)  C  -- end buoyancyR = Oceanic (z)
110        ENDIF        ENDIF
111    
# Line 123  C  -- P coordinate : Phi(R_low) is simpl Line 123  C  -- P coordinate : Phi(R_low) is simpl
123          ENDDO          ENDDO
124         ENDIF         ENDIF
125    
126    #ifdef NONLIN_FRSURF
127    c      IF ( select_rStar.GE.2 .AND. nonlinFreeSurf.GE.4 ) THEN
128           IF ( select_rStar.GT.0 .AND. nonlinFreeSurf.GE.4 ) THEN
129            DO j=jMin,jMax
130             DO i=iMin,iMax
131               phiHydLow(i,j,bi,bj) = phiHydLow(i,j,bi,bj)
132         &                          *rStarFacC(i,j,bi,bj)
133             ENDDO
134            ENDDO
135           ENDIF
136    #endif /* NONLIN_FRSURF */
137    
138          DO j=jMin,jMax          DO j=jMin,jMax
139           DO i=iMin,iMax           DO i=iMin,iMax
140             phiHydLow(i,j,bi,bj) = phiHydLow(i,j,bi,bj)             phiHydLow(i,j,bi,bj) = phiHydLow(i,j,bi,bj)
141       &            + Bo_surf(i,j,bi,bj)*etaN(i,j,bi,bj)       &            + Bo_surf(i,j,bi,bj)*etaN(i,j,bi,bj)
142       &            + phi0surf(i,j,bi,bj)       &            + phi0surf(i,j,bi,bj)
143           ENDDO           ENDDO
144          ENDDO          ENDDO
145    
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22