Changeset 4002 for trunk/WRF.COMMON

Timestamp:

Jan 14, 2026, 8:29:39 PM (29 hours ago)

Author:

aslmd

Message:

MESOSCALE: fix undefined last layer in module_lmd_driver

File:

• trunk/WRF.COMMON/INTERFACES/module_lmd_driver.F (modified) (3 diffs)

trunk/WRF.COMMON/INTERFACES/module_lmd_driver.F

@@ -419 +419 @@
 ! interm
 ALLOCATE(dz8w_prof(nlayer))
-ALLOCATE(p8w_prof(nlayer))
+ALLOCATE(p8w_prof(nlayer+1))
 ALLOCATE(p_prof(nlayer))
 ALLOCATE(t_prof(nlayer))
@@ -483 +483 @@
 !--------------------------------------!
 dz8w_prof(:) = dz8w(i,kps:kpe,j)   ! dz between full levels (m)   
-p8w_prof(:) = p8w(i,kps:kpe,j)     ! pressure full level (Pa) >> zplev_omp
+p8w_prof(:) = p8w(i,kps:kpe+1,j)     ! pressure full level (Pa) >> zplev_omp
 p_prof(:) = p(i,kps:kpe,j)         ! pressure half level (Pa) >> zplay_omp
 t_prof(:) = t(i,kps:kpe,j)         ! temperature half level (K) >> pt
@@ -542 +542 @@
 ! Dynamic fields for LMD physics !
 !--------------------------------!
-zplev_omp(subs,1:nlayer) = p8w_prof(1:nlayer)  !! NB: last level: no data
+zplev_omp(subs,1:nlayer+1) = p8w_prof(1:nlayer+1)
 zplay_omp(subs,:) = p_prof(:)
 ztfi_omp(subs,:) = t_prof(:)