Changeset 2079 for LMDZ5/trunk/libf


Ignore:
Timestamp:
Jul 7, 2014, 3:38:36 PM (10 years ago)
Author:
lguez
Message:

Protect against division by 0 in computation of proba_notrig in procedure physiq: initialize s2 at some low value instead of 0.

Location:
LMDZ5/trunk/libf/phylmd
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/cva_driver.F90

    r2007 r2079  
    882882                           asupmaxmin, cbmf, plfc, wbeff)
    883883
    884         PRINT *, 'cv3p1_closure-> plfc,wbeff ', plfc(1), wbeff(1)
     884        if (prt_level >= 10) &
     885             PRINT *, 'cv3p1_closure-> plfc,wbeff ', plfc(1), wbeff(1)
    885886      END IF
    886887    END IF ! iflag_con.eq.3
  • LMDZ5/trunk/libf/phylmd/thermcell_main.F90

    r1998 r2079  
    876876   depth(:)=0.
    877877   n2(:)=0.
    878    s2(:)=0.
     878   s2(:)=100. ! some low value, arbitrary
    879879   s_max(:)=0.
    880880
     
    894894!!fin jyg
    895895     else
    896        s2(ig)=0.
    897896       n2(ig)=0.
    898897       s_max(ig)=0.
Note: See TracChangeset for help on using the changeset viewer.