Changeset 532 for palm/trunk/SCRIPTS/NCL/spectra.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/spectra.ncl (modified) (6 diffs)

palm/trunk/SCRIPTS/NCL/spectra.ncl

@@ -34 +34 @@
 
    ;***************************************************
+   ; Retrieving the NCL version used
+   ;***************************************************
+   
+   ncl_version_ch = systemfunc("ncl -V") 
+   ncl_version    = stringtofloat(ncl_version_ch)
+
+   ;***************************************************
    ; Retrieving the double quote character
    ;***************************************************
@@ -52 +59 @@
    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'")
@@ -58 +65 @@
       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 (log_x .NE. 0 .AND. log_x .NE. 1)then
@@ -340 +356 @@
    res@gsnDraw                 = False
    res@gsnFrame                = False
-   res@gsnPaperOrientation     = "portrait"
-   res@gsnPaperWidth           = 8.27
-   res@gsnPaperHeight          = 11.69
-   res@gsnPaperMargin          = 0.79
    res@txFont                  = "helvetica"
    res@tiMainFont              = "helvetica"
@@ -417 +429 @@
    if ( dash .eq. 0 ) then
       res@xyMonoDashPattern = True 
+   end if
+
+   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
 
@@ -664 +684 @@
    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