Changeset 14908 for NEMO/branches/2021

Timestamp:

2021-05-26T17:43:13+02:00 (3 years ago)

Author:

mathiot

Message:

fix consistency between cav and par cases and sign issues

Location:

NEMO/branches/2021/ticket2669_isf_fluxes/src/OCE/ISF

Files:

• isfpar.F90 (modified) (1 diff)
• isfparmlt.F90 (modified) (2 diffs)

NEMO/branches/2021/ticket2669_isf_fluxes/src/OCE/ISF/isfpar.F90

@@ -88 +88 @@
       !
       ! set temperature content
-      ptsc(:,:,jp_tem) = zqh(:,:) * r1_rho0_rcp
+      ptsc(:,:,jp_tem) = - zqh(:,:) * r1_rho0_rcp
       !
       ! write restart variables (qoceisf, qhcisf, fwfisf for now and before)

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

@@ -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) 
-      pqoce(:,:) =   pqfwf(:,:) * rLfusisf             ! ocean/ice shelf flux assume to be equal to latent heat flux
-      pqhc (:,:) =   pqfwf(:,:) * ztfrz(:,:) * rcp     ! heat content flux 
+      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 )
       !
       CALL iom_put('isftfrz_par', ztfrz )
@@ -157 +157 @@
       !
       ! 2. ------------Net heat flux and fresh water flux due to the ice shelf
-      pqoce(:,:) =   rho0 * rcp * rn_gammat0 * risfLeff(:,:) * e1t(:,:) * ( ztavg(:,:) - ztfrz(:,:) ) * r1_e1e2t(:,:)
-      pqfwf(:,:) = - pqoce(:,:) / rLfusisf             ! derived from the latent heat flux
-      pqhc (:,:) =   pqfwf(:,:) * ztfrz(:,:) * rcp     ! heat content flux 
+      pqfwf(:,:) = - rho0 * rcp * rn_gammat0 * 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 )
       !
       ! 3. ------------BG03 output