Changeset 4000

Timestamp:

Jan 14, 2026, 8:29:27 PM (3 days ago)

Author:

aslmd

Message:

MESOSCALE: in module_initialize_real, replace 0 by 1E-30 in the enforcing positivity.

[kept commented for reference] MESOSCALE: sometimes tracers are coming slightly negative from the physics, enforce positivity. this seems to solve problems occurring in the physics which makes the model to crash in normal mode and to raise an exception in lwtt in debug mode

Location:

trunk/WRF.COMMON

Files:

• INTERFACES/module_lmd_driver.F (modified) (1 diff)
• WRFV2/dyn_em/module_initialize_real.F (modified) (2 diffs)

trunk/WRF.COMMON/INTERFACES/module_lmd_driver.F

@@ -816 +816 @@
 ENDDO
 ENDDO
+
+!!!!! KEPT COMMENTED FOR TESTING PURPOSES
+!!! sometimes tracers are coming slightly negative from the physics
+!!! enforce positivity here (rough patch admittedly)
+!where (scalar < 1.e-30)
+!        scalar = 1.e-30
+!endwhere
+
 CALL deallocate_interface
 DEALLOCATE(pdtheta)

trunk/WRF.COMMON/WRFV2/dyn_em/module_initialize_real.F

@@ -1699 +1699 @@
 !!! we want any scalar (i.e. tracer) to be positive
 !!! and because of interpolation it is possible that negative values occur...
-WHERE( scalar < 0. )  scalar = 0.
+WHERE( scalar < 1.e-30 )  scalar = 1.e-30
 
 !!!!!!****MARS****!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -6018 +6018 @@
       real p_pgcm(dimlevs)    ! high res to GCM pressure ratios
 
+
 ! 1. Coefficients
       rp=ps_hr/ps_gcm