Changeset 3026
- Timestamp:
- May 22, 2018 10:30:53 AM (7 years ago)
- Location:
- palm/trunk/SOURCE
- Files:
-
- 7 edited
-
diagnostic_quantities_mod.f90 (modified) (2 diffs)
-
land_surface_model_mod.f90 (modified) (7 diffs)
-
large_scale_forcing_nudging_mod.f90 (modified) (3 diffs)
-
microphysics_mod.f90 (modified) (4 diffs)
-
modules.f90 (modified) (6 diffs)
-
radiation_model_mod.f90 (modified) (2 diffs)
-
surface_mod.f90 (modified) (6 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/diagnostic_quantities_mod.f90
r2839 r3026 23 23 ! 24 24 ! Former revisions: 25 ! ----------------- !25 ! ----------------- 26 26 ! $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 27 31 ! Bugfix for Kessler microphysics 28 32 ! … … 94 98 e_s = magnus( t_l ) 95 99 ! 96 !-- Computation of saturation humidity:100 !-- Computation of saturation mixing ratio: 97 101 q_s = 0.622_wp * e_s / ( hyp(k) - e_s ) 98 102 ! -
palm/trunk/SOURCE/land_surface_model_mod.f90
r3014 r3026 25 25 ! ----------------- 26 26 ! $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 27 31 ! Bugfix: set some initial values 28 32 ! Bugfix: domain bounds of local_pf corrected … … 520 524 pavement_heat_capacity = 9999999.9_wp, & !< volumetric heat capacity of pavement (e.g. roads) (J/m3/K) 521 525 pavement_heat_conduct = 9999999.9_wp, & !< heat conductivity for pavements (e.g. roads) (W/m/K) 522 q_s = 0.0_wp, & !< saturation specific humidity526 q_s = 0.0_wp, & !< saturation water vapor mixing ratio 523 527 residual_moisture = 9999999.9_wp, & !< NAMELIST m_res 524 528 rho_cp, & !< rho_surface * cp … … 1947 1951 ENDIF 1948 1952 ! 1949 !-- Calculate saturation specific humidity1953 !-- Calculate saturation water vapor mixing ratio 1950 1954 q_s = 0.622_wp * e_s / ( surface_pressure - e_s ) 1951 1955 ! … … 2225 2229 2226 2230 ! 2227 !-- Calculate surface specific humidity2231 !-- Calculate surface water vapor mixing ratio 2228 2232 IF ( humidity ) THEN 2229 2233 CALL calc_q_surface … … 2238 2242 ! Description: 2239 2243 ! ------------ 2240 !> Calculation of specific humidityof the skin layer (surface). It is assumend2244 !> Calculation of mixing ratio of the skin layer (surface). It is assumend 2241 2245 !> that the skin is always saturated. 2242 2246 !------------------------------------------------------------------------------! … … 2259 2263 2260 2264 ! 2261 !-- Calculate specific humidityat saturation2265 !-- Calculate mixing ratio at saturation 2262 2266 q_s = 0.622_wp * e_s / ( surface_pressure - e_s ) 2263 2267 … … 2265 2269 2266 2270 ! 2267 !-- Calculate specific humidityat surface2271 !-- Calculate mixing ratio at surface 2268 2272 IF ( cloud_physics ) THEN 2269 2273 q(k+k_off,j+j_off,i+i_off) = resistance * q_s + & -
palm/trunk/SOURCE/large_scale_forcing_nudging_mod.f90
r2970 r3026 25 25 ! ----------------- 26 26 ! $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 27 31 ! Bugfix in old large-scale forcing mode 28 32 ! … … 101 105 102 106 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 humidityinterpolated to vertical grid (nudging)107 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: qnudge !< vertical profile of water vapor mixing ratio interpolated to vertical grid (nudging) 104 108 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: tnudge !< vertical profile of nudging time scale interpolated to vertical grid (nudging) 105 109 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 humiditytendency 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) 107 111 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 humiditytendency 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) 109 113 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: ug_vert !< vertical profile of geostrophic wind component in x-direction interpolated to vertical grid (large scale forcing) 110 114 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: unudge !< vertical profile of wind component in x-direction interpolated to vertical grid (nudging) … … 124 128 REAL(wp), DIMENSION(:), ALLOCATABLE :: pt_surf !< time-dependent surface temperature (large scale forcing) 125 129 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) 127 131 128 132 SAVE -
palm/trunk/SOURCE/microphysics_mod.f90
r2718 r3026 25 25 ! ----------------- 26 26 ! $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 27 31 ! Corrected "Former revisions" section 28 32 ! … … 2000 2004 ) 2001 2005 ! 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 ) 2004 2008 alpha = 0.622_wp * l_d_r * l_d_cp / ( t_l * t_l ) 2005 2009 q_s = q_s * ( 1.0_wp + alpha * q_1d(k) ) / & … … 2132 2136 ) 2133 2137 ! 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 ) 2136 2140 alpha = 0.622_wp * l_d_r * l_d_cp / ( t_l * t_l ) 2137 2141 q_s = q_s * ( 1.0_wp + alpha * q_1d(k) ) / & … … 2581 2585 ) 2582 2586 ! 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 ) 2585 2589 alpha = 0.622_wp * l_d_r * l_d_cp / ( t_l * t_l ) 2586 2590 q_s = q_s * ( 1.0_wp + alpha * q_1d(k) ) / ( 1.0_wp + alpha * q_s ) -
palm/trunk/SOURCE/modules.f90
r3014 r3026 25 25 ! ----------------- 26 26 ! $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 27 31 ! Added default values of u_max, v_max, and w_max to avoid floating invalid 28 32 ! during spinup … … 644 648 REAL(wp), DIMENSION(:), ALLOCATABLE :: ptdf_y !< damping factor for potential temperature in y-direction 645 649 REAL(wp), DIMENSION(:), ALLOCATABLE :: pt_init !< initial profile of potential temperature 646 REAL(wp), DIMENSION(:), ALLOCATABLE :: q_init !< initial profile of specific humidity650 REAL(wp), DIMENSION(:), ALLOCATABLE :: q_init !< initial profile of total water mixing ratio 647 651 !< (or total water content with active cloud physics) 648 652 REAL(wp), DIMENSION(:), ALLOCATABLE :: rdf !< rayleigh damping factor for velocity components … … 763 767 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: pt !< potential temperature 764 768 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: pt_p !< prognostic value of potential temperature 765 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: q !< specific humidity769 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: q !< mixing ratio 766 770 !< (or total water content with active cloud physics) 767 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: q_p !< prognostic value of specific humidity771 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: q_p !< prognostic value of mixing ratio 768 772 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: qc !< cloud water content 769 773 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: qc_p !< cloud water content … … 860 864 REAL(wp), DIMENSION(:,:,:), POINTER, CONTIGUOUS :: pt !< pointer: potential temperature 861 865 REAL(wp), DIMENSION(:,:,:), POINTER, CONTIGUOUS :: pt_p !< pointer: prognostic value of potential temperature 862 REAL(wp), DIMENSION(:,:,:), POINTER, CONTIGUOUS :: q !< pointer: specific humidity863 REAL(wp), DIMENSION(:,:,:), POINTER, CONTIGUOUS :: q_p !< pointer: prognostic value of specific humidity866 REAL(wp), DIMENSION(:,:,:), POINTER, CONTIGUOUS :: q !< pointer: mixing ratio 867 REAL(wp), DIMENSION(:,:,:), POINTER, CONTIGUOUS :: q_p !< pointer: prognostic value of mixing ratio 864 868 REAL(wp), DIMENSION(:,:,:), POINTER, CONTIGUOUS :: qc !< pointer: cloud water content 865 869 REAL(wp), DIMENSION(:,:,:), POINTER, CONTIGUOUS :: qc_p !< pointer: cloud water content … … 953 957 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: prr_av !< avg. precipitation rate 954 958 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: pt_av !< avg. potential temperature 955 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: q_av !< avg. specific humidity959 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: q_av !< avg. mixing ratio 956 960 !< (or total water content with active cloud physics) 957 961 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: qc_av !< avg. cloud water content … … 962 966 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: ql_vp_av !< avg. liquid water weighting factor 963 967 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) 965 969 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: rho_ocean_av !< avg. ocean density 966 970 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: s_av !< avg. passive scalar -
palm/trunk/SOURCE/radiation_model_mod.f90
r3016 r3026 28 28 ! ----------------- 29 29 ! $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 30 34 ! Revised structure of reading svf data according to PALM coding standard: 31 35 ! svf_code_field/len and fsvf removed, error messages PA0493 and PA0494 added, … … 590 594 591 595 REAL(wp), DIMENSION(:), ALLOCATABLE :: hyp_snd, & !< hypostatic pressure from sounding data (hPa) 592 q_snd, & !< specific humidityfrom sounding data (kg/kg) - dummy at the moment596 q_snd, & !< mixing ratio from sounding data (kg/kg) - dummy at the moment 593 597 rrtm_tsfc, & !< dummy array for storing surface temperature 594 598 t_snd !< actual temperature from sounding data (hPa) -
palm/trunk/SOURCE/surface_mod.f90
r2977 r3026 26 26 ! ----------------- 27 27 ! $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 28 32 ! Implement changes from branch radiation (r2948-2971) with minor modifications, 29 33 ! plus some formatting. … … 244 248 245 249 REAL(wp), DIMENSION(:), ALLOCATABLE :: pt1 !< Potential temperature at first grid level 246 REAL(wp), DIMENSION(:), ALLOCATABLE :: qv1 !< Specific humidityat first grid level250 REAL(wp), DIMENSION(:), ALLOCATABLE :: qv1 !< mixing ratio at first grid level 247 251 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: css !< scaling parameter chemical species 248 252 ! … … 1036 1040 ! 1037 1041 !-- Arrays for storing potential temperature and 1038 !-- specific humidityat first grid level1042 !-- mixing ratio at first grid level 1039 1043 DEALLOCATE ( surfaces%pt1 ) 1040 1044 DEALLOCATE ( surfaces%qv1 ) … … 1160 1164 ! 1161 1165 !-- Arrays for storing potential temperature and 1162 !-- specific humidityat first grid level1166 !-- mixing ratio at first grid level 1163 1167 ALLOCATE ( surfaces%pt1(1:surfaces%ns) ) 1164 1168 ALLOCATE ( surfaces%qv1(1:surfaces%ns) ) … … 1414 1418 ! 1415 1419 !-- Arrays for storing potential temperature and 1416 !-- specific humidityat first grid level1420 !-- mixing ratio at first grid level 1417 1421 DEALLOCATE ( surfaces%pt1 ) 1418 1422 DEALLOCATE ( surfaces%qv1 ) … … 1536 1540 ! 1537 1541 !-- Arrays for storing potential temperature and 1538 !-- specific humidityat first grid level1542 !-- mixing ratio at first grid level 1539 1543 ALLOCATE ( surfaces%pt1(1:surfaces%ns) ) 1540 1544 ALLOCATE ( surfaces%qv1(1:surfaces%ns) )
Note: See TracChangeset
for help on using the changeset viewer.
