Changeset 959 for palm/trunk/SCRIPTS/NCL/profiles.ncl

Timestamp:

Jul 24, 2012 1:13:41 PM (12 years ago)

Author:

hoffmann

Message:

bugfixes in palmplot and cross_profiles

File:

• palm/trunk/SCRIPTS/NCL/profiles.ncl (modified) (3 diffs)

palm/trunk/SCRIPTS/NCL/profiles.ncl

@@ -381 +381 @@
 
       ;************************************************************************
-      ; currently, more than 3 plots cannot be drawn together
+      ; currently, more than 3 plots cannot be drawn together. the profiles 
+      ; are split to a maximum of 3 plots per panel.
       ;***********************************************************************
 
@@ -387 +388 @@
          if(number_comb_all(com_i) .GT. 3) then
             print(" ")
-            print("Currently, it is not possible to plot more than three"+\
-                  "profiles together. Hence, they are split to a maximum"+\
-                  "of three profiles per coordinate plane.")
+            print("Currently, it is not possible to plot more than three "+\
+                  "profiles together ("+c_var_all(com_i)+"). Hence, they "+\
+                  "are split to a maximum of three profiles per "+\
+                  "coordinate plane.")
             print(" ") 
             c_var_all_temp = c_var_all
@@ -410 +412 @@
             end do
             number_comb_all(com_i+1) = max(dimsizes(c_var_long)) - 3 
-            c_var_all(com_i+2:max_com_i) = c_var_all_temp(com_i+1:max_com_i-1)
-            number_comb_all(com_i+2:max_com_i) = number_comb_all_temp(com_i+1:max_com_i-1)
+            if (com_i+2 .le. max_com_i) then
+               c_var_all(com_i+2:max_com_i) = c_var_all_temp(com_i+1:max_com_i-1)
+               number_comb_all(com_i+2:max_com_i) = number_comb_all_temp(com_i+1:max_com_i-1)
+            end if
             delete(c_var_all_temp)
             delete(number_comb_all_temp)