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/land_surface_model_mod.f90

    r3349 r3361  
    486486
    487487    USE basic_constants_and_equations_mod,                                     &
    488         ONLY:  c_p, g, lv_d_cp, l_v, magnus, rho_l, r_d, r_v
     488        ONLY:  c_p, g, lv_d_cp, l_v, magnus, rho_l, r_d, r_v, rd_d_rv
    489489
    490490    USE calc_mean_profile_mod,                                                 &
     
    19991999!
    20002000!--       Calculate vapour pressure
    2001           e  = surf%qv1(m) * surface_pressure / ( surf%qv1(m) + 0.622_wp )
     2001          e  = surf%qv1(m) * surface_pressure / ( surf%qv1(m) + rd_d_rv )
    20022002          f3 = EXP ( - surf%g_d(m) * (e_s - e) )
    20032003       ELSE
     
    20402040!
    20412041!--    Calculate saturation water vapor mixing ratio
    2042        q_s = 0.622_wp * e_s / ( surface_pressure - e_s )
     2042       q_s = rd_d_rv * e_s / ( surface_pressure - e_s )
    20432043!
    20442044!--    In case of dewfall, set evapotranspiration to zero
     
    20792079                       / ( surf_t_surface%var_1d(m) - 29.65_wp)**2 )
    20802080
    2081        dq_s_dt = 0.622_wp * e_s_dt / ( surface_pressure - e_s_dt )
     2081       dq_s_dt = rd_d_rv * e_s_dt / ( surface_pressure - e_s_dt )
    20822082!
    20832083!--    Calculate net radiation radiation without longwave outgoing flux because
     
    23502350!
    23512351!--       Calculate mixing ratio at saturation
    2352           q_s = 0.622_wp * e_s / ( surface_pressure - e_s )
     2352          q_s = rd_d_rv * e_s / ( surface_pressure - e_s )
    23532353
    23542354          resistance = surf%r_a(m) / ( surf%r_a(m) + surf%r_s(m) + 1E-5_wp )
Note: See TracChangeset for help on using the changeset viewer.