Changeset 5095 for LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/parasol.F90

Timestamp:

Jul 22, 2024, 9:46:57 AM (4 months ago)

Author:

abarral

Message:

Revert cosp*/ from the trunk, as it's external code
Add missing bits from FCM2 source

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/parasol.F90 (modified) (3 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/parasol.F90

@@ -81 +81 @@
     ! Lum_norm=f(PARASOL_SZA,tau_cloud) derived from adding-doubling calculations
     !        valid ONLY ABOVE OCEAN (albedo_sfce=5%)
-    !        valid only in one viewing direction (theta_v=30°, phi_s-phi_v=320°)
+    !        valid only in one viewing direction (theta_v=30�, phi_s-phi_v=320�)
     !        based on adding-doubling radiative transfer computation
     !        for PARASOL_TAU values (0 to 100) and for PARASOL_SZA values (0 to 80)
@@ -97 +97 @@
 
     ! Relative fraction of the opt. thick due to liquid or ice clouds
-    WHERE (tautot_S(1:npoints) > 0.)
+    WHERE (tautot_S(1:npoints) .gt. 0.)
        frac_taucol_liq(1:npoints) = tautot_S_liq(1:npoints) / tautot_S(1:npoints)
        frac_taucol_ice(1:npoints) = tautot_S_ice(1:npoints) / tautot_S(1:npoints)
@@ -118 +118 @@
     DO it=1,PARASOL_NREFL
        DO ny=1,PARASOL_NTAU-1
-          WHERE (tautot_S(1:npoints) >= PARASOL_TAU(ny).and. &
-                 tautot_S(1:npoints) <= PARASOL_TAU(ny+1))
+          WHERE (tautot_S(1:npoints) .ge. PARASOL_TAU(ny).and. &
+                 tautot_S(1:npoints) .le. PARASOL_TAU(ny+1))
              rlumA_mod(1:npoints,it) = aA(it,ny)*tautot_S(1:npoints) + bA(it,ny)
              rlumB_mod(1:npoints,it) = aB(it,ny)*tautot_S(1:npoints) + bB(it,ny)