46 |
_RL sumArea |
_RL sumArea |
47 |
_RL aC, aCw, aCs |
_RL aC, aCw, aCs |
48 |
|
|
|
C-- Initialise -Boyancy at surface level : Bo_surf |
|
|
C Bo_surf is defined as d/dr(Phi_surf) and set to g/rtoz (linear free surface) |
|
|
C with rtoz = conversion factor from r-unit to z-unit (=horiVertRatio) |
|
|
C an acurate formulation will include P_surf and T,S_surf effects on rho_surf: |
|
|
C z-ocean (rtoz=1) : Bo_surf = - Boyancy = gravity * rho_surf/rho_0 |
|
|
C p-atmos (rtoz=rho_c*g) : Bo_surf = (1/rho)_surf |
|
|
C-- |
|
|
IF ( buoyancyRelation .EQ. 'ATMOSPHERIC' ) THEN |
|
|
DO bj=myByLo(myThid),myByHi(myThid) |
|
|
DO bi=myBxLo(myThid),myBxHi(myThid) |
|
|
DO J=1-Oly,sNy+Oly |
|
|
DO I=1-Olx,sNx+Olx |
|
|
Bo_surf(I,J,bi,bj) = recip_rhoConst |
|
|
recip_Bo(I,J,bi,bj) = rhoConst |
|
|
ENDDO |
|
|
ENDDO |
|
|
ENDDO |
|
|
ENDDO |
|
|
ELSE |
|
|
C- gBaro = gravity (except for External mode test with reduced gravity) |
|
|
DO bj=myByLo(myThid),myByHi(myThid) |
|
|
DO bi=myBxLo(myThid),myBxHi(myThid) |
|
|
DO J=1-Oly,sNy+Oly |
|
|
DO I=1-Olx,sNx+Olx |
|
|
Bo_surf(I,J,bi,bj) = gBaro |
|
|
recip_Bo(I,J,bi,bj) = 1. _d 0 / gBaro |
|
|
ENDDO |
|
|
ENDDO |
|
|
ENDDO |
|
|
ENDDO |
|
|
ENDIF |
|
|
|
|
|
C-- Update overlap regions |
|
|
_EXCH_XY_R8(Bo_surf, myThid) |
|
|
_EXCH_XY_R8(recip_Bo, myThid) |
|
|
|
|
49 |
C-- Initialise laplace operator |
C-- Initialise laplace operator |
50 |
C aW2d: integral in Z Ax/dX |
C aW2d: integral in Z Ax/dX |
51 |
C aS2d: integral in Z Ay/dY |
C aS2d: integral in Z Ay/dY |