Changeset 1691 for palm/trunk/SOURCE/data_output_3d.f90

Timestamp:

Oct 26, 2015 4:17:44 PM (9 years ago)

Author:

maronga

Message:

various bugfixes and modifications of the atmosphere-land-surface-radiation interaction. Completely re-written routine to calculate surface fluxes (surface_layer_fluxes.f90) that replaces prandtl_fluxes. Minor formatting corrections and renamings

File:

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

palm/trunk/SOURCE/data_output_3d.f90

@@ -14 +14 @@
 ! PALM. If not, see <http://www.gnu.org/licenses/>.
 !
-! Copyright 1997-2014 Leibniz Universitaet Hannover
+! Copyright 1997-2015 Leibniz Universitaet Hannover
 !--------------------------------------------------------------------------------!
 !
 ! Current revisions:
 ! ------------------
-! 
+! Added output of radiative heating rates for RRTMG
 ! 
 ! Former revisions:
@@ -138 +138 @@
     USE radiation_model_mod,                                                   &
         ONLY:  rad_lw_in, rad_lw_in_av, rad_lw_out, rad_lw_out_av,             &
-               rad_sw_in, rad_sw_in_av, rad_sw_out, rad_sw_out_av
+               rad_lw_cs_hr, rad_lw_cs_hr_av, rad_lw_hr, rad_lw_hr_av,         &
+               rad_sw_in, rad_sw_in_av, rad_sw_out, rad_sw_out_av,             &
+               rad_sw_cs_hr, rad_sw_cs_hr_av, rad_sw_hr, rad_sw_hr_av
 
 
@@ -495 +497 @@
              ENDIF
 
+          CASE ( 'rad_sw_cs_hr' )
+             IF ( av == 0 )  THEN
+                to_be_resorted => rad_sw_cs_hr
+             ELSE
+                to_be_resorted => rad_sw_cs_hr_av
+             ENDIF
+
+          CASE ( 'rad_sw_hr' )
+             IF ( av == 0 )  THEN
+                to_be_resorted => rad_sw_hr
+             ELSE
+                to_be_resorted => rad_sw_hr_av
+             ENDIF
+
           CASE ( 'rad_lw_in' )
              IF ( av == 0 )  THEN
@@ -507 +523 @@
              ELSE
                 to_be_resorted => rad_lw_out_av
+             ENDIF
+
+          CASE ( 'rad_lw_cs_hr' )
+             IF ( av == 0 )  THEN
+                to_be_resorted => rad_lw_cs_hr
+             ELSE
+                to_be_resorted => rad_lw_cs_hr_av
+             ENDIF
+
+          CASE ( 'rad_lw_hr' )
+             IF ( av == 0 )  THEN
+                to_be_resorted => rad_lw_hr
+             ELSE
+                to_be_resorted => rad_lw_hr_av
              ENDIF