Changeset 3861

Timestamp:

Apr 4, 2019 6:27:41 AM (6 years ago)

Author:

maronga

Message:

bugfix in radiation model (rad_lw_out_change_0)

File:

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

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Bugfix: factor of 4.0 required instead of 3.0 in calculation of rad_lw_out_change_0
+! 
+! 3859 2019-04-03 20:30:31Z maronga
 ! Added some descriptions
 ! 
@@ -2940 +2943 @@
                              + surf%rad_lw_in - surf%rad_lw_out
 
-                surf%rad_lw_out_change_0 = 3.0_wp * emissivity_urb * sigma_sb  &
+                surf%rad_lw_out_change_0 = 4.0_wp * emissivity_urb * sigma_sb  &
                                            * (t_rad_urb)**3
 
@@ -2985 +2988 @@
                                       + surf%frac(ind_wat_win,m)   *           &
                                         surf%emissivity(ind_wat_win,m)         &
-                                      ) * 3.0_wp * sigma_sb                    &
+                                      ) * 4.0_wp * sigma_sb                    &
                                       * ( surf%pt_surface(m) * exner(nzb) )** 3
 
@@ -3114 +3117 @@
                                     * surf%rad_lw_in
 
-                surf%rad_lw_out_change_0 = 3.0_wp * emissivity_urb * sigma_sb  &
+                surf%rad_lw_out_change_0 = 4.0_wp * emissivity_urb * sigma_sb  &
                                            * t_rad_urb**3