Changeset 532 for palm/trunk/SCRIPTS/NCL/cross_sections.ncl

Timestamp:

Apr 21, 2010 1:33:38 PM (14 years ago)

Author:

heinze

Message:

NCL scripts allow the output of png files

File:

• palm/trunk/SCRIPTS/NCL/cross_sections.ncl (modified) (8 diffs)

palm/trunk/SCRIPTS/NCL/cross_sections.ncl

@@ -35 +35 @@
 begin
 
+   ;***************************************************
+   ; Retrieving the NCL version used
+   ;***************************************************
+   
+   ncl_version_ch = systemfunc("ncl -V") 
+   ncl_version    = stringtofloat(ncl_version_ch)
+
    ; ***************************************************
    ; Retrieving the double quote character
@@ -54 +61 @@
    end if
 
-   if (format_out .NE. "x11" .AND. format_out .NE. "pdf" .AND. format_out .NE. "eps" .AND. format_out .NE. "ps" .AND. format_out .NE. "epsi" .AND. format_out .NE. "ncgm")then
+   if (format_out .NE. "x11" .AND. format_out .NE. "pdf" .AND. format_out .NE. "eps" .AND. format_out .NE. "ps" .AND. format_out .NE. "epsi" .AND. format_out .NE. "ncgm" .AND. format_out .NE. "png")then
       print(" ")
       print("'format_out = "+format_out+"' is invalid and set to'x11'")
       print(" ")
       format_out="x11"
+   end if
+
+   if (ncl_version .LE. 5.1 .AND. format_out .EQ. "png") then
+      print(" ")
+      print("Output of png files not available")
+      print("png output is avaiable with NCL version 5.2.0 and higher ")
+      print("NCL version used: " + ncl_version_ch)
+      print(" ")
+      exit
    end if
 
@@ -336 +352 @@
    cs_res@gsnDraw                 = False
    cs_res@gsnFrame                = False
-   cs_res@gsnMaximize             = True 
+    cs_res@gsnMaximize                = True
    
    cs_res@tmXBLabelFontHeightF   = font_size
@@ -1446 +1462 @@
    ; open workstation(s) 
    ; ***************************************************
+
+   if ( format_out .EQ. "pdf" .OR. format_out .EQ. "ps" ) then
+      format_out@wkPaperSize = "A4"
+   end if
+   if ( format_out .EQ. "png" ) then
+      format_out@wkWidth  = 1000
+      format_out@wkHeight = 1000
+   end if
 
    wks_ps  = gsn_open_wks(format_out,file_out)
@@ -1807 +1831 @@
    ; merge plots onto one page 
    ; ***************************************************   
+
+   no_frames = 0
+
    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
@@ -1817 +1844 @@
             if ( np + no_rows*no_columns .gt. dim_plot-1) then 
                gsn_panel(wks_ps, plot(np:dim_plot-1),(/no_rows,no_columns/),cs_resP)
+               no_frames = no_frames + 1
             else
                gsn_panel(wks_ps, plot(np:np+no_rows*no_columns-1),(/no_rows,no_columns/),cs_resP)
+               no_frames = no_frames + 1
             end if
          end do
@@ -1832 +1861 @@
             if ( np + no_rows*no_columns .gt. dim_plot-1) then
                gsn_panel(wks_ps, plot(np:dim_plot-1),(/no_rows,no_columns/),cs_resP)
+               no_frames = no_frames + 1
             else
                gsn_panel(wks_ps, plot(np:np+no_rows*no_columns-1),(/no_rows,no_columns/),cs_resP)
+               no_frames = no_frames + 1
             end if
          end do
@@ -1839 +1870 @@
    end if
 
-   print(" ")
-   print("Output to: " + file_out +"."+ format_out)
-   print(" ")   
+   if (format_out .EQ. "png" ) then
+     png_output = new((/no_frames/), string)
+     j = 0
+     do i=0, no_frames-1
+       j = i + 1
+       if (j .LE. 9) then
+         png_output(i) = file_out+".00000"+j+".png"
+       end if
+       if (j .GT. 9 .AND. j .LE. 99) then
+         png_output(i) = file_out+".0000"+j+".png"
+       end if
+       if (j .GT. 99 .AND. j .LE. 999) then
+         png_output(i) = file_out+".000"+j+".png"
+       end if
+       if (j .GT. 999) then
+         png_output(i) = file_out+".00"+j+".png"
+       end if
+
+       ;using imagemagick's convert for reducing the white 
+       ;space around the plot
+       cmd = "convert -geometry 1000x1000 -density 300 -trim " +  \
+              png_output(i) + " " + png_output(i)
+       system(cmd)
+     end do
+
+     print(" ")
+     print("Output to: "+ png_output)
+     print(" ")
+   else
+     print(" ")
+     print("Output to: " + file_out +"."+ format_out)
+     print(" ")
+   end if
  
 end