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

    r3337 r3361  
    4545    USE arrays_3d,                                                             &
    4646        ONLY:  hyp, p, pt, q, u, v, w
     47
     48    USE basic_constants_and_equations_mod,                                     &
     49        ONLY:  rd_d_rv
    4750
    4851    USE kinds  !< to set precision of INTEGER and REAL arrays according to PALM
     
    231234                   CALL calculate_pet_static(                                      &
    232235                           pt(k,j,i) * (hyp(k) / 100000.0_wp )**0.286_wp + t_zero, &  !< Air temperature (°C)
    233                            q(k,j,i) * hyp(k) / ( q(k,j,i) + 0.622_wp ) / 100._wp,  &  !< Vapor pressure (hPa)
     236                           q(k,j,i) * hyp(k) / ( q(k,j,i) + rd_d_rv ) / 100._wp,   &  !< Vapor pressure (hPa)
    234237                           SQRT( MAX( ( ( u(k,j,i) + u(k,j,i+1) ) * 0.5_wp )**2 +  &
    235238                                      ( ( v(k,j,i) + v(k,j+1,i) ) * 0.5_wp )**2 +  &
     
    380383                    CALL calculate_pet_static(                                      &
    381384                       pt(k,j,i) * (hyp(k) / 100000.0_wp )**0.286_wp + t_zero, &  !< Air temperature (°C)
    382                        q(k,j,i) * hyp(k) / ( q(k,j,i) + 0.622_wp ) / 100.0_wp, &  !< Vapor pressure (hPa)
     385                       q(k,j,i) * hyp(k) / ( q(k,j,i) + rd_d_rv ) / 100.0_wp, &  !< Vapor pressure (hPa)
    383386                       SQRT( MAX( ( ( u(k,j,i) + u(k,j,i+1) ) * 0.5_wp )**2 +  &
    384387                                  ( ( v(k,j,i) + v(k,j+1,i) ) * 0.5_wp )**2 +  &
Note: See TracChangeset for help on using the changeset viewer.