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

Timestamp:

Mar 17, 2026, 4:22:44 PM (3 weeks ago)

Author:

jbclement

Message:

PEM:

• Move ice table variables from "ice_table" to the main program.
• Merge "job_id_mod" and "job_timelimit_mod" into "job" which is relocated to the PEM folder.
• Rename local variables in procedures to avoid masking variables in parent scope.
• Few cleanings to delete remaining PEM-external "include" and "use".

JBC

File:

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

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

@@ -408 +408 @@
 
 !=======================================================================
-SUBROUTINE compute_alt_coord(pa,preff,ap,bp,aps,bps)
+SUBROUTINE compute_alt_coord(pa,preff,ap_l,bp_l,aps,bps)
 !-----------------------------------------------------------------------
 ! NAME
@@ -449 +449 @@
 ! ---------
 real(dp)                                               :: pa, preff
-real(dp), dimension(nlayer + 1), intent(out)           :: ap, bp
+real(dp), dimension(nlayer + 1), intent(out)           :: ap_l, bp_l
 real(dp), dimension(nlayer),     intent(out), optional :: aps, bps
 
-! LOCAL VAIRABLES
+! LOCAL VARIABLES
 ! ---------------
 integer(di)                     :: i
@@ -471 +471 @@
         call compute_hybrid_sig(sig(i),pa,preff,newsig)
         bp(i) = exp(1._dp - 1._dp/(newsig**2))
-        ap(i) = pa*(newsig - bp(i))
+        ap_l(i) = pa*(newsig - bp(i))
     end do
-    ap(nlayer + 1) = 0._dp
+    ap_l(nlayer + 1) = 0._dp
     bp(nlayer + 1) = 0._dp
 else
     call print_msg('> Defining sigma altitude coordinates',LVL_NFO)
-    ap(:) = 0._dp
+    ap_l(:) = 0._dp
     bp(1:nlayer) = sig(1:nlayer)
     bp(nlayer + 1) = 0._dp
@@ -487 +487 @@
 
 do i = 1,nlayer - 1
-    aps(i) =  0.5_dp*(ap(i) + ap(i + 1))
+    aps(i) =  0.5_dp*(ap_l(i) + ap_l(i + 1))
     bps(i) =  0.5_dp*(bp(i) + bp(i + 1))
 end do