Changeset 10900 for branches/UKMO/dev_r5518_FOAMv14_output_heat_fluxes/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90
- Timestamp:
- 2019-04-26T13:51:00+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_FOAMv14_output_heat_fluxes/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90
r10879 r10900 392 392 & + sf(jp_snow)%fnow(:,:,1) * rn_pfac & ! add solid precip heat content at min(Tair,Tsnow) 393 393 & * ( MIN( sf(jp_tair)%fnow(:,:,1), rt0_snow ) - rt0 ) * cpic * tmask(:,:,1) 394 ! Compute and output the individual terms of non solar heat content 395 ! in order to check if qlw, qla, qsb are correctly evaluated 396 qlasprecip(:,:) = sf(jp_snow)%fnow(:,:,1) * rn_pfac * lfus ! latent melting heat for solid precip 397 qevapsst(:,:) = zevap(:,:) * pst(:,:) * rcp ! evap heat content at SST 398 qlprecip(:,:) = ( sf(jp_prec)%fnow(:,:,1) - sf(jp_snow)%fnow(:,:,1) ) * rn_pfac & ! liquid precip heat content at Tair 399 & * ( sf(jp_tair)%fnow(:,:,1) - rt0 ) * rcp 400 qsprecip(:,:) = sf(jp_snow)%fnow(:,:,1) * rn_pfac & ! solid precip heat content at min(Tair,Tsnow) 401 & * ( MIN( sf(jp_tair)%fnow(:,:,1), rt0_snow ) - rt0 ) * cpic * tmask(:,:,1) 394 402 ! 395 403 #if defined key_lim3
Note: See TracChangeset
for help on using the changeset viewer.