Changeset 758 for palm/trunk/SCRIPTS/NCL

Timestamp:

Sep 15, 2011 12:03:51 PM (13 years ago)

Author:

heinze

Message:

Bugfixes for NCL scripts timeseries.ncl and cross_sections.ncl concerning usage of norm_t and vector plot

Location:

palm/trunk/SCRIPTS/NCL

Files:

• cross_sections.ncl (modified) (9 diffs)
• timeseries.ncl (modified) (1 diff)

palm/trunk/SCRIPTS/NCL/cross_sections.ncl

@@ -1415 +1415 @@
    end if 
 
-   check = True
-   v1=0
-   v2=0
-   no_var=0
-   n=0
-   no_zu1=0
+   check  = True
+   v1     = 0
+   v2     = 0
+   no_var = 0
+   n      = 0
+   no_zu1 = 0
 
    ;****************************************************
@@ -1547 +1547 @@
    
    do varn=dim-1,0,1   
-    
-      ;if (vector .EQ. 1) then    
-      ;   check_vec1 = isStrSubset( vec1,","+vNam(varn)+",")
-      ;   check_vec2 = isStrSubset( vec2,","+vNam(varn)+",")
-      ;end if
            
       if ( vNam(varn) .eq. "zw" .or. vNam(varn) .eq. "yv" .or.          \
@@ -1630 +1625 @@
 
    if (vector .EQ. 1 .AND. plotvec .EQ. "plotvec") then
-      plot=new((/no_time*no_layer*no_var - no_time*(no_layer-1)*no_zu1 \
-                                         + no_time*no_layer/),graphic)
+      plot=new((/no_time*no_layer/),graphic)
    else
       plot=new((/no_time*no_layer*no_var - no_time*(no_layer-1)*no_zu1/)\
@@ -1839 +1833 @@
                      cs_res@gsnCenterString = "t=" + \
                        decimalPlaces(t_all(li)/3600,2,True) +"h  z"+level
-                     plot(n) = gsn_csm_contour(wks_ps,\
+
+                     if (vector .EQ. 1 .AND. plotvec .EQ. "plotvec") then
+                        ;nothing to be done
+                     else
+                         plot(n) = gsn_csm_contour(wks_ps,\
                                            data(varn,li,lo-los,:,:),cs_res) 
+                     end if
+
+
                      if (vector .EQ. 1 .AND. check_vecp) then
                         vecres                 = True           ; vector only
@@ -1897 +1898 @@
                      cs_res@gsnCenterString = "t=" + \
                            decimalPlaces(t_all(lo)/3600,2,True) + "h  z"+ level
-                     
-                     plot(n) = gsn_csm_contour(wks_ps,\
+
+                     if (vector .EQ. 1 .AND. plotvec .EQ. "plotvec") then
+                        ;nothing to be done
+                     else
+                        plot(n) = gsn_csm_contour(wks_ps,\
                                             data(varn,lo,li-lis,:,:),cs_res)
+                     end if
+
                      if (vector .EQ. 1 .AND. check_vecp) then
                         vecres                 = True           ; vector only
@@ -1941 +1947 @@
                      cs_res@gsnCenterString = "t=" + \
                            decimalPlaces(t_all(li)/3600,2,True) + "h  y"+ level
-                     plot(n) = gsn_csm_contour(wks_ps,\
+
+                     if (vector .EQ. 1 .AND. plotvec .EQ. "plotvec") then
+                        ;nothing to be done
+                     else
+                         plot(n) = gsn_csm_contour(wks_ps,\
                                               data(varn,li,:,lo-los,:),cs_res)
+                     end if
+
                      if (vector .EQ. 1 .AND. check_vecp) then
                         vecres                 = True           ; vector only 
@@ -1975 +1987 @@
                      cs_res@gsnCenterString = "t=" + \
                            decimalPlaces(t_all(lo)/3600,2,True) + "h  y"+ level
-                     plot(n) = gsn_csm_contour(wks_ps,\
-                                              data(varn,lo,:,li-lis,:),cs_res)
+
+                     if (vector .EQ. 1 .AND. plotvec .EQ. "plotvec") then
+                        ;nothing to be done
+                     else
+                         plot(n) = gsn_csm_contour(wks_ps,\
+                                              data(varn,lo,:,li-lis,:),cs_res)  
+                     end if
+
                      if (vector .EQ. 1 .AND. check_vecp) then
                         vecres                 = True           ; vector only 
@@ -2018 +2036 @@
                      cs_res@gsnCenterString = "t=" + \
                           decimalPlaces(t_all(li)/3600,2,True) + "h  x"+ level
-                     plot(n) = gsn_csm_contour(wks_ps,\
+
+                     if (vector .EQ. 1 .AND. plotvec .EQ. "plotvec") then
+                        ;nothing to be done
+                     else
+                        plot(n) = gsn_csm_contour(wks_ps,\
                                           data(varn,li,:,:,lo-los),cs_res)
+                     end if
+
                      if (vector .EQ. 1 .AND. check_vecp) then
                         vecres                 = True           ; vector only 
@@ -2052 +2076 @@
                      cs_res@gsnCenterString = "t=" + \
                            decimalPlaces(t_all(lo)/3600,2,True) + "h  x"+ level
-                     plot(n) = gsn_csm_contour(wks_ps,\
+
+                     if (vector .EQ. 1 .AND. plotvec .EQ. "plotvec") then
+                        ;nothing to be done
+                     else
+                        plot(n) = gsn_csm_contour(wks_ps,\
                                           data(varn,lo,:,:,li-lis),cs_res)
+                     end if
+
                      if (vector .EQ. 1 .AND. check_vecp)then
                         vecres                 = True           ; vector only 

palm/trunk/SCRIPTS/NCL/timeseries.ncl

@@ -328 +328 @@
    res@tmYLLabelFont           = "helvetica" 
    res@xyLineColors            = (/237/) 
-   res@trXMaxF                 = t_all(start_time_step)
-   res@trXMinF                 = t_all(end_time_step)
+   res@trXMaxF                 = t(start_time_step)
+   res@trXMinF                 = t(end_time_step)
    
    res@lgLabelFontHeightF     = 0.02