67 |
C rhoK, rhoKM1 - Density at current level, and level above |
C rhoK, rhoKM1 - Density at current level, and level above |
68 |
C phiHyd - Hydrostatic part of the potential phiHydi. |
C phiHyd - Hydrostatic part of the potential phiHydi. |
69 |
C In z coords phiHydiHyd is the hydrostatic |
C In z coords phiHydiHyd is the hydrostatic |
70 |
C pressure anomaly |
C Potential (=pressure/rho0) anomaly |
71 |
C In p coords phiHydiHyd is the geopotential |
C In p coords phiHydiHyd is the geopotential |
72 |
C surface height |
C surface height anomaly. |
|
C anomaly. |
|
73 |
C phiSurfX, - gradient of Surface potentiel (Pressure/rho, ocean) |
C phiSurfX, - gradient of Surface potentiel (Pressure/rho, ocean) |
74 |
C phiSurfY or geopotentiel (atmos) in X and Y direction |
C phiSurfY or geopotentiel (atmos) in X and Y direction |
75 |
C KappaRT, - Total diffusion in vertical for T and S. |
C KappaRT, - Total diffusion in vertical for T and S. |
531 |
jMin = 1-OLy+2 |
jMin = 1-OLy+2 |
532 |
jMax = sNy+OLy-1 |
jMax = sNy+OLy-1 |
533 |
|
|
534 |
C-- Explicit part of the Surface Pressure Gradient (add in TIMESTEP) |
C-- Explicit part of the Surface Potentiel Gradient (add in TIMESTEP) |
535 |
C (note: this loop will be replaced by CALL CALC_GRAD_ETA) |
C (note: this loop will be replaced by CALL CALC_GRAD_ETA) |
536 |
IF (implicSurfPress.NE.1.) THEN |
IF (implicSurfPress.NE.1.) THEN |
537 |
DO j=jMin,jMax |
CALL CALC_GRAD_PHI_SURF( |
538 |
DO i=iMin,iMax |
I bi,bj,iMin,iMax,jMin,jMax, |
539 |
phiSurfX(i,j) = _recip_dxC(i,j,bi,bj)*gBaro |
I etaN, |
540 |
& *(etaN(i,j,bi,bj)-etaN(i-1,j,bi,bj)) |
O phiSurfX,phiSurfY, |
541 |
phiSurfY(i,j) = _recip_dyC(i,j,bi,bj)*gBaro |
I myThid ) |
|
& *(etaN(i,j,bi,bj)-etaN(i,j-1,bi,bj)) |
|
|
ENDDO |
|
|
ENDDO |
|
542 |
ENDIF |
ENDIF |
543 |
|
|
544 |
C-- Start of dynamics loop |
C-- Start of dynamics loop |