Changeset 183 for palm/trunk/SCRIPTS/NCL


Ignore:
Timestamp:
Aug 4, 2008 3:39:12 PM (13 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:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SCRIPTS/NCL/spectra.ncl

    r178 r183  
    1313      delete(parameter@_FillValue)
    1414   else
    15       if (isfilepresent("~/palm/current_version/trunk/SRIPTS/NCL/.ncl_preferences")) then
    16          parameter = asciiread("~/palm/current_version/trunk/SRIPTS/NCL/.ncl_preferences",129,"string")
     15      if (isfilepresent("~/palm/current_version/trunk/SCRIPTS/NCL/.ncl_preferences")) then
     16         parameter = asciiread("~/palm/current_version/trunk/SCRIPTS/NCL/.ncl_preferences",129,"string")
    1717         delete(parameter@_FillValue)
    1818      else
    1919         print(" ")
    20          print("'.ncl_preferences' does not exist in '~/palm/current_version/trunk/SRIPTS/NCL/'")
     20         print("'.ncl_preferences' does not exist in '~/palm/current_version/trunk/SCRIPTS/NCL/'")
    2121         print(" ")
    2222         exit
     
    372372   if (sort .EQ. "time")
    373373      plot  = new(dim*dimz,graphic)
    374       plot_ = new(dim*dimz,graphic)
    375374      np=dimt
    376375      res@lgTitleString = "Time [h]"
    377376   else
    378377      plot  = new(dim*dimt,graphic)
    379       plot_ = new(dim*dimt,graphic)
    380378      np=dimz
    381379      res@lgTitleString = "Height [m]" 
     
    455453               res@trXMaxF = max(x_axis)
    456454               res@gsnLeftString      = vNam(varn)
    457                res@gsnRightString     = "Height: "+z(p)+"m"
     455               res@gsnRightString     = "Height = "+z(p)+"m"
    458456               if (norm .NE. 1)then
    459457                  res@tiYAxisString      = vNam(varn)+" / "+norm
     
    482480                  res@trXMaxF = max(x_axis)
    483481                  res@gsnLeftString      = vNam(varn)
    484                   res@gsnRightString     = "Time: "+legend_label(p)+"h"
     482                  res@gsnRightString     = "Time = "+legend_label(p)+"h"
    485483                  if (norm .NE. 1)then
    486484                     res@tiYAxisString      = vNam(varn)+" / "+norm
     
    521519   resP@txFontHeightF          = 0.014
    522520
    523    do m=0,n-1
    524       plot_(m)=plot(n-1-m)
    525    end do
    526 
    527521   if (format_out .EQ. "eps" .OR. format_out .EQ. "epsi") then
    528       gsn_panel(wks,plot_,(/n,1/),resP)
     522      gsn_panel(wks,plot,(/n,1/),resP)
    529523   else   
    530524      do i = 0,n-1, no_lines*no_columns
    531525         if( (i+no_lines*no_columns) .gt. (n-1)) then
    532             gsn_panel(wks,plot_(i:n-1),(/no_lines,no_columns/),resP)
     526            gsn_panel(wks,plot(i:n-1),(/no_lines,no_columns/),resP)
    533527         else
    534            gsn_panel(wks,plot_(i:i+no_lines*no_columns-1),(/no_lines,no_columns/),resP)
     528           gsn_panel(wks,plot(i:i+no_lines*no_columns-1),(/no_lines,no_columns/),resP)
    535529         end if
    536530      end do
Note: See TracChangeset for help on using the changeset viewer.