Changeset 4021 for palm/trunk/SCRIPTS/palm_csd_files/palm_csd_canopy_generator.py

Timestamp:

Jun 12, 2019 6:33:33 AM (5 years ago)

Author:

maronga

Message:

minor changes in palm_csd

File:

• palm/trunk/SCRIPTS/palm_csd_files/palm_csd_canopy_generator.py (modified) (6 diffs)

palm/trunk/SCRIPTS/palm_csd_files/palm_csd_canopy_generator.py

@@ -20 +20 @@
 # Current revisions:
 # -----------------
-# 
+# Added support for tree shape
 # 
 # Former revisions:
@@ -79 +79 @@
                tree_id_counter = tree_id_counter + 1
 #               print("   Processing tree No " +  str(tree_id_counter) + " ...", end="")
-               lad_loc, bad_loc, x_loc, y_loc, z_loc, status = process_single_tree(dx,dz,tree_type[j,i],tree_height[j,i],tree_lai[j,i],tree_dia[j,i],trunk_dia[j,i],season)
+               lad_loc, bad_loc, x_loc, y_loc, z_loc, status = process_single_tree(dx,dz,tree_type[j,i],fill,tree_height[j,i],tree_lai[j,i],tree_dia[j,i],trunk_dia[j,i],season)
                if ( np.any(lad_loc) != fill ):
  
@@ -116 +116 @@
 
 
-def process_single_tree(dx,dz,tree_type,tree_height,tree_lai,tree_dia,trunk_dia,season):
+def process_single_tree(dx,dz,tree_type,tree_shape,tree_height,tree_lai,tree_dia,trunk_dia,season):
 
 #  Set some parameters
@@ -236 +236 @@
       tree_type = int(tree_type)
 
+   if ( tree_shape == fillvalues["tree_data"] ):
+      tree_shape = default_trees[tree_type].shape
+
    if ( tree_height == fillvalues["tree_data"] ):
       tree_height = default_trees[tree_type].height
@@ -254 +257 @@
  
 #  Assign values that are not defined as user input from lookup table
-   tree_shape          = default_trees[tree_type].shape
    tree_ratio          = default_trees[tree_type].ratio
    lad_max_height      = default_trees[tree_type].lad_max_height
@@ -453 +455 @@
 
 
-#  Dreate BAD array and populate
-   bad_loc = np.where(lad_loc != fillvalues["tree_data"],lad_loc*0.01,lad_loc)
+#  Create BAD array and populate. TODO: revise as low LAD inside the foliage does not result in low BAD values.
+   bad_loc = np.where(lad_loc != fillvalues["tree_data"],(1.0 - lad_loc)*0.1,lad_loc)