Changeset 1571 for palm/trunk/SOURCE/radiation_model.f90
- Timestamp:
- Mar 12, 2015 4:12:49 PM (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/radiation_model.f90
r1552 r1571 20 20 ! Current revisions: 21 21 ! ----------------- 22 ! 22 ! Added missing KIND attribute. Removed upper-case variable names 23 23 ! 24 24 ! Former revisions: … … 68 68 LOGICAL :: radiation = .FALSE. !: flag parameter indicating wheather the radiation model is used 69 69 70 REAL(wp), PARAMETER :: SW_0 = 1368.0, &!: solar constant70 REAL(wp), PARAMETER :: sw_0 = 1368.0_wp, & !: solar constant 71 71 pi = 3.14159265358979323_wp, & 72 sigma_ SB= 5.67E-8_wp !: Stefan-Boltzmann constant72 sigma_sb = 5.67E-8_wp !: Stefan-Boltzmann constant 73 73 74 74 REAL(wp) :: albedo = 0.2_wp, & !: NAMELIST alpha … … 126 126 irad_scheme, lambda, net_radiation, rad_net, rad_net_av, radiation, & 127 127 radiation_clearsky, radiation_constant, radiation_rrtm, & 128 radiation_scheme, rad_sw_in, rad_sw_in_av, sigma_ SB, &128 radiation_scheme, rad_sw_in, rad_sw_in_av, sigma_sb, & 129 129 time_radiation, time_utc_init 130 130 … … 231 231 232 232 k = nzb_s_inner(j,i) 233 rad_sw_in(j,i) = SW_0 * sky_trans * zenith233 rad_sw_in(j,i) = sw_0 * sky_trans * zenith 234 234 rad_sw_out(j,i) = - alpha(j,i) * rad_sw_in(j,i) 235 rad_lw_out(j,i) = - sigma_ SB* (pt(k,j,i) * exn)**4236 rad_lw_in(j,i) = 0.8_wp * sigma_ SB* (pt(k+1,j,i) * exn)**4235 rad_lw_out(j,i) = - sigma_sb * (pt(k,j,i) * exn)**4 236 rad_lw_in(j,i) = 0.8_wp * sigma_sb * (pt(k+1,j,i) * exn)**4 237 237 rad_net(j,i) = rad_sw_in(j,i) + rad_sw_out(j,i) & 238 238 + rad_lw_in(j,i) + rad_lw_out(j,i)
Note: See TracChangeset
for help on using the changeset viewer.