Changeset 15001 for NEMO/branches/2021/ticket2690_consign_isf/src/OCE/ISF/isfparmlt.F90

Timestamp:

2021-06-16T11:30:56+02:00 (3 years ago)

Author:

mathiot

Message:

ticket #2690: change sign convention of isf fluxes, long name of output and code comments. No changes in WED025 and ISOMIP+ results (output changed because of convention changed, trunk + branch output = 0.0)

File:

• NEMO/branches/2021/ticket2690_consign_isf/src/OCE/ISF/isfparmlt.F90 (modified) (5 diffs)

NEMO/branches/2021/ticket2690_consign_isf/src/OCE/ISF/isfparmlt.F90

@@ -102 +102 @@
       !!--------------------------------------------------------------------
       !
-      ! 0. ------------Read specified runoff
+      ! 0. ------------Read specified fwf from isf to oce
       CALL fld_read ( kt, 1, sf_isfpar_fwf   )
       !
@@ -112 +112 @@
       CALL isf_tbl(Kmm, ztfrz3d, ztfrz, 'T', misfkt_par, rhisf_tbl_par, misfkb_par, rfrac_tbl_par )
       !
-      pqfwf(:,:) = - sf_isfpar_fwf(1)%fnow(:,:,1)      ! fresh water flux from the isf (fwfisf <0 mean melting)       ( > 0 out )
-      pqoce(:,:) = - pqfwf(:,:) * rLfusisf             ! ocean/ice shelf flux assume to be equal to latent heat flux  ( > 0 out )
-      pqhc (:,:) =   pqfwf(:,:) * ztfrz(:,:) * rcp     ! heat content flux                                            ( > 0 out )
+      pqfwf(:,:) =   sf_isfpar_fwf(1)%fnow(:,:,1)      ! fresh water flux from the isf (fwfisf <0 mean melting)       ( > 0 from isf to oce)
+      pqoce(:,:) = - pqfwf(:,:) * rLfusisf             ! ocean/ice shelf flux assume to be equal to latent heat flux  ( > 0 from isf to oce)
+      pqhc (:,:) =   pqfwf(:,:) * ztfrz(:,:) * rcp     ! heat content flux                                            ( > 0 from isf to oce)
       !
       CALL iom_put('isftfrz_par', ztfrz(:,:) * mskisf_par(:,:) )
@@ -157 +157 @@
       !
       ! 2. ------------Net heat flux and fresh water flux due to the ice shelf
-      pqfwf(:,:) = - rho0 * rcp * rn_isfpar_bg03_gt0 * risfLeff(:,:) * e1t(:,:) * (ztavg(:,:) - ztfrz(:,:) ) * r1_e1e2t(:,:) / rLfusisf  ! ( > 0 out )
-      pqoce(:,:) = - pqfwf(:,:) * rLfusisf             ! ocean/ice shelf flux assume to be equal to latent heat flux  ( > 0 out )
-      pqhc (:,:) =   pqfwf(:,:) * ztfrz(:,:) * rcp     ! heat content flux                                            ( > 0 out )
+      pqfwf(:,:) =   rho0 * rcp * rn_isfpar_bg03_gt0 * risfLeff(:,:) * e1t(:,:) * (ztavg(:,:) - ztfrz(:,:) ) * r1_e1e2t(:,:) / rLfusisf  ! ( > 0 from isf to oce)
+      pqoce(:,:) = - pqfwf(:,:) * rLfusisf             ! ocean/ice shelf flux assume to be equal to latent heat flux  ( > 0 from isf to oce)
+      pqhc (:,:) =   pqfwf(:,:) * ztfrz(:,:) * rcp     ! heat content flux                                            ( > 0 from isf to oce)
       !
       ! 3. ------------BG03 output
@@ -207 +207 @@
       !
       ! 2. ------------Scale isf melt pattern with total amount from oasis
-      ! ice shelf 2d map
-      zfwf(:,:) = - sf_isfpar_fwf(1)%fnow(:,:,1)
+      ! ice shelf 2d map of fwf from isf to oce
+      zfwf(:,:) = sf_isfpar_fwf(1)%fnow(:,:,1)
       !
       ! compute glob sum from input file
@@ -223 +223 @@
       ! 3. -----------Define fwf and qoce
       ! ocean heat flux is assume to be equal to the latent heat
-      pqfwf(:,:) =   zfwf(:,:)                         ! fwf                ( >0 out )
-      pqoce(:,:) = - pqfwf(:,:) * rLfusisf             ! ocean heat flux    ( >0 out ) (assumed to be the latent heat flux)
-      pqhc (:,:) =   pqfwf(:,:) * ztfrz(:,:) * rcp     ! heat content flux  ( >0 out )
+      pqfwf(:,:) =   zfwf(:,:)                         ! fwf                ( > 0 from isf to oce)
+      pqoce(:,:) = - pqfwf(:,:) * rLfusisf             ! ocean heat flux    ( > 0 from isf to oce) (assumed to be the latent heat flux)
+      pqhc (:,:) =   pqfwf(:,:) * ztfrz(:,:) * rcp     ! heat content flux  ( > 0 from isf to oce)
       !
       CALL iom_put('isftfrz_par', ztfrz )