From eb8520994361dfdc765e494fd2f14dff13a7882f Mon Sep 17 00:00:00 2001 From: "yuanjian.z" Date: Mon, 30 Sep 2024 21:48:52 -0500 Subject: [PATCH] Raw mass flux multiplied by 1-SPHU to generate dry mass flux --- .../GCHPctmEnv_GridComp/GCHPctmEnv_GridCompMod.F90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/GCHP_GridComp/GCHPctmEnv_GridComp/GCHPctmEnv_GridCompMod.F90 b/src/GCHP_GridComp/GCHPctmEnv_GridComp/GCHPctmEnv_GridCompMod.F90 index a6b47ff27..d7f424d59 100644 --- a/src/GCHP_GridComp/GCHPctmEnv_GridComp/GCHPctmEnv_GridCompMod.F90 +++ b/src/GCHP_GridComp/GCHPctmEnv_GridComp/GCHPctmEnv_GridCompMod.F90 @@ -871,7 +871,7 @@ subroutine prepare_massflux_exports(IMPORT, EXPORT, PLE, dt, RC) call MAPL_GetPointer(IMPORT, temp3d_r4, 'MFXC', RC=STATUS) _VERIFY(STATUS) if ( correct_mass_flux_for_humidity > 0 ) then - MFX_EXPORT = dble(temp3d_r4) / ( 1.d0 - SPHU0_EXPORT ) + MFX_EXPORT = dble(temp3d_r4) * ( 1.d0 - SPHU0_EXPORT ) else MFX_EXPORT = dble(temp3d_r4) endif @@ -879,7 +879,7 @@ subroutine prepare_massflux_exports(IMPORT, EXPORT, PLE, dt, RC) call MAPL_GetPointer(IMPORT, temp3d_r4, 'MFYC', RC=STATUS) _VERIFY(STATUS) if ( correct_mass_flux_for_humidity > 0 ) then - MFY_EXPORT = dble(temp3d_r4) / ( 1.d0 - SPHU0_EXPORT ) + MFY_EXPORT = dble(temp3d_r4) * ( 1.d0 - SPHU0_EXPORT ) else MFY_EXPORT = dble(temp3d_r4) endif