Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/rrtm/swni.F90

    r1990 r2401  
    419419        ZRRJ=ZRJ(JL,JN,JK) / ZRJ(JL,JN2J,JK)
    420420        ZRRK=ZRK(JL,JN,JK) / ZRK(JL,JN2J,JK)
    421         ZW2(JL,1) = LOG( ZRRJ ) * ZRR
    422         ZW2(JL,2) = LOG( ZRRK ) * ZRR
     421!        ZW2(JL,1) = LOG( ZRRJ ) * ZRR
     422!        ZW2(JL,2) = LOG( ZRRK ) * ZRR
     423!--correction Olivier Boucher based on ECMWF code
     424        ZW2(JL,1) = LOG( MAX(1.0_JPRB,ZRRJ) ) * ZRR
     425        ZW2(JL,2) = LOG( MAX(1.0_JPRB,ZRRK) ) * ZRR
    423426      ENDDO
    424427
Note: See TracChangeset for help on using the changeset viewer.