9 |
SUBROUTINE SUFLUX_OCEAN( |
SUBROUTINE SUFLUX_OCEAN( |
10 |
I PSA, FMASK, |
I PSA, FMASK, |
11 |
I Tsurf, SSR, SLRD, |
I Tsurf, SSR, SLRD, |
12 |
I T0, Q0, CDENVV, |
I T1, T0, Q0, DENVV, |
13 |
O SHF, EVAP, SLRU, |
O SHF, EVAP, SLRU, |
14 |
I bi,bj,myThid) |
I bi,bj,myThid) |
15 |
|
|
43 |
C Tsurf :: surface temperature (2-dim) |
C Tsurf :: surface temperature (2-dim) |
44 |
C SSR :: sfc sw radiation (net flux) (2-dim) |
C SSR :: sfc sw radiation (net flux) (2-dim) |
45 |
C SLRD :: sfc lw radiation (downward flux)(2-dim) |
C SLRD :: sfc lw radiation (downward flux)(2-dim) |
46 |
|
C T1 :: near-surface air temperature (from Pot.temp) |
47 |
C T0 :: near-surface air temperature (2-dim) |
C T0 :: near-surface air temperature (2-dim) |
48 |
C Q0 :: near-surface sp. humidity [g/kg](2-dim) |
C Q0 :: near-surface sp. humidity [g/kg](2-dim) |
49 |
C CDENVV :: sensible heat flux coefficient (2-dim) |
C DENVV :: surface flux (sens,lat.) coeff. (=Rho*|V|) [kg/m2/s] |
50 |
C-- Output: |
C-- Output: |
51 |
C SHF :: sensible heat flux (2-dim) |
C SHF :: sensible heat flux (2-dim) |
52 |
C EVAP :: evaporation [g/(m^2 s)] (2-dim) |
C EVAP :: evaporation [g/(m^2 s)] (2-dim) |
57 |
C-- |
C-- |
58 |
_RL PSA(NGP) |
_RL PSA(NGP) |
59 |
_RL FMASK(NGP), Tsurf(NGP) |
_RL FMASK(NGP), Tsurf(NGP) |
60 |
_RL SSR(NGP), SLRD(NGP), T0(NGP), Q0(NGP), CDENVV(NGP) |
_RL SSR(NGP), SLRD(NGP) |
61 |
|
_RL T1(NGP), T0(NGP), Q0(NGP), DENVV(NGP) |
62 |
|
|
63 |
_RL SHF(NGP), EVAP(NGP), SLRU(NGP) |
_RL SHF(NGP), EVAP(NGP), SLRU(NGP) |
64 |
|
|
68 |
#ifdef ALLOW_AIM |
#ifdef ALLOW_AIM |
69 |
|
|
70 |
C-- Local variables: |
C-- Local variables: |
71 |
|
C CDENVV :: surf. heat flux (sens.,lat.) coeff including stability effect |
72 |
|
_RL CDENVV(NGP), RDTH, FSSEA |
73 |
_RL QSAT0(NGP) |
_RL QSAT0(NGP) |
74 |
_RL QDUMMY(1), RDUMMY(1), TS4 |
_RL QDUMMY(1), RDUMMY(1), TS4 |
75 |
INTEGER J |
INTEGER J |
82 |
|
|
83 |
C 2.2 Sensible heat flux (from clim. TS over land) |
C 2.2 Sensible heat flux (from clim. TS over land) |
84 |
|
|
85 |
|
C Stability correction |
86 |
|
|
87 |
|
RDTH = FSTAB/DTHETA |
88 |
|
|
89 |
|
DO J=1,NGP |
90 |
|
FSSEA =1.+MIN(DTHETA,MAX(-DTHETA, Tsurf(J)-T1(J)))*RDTH |
91 |
|
CDENVV(J)=CHS*DENVV(J)*FSSEA |
92 |
|
ENDDO |
93 |
|
|
94 |
DO J=1,NGP |
DO J=1,NGP |
95 |
SHF(J) = CDENVV(J)*CP*(Tsurf(J)-T0(J)) |
SHF(J) = CDENVV(J)*CP*(Tsurf(J)-T0(J)) |
96 |
ENDDO |
ENDDO |