Changeset 4150 for trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90

Timestamp:

Mar 24, 2026, 5:27:57 PM (2 weeks ago)

Author:

yluo

Message:

Mars PCM:
Refactored the tracer mass-conservation fixer for dynamics introduced in r4142:
Added a runtime flag to enable/disable the correction of tracer mass non-conservation in dynamics.
Reorganized the implementation into a modular structure.
Restricted the correction to a subset of tracers: CO, O2, H2, HO2, H2O2, N2, Ar, and He.
Enabled correction for dynamics timesteps preceding the first physics timestep of a simulation.
YCL

File:

• trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90

@@ -26 +26 @@
 use nonoro_gwd_mix_mod,  only: du_eddymix_gwd, dv_eddymix_gwd, de_eddymix_rto, &
                                df_eddymix_flx, dh_eddymix_gwd, dq_eddymix_gwd
+use tracer_mass_fixer_dyn_mod, only: mass_predyn_co, mass_predyn_o2, mass_predyn_h2, mass_predyn_ho2,  &
+                                     mass_predyn_h2o2, mass_predyn_n2, mass_predyn_ar, mass_predyn_he, &
+                                     found_startfi_co, found_startfi_o2, found_startfi_h2, found_startfi_ho2, &
+                                     found_startfi_h2o2, found_startfi_n2, found_startfi_ar, found_startfi_he, &
+                                     call_mass_fixer_dyn
 use compute_dtau_mod,    only: dtau
 use dust_rad_adjust_mod, only: dust_rad_adjust_prev,dust_rad_adjust_next
@@ -718 +723 @@
              minval(df_eddymix_flx ), maxval(df_eddymix_flx )
 
+! Correction to tracer mass non-conservation in dynamics
+if (startphy_file .and. call_mass_fixer_dyn) then
+   call get_var("mass_predyn_co", mass_predyn_co, found)
+   if (found) then
+      found_startfi_co = .true.
+   else
+      write(*,*) "phyetat0: <mass_predyn_co> not in file"
+   endif
+
+   call get_var("mass_predyn_o2", mass_predyn_o2, found)
+   if (found) then
+      found_startfi_o2 = .true.
+   else
+      write(*,*) "phyetat0: <mass_predyn_o2> not in file"
+   endif
+
+   call get_var("mass_predyn_h2", mass_predyn_h2, found)
+   if (found) then
+      found_startfi_h2 = .true.
+   else
+      write(*,*) "phyetat0: <mass_predyn_h2> not in file"
+   endif
+
+   call get_var("mass_predyn_ho2", mass_predyn_ho2, found)
+   if (found) then
+      found_startfi_ho2 = .true.
+   else
+      write(*,*) "phyetat0: <mass_predyn_ho2> not in file"
+   endif
+
+   call get_var("mass_predyn_h2o2", mass_predyn_h2o2, found)
+   if (found) then
+      found_startfi_h2o2 = .true.
+   else
+      write(*,*) "phyetat0: <mass_predyn_h2o2> not in file"
+   endif
+
+   call get_var("mass_predyn_n2", mass_predyn_n2, found)
+   if (found) then
+      found_startfi_n2 = .true.
+   else
+      write(*,*) "phyetat0: <mass_predyn_n2> not in file"
+   endif
+
+   call get_var("mass_predyn_ar", mass_predyn_ar, found)
+   if (found) then
+      found_startfi_ar = .true.
+   else
+      write(*,*) "phyetat0: <mass_predyn_ar> not in file"
+   endif
+
+   call get_var("mass_predyn_he", mass_predyn_he, found)
+   if (found) then
+      found_startfi_he = .true.
+   else
+      write(*,*) "phyetat0: <mass_predyn_he> not in file"
+   endif
+endif ! if (startphy_file .and. call_mass_fixer_dyn)
+
 ! tracer on surface
 if (nq.ge.1) then