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

Timestamp:

Mar 17, 2026, 4:22:44 PM (12 days 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/maths.F90 (modified) (2 diffs)

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

@@ -301 +301 @@
 integer(di)            :: i
 real(dp)               :: m
-real(dp), dimension(n) :: cp, dp
+real(dp), dimension(n) :: c_p, d_p
 
 ! CODE
@@ -323 +323 @@
 
 ! Initialization
-cp(1) = c(1)/b(1)
-dp(1) = d(1)/b(1)
+c_p(1) = c(1)/b(1)
+d_p(1) = d(1)/b(1)
 
 ! Forward sweep
 do i = 2,n - 1
-    m = b(i) - a(i - 1)*cp(i - 1)
-    cp(i) = c(i)/m
-    dp(i) = (d(i) - a(i - 1)*dp(i - 1))/m
-end do
-m = b(n) - a(n - 1)*cp(n - 1)
-dp(n) = (d(n) - a(n - 1)*dp(n - 1))/m
+    m = b(i) - a(i - 1)*c_p(i - 1)
+    c_p(i) = c(i)/m
+    d_p(i) = (d(i) - a(i - 1)*d_p(i - 1))/m
+end do
+m = b(n) - a(n - 1)*c_p(n - 1)
+d_p(n) = (d(n) - a(n - 1)*d_p(n - 1))/m
 
 ! Backward substitution
-x(n) = dp(n)
+x(n) = d_p(n)
 do i = n - 1,1,-1
-    x(i) = dp(i) - cp(i)*x(i + 1)
+    x(i) = d_p(i) - c_p(i)*x(i + 1)
 end do