Raw mass flux multiplied by 1-SPHU to generate dry mass flux

geoschem · Oct 1, 2024 · eb85209 · eb85209
1 parent e499969
commit eb85209
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/GCHP_GridComp/GCHPctmEnv_GridComp/GCHPctmEnv_GridCompMod.F90 b/src/GCHP_GridComp/GCHPctmEnv_GridComp/GCHPctmEnv_GridCompMod.F90
@@ -871,15 +871,15 @@ 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
 
          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