Changeset 10879 for branches/UKMO/dev_r5518_FOAMv14_output_heat_fluxes/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90

Timestamp:

2019-04-16T14:58:06+02:00 (5 years ago)

Author:

anaguiar

Message:

Changes to reinstate heat fluxes, see ​https://code.metoffice.gov.uk/trac/gmed/ticket/454

File:

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

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

@@ -668 +668 @@
 ! Solar penetrative radiation and non solar surface heat flux
 
-! Scale qsr and qns according to ice fraction (bulk formulae only)
+! Scale qsr, qns, qlw, qsb, qla according to ice fraction (bulk formulae only)
 
       IF (ksbc == jp_core) THEN
          qsr(:,:)=qsr(:,:)*(1.0-fr_i(:,:))
          qns(:,:)=qns(:,:)*(1.0-fr_i(:,:))
+         qlw(:,:)=qlw(:,:)*(1.0-fr_i(:,:))
+         qsb(:,:)=qsb(:,:)*(1.0-fr_i(:,:))
+         qla(:,:)=qla(:,:)*(1.0-fr_i(:,:))
       ENDIF
 ! Take into account snow melting except for fully coupled when already in qns_tot
@@ -680 +683 @@
       ELSE
          qns(:,:)= qns(:,:)-sprecip(:,:)*Lfresh*(1.0-fr_i(:,:))
+         qla(:,:)= qla(:,:)+sprecip(:,:)*Lfresh*(1.0-fr_i(:,:))
       ENDIF