New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 10900 for branches/UKMO/dev_r5518_FOAMv14_output_heat_fluxes/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90 – NEMO

Ignore:
Timestamp:
2019-04-26T13:51:00+02:00 (5 years ago)
Author:
anaguiar
Message:

Changes to allow checking qns balance and if all terms are correctly evaluated - Test 4 described in https://code.metoffice.gov.uk/trac/gmed/ticket/454#comment:2

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  
    392392         &     + sf(jp_snow)%fnow(:,:,1) * rn_pfac                                &   ! add solid  precip heat content at min(Tair,Tsnow) 
    393393         &     * ( 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) 
    394402      ! 
    395403#if defined key_lim3 
Note: See TracChangeset for help on using the changeset viewer.