Changeset 1571 for palm/trunk/SOURCE/radiation_model.f90

Timestamp:

Mar 12, 2015 4:12:49 PM (9 years ago)

Author:

maronga

Message:

minor bugfixes in lasm/radiation

File:

• palm/trunk/SOURCE/radiation_model.f90 (modified) (4 diffs)

palm/trunk/SOURCE/radiation_model.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Added missing KIND attribute. Removed upper-case variable names
 ! 
 ! Former revisions:
@@ -68 +68 @@
     LOGICAL :: radiation = .FALSE.  !: flag parameter indicating wheather the radiation model is used
 
-    REAL(wp), PARAMETER :: SW_0 = 1368.0, &       !: solar constant  
+    REAL(wp), PARAMETER :: sw_0 = 1368.0_wp, &    !: solar constant  
                            pi = 3.14159265358979323_wp, &
-                           sigma_SB  = 5.67E-8_wp !: Stefan-Boltzmann constant
+                           sigma_sb  = 5.67E-8_wp !: Stefan-Boltzmann constant
  
     REAL(wp) :: albedo = 0.2_wp,             & !: NAMELIST alpha
@@ -126 +126 @@
            irad_scheme, lambda, net_radiation, rad_net, rad_net_av, radiation, &
            radiation_clearsky, radiation_constant, radiation_rrtm,             &
-           radiation_scheme, rad_sw_in, rad_sw_in_av, sigma_SB,                &
+           radiation_scheme, rad_sw_in, rad_sw_in_av, sigma_sb,                &
            time_radiation, time_utc_init 
 
@@ -231 +231 @@
 
              k = nzb_s_inner(j,i)
-             rad_sw_in(j,i)  = SW_0 * sky_trans * zenith
+             rad_sw_in(j,i)  = sw_0 * sky_trans * zenith
              rad_sw_out(j,i) = - alpha(j,i) * rad_sw_in(j,i)
-             rad_lw_out(j,i) = - sigma_SB * (pt(k,j,i) * exn)**4
-             rad_lw_in(j,i)  = 0.8_wp * sigma_SB * (pt(k+1,j,i) * exn)**4
+             rad_lw_out(j,i) = - sigma_sb * (pt(k,j,i) * exn)**4
+             rad_lw_in(j,i)  = 0.8_wp * sigma_sb * (pt(k+1,j,i) * exn)**4
              rad_net(j,i)    = rad_sw_in(j,i) + rad_sw_out(j,i)                &
                                 + rad_lw_in(j,i) + rad_lw_out(j,i)