Changeset 4555 for palm

Timestamp:

Jun 5, 2020 9:52:00 PM (4 years ago)

Author:

moh.hefny

Message:

Bugfix in averaging PC and MRT related quantities

File:

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

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Bugfix in averaging PC and MRT related quantities
+! 
+! 4552 2020-06-02 20:33:29Z moh.hefny
 ! Bugfix in IF statement in the emissivity coupling parameter for radiation-RTM
 ! 
@@ -9863 +9866 @@
 
     INTEGER(iwp) ::  i !<
+    INTEGER(iwp) ::  imrt !< index of MRT
     INTEGER(iwp) ::  j !<
     INTEGER(iwp) ::  k !<
     INTEGER(iwp) ::  l, m !< index of current surface element
-
+    
     INTEGER(iwp)                                       :: ids, idsint_u, idsint_l, isurf
     CHARACTER(LEN=varnamelength)                       :: var
@@ -10719 +10723 @@
           CASE ( 'rtm_rad_pc_inlw' )
               DO l = 1, npcbl
-                 pcbinlw_av(:) = pcbinlw_av(:) / REAL( average_count_3d, kind=wp )
+                 pcbinlw_av(l) = pcbinlw_av(l) / REAL( average_count_3d, kind=wp )
               ENDDO
 
           CASE ( 'rtm_rad_pc_insw' )
               DO l = 1, npcbl
-                 pcbinsw_av(:) = pcbinsw_av(:) / REAL( average_count_3d, kind=wp )
+                 pcbinsw_av(l) = pcbinsw_av(l) / REAL( average_count_3d, kind=wp )
               ENDDO
 
           CASE ( 'rtm_rad_pc_inswdir' )
               DO l = 1, npcbl
-                 pcbinswdir_av(:) = pcbinswdir_av(:) / REAL( average_count_3d, kind=wp )
+                 pcbinswdir_av(l) = pcbinswdir_av(l) / REAL( average_count_3d, kind=wp )
               ENDDO
 
           CASE ( 'rtm_rad_pc_inswdif' )
               DO l = 1, npcbl
-                 pcbinswdif_av(:) = pcbinswdif_av(:) / REAL( average_count_3d, kind=wp )
+                 pcbinswdif_av(l) = pcbinswdif_av(l) / REAL( average_count_3d, kind=wp )
               ENDDO
 
           CASE ( 'rtm_rad_pc_inswref' )
               DO l = 1, npcbl
-                 pcbinswref_av(:) = pcbinswref_av(:) / REAL( average_count_3d, kind=wp )
+                 pcbinswref_av(l) = pcbinswref_av(l) / REAL( average_count_3d, kind=wp )
               ENDDO
 
           CASE ( 'rad_mrt_lw' )
              IF ( ALLOCATED( mrtinlw_av ) )  THEN
-                mrtinlw_av(:) = mrtinlw_av(:) / REAL( average_count_3d, KIND=wp )
+                DO imrt = 1, nmrtbl
+                   mrtinlw_av(imrt) = mrtinlw_av(imrt) / REAL( average_count_3d, KIND=wp )
+                ENDDO
              ENDIF
 
           CASE ( 'rad_mrt' )
              IF ( ALLOCATED( mrt_av ) )  THEN
-                mrt_av(:) = mrt_av(:) / REAL( average_count_3d, KIND=wp )
+                DO imrt = 1, nmrtbl
+                   mrt_av(imrt) = mrt_av(imrt) / REAL( average_count_3d, KIND=wp )
+                ENDDO
              ENDIF