Changeset 4749 for palm

Timestamp:

Oct 16, 2020 11:21:27 AM (4 years ago)

Author:

maronga

Message:

bugfix in palm_csd

File:

• palm/trunk/SCRIPTS/palm_csd (modified) (2 diffs)

palm/trunk/SCRIPTS/palm_csd

@@ -25 +25 @@
 # -----------------
 # $Id$
+# Bugfix: processing of patch height faulty for multiple domains
+#
+# 4746 2020-10-15 17:16:30Z maronga
 # Add support for backyard high vegetation, which is commonly classified as 
 # short grass
@@ -1102 +1105 @@
 
 #        2D loop in order to avoid memory problems with large arrays
-         for ii in range(0,domain_nx[i]):
+         for iii in range(0,domain_nx[i]):
             for jj in range(0,domain_ny[i]):            
 
-               lad[0:patch_nz+1,jj,ii] = np.where( (lad[0:patch_nz+1,jj,ii] == fillvalues["tree_data"]),lad_patch[0:patch_nz+1,jj,ii], lad[0:patch_nz+1,jj,ii] )
+               lad[0:patch_nz+1,jj,iii] = np.where( (lad[0:patch_nz+1,jj,iii] == fillvalues["tree_data"]),lad_patch[0:patch_nz+1,jj,iii], lad[0:patch_nz+1,jj,iii] )
 
 #     Remove high vegetation wherever it is replaced by a leaf area density. This should effectively remove all high vegetation pixels