Changeset 3999

Timestamp:

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

Author:

aslmd

Message:

MESOSCALE: workaround to a problem in initialization of vertical levels arising in debug mode

File:

• trunk/WRF.COMMON/WRFV2/dyn_em/module_initialize_real.F (modified) (1 diff)

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

@@ -5320 +5320 @@
                DO k = fixedpbl, kte-1
                   pb = (znw(k)+znw(k+1))*0.5 * (p_surf - p_top) + p_top
+!! sometimes pb could be negative in the iteration process
+!! because one of znw(k) znw(k+1) is negative
+!! this triggers error in the LOG below in debug mode
+!! so if the case arises, we set pb to p_top
+!! in the end with that patch convergence looks OK
+IF (pb<0.) THEN
+        print*,"AAAA",pb,znw(k), znw(k+1), loop, loop1
+        pb = p_top
+ENDIF
 !!                 temp = MAX ( 200., t00 + A*LOG(pb/p00) )
 !                  temp =             t00 + A*LOG(pb/p00)