Changeset 183 for palm/trunk/SCRIPTS/NCL

Timestamp:

Aug 4, 2008 3:39:12 PM (16 years ago)

Author:

letzel

Message:

• mrun: (optionally) link local input
• diffusion_s: bugfix for calculation of horizontal fluxes at vertical walls
• NCL scripts in trunk/SCRIPTS/NCL updated

File:

• palm/trunk/SCRIPTS/NCL/spectra.ncl (modified) (5 diffs)

palm/trunk/SCRIPTS/NCL/spectra.ncl

@@ -13 +13 @@
       delete(parameter@_FillValue)
    else
-      if (isfilepresent("~/palm/current_version/trunk/SRIPTS/NCL/.ncl_preferences")) then
-         parameter = asciiread("~/palm/current_version/trunk/SRIPTS/NCL/.ncl_preferences",129,"string")
+      if (isfilepresent("~/palm/current_version/trunk/SCRIPTS/NCL/.ncl_preferences")) then
+         parameter = asciiread("~/palm/current_version/trunk/SCRIPTS/NCL/.ncl_preferences",129,"string")
          delete(parameter@_FillValue)
       else
          print(" ")
-         print("'.ncl_preferences' does not exist in '~/palm/current_version/trunk/SRIPTS/NCL/'")
+         print("'.ncl_preferences' does not exist in '~/palm/current_version/trunk/SCRIPTS/NCL/'")
          print(" ")
          exit
@@ -372 +372 @@
    if (sort .EQ. "time")
       plot  = new(dim*dimz,graphic)
-      plot_ = new(dim*dimz,graphic)
       np=dimt
       res@lgTitleString = "Time [h]"
    else
       plot  = new(dim*dimt,graphic)
-      plot_ = new(dim*dimt,graphic)
       np=dimz
       res@lgTitleString = "Height [m]"  
@@ -455 +453 @@
                res@trXMaxF = max(x_axis)
                res@gsnLeftString      = vNam(varn)
-               res@gsnRightString     = "Height: "+z(p)+"m"
+               res@gsnRightString     = "Height = "+z(p)+"m"
                if (norm .NE. 1)then
                   res@tiYAxisString      = vNam(varn)+" / "+norm
@@ -482 +480 @@
                   res@trXMaxF = max(x_axis)
                   res@gsnLeftString      = vNam(varn)
-                  res@gsnRightString     = "Time: "+legend_label(p)+"h"
+                  res@gsnRightString     = "Time = "+legend_label(p)+"h"
                   if (norm .NE. 1)then
                      res@tiYAxisString      = vNam(varn)+" / "+norm
@@ -521 +519 @@
    resP@txFontHeightF          = 0.014
 
-   do m=0,n-1
-      plot_(m)=plot(n-1-m)
-   end do
-
    if (format_out .EQ. "eps" .OR. format_out .EQ. "epsi") then
-      gsn_panel(wks,plot_,(/n,1/),resP)
+      gsn_panel(wks,plot,(/n,1/),resP)
    else    
       do i = 0,n-1, no_lines*no_columns
          if( (i+no_lines*no_columns) .gt. (n-1)) then
-            gsn_panel(wks,plot_(i:n-1),(/no_lines,no_columns/),resP)
+            gsn_panel(wks,plot(i:n-1),(/no_lines,no_columns/),resP)
          else
-           gsn_panel(wks,plot_(i:i+no_lines*no_columns-1),(/no_lines,no_columns/),resP)
+           gsn_panel(wks,plot(i:i+no_lines*no_columns-1),(/no_lines,no_columns/),resP)
          end if
       end do