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

Timestamp:

Jul 20, 2012 10:57:49 AM (11 years ago)

Author:

heinze

Message:

png output adjusted for NCL versions higher than 5.2.1

File:

• palm/trunk/SCRIPTS/NCL/spectra.ncl (modified) (1 diff)

palm/trunk/SCRIPTS/NCL/spectra.ncl

@@ -751 +751 @@
      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)
+
+     if (no_frames .eq. 1) then
+        if (ncl_version .GE. 6 ) then
+           png_output(0) = file_out+".png"
+        else
+           png_output(0) = file_out+".00000"+1+".png"
+        end if
+        ;using imagemagick's convert for reducing the white 
+        ;space around the plot
+        cmd = "convert -geometry 1000x1000 -density 300 -trim " +  \
+              png_output(0) + " " + png_output(0)
        system(cmd)
-     end do
+     else
+
+       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
+     end if
 
      print(" ")