Ignore:
Timestamp:
Sep 19, 2023, 5:31:25 PM (15 months ago)
Author:
jleconte
Message:

fixed some bugs and errors

Location:
trunk/LMDZ.GENERIC/libf/phystd
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.GENERIC/libf/phystd/calc_rayleigh.F90

    r2890 r3043  
    174174                  tauvarvar=tauvarvar+tauconsti(igas)*tauvari(igas)
    175175                  tauvar=tauvar+0.0*0.0*gfrac(igas)
    176                elseif(igas.eq.igas_CH4) then !!!!!!! methane rayleigh verification
    177                   tauvari(igas_CH4) = (4.6662E-4+4.02E-6/wl**2)**2 / wl**4
    178                   tauvarvar=tauvarvar+tauconsti(igas_CH4)*tauvari(igas_CH4)*gfrac(igas_CH4)
    179                   tauvar=tauvar+0.0*0.0*gfrac(igas)
    180176               else
    181177                  tauvar=tauvar+tauconsti(igas)*tauvari(igas)*gfrac(igas)
  • trunk/LMDZ.GENERIC/libf/phystd/rain_generic.F90

    r3041 r3043  
    240240         precip_rate_tmp(1:ngrid) = 0.0
    241241
     242         dq_rain_generic_vap(1:ngrid,1:nlayer,1:nq) = 0.0
     243         dq_rain_generic_cld(1:ngrid,1:nlayer,1:nq) = 0.0
     244         dqsrain_generic(1:ngrid,1:nq) = 0.0
     245         dqssnow_generic(1:ngrid,1:nq) = 0.0
     246
    242247         ! calculate saturation mixing ratio
    243248         do k = 1, nlayer
Note: See TracChangeset for help on using the changeset viewer.