Changeset 3726 for palm/trunk/SCRIPTS/palm_csd_files

Timestamp:

Feb 7, 2019 6:22:49 PM (6 years ago)

Author:

maronga

Message:

bugfixes in palm_csd

File:

• palm/trunk/SCRIPTS/palm_csd_files/palm_csd_tools.py (modified) (3 diffs)

palm/trunk/SCRIPTS/palm_csd_files/palm_csd_tools.py

@@ -25 +25 @@
 # -----------------
 # $Id$
+# Index bound bugfix
+# 
+# 3668 2019-01-14 12:49:24Z maronga
 # Minor changes
 # 
@@ -90 +93 @@
       for l in range(0,len(x)):
          for m in range(0,len(y)):
-            for k in range(0,len(k_tmp+1)):
+            for k in range(1,len(k_tmp+1)):
                if k_tmp[k] > array[m,l]:
-                  array[m,l] = k_tmp[k]-0.5*dz
+                  array[m,l] = k_tmp[k-1]
                   break
 
@@ -105 +108 @@
       array_3d = np.ones((len(k_tmp),len(y),len(x)))
       
-      for l in range(0,len(x)-1):
-         for m in range(0,len(y)-1):
-            for k in range(0,len(k_tmp)-1):
+      for l in range(0,len(x)):
+         for m in range(0,len(y)):
+            for k in range(0,len(k_tmp)):
                if k_tmp[k] > array_2d[m,l]:
                   array_3d[k,m,l] = 0