Changeset 5158 for LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad.v1.5.1/radiation_mcica_lw.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_mcica_lw.F90 (modified) (7 diffs)

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

@@ -150 +150 @@
 
     ! Loop through columns
-    do jcol = istartcol,iendcol
+    DO jcol = istartcol,iendcol
 
       ! Clear-sky calculation
@@ -156 +156 @@
         ! Scattering case: first compute clear-sky reflectance,
         ! transmittance etc at each model level
-        do jlev = 1,nlev
+        DO jlev = 1,nlev
           call calc_two_stream_gammas_lw(ng, ssa(:,jlev,jcol), g(:,jlev,jcol), &
                &  gamma1, gamma2)
@@ -174 +174 @@
         ! Non-scattering case: use simpler functions for
         ! transmission and emission
-        do jlev = 1,nlev
+        DO jlev = 1,nlev
           call calc_no_scattering_transmittance_lw(ng, od(:,jlev,jcol), &
                &  planck_hl(:,jlev,jcol), planck_hl(:,jlev+1, jcol), &
@@ -215 +215 @@
         is_clear_sky_layer = .true.
         i_cloud_top = nlev+1
-        do jlev = 1,nlev
+        DO jlev = 1,nlev
           ! Compute combined gas+aerosol+cloud optical properties
           if (cloud%fraction(jcol,jlev) >= config%cloud_fraction_threshold) then
@@ -224 +224 @@
             end if
 
-            do jg = 1,ng
+            DO jg = 1,ng
               od_cloud_new(jg) = od_scaling(jg,jlev) &
                  &  * od_cloud(config%i_band_from_reordered_g_lw(jg),jlev,jcol)
@@ -240 +240 @@
                 ! case that od_total > 0.0 and ssa_total > 0.0 but
                 ! od_total*ssa_total == 0 due to underflow
-                do jg = 1,ng
+                DO jg = 1,ng
                   if (od_total(jg) > 0.0_jprb) then
                     scat_od_total(jg) = ssa(jg,jlev,jcol)*od(jg,jlev,jcol) &
@@ -259 +259 @@
               else
 
-                do jg = 1,ng
+                DO jg = 1,ng
                   if (od_total(jg) > 0.0_jprb) then
                     scat_od = ssa_cloud(config%i_band_from_reordered_g_lw(jg),jlev,jcol) &