Changeset 5185 for LMDZ6/branches/Amaury_dev/libf/phylmd/ecrad.v1.5.1/radiation_cloud_generator.F90

Timestamp:

Sep 11, 2024, 4:27:07 PM (5 months ago)

Author:

abarral

Message:

Replace REPROBUS CPP KEY by logical using handmade wonky wrapper

File:

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

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

@@ -213 +213 @@
           trigger = rand_top(jg) * total_cloud_cover
           jlev = ibegin
-          DO while (trigger > cum_cloud_cover(jlev) .and. jlev < iend)
+          DO while (trigger > cum_cloud_cover(jlev) .AND. jlev < iend)
             jlev = jlev + 1
           end do
@@ -689 +689 @@
 
           ! For each spectral interval, has the first cloud appeared at this level?
-          first_cloud(jg) = (trigger(jg) <= cum_cloud_cover(jlev) .and. .not. found_cloud(jg))
+          first_cloud(jg) = (trigger(jg) <= cum_cloud_cover(jlev) .AND. .not. found_cloud(jg))
 
           ! ...if so, add to found_cloud
@@ -699 +699 @@
           ! prev_cloud)
           is_cloud(jg) = first_cloud(jg) &
-               &  .or. found_cloud(jg) .and. merge(rand_cloud(jg,jlev)*frac(jlev-1) &
+               &  .or. found_cloud(jg) .AND. merge(rand_cloud(jg,jlev)*frac(jlev-1) &
                &               < frac(jlev)+frac(jlev-1)-pair_cloud_cover(jlev-1), &
                &             rand_cloud(jg,jlev)*(cum_cloud_cover(jlev-1) - frac(jlev-1)) &
@@ -712 +712 @@
           rand_inhom(jg,jlev) = merge(merge(rand_inhom(jg,jlev-1), rand_inhom(jg,jlev), &
                &                           rand_inhom2(jg,jlev) < overlap_param_inhom(jlev-1) &
-               &                           .and. prev_cloud(jg)), &
+               &                           .AND. prev_cloud(jg)), &
                &                     0.0_jprb, is_cloud(jg))
         end do