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

palm/trunk/SCRIPTS/NCL/profiles.ncl

@@ -33 +33 @@
     
 begin
+
+   ;***************************************************
+   ; Retrieving the NCL version used
+   ;***************************************************
+   
+   ncl_version_ch = systemfunc("ncl -V") 
+   ncl_version    = stringtofloat(ncl_version_ch)
 
    ;***************************************************
@@ -241 +248 @@
 
    
-   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'")
@@ -247 +254 @@
       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
@@ -1088 +1104 @@
    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"
@@ -1164 +1176 @@
    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
+
    wks=gsn_open_wks(format_out,file_out)
    gsn_define_colormap(wks,"rainbow+white") 
@@ -4092 +4112 @@
    end do
 
+   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)
@@ -4101 +4123 @@
          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)
+            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