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

Timestamp:

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

Author:

heinze

Message:

NCL scripts allow the output of png files

File:

• palm/trunk/SCRIPTS/NCL/timeseries.ncl (modified) (7 diffs)

palm/trunk/SCRIPTS/NCL/timeseries.ncl

@@ -35 +35 @@
 
    ;***************************************************
+   ; Retrieving the NCL version used
+   ;***************************************************
+   
+   ncl_version_ch = systemfunc("ncl -V") 
+   ncl_version    = stringtofloat(ncl_version_ch)
+
+   ;***************************************************
    ; Retrieving the double quote character
    ;***************************************************
@@ -53 +60 @@
    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'")
@@ -59 +66 @@
       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
 
    if (over .NE. 0 .AND. over .NE. 1) then
@@ -275 +291 @@
    res@gsnDraw                 = False
    res@gsnFrame                = False
-   res@gsnPaperOrientation     = "portrait"
-   res@gsnPaperWidth           = 8.27
-   res@gsnPaperHeight          = 11.69
-   res@gsnPaperMargin          = 0.79
    res@tmXBMode                = True
    res@tmYLMode                = True
@@ -308 +320 @@
    ; ***************************************************
    
+   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)           
    gsn_define_colormap(wks_ps,"rainbow+white") 
@@ -1207 +1227 @@
    ; merge plots onto one page 
    ; ***************************************************
+
+   no_frames = 0
  
    if ((format_out .EQ. "eps" .OR. format_out .EQ. "epsi") .AND. n .gt. no_rows*no_columns) then
@@ -1217 +1239 @@
          if ( np + no_rows*no_columns .gt. n) then   
             gsn_panel(wks_ps, plot_ps(np:n),(/no_rows,no_columns/),resP)
+            no_frames = no_frames + 1            
          else
             gsn_panel(wks_ps, plot_ps(np:np+no_rows*no_columns-1),(/no_rows,no_columns/),resP)
+            no_frames = no_frames + 1
          end if
       end do
    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