64 |
C dEvp :: evaporation derivative relative to surf. temp |
C dEvp :: evaporation derivative relative to surf. temp |
65 |
C Slr0 :: upward long wave radiation over freezing surf. |
C Slr0 :: upward long wave radiation over freezing surf. |
66 |
C dSlr :: upward long wave rad. derivative relative to surf. temp |
C dSlr :: upward long wave rad. derivative relative to surf. temp |
67 |
C sFlx :: net surface flux (+=down) function of surf. temp Ts: |
C sFlx :: net heat flux (+=down) except SW, function of surf. temp Ts: |
68 |
C 0: Flux(Ts=0.oC) ; 1: Flux(Ts^n) ; 2: d.Flux/d.Ts(Ts^n) |
C 0: Flux(Ts=0.oC) ; 1: Flux(Ts^n) ; 2: d.Flux/d.Ts(Ts^n) |
69 |
C TSFC :: surface temperature (clim.) (2-dim) |
C TSFC :: surface temperature (clim.) (2-dim) |
70 |
C TSKIN :: skin surface temperature (2-dim) |
C TSKIN :: skin surface temperature (2-dim) |
143 |
C-- Compute net surface heat flux and its derivative ./. surf. temp. |
C-- Compute net surface heat flux and its derivative ./. surf. temp. |
144 |
DO J=1,NGP |
DO J=1,NGP |
145 |
sFlx(J,0)= sFlx(J,0) |
sFlx(J,0)= sFlx(J,0) |
146 |
& - ALHC*Evp0(J) - EMISloc*Slr0(J) + SLRD(J) + SSR(J) |
& - ALHC*Evp0(J) - EMISloc*Slr0(J) + SLRD(J) |
147 |
sFlx(J,1)= sFlx(J,1) |
sFlx(J,1)= sFlx(J,1) |
148 |
& - ALHC*EVAP(J) - EMISloc*SLRU(J) + SLRD(J) + SSR(J) |
& - ALHC*EVAP(J) - EMISloc*SLRU(J) + SLRD(J) |
149 |
sFlx(J,2)= sFlx(J,2) |
sFlx(J,2)= sFlx(J,2) |
150 |
& - ALHC*dEvp(J) - EMISloc*dSlr(J) |
& - ALHC*dEvp(J) - EMISloc*dSlr(J) |
151 |
ENDDO |
ENDDO |
152 |
IF ( aim_energPrecip ) THEN |
IF ( aim_energPrecip ) THEN |
153 |
C- Evap of snow: substract Latent Heat of freezing from heatFlux |
C- Evap of snow/ice: substract Latent Heat of freezing from heatFlux |
154 |
DO J=1,NGP |
DO J=1,NGP |
|
IF ( EnPrec(J) .LT. 0. ) THEN |
|
155 |
sFlx(J,0) = sFlx(J,0) - ALHF*Evp0(J) |
sFlx(J,0) = sFlx(J,0) - ALHF*Evp0(J) |
156 |
sFlx(J,1) = sFlx(J,1) - ALHF*EVAP(J) |
sFlx(J,1) = sFlx(J,1) - ALHF*EVAP(J) |
157 |
sFlx(J,2) = sFlx(J,2) - ALHF*dEvp(J) |
sFlx(J,2) = sFlx(J,2) - ALHF*dEvp(J) |
|
ENDIF |
|
158 |
ENDDO |
ENDDO |
159 |
ENDIF |
ENDIF |
160 |
|
|