Changeset 4441 for palm/trunk/SOURCE/data_output_2d.f90

Timestamp:

Mar 4, 2020 7:20:35 PM (4 years ago)

Author:

suehring

Message:

Change order of dimension in surface arrays %frac, %emissivity and %albedo to allow for better vectorization in the radiation interactions; Set back turbulent length scale to 8 x grid spacing in the parametrized mode for the synthetic turbulence generator (was accidentally changed in last commit)

File:

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

palm/trunk/SOURCE/data_output_2d.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
+! Change order of dimension in surface array %frac to allow for better 
+! vectorization.
 ! 
 ! Former revisions:
@@ -426 +427 @@
                       i                   = surf_usm_h%i(m)            
                       j                   = surf_usm_h%j(m)
-                      local_pf(i,j,nzb+1) = surf_usm_h%frac(ind_veg_wall,m)  *  &
+                      local_pf(i,j,nzb+1) = surf_usm_h%frac(m,ind_veg_wall)  *  &
                                             surf_usm_h%wghf_eb(m)        +      &
-                                            surf_usm_h%frac(ind_pav_green,m) *  &
+                                            surf_usm_h%frac(m,ind_pav_green) *  &
                                             surf_usm_h%wghf_eb_green(m)  +      &
-                                            surf_usm_h%frac(ind_wat_win,m)   *  &
+                                            surf_usm_h%frac(m,ind_wat_win)   *  &
                                             surf_usm_h%wghf_eb_window(m)
                    ENDDO
@@ -755 +756 @@
                       j   = surf_usm_h%j(m)
                       local_pf(i,j,nzb+1) =                                          &
-                                 ( surf_usm_h%frac(ind_veg_wall,m)  *                &
+                                 ( surf_usm_h%frac(m,ind_veg_wall)  *                &
                                    surf_usm_h%r_a(m)       +                         & 
-                                   surf_usm_h%frac(ind_pav_green,m) *                &
+                                   surf_usm_h%frac(m,ind_pav_green) *                &
                                    surf_usm_h%r_a_green(m) +                         & 
-                                   surf_usm_h%frac(ind_wat_win,m)   *                &
+                                   surf_usm_h%frac(m,ind_wat_win)   *                &
                                    surf_usm_h%r_a_window(m) )
                    ENDDO