94 |
phiRef(K) = phiRef(K)*recip_gravity |
phiRef(K) = phiRef(K)*recip_gravity |
95 |
ENDDO |
ENDDO |
96 |
|
|
97 |
|
_BEGIN_MASTER( myThid ) |
98 |
C- Write to check : |
C- Write to check : |
99 |
write(6,'(2A)') 'INI_P_GROUND: PhiRef/g [m] at Center ', |
WRITE(standardMessageUnit,'(2A)') |
100 |
& '(integer) and Interface (half-int.) levels:' |
& 'INI_P_GROUND: PhiRef/g [m] at Center (integer) ', |
101 |
|
& 'and Interface (half-int.) levels:' |
102 |
DO K=1,2*Nr+1 |
DO K=1,2*Nr+1 |
103 |
write(6,'(A,F5.1,A,F10.1,A,F12.3)') ' K=', K*0.5, |
WRITE(standardMessageUnit,'(A,F5.1,A,F10.1,A,F12.3)') |
104 |
& ' ; r=', rHalf(K), ' ; phiRef/g=', phiRef(K) |
& ' K=',K*0.5,' ; r=',rHalf(K),' ; phiRef/g=', phiRef(K) |
105 |
ENDDO |
ENDDO |
106 |
|
_END_MASTER( myThid ) |
107 |
|
|
108 |
C- Find Po_surf : Linear between 2 half levels : |
C- Find Po_surf : Linear between 2 half levels : |
109 |
DO bj = myByLo(myThid), myByHi(myThid) |
DO bj = myByLo(myThid), myByHi(myThid) |