Changeset 4000 for trunk/WRF.COMMON/WRFV2

Timestamp:

Jan 14, 2026, 8:29:27 PM (10 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

File:

• trunk/WRF.COMMON/WRFV2/dyn_em/module_initialize_real.F (modified) (2 diffs)

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