Changeset 4671 for palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90

Timestamp:

Sep 9, 2020 8:27:58 PM (4 years ago)

Author:

pavelkrc

Message:

Radiative transfer model RTM version 4.1

File:

• palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90 (modified) (2 diffs)

palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Implementation of downward facing USM and LSM surfaces
+!
+! 4647 2020-08-24 16:36:18Z suehring
 ! Change default value of synthetic turbulence adjustment as well as compute_velocity_seeds_local 
 ! By default, the random-seed computation is now distributed among several cores. Especially for 
@@ -2410 +2413 @@
        pt_surf_mean_l = pt_surf_mean_l  + pt(k,j,i)
     ENDDO
-    DO  m = 1, surf_lsm_h%ns
-       i = surf_lsm_h%i(m)
-       j = surf_lsm_h%j(m)
-       k = surf_lsm_h%k(m)
-       friction_vel_l = friction_vel_l  + surf_lsm_h%us(m)
-       shf_mean_l     = shf_mean_l      + surf_lsm_h%shf(m) * drho_air(k)
-       scale_l_l      = scale_l_l       + surf_lsm_h%ol(m)
+    DO  m = 1, surf_lsm_h(0)%ns
+       i = surf_lsm_h(0)%i(m)
+       j = surf_lsm_h(0)%j(m)
+       k = surf_lsm_h(0)%k(m)
+       friction_vel_l = friction_vel_l  + surf_lsm_h(0)%us(m)
+       shf_mean_l     = shf_mean_l      + surf_lsm_h(0)%shf(m) * drho_air(k)
+       scale_l_l      = scale_l_l       + surf_lsm_h(0)%ol(m)
        pt_surf_mean_l = pt_surf_mean_l  + pt(k,j,i)
     ENDDO
-    DO  m = 1, surf_usm_h%ns
-       i = surf_usm_h%i(m)
-       j = surf_usm_h%j(m)
-       k = surf_usm_h%k(m)
-       friction_vel_l = friction_vel_l  + surf_usm_h%us(m)
-       shf_mean_l     = shf_mean_l      + surf_usm_h%shf(m) * drho_air(k)
-       scale_l_l      = scale_l_l       + surf_usm_h%ol(m)
+    DO  m = 1, surf_usm_h(0)%ns
+       i = surf_usm_h(0)%i(m)
+       j = surf_usm_h(0)%j(m)
+       k = surf_usm_h(0)%k(m)
+       friction_vel_l = friction_vel_l  + surf_usm_h(0)%us(m)
+       shf_mean_l     = shf_mean_l      + surf_usm_h(0)%shf(m) * drho_air(k)
+       scale_l_l      = scale_l_l       + surf_usm_h(0)%ol(m)
        pt_surf_mean_l = pt_surf_mean_l  + pt(k,j,i)
     ENDDO