Changeset 4074 for trunk/LMDZ.COMMON/libf/evolution/atmosphere.F90

Timestamp:

Feb 17, 2026, 2:45:53 PM (11 days ago)

Author:

jbclement

Message:

PEM:

• Correct management of H2O ice tendency in 1D when there is not enough ice anymore.
• Clean initialization of allocatable module arrays (especially needed when no slope)
• One more renaming for consistency + few small updates thoughout the code.

JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/atmosphere.F90 (modified) (2 diffs)

trunk/LMDZ.COMMON/libf/evolution/atmosphere.F90

@@ -69 +69 @@
 if (.not. allocated(u_PCM)) allocate(u_PCM(ngrid,nlayer))
 if (.not. allocated(v_PCM)) allocate(v_PCM(ngrid,nlayer))
+ap(:) = 0._dp
+bp(:) = 0._dp
+ps_PCM(:) = 0._dp
+teta_PCM(:,:) = 0._dp
+u_PCM(:,:) = 0._dp
+v_PCM(:,:) = 0._dp
 
 END SUBROUTINE ini_atmosphere
@@ -662 +668 @@
 ! The pressure deviation is rescaled to avoid disproportionate oscillations in case of huge change of average pressure during the PEM run
 ps4PCM(:) = ps_avg(:) + ps_dev(:)*ps_avg_glob/ps_avg_glob_ini
-pa4PCM = ps_avg_glob_ini/30.  ! For now the altitude grid is not changed
-preff4PCM = ps_avg_glob_ini   ! For now the altitude grid is not changed
-
+pa4PCM = ps_avg_glob_ini/30. ! For now the altitude grid is not changed
+preff4PCM = ps_avg_glob_ini  ! For now the altitude grid is not changed
+
+! Correction on teta due to surface pressure changes
 call print_msg('> Building potential temperature for the PCM')
 do l = 1,nlayer
-    ! Correction on teta due to surface pressure changes
     teta4PCM(:,l) = teta_PCM(:,l)*ps4PCM(:)**rcp
 end do
 
+! Compute atmospheric pressure
 call print_msg('> Building air mass for the PCM')
-! Compute atmospheric pressure
 do l = 1,nlayer + 1
     p(:,l) = ap(l) + bp(l)*ps4PCM(:)