Changeset 3026 for palm


Ignore:
Timestamp:
May 22, 2018 10:30:53 AM (6 years ago)
Author:
schwenkel
Message:

Changed the name specific humidity to mixing ratio

Location:
palm/trunk/SOURCE
Files:
7 edited

Legend:

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

    r2839 r3026  
    2323!
    2424! Former revisions:
    25 ! -----------------!
     25! -----------------
    2626! $Id$
     27! Changed the name specific humidity to mixing ratio, since we are computing
     28! mixing ratios.
     29!
     30! 2839 2018-02-27 09:49:06Z schwenkel
    2731! Bugfix for Kessler microphysics
    2832!
     
    9498       e_s = magnus( t_l )
    9599!
    96 !--    Computation of saturation humidity:
     100!--    Computation of saturation mixing ratio:
    97101       q_s   = 0.622_wp * e_s / ( hyp(k) - e_s )
    98102!
  • palm/trunk/SOURCE/land_surface_model_mod.f90

    r3014 r3026  
    2525! -----------------
    2626! $Id$
     27! Changed the name specific humidity to mixing ratio, since we are computing
     28! mixing ratios.
     29!
     30! 3014 2018-05-09 08:42:38Z maronga
    2731! Bugfix: set some initial values
    2832! Bugfix: domain bounds of local_pf corrected
     
    520524                pavement_heat_capacity = 9999999.9_wp,  & !< volumetric heat capacity of pavement (e.g. roads) (J/m3/K)
    521525                pavement_heat_conduct  = 9999999.9_wp,  & !< heat conductivity for pavements (e.g. roads) (W/m/K)
    522                 q_s = 0.0_wp,                           & !< saturation specific humidity
     526                q_s = 0.0_wp,                           & !< saturation water vapor mixing ratio
    523527                residual_moisture = 9999999.9_wp,       & !< NAMELIST m_res
    524528                rho_cp,                                 & !< rho_surface * cp
     
    19471951       ENDIF
    19481952!
    1949 !--    Calculate saturation specific humidity
     1953!--    Calculate saturation water vapor mixing ratio
    19501954       q_s = 0.622_wp * e_s / ( surface_pressure - e_s )
    19511955!
     
    22252229
    22262230!
    2227 !-- Calculate surface specific humidity
     2231!-- Calculate surface water vapor mixing ratio
    22282232    IF ( humidity )  THEN
    22292233       CALL calc_q_surface
     
    22382242! Description:
    22392243! ------------
    2240 !> Calculation of specific humidity of the skin layer (surface). It is assumend
     2244!> Calculation of mixing ratio of the skin layer (surface). It is assumend
    22412245!> that the skin is always saturated.
    22422246!------------------------------------------------------------------------------!
     
    22592263
    22602264!
    2261 !--       Calculate specific humidity at saturation
     2265!--       Calculate mixing ratio at saturation
    22622266          q_s = 0.622_wp * e_s / ( surface_pressure - e_s )
    22632267
     
    22652269
    22662270!
    2267 !--       Calculate specific humidity at surface
     2271!--       Calculate mixing ratio at surface
    22682272          IF ( cloud_physics )  THEN
    22692273             q(k+k_off,j+j_off,i+i_off) = resistance * q_s +                   &
  • palm/trunk/SOURCE/large_scale_forcing_nudging_mod.f90

    r2970 r3026  
    2525! -----------------
    2626! $Id$
     27! Changed the name specific humidity to mixing ratio, since we are computing
     28! mixing ratios.
     29!
     30! 2970 2018-04-13 15:09:23Z suehring
    2731! Bugfix in old large-scale forcing mode
    2832!
     
    101105
    102106    REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  ptnudge     !< vertical profile of pot. temperature interpolated to vertical grid (nudging)
    103     REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  qnudge      !< vertical profile of specific humidity interpolated to vertical grid (nudging)
     107    REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  qnudge      !< vertical profile of water vapor mixing ratio interpolated to vertical grid (nudging)
    104108    REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  tnudge      !< vertical profile of nudging time scale interpolated to vertical grid (nudging) 
    105109    REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  td_lsa_lpt  !< temperature tendency due to large scale advection (large scale forcing)
    106     REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  td_lsa_q    !< specific humidity tendency due to large scale advection (large scale forcing)
     110    REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  td_lsa_q    !< water vapor mixing ratio tendency due to large scale advection (large scale forcing)
    107111    REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  td_sub_lpt  !< temperature tendency due to subsidence/ascent (large scale forcing)
    108     REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  td_sub_q    !< specific humidity tendency due to subsidence/ascent (large scale forcing)
     112    REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  td_sub_q    !< water vapor mixing ratio tendency due to subsidence/ascent (large scale forcing)
    109113    REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  ug_vert     !< vertical profile of geostrophic wind component in x-direction interpolated to vertical grid (large scale forcing)
    110114    REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  unudge      !< vertical profile of wind component in x-direction interpolated to vertical grid (nudging)
     
    124128    REAL(wp), DIMENSION(:), ALLOCATABLE ::  pt_surf       !< time-dependent surface temperature (large scale forcing)
    125129    REAL(wp), DIMENSION(:), ALLOCATABLE ::  qsws_surf     !< time-dependent surface latent heat flux (large scale forcing)
    126     REAL(wp), DIMENSION(:), ALLOCATABLE ::  q_surf        !< time-dependent surface specific humidity (large scale forcing)
     130    REAL(wp), DIMENSION(:), ALLOCATABLE ::  q_surf        !< time-dependent surface water vapor mixing ratio (large scale forcing)
    127131
    128132    SAVE
  • palm/trunk/SOURCE/microphysics_mod.f90

    r2718 r3026  
    2525! -----------------
    2626! $Id$
     27! Changed the name specific humidity to mixing ratio, since we are computing
     28! mixing ratios.
     29!
     30! 2718 2018-01-02 08:49:38Z maronga
    2731! Corrected "Former revisions" section
    2832!
     
    20002004                                 )
    20012005!
    2002 !--       Computation of saturation humidity:
    2003           q_s   = 0.622_wp * e_s / ( hyp(k) - 0.378_wp * e_s )
     2006!--       Computation of saturation mixing ratio:
     2007          q_s   = 0.622_wp * e_s / ( hyp(k) - e_s )
    20042008          alpha = 0.622_wp * l_d_r * l_d_cp / ( t_l * t_l )
    20052009          q_s   = q_s * ( 1.0_wp + alpha * q_1d(k) ) /                         &
     
    21322136                                 )
    21332137!
    2134 !--       Computation of saturation humidity:
    2135           q_s   = 0.622_wp * e_s / ( hyp(k) - 0.378_wp * e_s )
     2138!--       Computation of saturation mixing ratio:
     2139          q_s   = 0.622_wp * e_s / ( hyp(k) - e_s )
    21362140          alpha = 0.622_wp * l_d_r * l_d_cp / ( t_l * t_l )
    21372141          q_s   = q_s * ( 1.0_wp + alpha * q_1d(k) ) /               &
     
    25812585                                  )
    25822586!
    2583 !--          Computation of saturation humidity:
    2584              q_s   = 0.622_wp * e_s / ( hyp(k) - 0.378_wp * e_s )
     2587!--          Computation of saturation mixing ratio:
     2588             q_s   = 0.622_wp * e_s / ( hyp(k) - e_s )
    25852589             alpha = 0.622_wp * l_d_r * l_d_cp / ( t_l * t_l )
    25862590             q_s   = q_s * ( 1.0_wp + alpha * q_1d(k) ) / ( 1.0_wp + alpha * q_s )
  • palm/trunk/SOURCE/modules.f90

    r3014 r3026  
    2525! -----------------
    2626! $Id$
     27! Changed the name specific humidity to mixing ratio, since we are computing
     28! mixing ratios.
     29!
     30! 3014 2018-05-09 08:42:38Z maronga
    2731! Added default values of u_max, v_max, and w_max to avoid floating invalid
    2832! during spinup
     
    644648    REAL(wp), DIMENSION(:), ALLOCATABLE ::  ptdf_y                 !< damping factor for potential temperature in y-direction
    645649    REAL(wp), DIMENSION(:), ALLOCATABLE ::  pt_init                !< initial profile of potential temperature
    646     REAL(wp), DIMENSION(:), ALLOCATABLE ::  q_init                 !< initial profile of specific humidity
     650    REAL(wp), DIMENSION(:), ALLOCATABLE ::  q_init                 !< initial profile of total water mixing ratio
    647651                                                                   !< (or total water content with active cloud physics)
    648652    REAL(wp), DIMENSION(:), ALLOCATABLE ::  rdf                    !< rayleigh damping factor for velocity components
     
    763767    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET ::  pt         !< potential temperature
    764768    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET ::  pt_p       !< prognostic value of potential temperature
    765     REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET ::  q          !< specific humidity
     769    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET ::  q          !< mixing ratio 
    766770                                                                   !< (or total water content with active cloud physics)
    767     REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET ::  q_p        !< prognostic value of specific humidity
     771    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET ::  q_p        !< prognostic value of mixing ratio
    768772    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET ::  qc         !< cloud water content
    769773    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET ::  qc_p       !< cloud water content
     
    860864    REAL(wp), DIMENSION(:,:,:), POINTER, CONTIGUOUS ::  pt         !< pointer: potential temperature
    861865    REAL(wp), DIMENSION(:,:,:), POINTER, CONTIGUOUS ::  pt_p       !< pointer: prognostic value of potential temperature
    862     REAL(wp), DIMENSION(:,:,:), POINTER, CONTIGUOUS ::  q          !< pointer: specific humidity
    863     REAL(wp), DIMENSION(:,:,:), POINTER, CONTIGUOUS ::  q_p        !< pointer: prognostic value of specific humidity
     866    REAL(wp), DIMENSION(:,:,:), POINTER, CONTIGUOUS ::  q          !< pointer: mixing ratio
     867    REAL(wp), DIMENSION(:,:,:), POINTER, CONTIGUOUS ::  q_p        !< pointer: prognostic value of mixing ratio
    864868    REAL(wp), DIMENSION(:,:,:), POINTER, CONTIGUOUS ::  qc         !< pointer: cloud water content
    865869    REAL(wp), DIMENSION(:,:,:), POINTER, CONTIGUOUS ::  qc_p       !< pointer: cloud water content
     
    953957    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET ::  prr_av        !< avg. precipitation rate
    954958    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET ::  pt_av         !< avg. potential temperature
    955     REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET ::  q_av          !< avg. specific humidity
     959    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET ::  q_av          !< avg. mixing ratio
    956960                                                                      !< (or total water content with active cloud physics)
    957961    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET ::  qc_av         !< avg. cloud water content
     
    962966    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET ::  ql_vp_av      !< avg. liquid water weighting factor
    963967    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET ::  qr_av         !< avg. rain water content
    964     REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET ::  qv_av         !< avg. water vapor content (specific humidity)
     968    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET ::  qv_av         !< avg. water vapor content (mixing ratio)
    965969    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET ::  rho_ocean_av  !< avg. ocean density
    966970    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET ::  s_av          !< avg. passive scalar
  • palm/trunk/SOURCE/radiation_model_mod.f90

    r3016 r3026  
    2828! -----------------
    2929! $Id$
     30! Changed the name specific humidity to mixing ratio, since we are computing
     31! mixing ratios.
     32!
     33! 3016 2018-05-09 10:53:37Z Giersch
    3034! Revised structure of reading svf data according to PALM coding standard:
    3135! svf_code_field/len and fsvf removed, error messages PA0493 and PA0494 added,
     
    590594
    591595    REAL(wp), DIMENSION(:), ALLOCATABLE :: hyp_snd,     & !< hypostatic pressure from sounding data (hPa)
    592                                            q_snd,       & !< specific humidity from sounding data (kg/kg) - dummy at the moment
     596                                           q_snd,       & !< mixing ratio from sounding data (kg/kg) - dummy at the moment
    593597                                           rrtm_tsfc,   & !< dummy array for storing surface temperature
    594598                                           t_snd          !< actual temperature from sounding data (hPa)
  • palm/trunk/SOURCE/surface_mod.f90

    r2977 r3026  
    2626! -----------------
    2727! $Id$
     28! Changed the name specific humidity to mixing ratio, since we are computing
     29! mixing ratios.
     30!
     31! 2977 2018-04-17 10:27:57Z kanani
    2832! Implement changes from branch radiation (r2948-2971) with minor modifications,
    2933! plus some formatting.
     
    244248
    245249       REAL(wp), DIMENSION(:), ALLOCATABLE ::  pt1       !< Potential temperature at first grid level
    246        REAL(wp), DIMENSION(:), ALLOCATABLE ::  qv1       !< Specific humidity at first grid level
     250       REAL(wp), DIMENSION(:), ALLOCATABLE ::  qv1       !< mixing ratio at first grid level
    247251       REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  css     !< scaling parameter chemical species
    248252!
     
    10361040!
    10371041!--    Arrays for storing potential temperature and
    1038 !--    specific humidity at first grid level
     1042!--    mixing ratio at first grid level
    10391043       DEALLOCATE ( surfaces%pt1 )
    10401044       DEALLOCATE ( surfaces%qv1 )
     
    11601164!
    11611165!--    Arrays for storing potential temperature and
    1162 !--    specific humidity at first grid level
     1166!--    mixing ratio at first grid level
    11631167       ALLOCATE ( surfaces%pt1(1:surfaces%ns) )
    11641168       ALLOCATE ( surfaces%qv1(1:surfaces%ns) )
     
    14141418!
    14151419!--    Arrays for storing potential temperature and
    1416 !--    specific humidity at first grid level
     1420!--    mixing ratio at first grid level
    14171421       DEALLOCATE ( surfaces%pt1 )
    14181422       DEALLOCATE ( surfaces%qv1 )
     
    15361540!
    15371541!--    Arrays for storing potential temperature and
    1538 !--    specific humidity at first grid level
     1542!--    mixing ratio at first grid level
    15391543       ALLOCATE ( surfaces%pt1(1:surfaces%ns) )
    15401544       ALLOCATE ( surfaces%qv1(1:surfaces%ns) )
Note: See TracChangeset for help on using the changeset viewer.