Changeset 4746

Timestamp:

Oct 15, 2020 5:16:30 PM (4 years ago)

Author:

maronga

Message:

support for backyard trees in palm_csd

File:

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

palm/trunk/SCRIPTS/palm_csd

@@ -20 +20 @@
 # Current revisions:
 # -----------------
-# 
-# 
+#
+#
 # Former revisions:
 # -----------------
 # $Id$
+# Add support for backyard high vegetation, which is commonly classified as 
+# short grass
+# 
+# 4490 2020-04-10 05:32:45Z maronga
 # Bugfix: call of canopy generator was not made if tree heights were only
 # available via patch_height, but not in tree_height
@@ -1066 +1070 @@
       lai_high = np.where( (lad[0,:,:] == fillvalues["tree_data"]) & ( ( (vegetation_type == 4) | (vegetation_type == 5) | (vegetation_type == 6) |(vegetation_type == 7) | (vegetation_type == 17) | (vegetation_type == 18) ) & ( (patch_height == fillvalues["tree_data"]) | (patch_height >= domain_dz[i])) ),1.0,fillvalues["tree_data"])
 
+
+#     Treat all pixels where short grass is defined, but where a patch_height >= dz is found, as high vegetation (often the case in backyards)
+      lai_high = np.where( (lai_high == fillvalues["tree_data"]) & (patch_height >= domain_dz[i]) & (vegetation_type == 3), 1.0, lai_high)
+                          
 #     Now, assign either the default LAI for high vegetation or keep 1.0 from the lai_high array.
       lai_high = np.where( (lai_high != fillvalues["tree_data"]) & (lai == fillvalues["tree_data"]), settings_lai_high_default, lai_high)
 
-#     If lai values are available in the lai array, write them on the lai_high array
+#     If LAI values are available in the LAI array, write them on the lai_high array
       lai_high = np.where( (lai_high != fillvalues["tree_data"]) & (lai != fillvalues["tree_data"]), lai, lai_high)
+
 
 #     Define a patch height wherever it is missing, but where a high vegetation LAI was set
       patch_height = np.where ( (lai_high != fillvalues["tree_data"]) & (patch_height == fillvalues["tree_data"]), settings_patch_height_default, patch_height)
+
 
 #     Remove pixels where street trees were already set
@@ -1091 +1101 @@
          lad_patch, patch_nz, status = process_patch(domain_dz[i],patch_height,max(zlad),lai_high,settings_lai_alpha,settings_lai_beta)
 
-         lad[0:patch_nz+1,:,:] = np.where( (lad[0:patch_nz+1,:,:] == fillvalues["tree_data"]),lad_patch[0:patch_nz+1,:,:], lad[0:patch_nz+1,:,:] )
+#        2D loop in order to avoid memory problems with large arrays
+         for ii 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] )
 
 #     Remove high vegetation wherever it is replaced by a leaf area density. This should effectively remove all high vegetation pixels