Ignore:
Timestamp:
Oct 16, 2018 8:39:37 PM (3 years ago)
Author:
knoop
Message:

Introduced global constant rd_d_rv=0.622

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/bulk_cloud_model_mod.f90

    r3274 r3361  
    186186               molecular_weight_of_water, pi, rho_l, rho_s, r_d, r_v, vanthoff,&
    187187               exner_function, exner_function_invers, ideal_gas_law_rho,       &
    188                ideal_gas_law_rho_pt, barometric_formula
     188               ideal_gas_law_rho_pt, barometric_formula, rd_d_rv
    189189
    190190    USE control_parameters,                                                    &
     
    33083308!
    33093309!--       Computation of saturation mixing ratio:
    3310           q_s   = 0.622_wp * e_s / ( hyp(k) -  e_s )
    3311           alpha = 0.622_wp * lv_d_rd * lv_d_cp / ( t_l * t_l )
     3310          q_s   = rd_d_rv * e_s / ( hyp(k) -  e_s )
     3311          alpha = rd_d_rv * lv_d_rd * lv_d_cp / ( t_l * t_l )
    33123312          q_s   = q_s * ( 1.0_wp + alpha * q(k,j,i) ) /                        &
    33133313                        ( 1.0_wp + alpha * q_s )
     
    34233423!
    34243424!--       Computation of saturation mixing ratio:
    3425           q_s   = 0.622_wp * e_s / ( hyp(k) - e_s )
    3426           alpha = 0.622_wp * lv_d_rd * lv_d_cp / ( t_l * t_l )
     3425          q_s   = rd_d_rv * e_s / ( hyp(k) - e_s )
     3426          alpha = rd_d_rv * lv_d_rd * lv_d_cp / ( t_l * t_l )
    34273427          q_s   = q_s * ( 1.0_wp + alpha * q(k,j,i) ) /               &
    34283428                        ( 1.0_wp + alpha * q_s )
     
    38033803!
    38043804!--          Computation of saturation mixing ratio:
    3805              q_s   = 0.622_wp * e_s / ( hyp(k) - e_s )
    3806              alpha = 0.622_wp * lv_d_rd * lv_d_cp / ( t_l * t_l )
     3805             q_s   = rd_d_rv * e_s / ( hyp(k) - e_s )
     3806             alpha = rd_d_rv * lv_d_rd * lv_d_cp / ( t_l * t_l )
    38073807             q_s   = q_s * ( 1.0_wp + alpha * q(k,j,i) ) / ( 1.0_wp + alpha * q_s )
    38083808!
     
    42474247!
    42484248!--    Computation of saturation mixing ratio:
    4249        q_s   = 0.622_wp * e_s / ( hyp(k) - e_s )
     4249       q_s   = rd_d_rv * e_s / ( hyp(k) - e_s )
    42504250!
    42514251!--    Correction factor
    4252        alpha = 0.622_wp * lv_d_rd * lv_d_cp / ( t_l * t_l )
     4252       alpha = rd_d_rv * lv_d_rd * lv_d_cp / ( t_l * t_l )
    42534253!
    42544254!--    Correction of the approximated value
Note: See TracChangeset for help on using the changeset viewer.