Changeset 534 for palm/trunk/SCRIPTS/NCL/spectra.ncl

Timestamp:

Apr 29, 2010 1:58:28 PM (14 years ago)

Author:

heinze

Message:

Bugfix in spectra.ncl concerning output in png and small changes in layout

File:

• palm/trunk/SCRIPTS/NCL/spectra.ncl (modified) (9 diffs)

palm/trunk/SCRIPTS/NCL/spectra.ncl

@@ -401 +401 @@
       plot  = new(dim*dimz,graphic)
       np=dimt
-      res@lgTitleString = "Time [h]"
+      res@lgTitleString = "Time (h)"
    else
       plot  = new(dim*dimt,graphic)
@@ -463 +463 @@
          a=getvardims(temp_att)
          b=dimsizes(a)
+         delete(temp_att)
 
          if (height_level(0) .NE. -1)then
@@ -506 +507 @@
             end do
             if (norm_x .NE. 1.)then
-               res@tiXAxisString = "k~B~x~N~ ["+unit_x+"]"
+               res@tiXAxisString = "k~B~x~N~ ("+unit_x+")"
             else
                if (norm_height .EQ. 1)then
-                  res@tiXAxisString = "k~B~x~N~ x z [1]"
+                  res@tiXAxisString = "k~B~x~N~ x z (1"
                else
-                  res@tiXAxisString = "k~B~x~N~ [1/m]"
+                  res@tiXAxisString = "k~B~x~N~ (1/m)"
                end if
             end if
@@ -522 +523 @@
             end do
             if (norm_x .NE. 1.)then
-               res@tiXAxisString = "k~B~x~N~ ["+unit_x+"]"
+               res@tiXAxisString = "k~B~x~N~ ("+unit_x+")"
             else
                if (norm_height .EQ. 1)then
-                  res@tiXAxisString = "k~B~x~N~ x z [1]"
+                  res@tiXAxisString = "k~B~x~N~ x z (1)"
                else
-                  res@tiXAxisString = "k~B~x~N~ [1/m]"
+                  res@tiXAxisString = "k~B~x~N~ (1/m)"
                end if 
             end if
@@ -557 +558 @@
                res@gsnLeftString      = vNam(varn)
                res@gsnRightString     = "Height = "+z(p)+"m"               
-               res@tiYAxisString      = "["+unit_y+"]"           
+               res@tiYAxisString      = "("+unit_y+")"           
                res@xyExplicitLegendLabels  = legend_label              
                if (norm_height .EQ. 1)then
@@ -596 +597 @@
                      end if
                      if (q .EQ. 0)then
-                        res@tiYAxisString      = "["+unit_y+"]"                 
+                        res@tiYAxisString      = "("+unit_y+")"                 
                         res@gsnLeftString      = vNam(varn)
                         res@gsnRightString     = "Time = "+legend_label(p)+"h"
@@ -617 +618 @@
                            lgres@lgLineColors = color
                         end if
-                        lgres@lgTitleString = "Height [m]"  
+                        lgres@lgTitleString = "Height (m)"  
                         lgres@lgLabelFont        = "helvetica"
-                        lgres@lgLabelFontHeightF     = font_size_legend 
+                        lgres@lgLabelFontHeightF     = font_size_legend*6 
                         lgres@lgTitleFontHeightF     = font_size      
                         lgres@vpWidthF           = 0.12           
-                        lgres@vpHeightF          = font_size_legend/5*dimz           
+                        lgres@vpHeightF          = font_size_legend*(dimz+3)           
   
                         lbid = gsn_create_legend(wks,dimz,legend_label_z,lgres)       
@@ -669 +670 @@
    resP@txFontHeightF          = 0.015
 
+   no_frames = 0
+
    if ((format_out .EQ. "eps" .OR. format_out .EQ. "epsi") .AND. n .gt. no_rows*no_columns) then
       gsn_panel(wks,plot,(/n,1/),resP)
@@ -678 +681 @@
          if( (i+no_rows*no_columns) .gt. (n-1)) then
             gsn_panel(wks,plot(i:n-1),(/no_rows,no_columns/),resP)
+            no_frames = no_frames + 1
          else
-           gsn_panel(wks,plot(i:i+no_rows*no_columns-1),(/no_rows,no_columns/),resP)
+            gsn_panel(wks,plot(i:i+no_rows*no_columns-1),(/no_rows,no_columns/),resP)
+            no_frames = no_frames + 1   
          end if
       end do