/[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.5 by jmc, Fri Jan 15 19:25:21 2010 UTC revision 1.6 by jmc, Fri Jan 15 22:39:54 2010 UTC
# Line 64  C     i,j :: Loop counters Line 64  C     i,j :: Loop counters
64        PARAMETER ( zero= 0. _d 0 , one= 1. _d 0 , half= .5 _d 0 )        PARAMETER ( zero= 0. _d 0 , one= 1. _d 0 , half= .5 _d 0 )
65  CEOP  CEOP
66    
67        IF ( buoyancyRelation .EQ. 'OCEANIC' ) THEN        IF ( usingZCoords ) 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    C--    Initialise to zero (otherwise phi0surf accumulate over land)
74           IF ( k.EQ.1 ) THEN
75             DO j=1-Oly,sNy+Oly
76              DO i=1-Olx,sNx+Olx
77                phiHydLow(i,j,bi,bj) = 0. _d 0
78              ENDDO
79             ENDDO
80           ENDIF
81    
82         IF (integr_GeoPot.EQ.1) THEN         IF (integr_GeoPot.EQ.1) THEN
83  C  --  Finite Volume Form  C  --  Finite Volume Form
84    
# Line 106  C  --  Finite Difference Form Line 115  C  --  Finite Difference Form
115  C  --  end if integr_GeoPot = ...  C  --  end if integr_GeoPot = ...
116         ENDIF         ENDIF
117    
118  C  -- end buoyancyR = Oceanic (z)  C  -- end usingZCoords
119        ENDIF        ENDIF
120    
121        IF (k.EQ.Nr) THEN        IF ( k.EQ.Nr ) THEN
122  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
123  C  --  last level (bottom): rescale (r*) and add surface contribution  C  --  last level (bottom): rescale (r*) and add surface contribution
124    
125         IF ( buoyancyRelation .EQ. 'OCEANICP' .OR.         IF ( usingPCoords ) THEN
      &      buoyancyRelation .EQ. 'ATMOSPHERIC' ) THEN  
126  C  -- P coordinate : Phi(R_low) is simply at the top :  C  -- P coordinate : Phi(R_low) is simply at the top :
127          DO j=jMin,jMax          DO j=jMin,jMax
128           DO i=iMin,iMax           DO i=iMin,iMax
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22