Changeset 10900 for branches/UKMO/dev_r5518_FOAMv14_output_heat_fluxes/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90

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:

• branches/UKMO/dev_r5518_FOAMv14_output_heat_fluxes/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90 (modified) (2 diffs)

branches/UKMO/dev_r5518_FOAMv14_output_heat_fluxes/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90

@@ -474 +474 @@
          CALL iom_put( "fmmflx", fmmflx  )                      ! Freezing-melting water flux
          CALL iom_put( "qt"    , qns  + qsr )                   ! total heat flux 
-         CALL iom_put( "qns"   , qns        )                   ! solar heat flux
+         CALL iom_put( "qns"   , qns        )                   ! non solar heat flux
          CALL iom_put( "qsr"   ,       qsr  )                   ! solar heat flux
          IF( nn_ice > 0 .OR. nn_components == jp_iam_opa )   CALL iom_put( "ice_cover", fr_i )   ! ice fraction 
@@ -483 +483 @@
             CALL iom_put( "qsb_oce"  , - qsb   )                ! output surface_downward_sensible_heat_flux
             CALL iom_put( "qla_oce"  , - qla   )                ! output surface_downward_latent_heat_flux
+            !Output individual terms to check if their sum matches qns_oce
+            !To be consistent with qns in sbcice_cice.F90, evaluate over ocean only *(1.0-fr_i) 
+            CALL iom_put( "qlasprecip" , - qlasprecip*(1.0-fr_i) )  ! latent heat due to melting of solid precip non solar
+            CALL iom_put( "qevapsst"   , - qevapsst  *(1.0-fr_i) )    ! evap heat content at SST non solar
+            CALL iom_put( "qlprecip"   ,   qlprecip  *(1.0-fr_i) )    ! liquid precip heat content at Tair non solar
+            CALL iom_put( "qsprecip"   ,   qsprecip  *(1.0-fr_i) )    ! solid precip heat content at min(Tair,Tsnow) non solar
          ENDIF
       ENDIF