Changeset 2753 for palm/trunk/SOURCE/flow_statistics.f90

Timestamp:

Jan 16, 2018 2:16:49 PM (6 years ago)

Author:

suehring

Message:

Tile approach for spectral albedo implemented.

File:

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

palm/trunk/SOURCE/flow_statistics.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Tile approach for spectral albedo implemented. 
+! 
+! 2718 2018-01-02 08:49:38Z maronga
 ! Corrected "Former revisions" section
 ! 
@@ -1165 +1168 @@
                       m = surf_def_h(0)%start_index(j,i)
                       sums_l(nzb,108,tn)  = sums_l(nzb,108,tn)  +              &
-                                           surf_def_h(0)%rrtm_aldif(m) * rmask(j,i,sr)
+                                   surf_def_h(0)%rrtm_aldif(0,m) * rmask(j,i,sr)
                       sums_l(nzb,109,tn) = sums_l(nzb,109,tn)  +               &
-                                           surf_def_h(0)%rrtm_aldir(m) * rmask(j,i,sr)
+                                   surf_def_h(0)%rrtm_aldir(0,m) * rmask(j,i,sr)
                       sums_l(nzb,110,tn) = sums_l(nzb,110,tn)  +               &
-                                           surf_def_h(0)%rrtm_asdif(m) * rmask(j,i,sr)
+                                   surf_def_h(0)%rrtm_asdif(0,m) * rmask(j,i,sr)
                       sums_l(nzb,111,tn) = sums_l(nzb,111,tn)  +               &
-                                           surf_def_h(0)%rrtm_asdir(m) * rmask(j,i,sr)
+                                   surf_def_h(0)%rrtm_asdir(0,m) * rmask(j,i,sr)
                    ENDIF
                    IF ( surf_lsm_h%end_index(j,i) >=                           &
@@ -1177 +1180 @@
                       m = surf_lsm_h%start_index(j,i)
                       sums_l(nzb,108,tn)  = sums_l(nzb,108,tn)  +              &
-                                           surf_lsm_h%rrtm_aldif(m) * rmask(j,i,sr)
+                               SUM( surf_lsm_h%frac(:,m) *                     &
+                                    surf_lsm_h%rrtm_aldif(:,m) ) * rmask(j,i,sr)
                       sums_l(nzb,109,tn) = sums_l(nzb,109,tn)  +               &
-                                           surf_lsm_h%rrtm_aldir(m) * rmask(j,i,sr)
+                               SUM( surf_lsm_h%frac(:,m) *                     &
+                                    surf_lsm_h%rrtm_aldir(:,m) ) * rmask(j,i,sr)
                       sums_l(nzb,110,tn) = sums_l(nzb,110,tn)  +               &
-                                           surf_lsm_h%rrtm_asdif(m) * rmask(j,i,sr)
+                               SUM( surf_lsm_h%frac(:,m) *                     &
+                                    surf_lsm_h%rrtm_asdif(:,m) ) * rmask(j,i,sr)
                       sums_l(nzb,111,tn) = sums_l(nzb,111,tn)  +               &
-                                           surf_lsm_h%rrtm_asdir(m) * rmask(j,i,sr)
+                               SUM( surf_lsm_h%frac(:,m) *                     &
+                                    surf_lsm_h%rrtm_asdir(:,m) ) * rmask(j,i,sr)
                    ENDIF
                    IF ( surf_usm_h%end_index(j,i) >=                           &
@@ -1189 +1196 @@
                       m = surf_usm_h%start_index(j,i)
                       sums_l(nzb,108,tn)  = sums_l(nzb,108,tn)  +              &
-                                           surf_usm_h%rrtm_aldif(m) * rmask(j,i,sr)
+                               SUM( surf_usm_h%frac(:,m) *                     &
+                                    surf_usm_h%rrtm_aldif(:,m) ) * rmask(j,i,sr)
                       sums_l(nzb,109,tn) = sums_l(nzb,109,tn)  +               &
-                                           surf_usm_h%rrtm_aldir(m) * rmask(j,i,sr)
+                               SUM( surf_usm_h%frac(:,m) *                     &
+                                    surf_usm_h%rrtm_aldir(:,m) ) * rmask(j,i,sr)
                       sums_l(nzb,110,tn) = sums_l(nzb,110,tn)  +               &
-                                           surf_usm_h%rrtm_asdif(m) * rmask(j,i,sr)
+                               SUM( surf_usm_h%frac(:,m) *                     &
+                                    surf_usm_h%rrtm_asdif(:,m) ) * rmask(j,i,sr)
                       sums_l(nzb,111,tn) = sums_l(nzb,111,tn)  +               &
-                                           surf_usm_h%rrtm_asdir(m) * rmask(j,i,sr)
+                               SUM( surf_usm_h%frac(:,m) *                     &
+                                    surf_usm_h%rrtm_asdir(:,m) ) * rmask(j,i,sr)
                    ENDIF