Changeset 4587 for palm/trunk/SOURCE/radiation_model_mod.f90

Timestamp:

Jul 6, 2020 8:53:45 AM (4 years ago)

Author:

pavelkrc

Message:

Marked RTM version 3.1, bugfixes

File:

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

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -23 +23 @@
 ! Current revisions:
 ! ------------------
-! 
-! 
+!
+!
 ! Former revisions:
 ! -----------------
 ! $Id$
+! RTM version 3.1 (see also previous commits):
+! - rotation_angle supported
+! - plant canopy box count minimized
+! - multiple enhancements and bugfixes
+! 
+! 4584 2020-06-29 13:16:14Z pavelkrc
 ! Consider only boxes with LAD>0 as plant canopy (credit: S. Schubert)
 ! 
@@ -10502 +10508 @@
               ENDDO
 
-          CASE ( 'rad_mrt_sw' )
+          CASE ( 'rtm_mrt_sw' )
              IF ( ALLOCATED( mrtinsw_av ) )  THEN
                 mrtinsw_av(:) = mrtinsw_av(:) + mrtinsw(:)
              ENDIF
 
-          CASE ( 'rad_mrt_lw' )
+          CASE ( 'rtm_mrt_lw' )
              IF ( ALLOCATED( mrtinlw_av ) )  THEN
                 mrtinlw_av(:) = mrtinlw_av(:) + mrtinlw(:)
              ENDIF
 
-          CASE ( 'rad_mrt' )
+          CASE ( 'rtm_mrt' )
              IF ( ALLOCATED( mrt_av ) )  THEN
                 mrt_av(:) = mrt_av(:) + mrt(:)
@@ -10794 +10800 @@
               ENDDO
 
-          CASE ( 'rad_mrt_lw' )
+          CASE ( 'rtm_mrt_sw' )
+             IF ( ALLOCATED( mrtinsw_av ) )  THEN
+                DO imrt = 1, nmrtbl
+                   mrtinsw_av(imrt) = mrtinsw_av(imrt) / REAL( average_count_3d, KIND=wp )
+                ENDDO
+             ENDIF
+
+          CASE ( 'rtm_mrt_lw' )
              IF ( ALLOCATED( mrtinlw_av ) )  THEN
                 DO imrt = 1, nmrtbl
@@ -10801 +10814 @@
              ENDIF
 
-          CASE ( 'rad_mrt' )
+          CASE ( 'rtm_mrt' )
              IF ( ALLOCATED( mrt_av ) )  THEN
                 DO imrt = 1, nmrtbl
@@ -10866 +10879 @@
                  'rad_sw_hr_xy', 'rad_lw_cs_hr_xz', 'rad_lw_hr_xz',               &
                  'rad_sw_cs_hr_xz', 'rad_sw_hr_xz', 'rad_lw_cs_hr_yz',            &
-                 'rad_lw_hr_yz', 'rad_sw_cs_hr_yz', 'rad_sw_hr_yz',               &
-                 'rad_mrt', 'rad_mrt_sw', 'rad_mrt_lw' )
+                 'rad_lw_hr_yz', 'rad_sw_cs_hr_yz', 'rad_sw_hr_yz' )
              grid_x = 'x'
              grid_y = 'y'