Changeset 267 for palm/trunk/SCRIPTS/NCL


Ignore:
Timestamp:
Mar 25, 2009 10:25:50 AM (13 years ago)
Author:
weinreis
Message:

NCL scripts profiles cross_sections ncl_preference and UserGuide? modified

Location:
palm/trunk/SCRIPTS/NCL
Files:
4 edited

Legend:

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

    r218 r267  
    11431143      end if
    11441144      if (xyc .EQ. 1)then
    1145          cs_res@tiXAxisString = "x [gridpoints]"
    1146          cs_res@tiYAxisString = "y [gridpoints]"
     1145         cs_res@tiXAxisString = "x [m]"
     1146         cs_res@tiYAxisString = "y [m]"
    11471147      end if
    11481148      if (xzc .EQ. 1)then
    1149          cs_res@tiXAxisString = "x [gridpoints]"
    1150          cs_res@tiYAxisString = "z [gridpoints]"
     1149         cs_res@tiXAxisString = "x [m]"
     1150         cs_res@tiYAxisString = "z [m]"
    11511151      end if
    11521152      if (yzc .EQ. 1)then
    1153          cs_res@tiXAxisString = "y [gridpoints]"
    1154          cs_res@tiYAxisString = "z [gridpoints]"
     1153         cs_res@tiXAxisString = "y [m]"
     1154         cs_res@tiYAxisString = "z [m]"
    11551155      end if
    11561156   end if
     
    14921492                 
    14931493                  if ( sort .eq. "time" ) then
    1494                      if ( z_d(lo) .eq. -1)then
     1494                     if ( z_d(zs+lo) .eq. -1)then
    14951495                        if (delta_z .EQ. -1) then
    14961496                           level = "-average"
    14971497                        else
    1498                            level = "=" + z_d(lo) + "m"
     1498                           level = "=" + z_d(zs+lo) + "m"
    14991499                        end if
    15001500                     else
    1501                         level = "=" + z_d(lo) + "m"
     1501                        level = "=" + z_d(zs+lo) + "m"
    15021502                     end if
    15031503                     cs_res@gsnCenterString = "t=" + decimalPlaces(t_all(li)/3600,2,True) +"h  z"+level             
     
    15191519                 
    15201520                  if ( sort .eq. "layer" ) then
    1521                      if (z_d(li) .eq. -1) then
     1521                     if (z_d(zs+li) .eq. -1) then
    15221522                        if (delta_z .EQ. -1) then
    15231523                           level = "-average"
    15241524                        else
    1525                            level = "=" + z_d(li) + "m"
     1525                           level = "=" + z_d(zs+li) + "m"
    15261526                        end if
    15271527                     else
    1528                         level = "=" + z_d(li) + "m"
     1528                        level = "=" + z_d(zs+li) + "m"
    15291529                     end if
    15301530                     cs_res@gsnCenterString = "t=" + decimalPlaces(t_all(lo)/3600,2,True) + "h  z"+ level               
     
    15551555               
    15561556                  if ( sort .eq. "time" ) then
    1557                      if ( y_d(lo) .eq. -1 ) then
     1557                     if ( y_d(ys+lo) .eq. -1 ) then
    15581558                        level = "-average"
    15591559                     else
    1560                         level = "=" + y_d(lo) + "m"
     1560                        level = "=" + y_d(ys+lo) + "m"
    15611561                     end if
    15621562                     cs_res@gsnCenterString = "t=" + decimalPlaces(t_all(li)/3600,2,True) + "h  y"+ level
     
    15781578
    15791579                  if ( sort .eq. "layer" ) then
    1580                      if ( y_d(li) .eq. -1 ) then
     1580                     if ( y_d(ys+li) .eq. -1 ) then
    15811581                        level = "-average"
    15821582                     else
    1583                         level = "=" + y_d(li) + "m"
     1583                        level = "=" + y_d(ys+li) + "m"
    15841584                     end if
    15851585                     cs_res@gsnCenterString = "t=" + decimalPlaces(t_all(lo)/3600,2,True) + "h  y"+ level
     
    16101610
    16111611                  if ( sort .eq. "time" ) then
    1612                      if ( x_d(lo) .eq. -1 ) then
     1612                     if ( x_d(xs+lo) .eq. -1 ) then
    16131613                        level = "-average"
    16141614                     else
    1615                         level = "=" + x_d(lo) + "m"
     1615                        level = "=" + x_d(xs+lo) + "m"
    16161616                     end if
    16171617                     cs_res@gsnCenterString = "t=" + decimalPlaces(t_all(li)/3600,2,True) + "h  x"+ level
     
    16331633
    16341634                  if ( sort .eq. "layer" ) then
    1635                      if ( x_d(li) .eq. -1 ) then
     1635                     if ( x_d(xs+li) .eq. -1 ) then
    16361636                        level = "-average"
    16371637                     else
    1638                         level = "=" + x_d(li) + "m"
     1638                        level = "=" + x_d(xs+li) + "m"
    16391639                     end if
    16401640                     cs_res@gsnCenterString = "t=" + decimalPlaces(t_all(lo)/3600,2,True) + "h  x"+ level
     
    16801680         end do
    16811681      end if
    1682    else        
     1682   else       
    16831683      if (format_out .EQ. "eps" .OR. format_out .EQ. "epsi") then
    16841684         gsn_panel(wks_ps,plot(0:(no_time*no_layer*no_var)-1),(/no_var,no_layer*no_time/),cs_resP)
     
    16871687         print(" ")
    16881688      else
    1689          do np = 0,no_layer*no_time*no_var-1,no_rows*no_columns   
     1689         do np = 0,no_layer*no_time*no_var-1,no_rows*no_columns 
    16901690            if ( np + no_rows*no_columns .gt. (no_layer*no_time*no_var)-1) then
    16911691               gsn_panel(wks_ps, plot(np:(no_layer*no_time*no_var)-1),(/no_rows,no_columns/),cs_resP)
  • palm/trunk/SCRIPTS/NCL/ncl_preferences.ncl

    r236 r267  
    3030
    3131
    32          file_1 = "File in" 
     32         file_1 = "File in"
    3333         
    3434         
     
    677677
    678678
    679          file_1 = "File in" 
     679         file_1 = "File in"
    680680         
    681681         
     
    717717         
    718718      end if
    719       ;***************************************************
    720       ; REQUIRED --prof3d-- 3D OR PROFILE DATA
    721       ;
    722       ; must be switched on [=1] if 3D data are used or switched off [=0] if profile data are used for profiles
    723       ;
    724       ; data type: integer
    725       ;
    726       ; default:   -1
    727       ;***************************************************
    728       if (.not. isvar("prof3d"))then 
    729      
    730          
    731          prof3d = -1
    732          
    733          
    734       end if     
    735719      ;***************************************************
    736720      ; OPTIONAL --format_out-- OUTPUT FORMAT
     
    12671251     
    12681252         
    1269          end_time_step = -1.d     
     1253         end_time_step = -1.d   
    12701254         
    12711255         
     
    20382022
    20392023
    2040          file_1 = "File in" 
     2024         file_1 = "File in"
    20412025         
    20422026         
     
    21462130     
    21472131         
    2148          no_rows = 2
     2132         no_rows = 1
    21492133         
    21502134         
     
    22042188         end_time_step = -1.d
    22052189         
    2206          
     2190
    22072191      end if     
    22082192      ;***************************************************
  • palm/trunk/SCRIPTS/NCL/profiles.ncl

    r250 r267  
    7373   start_f(0) = start_f_1
    7474   end_f(0) = end_f_1   
    75      
    76    if (prof3d .EQ. -1)then
    77       print(" ")
    78       print("Set 'prof3d' to 0 or 1")
    79       print(" ")
    80       exit
    81    else
    82       if (prof3d .NE. 0 .AND. prof3d .NE. 1)then
    83          print(" ")
    84          print("'prof3d'= "+prof3d+" is invalid; set 'prof3d' to 0 or 1")
    85          print(" ")
    86          exit
    87       end if
    88    end if
    8975 
    9076   if (no_files .GT. 1) then
     
    372358   
    373359   vNam = getfilevarnames(f_att)
    374 
     360   
    375361   if (nof .EQ. 0)then
    376362      vNam0=vNam
     
    407393      print(" ")
    408394   end if
     395   
     396   prof3d = 0
     397   do varn = dim-1,0,1
     398      if ( isStrSubset( vNam(varn), "zu_3d") .OR. isStrSubset( vNam(varn), "zw_3d")) then
     399         prof3d = 1
     400         break
     401      end if
     402   end do
    409403
    410404   ;-----------------------------------------------------below steps only for first file -> nof=0
Note: See TracChangeset for help on using the changeset viewer.