Changeset 529 for palm/trunk/SCRIPTS/NCL

Timestamp:

Mar 31, 2010 10:36:23 AM (15 years ago)

Author:

heinze

Message:

palmplot expects lists which are separated by blanks. The lists have to be enclosed in single quotes in the prompt.

Location:

palm/trunk/SCRIPTS/NCL

Files:

• UserGuide_for_NCLscripts.odt (modified) (previous)
• UserGuide_for_NCLscripts.pdf (modified) (previous)
• cross_sections.ncl (modified) (6 diffs)

palm/trunk/SCRIPTS/NCL/cross_sections.ncl

@@ -1533 +1533 @@
                vecres@gsnLeftString = "t=" + decimalPlaces(t_all(lo)/3600,2,True) +"h   "+level 
             end if
-            vecres@tiXAxisString    = " " 
-            if (xyc .EQ. 1)then  
+            if (xyc .EQ. 1) then  
+               vecres@tiXAxisString = "x [m]"
+               vecres@tiYAxisString = "y [m]"   
                if (sort .EQ. "time")then                                          
                   plot(n) = gsn_csm_vector(wks_ps,vect1(li,lo-los,:,:),vect2(li,lo-los,:,:),vecres)
@@ -1542 +1543 @@
             end if
             if (xzc .EQ. 1) then
+               vecres@tiXAxisString = "x [m]"
+               vecres@tiYAxisString = "z [m]" 
                if (sort .EQ. "time")then 
                   plot(n) = gsn_csm_vector(wks_ps,vect1(li,:,lo-los,:),vect2(li,:,lo-los,:),vecres)
@@ -1549 +1552 @@
             end if
             if (yzc .EQ. 1) then
+               vecres@tiXAxisString = "y [m]"
+               vecres@tiYAxisString = "z [m]" 
                if (sort .EQ. "time")then 
                   plot(n) = gsn_csm_vector(wks_ps,vect1(li,:,:,lo-los),vect2(li,:,:,lo-los),vecres)
@@ -1562 +1567 @@
    do varn=dim-1,0,1 
 
-      if (vector .EQ. 1) then    
+      if (vector .EQ. 1 ) then    
          check_vecp = isStrSubset( plotvec,","+vNam(varn)+",") 
       end if 
@@ -1678 +1683 @@
                         vecres@gsnLeftString    = " "             ; turn off left string
                         vecres@tiXAxisString    = " "   
-                        plot_vec=gsn_csm_vector(wks_ps,vect1(lo,li,:,:),vect2(lo,li,:,:),vecres)
+                        plot_vec=gsn_csm_vector(wks_ps,vect1(lo,li-lis,:,:),vect2(lo,li-lis,:,:),vecres)
                         overlay(plot(n), plot_vec)
                      end if
@@ -1802 +1807 @@
    ; merge plots onto one page 
    ; ***************************************************   
-   if (vector .EQ. 1 .AND. plotvec .EQ. "plotvec")then 
+   if (vector .EQ. 1 .AND. plotvec .EQ. "plotvec") then 
       if ((format_out .EQ. "eps" .OR. format_out .EQ. "epsi") .AND. no_time*no_layer*(no_var+1) .gt. no_rows*no_columns) then
          gsn_panel(wks_ps,plot(0:(no_time*no_layer*(no_var+1))-1),(/no_var+1,no_layer*no_time/),cs_resP)