Changeset 2401 for LMDZ5

Timestamp:

Nov 22, 2015, 5:14:24 PM (9 years ago)

Author:

oboucher

Message:

Correction of a bug in swni.F90
should enable new SW radiation to run smoothly

File:

• LMDZ5/trunk/libf/phylmd/rrtm/swni.F90 (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/rrtm/swni.F90

@@ -419 +419 @@
         ZRRJ=ZRJ(JL,JN,JK) / ZRJ(JL,JN2J,JK)
         ZRRK=ZRK(JL,JN,JK) / ZRK(JL,JN2J,JK)
-        ZW2(JL,1) = LOG( ZRRJ ) * ZRR
-        ZW2(JL,2) = LOG( ZRRK ) * ZRR
+!        ZW2(JL,1) = LOG( ZRRJ ) * ZRR
+!        ZW2(JL,2) = LOG( ZRRK ) * ZRR
+!--correction Olivier Boucher based on ECMWF code 
+        ZW2(JL,1) = LOG( MAX(1.0_JPRB,ZRRJ) ) * ZRR
+        ZW2(JL,2) = LOG( MAX(1.0_JPRB,ZRRK) ) * ZRR
       ENDDO