10 |
I myThid) |
I myThid) |
11 |
C /==========================================================\ |
C /==========================================================\ |
12 |
C | SUBROUTINE CALC_PHI_HYD | |
C | SUBROUTINE CALC_PHI_HYD | |
13 |
C | o Integrate the hydrostatic relation to find phiHyd. | |
C | o Integrate the hydrostatic relation to find the Hydros. | |
14 |
C | | |
C | Potential (ocean: Pressure/rho ; atmos = geopotential)| |
15 |
C | On entry: | |
C | On entry: | |
16 |
C | theta,salt are the current thermodynamics quantities| |
C | theta,salt are the current thermodynamics quantities| |
17 |
C | (unchanged on exit) | |
C | (unchanged on exit) | |
18 |
C | phiHyd(i,j,1:k-1) is the hydrostatic pressure/geopot. | |
C | phiHyd(i,j,1:k-1) is the hydrostatic Potential | |
19 |
C | at cell centers (tracer points) | |
C | at cell centers (tracer points) | |
20 |
C | - 1:k-1 layers are valid | |
C | - 1:k-1 layers are valid | |
21 |
C | - k:Nr layers are invalid | |
C | - k:Nr layers are invalid | |
22 |
C | phiHyd(i,j,k) is the hydrostatic pressure/geop. | |
C | phiHyd(i,j,k) is the hydrostatic Potential | |
23 |
C | at cell the interface k (w point above) | |
C | at cell the interface k (w point above) | |
24 |
C | On exit: | |
C | On exit: | |
25 |
C | phiHyd(i,j,1:k) is the hydrostatic pressure/geopot. | |
C | phiHyd(i,j,1:k) is the hydrostatic Potential | |
26 |
C | at cell centers (tracer points) | |
C | at cell centers (tracer points) | |
27 |
C | - 1:k layers are valid | |
C | - 1:k layers are valid | |
28 |
C | - k+1:Nr layers are invalid | |
C | - k+1:Nr layers are invalid | |
29 |
C | phiHyd(i,j,k+1) is the hydrostatic pressure/geop. | |
C | phiHyd(i,j,k+1) is the hydrostatic Potential (P/rho) | |
30 |
C | at cell the interface k+1 (w point below)| |
C | at cell the interface k+1 (w point below)| |
31 |
C | | |
C | | |
32 |
C \==========================================================/ |
C \==========================================================/ |
94 |
C conserve KE+PE exactly even though it is more natural |
C conserve KE+PE exactly even though it is more natural |
95 |
C |
C |
96 |
c IF (k.LT.Nr) phiHyd(i,j,k+1)=phiHyd(i,j,k)+ |
c IF (k.LT.Nr) phiHyd(i,j,k+1)=phiHyd(i,j,k)+ |
97 |
c & drF(K)*gravity*alphaRho(i,j) |
c & drF(K)*gravity*alphaRho(i,j)*recip_rhoConst |
98 |
c phiHyd(i,j,k)=phiHyd(i,j,k)+ |
c phiHyd(i,j,k)=phiHyd(i,j,k)+ |
99 |
c & 0.5*drF(K)*gravity*alphaRho(i,j) |
c & 0.5*drF(K)*gravity*alphaRho(i,j)*recip_rhoConst |
100 |
C----------------------------------------------------------------------- |
C----------------------------------------------------------------------- |
101 |
|
|
102 |
C---------- This discretization is the "energy conserving" form |
C---------- This discretization is the "energy conserving" form |
103 |
C which has been used since at least Adcroft et al., MWR 1997 |
C which has been used since at least Adcroft et al., MWR 1997 |
104 |
C |
C |
105 |
phiHyd(i,j,k)=phiHyd(i,j,k)+ |
phiHyd(i,j,k)=phiHyd(i,j,k)+ |
106 |
& 0.5*dRloc*gravity*alphaRho(i,j) |
& 0.5*dRloc*gravity*alphaRho(i,j)*recip_rhoConst |
107 |
IF (k.LT.Nr) phiHyd(i,j,k+1)=phiHyd(i,j,k)+ |
IF (k.LT.Nr) phiHyd(i,j,k+1)=phiHyd(i,j,k)+ |
108 |
& 0.5*dRlocKp1*gravity*alphaRho(i,j) |
& 0.5*dRlocKp1*gravity*alphaRho(i,j)*recip_rhoConst |
109 |
C----------------------------------------------------------------------- |
C----------------------------------------------------------------------- |
110 |
ENDDO |
ENDDO |
111 |
ENDDO |
ENDDO |
112 |
|
|
113 |
|
|
114 |
|
|
|
|
|
115 |
ELSEIF ( buoyancyRelation .eq. 'ATMOSPHERIC' ) THEN |
ELSEIF ( buoyancyRelation .eq. 'ATMOSPHERIC' ) THEN |
116 |
C This is the hydrostatic geopotential calculation for the Atmosphere |
C This is the hydrostatic geopotential calculation for the Atmosphere |
117 |
C The ideal gas law is used implicitly here rather than calculating |
C The ideal gas law is used implicitly here rather than calculating |
184 |
|
|
185 |
#endif |
#endif |
186 |
|
|
187 |
return |
RETURN |
188 |
end |
END |