Changeset 4151 for trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

Timestamp:

Mar 24, 2026, 5:50:52 PM (13 days ago)

Author:

yluo

Message:

Mars PCM:
Add missing updates to physiq_mod.F.
The previous commit only reverted the changes introduced in r4142 and omitted these updates.
YCL

File:

• trunk/LMDZ.MARS/libf/phymars/physiq_mod.F (modified) (4 diffs)

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -98 +98 @@
       use nonoro_gwd_ran_mod, only: nonoro_gwd_ran
       use nonoro_gwd_mix_mod, only: nonoro_gwd_mix, calljliu_gwimix
+      use tracer_mass_fixer_dyn_mod, only: call_mass_fixer_dyn,
+     &                                   mass_predyn,
+     &                                   tracer_mass_fixer_dyn, 
+     &                                   compute_tracer_mass_global,
+     &                                   set_mass_predyn_from_startfi
       use check_fields_mod, only: check_physics_fields
       use surfini_mod, only: surfini
@@ -584 +589 @@
       REAL :: viscoh                              ! kinematic molecular viscosity for heat
 
+! Variables for correcting tracer mass non-conservation in dynamics
+      real pdq_corrdyn(ngrid, nlayer, nq) ! tracer mass mixing ratio tendencies in order to correct
+                                          ! mass non-conservations in dynamics
+
 c=======================================================================
       pdq(:,:,:) = 0.
@@ -878 +887 @@
       zplay(:,:) = pplay(:,:)
       ps(:) = pplev(:,1)
+
+! Correct mass non-conservation in dynamics
+! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+      if (call_mass_fixer_dyn) then
+         if (firstcall) then
+            call set_mass_predyn_from_startfi
+         endif ! if (firstcall)
+
+         call tracer_mass_fixer_dyn(ngrid, nlayer, nq, ptimestep, 
+     &                              pq, zplev, firstcall, pdq_corrdyn)
+         pdq = pdq + pdq_corrdyn
+      endif ! if (call_mass_fixer_dyn)
 
 !#ifndef MESOSCALE
@@ -2568 +2589 @@
         ENDDO
       ENDDO
+
+      if (call_mass_fixer_dyn) then      
+          call compute_tracer_mass_global(ngrid, nlayer, nq, zq, 
+     &                                    zplev, mass_predyn)
+      end if
 
       ! Density