Changeset 5158 for LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad.v1.5.1/radiation_overlap.F90

Timestamp:

Aug 2, 2024, 2:12:03 PM (3 months ago)

Author:

abarral

Message:

Add missing klon on strataer_emiss_mod.F90
Correct various missing explicit declarations
Replace tabs by spaces (tabs are not part of the fortran charset)
Continue cleaning modules
Removed unused arguments and variables

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad.v1.5.1/radiation_overlap.F90 (modified) (6 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad.v1.5.1/radiation_overlap.F90

@@ -94 +94 @@
     ! multiply by "factor"
     denominator = 1.0_jprb
-    do jreg = 1,nreg
+    DO jreg = 1,nreg
       op_x_frac_min(jreg) = op(jreg) &
            &  * min(frac_upper(jreg), frac_lower(jreg))
@@ -103 +103 @@
       factor = 1.0_jprb / denominator
       ! Create the random part of the overlap matrix
-      do jupper = 1,nreg
-        do jlower = 1,nreg
+      DO jupper = 1,nreg
+        DO jlower = 1,nreg
           overlap_matrix(jupper,jlower) = factor &
                &  * (frac_lower(jlower)-op_x_frac_min(jlower)) &
@@ -115 +115 @@
     
     ! Add on the maximum part of the overlap matrix
-    do jreg = 1,nreg
+    DO jreg = 1,nreg
       overlap_matrix(jreg,jreg) = overlap_matrix(jreg,jreg) &
            &  + op_x_frac_min(jreg)
@@ -368 +368 @@
 
     ! Loop through each atmospheric column
-    do jcol = istartcol, iendcol
+    DO jcol = istartcol, iendcol
       ! For this column, outer space is treated as one clear-sky
       ! region, so the fractions are assigned as such
@@ -381 +381 @@
       ! half-level starting at 1 for the top-of-atmosphere, as well
       ! as indexing each level starting at 1 for the top-most level.
-      do jlev = 1,nlev+1
+      DO jlev = 1,nlev+1
         ! Fraction of each region just below the interface
         if (jlev > nlev) then
@@ -426 +426 @@
 
         ! Convert to directional overlap matrices
-        do jupper = 1,nreg
-          do jlower = 1,nreg
+        DO jupper = 1,nreg
+          DO jlower = 1,nreg
             if (frac_lower(jlower) >= frac_threshold) then
               u_matrix(jupper,jlower,jlev,jcol) = overlap_matrix(jupper,jlower) &