Changeset 4045

Timestamp:

Dec 7, 2021, 9:14:49 AM (2 years ago)

Author:

idelkadi

Message:

Correction in radlwsw_m.F90:
Addition of machine threshold to avoid a division by 0

File:

• LMDZ6/trunk/libf/phylmd/radlwsw_m.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/radlwsw_m.F90

@@ -335 +335 @@
   REAL(KIND=8) ZSWDIFFUSEBAND(klon,NSW), &  ! SW DN flux in diffuse albedo band 
                ZSWDIRECTBAND(klon,NSW)      ! SW DN flux in direct albedo band
+  REAL(KIND=8) seuilmach
 #endif
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -506 +507 @@
   ZAEROSOL_OLD = 0.0
   ZAEROSOL = 0.0
+  seuilmach=tiny(seuilmach)
 #endif
 
@@ -1453 +1455 @@
       !--fraction of diffuse radiation in surface SW downward radiation
       DO i = 1, kdlon
-       IF (fract(i).GT.0.0) THEN
          zdir=SUM(PSFSWDIR(i,:))
          zdif=SUM(PSFSWDIF(i,:))
+       IF (fract(i).GT.0.0.and.(zdir+zdif).gt.seuilmach) THEN
          zsolswfdiff(i) = zdif/(zdir+zdif)
        ELSE  !--night