Ignore:
Timestamp:
Apr 14, 2025, 11:29:13 PM (3 months ago)
Author:
aborella
Message:

Comments + small bugfixes

File:
1 edited

Legend:

Unmodified
Added
Removed
  • TabularUnified LMDZ6/branches/contrails/libf/phylmd/lmdz_aviation.f90

    r5613 r5615  
    179179        pdf_q_above_qcritcont = 0.
    180180      ELSEIF ( pdf_y .LT. -10. ) THEN
    181         pdf_fra_above_qcritcont = 1.
    182         pdf_q_above_qcritcont = qclr(i) / clrfra(i)
     181        pdf_fra_above_qcritcont = clrfra(i)
     182        pdf_q_above_qcritcont = qclr(i)
    183183      ELSE
    184184        pdf_y = EXP( pdf_y )
     
    196196        pdf_q_above_qsat = 0.
    197197      ELSEIF ( pdf_y .LT. -10. ) THEN
    198         pdf_fra_above_qsat = 1.
    199         pdf_q_above_qsat = qclr(i) / clrfra(i)
     198        pdf_fra_above_qsat = clrfra(i)
     199        pdf_q_above_qsat = qclr(i)
    200200      ELSE
    201201        pdf_y = EXP( pdf_y )
Note: See TracChangeset for help on using the changeset viewer.