Changeset 14917 for NEMO/branches/2021

Timestamp:

2021-05-27T19:30:09+02:00 (3 years ago)

Author:

mathiot

Message:

ticket #2669: add masking of isftfrz_* variables, fix issue on vertical mean in bg03 case, fix sign of qlatisf output

Location:

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

Files:

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

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

@@ -124 +124 @@
       !
       ! output freezing point at the interface
-      CALL iom_put('isftfrz_cav', ztfrz )
+      CALL iom_put('isftfrz_cav', ztfrz(:,:) * mskisf_cav(:,:) )
       !
    END SUBROUTINE isfcav_mlt_spe
@@ -169 +169 @@
       ! output thermal driving and freezinpoint at the ice shelf interface
       CALL iom_put('isfthermald_cav', zthd )
-      CALL iom_put('isftfrz_cav'    , ztfrz )
+      CALL iom_put('isftfrz_cav'    , ztfrz(:,:) * mskisf_cav(:,:) )
       !
    END SUBROUTINE isfcav_mlt_2eq

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

@@ -75 +75 @@
       zqhc (:,:) = zqhc(:,:)  * mskisf_par(:,:)
       !
-      ! compute heat content flux ( > 0 out )
-      zqlat(:,:) = pqfwf(:,:) * rLfusisf    ! 2d latent heat flux (W/m2)
+      ! compute latent heat flux ( > 0 out )
+      zqlat(:,:) = - pqfwf(:,:) * rLfusisf    ! 2d latent heat flux (W/m2)
       !
       ! total heat flux ( > 0 out )

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

@@ -116 +116 @@
       pqhc (:,:) =   pqfwf(:,:) * ztfrz(:,:) * rcp     ! heat content flux                                            ( > 0 out )
       !
-      CALL iom_put('isftfrz_par', ztfrz )
+      CALL iom_put('isftfrz_par', ztfrz(:,:) * mskisf_par(:,:) )
       !
    END SUBROUTINE isfpar_mlt_spe
@@ -154 +154 @@
       !
       ! 1. ------------Mean temperature
-      CALL isf_tbl(Kmm, ts(:,:,jk,jp_tem,Kmm), ztavg, 'T', misfkt_par, rhisf_tbl_par, misfkb_par, rfrac_tbl_par )
+      CALL isf_tbl(Kmm, ts(:,:,:,jp_tem,Kmm), ztavg, 'T', misfkt_par, rhisf_tbl_par, misfkb_par, rfrac_tbl_par )
       !
       ! 2. ------------Net heat flux and fresh water flux due to the ice shelf
@@ -169 +169 @@
       !
       ! output freezing point used to define the thermal driving and heat content fluxes
-      CALL iom_put('isftfrz_par', ztfrz )
+      CALL iom_put('isftfrz_par', ztfrz(:,:) * mskisf_par(:,:) )
       !
    END SUBROUTINE isfpar_mlt_bg03