Ignore:
Timestamp:
Dec 10, 2008 9:14:34 AM (15 years ago)
Author:
letzel
Message:
  • User can check user parameters and deduce further quantities in user_check_parameters
  • Topography definition according to new user parameter topography_grid_convention (init_grid, modules, user_header, user_init_grid, user_parin)
  • New subdirectory trunk/EXAMPLES with two examples highlighting the topography_grid_convention
File:
1 edited

Legend:

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

    r194 r218  
    1616      ; REQUIRED --file_1-- INPUT FILE
    1717      ;
    18       ; input file produced by PALM; please declare path and name
    19       ;
    20       ; data type: string
     18      ; input file produced by PALM; if there is more than one file for all timesteps of a job chain,
     19      ; declare path and file name without cycle number and ".nc" (first and last cycle numbers are given with the parameters 'start_f' and 'end_f');
     20      ; otherwise declare path and full file name
     21      ;
     22      ; data type: string
     23      ;
     24      ; example: file_1 = "/path_to_file/file_xy_av"      if you use more than one cycle number
     25      ;          file_1 = "/path_to_file/file_xy_av.4.nc" if you use one file (e.g. with cycle number 4)
    2126      ;
    2227      ; default:   "File in"
     
    2530
    2631
    27          file_1 = "File in"
    28          
    29          
    30       end if
    31       ;***************************************************
    32       ; REQUIRED --xyc, xzc, yzc-- SECTION of XY or XZ or YZ
     32         file_1 = "file_1"
     33         
     34         
     35      end if
     36      ;***************************************************
     37      ; REQUIRED if more than one file of a job chain --start_f-- START OF CYCLIC NUMBER
     38      ;
     39      ; if there is more than one file for all timesteps of a job chain, declare the first cyclicnumber which contains timesteps you'd like to use
     40      ;
     41      ; data type: integer
     42      ;
     43      ; example:   start_f = 3 (if first file: example_ts.3.nc)
     44      ;
     45      ; default:   -1 
     46      ;***************************************************
     47      if (.not. isvar("start_f"))then 
     48     
     49         
     50         start_f = -1
     51         
     52         
     53      end if
     54      ;***************************************************
     55      ; REQUIRED if more than one file of a job chain --end_f-- END OF CYCLIC NUMBER
     56      ;
     57      ; if there is more than one file for all timesteps of a job chain, declare the last cyclicnumber which contains timesteps you'd like to use
     58      ;
     59      ; data type: integer
     60      ;
     61      ; example:   end_f = 48 (if last file: example_ts.48.nc)
     62      ;
     63      ; default:   -1 
     64      ;***************************************************
     65      if (.not. isvar("end_f"))then 
     66     
     67         
     68         end_f = -1
     69         
     70         
     71      end if
     72      ;***************************************************
     73      ; REQUIRED --xyc--xzc--yzc-- SELECTION of XY or XZ or YZ CROSS SECTIONS
    3374      ;
    3475      ; desired section has to be set to 1; e.g.: xyc=1, xzc=0 and yzc=0 for xy-section
    3576      ;
    3677      ; data type: integer
     78      ;
     79      ; example: xyc = 1
     80      ;          xzc = 0
     81      ;          yzc = 0
    3782      ;
    3883      ; default:   0
     
    56101      ; OPTIONAL --format_out-- OUTPUT FORMAT
    57102      ;
    58       ; supported file formats: x11, pdf, ps, eps, epsi or ncgm
    59       ;
    60       ; data type: string
     103      ; format of the output file; supported file formats: x11, pdf, ps, eps, epsi or ncgm
     104      ;
     105      ; data type: string
     106      ;
     107      ; example: format_out = "pdf"
    61108      ;
    62109      ; default:   "x11"
     
    72119      ; OPTIONAL --file_out-- OUTPUT FILE
    73120      ;
    74       ; please declare path and name
    75       ;
    76       ; data type: string
     121      ; path and name of the output file
     122      ;
     123      ; data type: string
     124      ;
     125      ; example: file_out = "/home/path_to_output/test_cs"
    77126      ;
    78127      ; default:   "~/test_cs"
     
    86135      end if     
    87136      ;***************************************************
    88       ; OPTIONAL --no_columns-- NUMBER OF ROWS
     137      ; OPTIONAL --no_columns-- NUMBER OF COLUMNS
    89138      ;
    90139      ; number of plots in one row
    91140      ;
    92141      ; data type: integer
     142      ;
     143      ; example: no_columns = 1
    93144      ;
    94145      ; default:   1
     
    102153      end if     
    103154      ;***************************************************
    104       ; OPTIONAL --no_lines-- NUMBER OF COLUMNS
     155      ; OPTIONAL --no_rows-- NUMBER OF ROWS
    105156      ;
    106157      ; number of plots in one column
    107158      ;
    108159      ; data type: integer
     160      ;
     161      ; example: no_rows = 2
    109162      ;
    110163      ; default:   2
    111164      ;***************************************************
    112       if(.not. isvar("no_lines"))then   
    113      
    114      
    115          no_lines = 2
     165      if(.not. isvar("no_rows"))then     
     166     
     167     
     168         no_rows = 2
    116169         
    117170         
     
    124177      ; data type: string
    125178      ;
     179      ; example: sort = "layer"
     180      ;
    126181      ; default:   "layer"
    127182      ;***************************************************
     
    136191      ; OPTIONAL --var-- OUTPUT VARIABLES
    137192      ;
    138       ; names of output variables; please inidcate them with one comma before and after each variable;
    139       ; e.g.: var=",u,v,w," (the windcomponents u, v and w will be plotted)
     193      ; names of output variables; inidcate them with one comma before and after each variable
    140194      ;
    141195      ; data type: string
    142       ;
    143       ; default:   "all" (all available variables will be plotted)
     196      ;
     197      ; example: var = ",u,v,w"
     198      ;
     199      ; default:   "all" (all available variables are plotted)
    144200      ;***************************************************
    145201      if(.not. isvar("var"))then
     
    153209      ; OPTIONAL --start_time_step-- FIRST TIME STEP
    154210      ;
    155       ; first time step of plot; you don't have to specify a precise time value from the input file,
    156       ; it will be rounded to the next existent time step
     211      ; first time step of plot; there is no need to specify a precise time value from the input file,
     212      ; 'start_time_step' is rounded to the next existent time step
    157213      ;
    158214      ; data type: double; dimension: [hour]
    159215      ;
     216      ; example:   start_time_step = 0.5
     217      ;
    160218      ; default:   -1.d (first available time step)
    161219      ;***************************************************
     
    170228      ; OPTIONAL --end_time_step-- LAST TIME STEP
    171229      ;
    172       ; last time step of plot; you don't have to specify a precise time value from the input file,
    173       ; it will be rounded to the next existent time step
     230      ; last time step of plot; there is no need to specify a precise time value from the input file,
     231      ; 'end_time_step' is rounded to the next existent time step
    174232      ;
    175233      ; data type: double; dimension: [hour]
    176234      ;
     235      ; example:   end_time_step = 2.0
     236      ;
    177237      ; default:   -1.d (last available time step)
    178238      ;***************************************************
     
    188248      ;
    189249      ; value for minimum of x-axis, selectable if there are NO preseted layers for x (as for XY or XZ 2D-DATA);
    190       ; you don't have to specify a precise meter value from the input file, it will be rounded to the next
    191       ; existent value
     250      ; there is no need to specify a precise meter value from the input file, it is rounded to the next existent value
    192251      ;
    193252      ; data type: double; dimension: [meter]
    194253      ;
     254      ; example: xs = 150.0
     255      ;
    195256      ; default:   -1.d (minimum of range)
    196257      ;***************************************************
     
    206267      ;
    207268      ; value for maximum of x-axis, selectable if there are NO preseted layers for x (as for XY or XZ 2D-DATA);
    208       ; you don't have to specify a precise meter value from the input file, it will be rounded to the next
    209       ; existent value
     269      ; there is no need to specify a precise meter value from the input file, it is rounded to the next existent value
    210270      ;
    211271      ; data type: double; dimension: [meter]
    212272      ;
     273      ; example: xe = 300.0
     274      ;
    213275      ; default:   -1.d (maximum of range)
    214276      ;***************************************************
     
    224286      ;
    225287      ; value for minimum of y-axis, selectable if there are NO preseted layers for y (as for XY or YZ 2D-DATA);
    226       ; you don't have to specify a precise meter value from the input file, it will be rounded to the next
    227       ; existent value
     288      ; there is no need to specify a precise meter value from the input file, it is rounded to the next existent value
    228289      ;
    229290      ; data type: double; dimension: [meter]
    230291      ;
     292      ; example: ys = 150.0
     293      ;
    231294      ; default:   -1.d (minimum of range)
    232295      ;***************************************************
     
    242305      ;
    243306      ; value for maximum of y-axis, selectable if there are NO preseted layers for y (as for XY or YZ 2D-DATA);
    244       ; you don't have to specify a precise meter value from the input file, it will be rounded to the next
    245       ; existent value
     307      ; there is no need to specify a precise meter value from the input file, it is rounded to the next existent value
    246308      ;
    247309      ; data type: double; dimension: [meter]
     310      ;
     311      ; example: ye = 300.0
    248312      ;
    249313      ; default:   -1.d (maximum of range)
     
    264328      ; data type: integer
    265329      ;
     330      ; example: zs = 0
     331      ;
    266332      ; default:   -1 (minimum of range)
    267333      ;***************************************************
     
    281347      ; data type: integer
    282348      ;
     349      ; example: ze = 15
     350      ;
    283351      ; default:   -1 (maximum of range)
    284352      ;***************************************************
     
    297365      ; data type: string
    298366      ;
     367      ; example: mode = "Fill"
     368      ;
    299369      ; default:   "Fill"
    300370      ;***************************************************
     
    313383      ; data type: string
    314384      ;
     385      ; example: fill_mode = "AreaFill"
     386      ;
    315387      ; default:   "AreaFill"
    316388      ;***************************************************
     
    325397      ; OPTIONAL --shape-- ASPECT RATIO
    326398      ;
    327       ; aspect ratio of axis is kept [=1] or  not[=0]
    328       ;
    329       ; data type: integer
     399      ; aspect ratio of axis is kept [=1] or not[=0]
     400      ;
     401      ; data type: integer
     402      ;
     403      ; example: shape = 1
    330404      ;
    331405      ; default:   1
     
    337411         
    338412         
     413      end if
     414      ;***************************************************
     415      ; OPTIONAL --font_size-- FONT SIZE OF STRINGS
     416      ;
     417      ; font size can be changed to fit to the current plot
     418      ;
     419      ; data type: float
     420      ;
     421      ; example: font_size = 0.04
     422      ;
     423      ; default:   0.04
     424      ;***************************************************
     425      if(.not. isvar("font_size"))then
     426     
     427         
     428         font_size = 0.04
     429         
     430         
     431      end if
     432      ;***************************************************
     433      ; OPTIONAL --font_size_legend-- FONT SIZE OF LEGEND STRINGS
     434      ;
     435      ; font size of the legend can be changed to fit to the current plot
     436      ;
     437      ; data type: float
     438      ;
     439      ; example: font_size_legend = 0.02
     440      ;
     441      ; default:   0.02
     442      ;***************************************************
     443      if(.not. isvar("font_size_legend"))then
     444     
     445         
     446         font_size_legend = 0.02
     447         
     448         
     449      end if   
     450      ;***************************************************
     451      ; OPTIONAL --legend_label_stride-- REDUCTION OF THE NUMBER OF LABELS IN THE LEGEND
     452      ;
     453      ; if legend_label_stride is set to 1 every label is shown; if set to 2 every second,...
     454      ;
     455      ; data type: integer
     456      ;
     457      ; example: legend_label_stride = 2
     458      ;
     459      ; default:   2
     460      ;***************************************************
     461      if(.not. isvar("legend_label_reduction"))then
     462     
     463         
     464         legend_label_stride = 2
     465         
     466         
     467      end if
     468      ;***************************************************
     469      ; OPTIONAL --axes_explicit-- EXPLICIT SETTINGS OF X- AND Y-AXES ARE SWITCHED ON
     470      ;
     471      ; to use several parameters to set up the axes, 'axes_explicit' must be switched [=1]; if switched off [=0] units of axes are [gridpoints]
     472      ;
     473      ; data type: integer
     474      ;
     475      ; example: axes_explicit = 0
     476      ;
     477      ; default:   0
     478      ;***************************************************
     479      if(.not. isvar("axes_explicit"))then
     480     
     481         
     482         axes_explicit = 0
     483         
     484         
     485      end if
     486      ;***************************************************
     487      ; OPTIONAL if axes_explicit = 1 --major_ticks_x-- NUMBER OF MAJOR TICK MARKS AT X-AXIS
     488      ;
     489      ; major_ticks_x can be set to any integer, best segmentation of x-axis can be achieved, if
     490      ; (number of gripdpoints in x-direction)/(major_ticks_x-1)=INTEGER, e.g. 128/8 = 16 -> major_ticks_x should be set to 9
     491      ;
     492      ; data type: integer
     493      ;
     494      ; example: major_ticks_x = 6
     495      ;
     496      ; default:   6
     497      ;***************************************************
     498      if(.not. isvar("major_ticks_x"))then
     499     
     500         
     501         major_ticks_x = 6
     502         
     503         
     504      end if
     505      ;***************************************************
     506      ; OPTIONAL if axes_explicit = 1 --major_ticks_y-- NUMBER OF MAJOR TICK MARKS AT Y-AXIS
     507      ;
     508      ; major_ticks_y can be set to any integer, best segmentation of y-axis can be achieved, if
     509      ; (number of gripdpoints in y-direction)/(major_ticks_y-1)=INTEGER, e.g. 128/8 = 16 -> major_ticks_y should be set to 9
     510      ;
     511      ; data type: integer
     512      ;
     513      ; example: major_ticks_y = 6
     514      ;
     515      ; default:   6
     516      ;***************************************************
     517      if(.not. isvar("major_ticks_y"))then
     518     
     519         
     520         major_ticks_y = 6
     521         
     522         
     523      end if
     524      ;***************************************************
     525      ; OPTIONAL if axes_explicit = 1 --norm_x--norm_y--norm_z-- NORMALISE X-, Y- and Z-AXES
     526      ;
     527      ; normalising the axes with the given value; units can explicitly set; by default no normalising is used and units are [meter]
     528      ;
     529      ; data type: float
     530      ;
     531      ; default:   1.
     532      ;***************************************************
     533      if(.not. isvar("norm_x"))then     
     534         
     535         norm_x = 1.     
     536         
     537      end if
     538      if(.not. isvar("norm_y"))then     
     539         
     540         norm_y = 1.     
     541         
     542      end if
     543      if(.not. isvar("norm_z"))then     
     544         
     545         norm_z = 1.     
     546         
     547      end if
     548      ;***************************************************
     549      ; OPTIONAL if axes_explicit = 1 --unit_x--unit_y--unit_z-- UNITS OF X-, Y- and Z-AXES
     550      ;
     551      ; units of the axes can be set; by default no units are given, if norm_x, norm_y or norm_z are unequal 1., otherwise units are [meter]
     552      ;
     553      ; data type: string,
     554      ;
     555      ; example: unit_x = "m/s"
     556      ;
     557      ; default:   " "
     558      ;***************************************************
     559      if(.not. isvar("unit_x"))then     
     560         
     561         unit_x = " "   
     562         
     563      end if
     564      if(.not. isvar("unit_y"))then     
     565         
     566         unit_y = " "   
     567         
     568      end if
     569      if(.not. isvar("unit_z"))then     
     570         
     571         unit_z = " "   
     572         
     573      end if
     574      ;***************************************************
     575      ; OPTIONAL --vector-- VECTOR PLOT
     576      ;
     577      ; a vector plot for one defined vector is switched on [=1] or off [=0]
     578      ;
     579      ; data type: integer
     580      ;
     581      ; example: vector = 0
     582      ;
     583      ; default:   0
     584      ;***************************************************
     585      if(.not. isvar("vector"))then
     586     
     587         
     588         vector = 0
     589         
     590         
    339591      end if     
    340592      ;***************************************************
    341       ; OPTIONAL --vector-- VECTOR PLOT
    342       ;
    343       ; a vector plot for one defined vector will be switched on [=1] or off [=0]
    344       ;
    345       ; data type: integer
    346       ;
    347       ; default:   0
    348       ;***************************************************
    349       if(.not. isvar("vector"))then
    350      
    351          
    352          vector = 0
    353          
     593      ; REQUIRED IF vector=1 --vec1, vec2-- VECTOR FOR VECTOR PLOT
     594      ;
     595      ; variables for the vector, vec1 hast to be set to component 1 and vec2 to component 2;
     596      ; inidcate them with one comma before and after each variable;
     597      ;
     598      ; data type: string
     599      ;
     600      ; example: vec1=",u,"
     601      ;          vec2=",v,"
     602      ;
     603      ; default:   vec1 = "component1", vec2 = "component2"
     604      ;***************************************************
     605      if(.not. isvar("vec1"))then     
     606         
     607         vec1 = "component1"
     608         
     609      end if
     610      if (.not. isvar("vec2"))then
     611     
     612         vec2 = "component2"     
    354613         
    355614      end if     
    356615      ;***************************************************
    357       ; REQUIRED IF vector=1 --vec1, vec2-- VECTOR FOR VECTOR PLOT
    358       ;
    359       ; variables for the vector, vec1 hast to be set to component 1 and vec2 to component 2;
    360       ; please inidcate them with one comma before and after each variable;
    361       ; e.g.: vec1=",u," and vec2=",v," for plotting the horizontal windvector
    362       ;
    363       ; data type: string
    364       ;
    365       ; default:   vec1 = "component1", vec2 = "component2"
    366       ;***************************************************
    367       if(.not. isvar("vec1"))then     
    368          
    369          vec1 = "component1"
    370          
    371       end if
    372       if (.not. isvar("vec2"))then
    373      
    374          vec2 = "component2"     
     616      ; OPTIONAL --plotvec-- VARIABLE FOR EXTRA VECTOR PLOT
     617      ;
     618      ; the vector plot can be combined with other plots of variables; if this is desired inidcate the variables,
     619      ; where a vector plot shall overlay; inidcate them with one comma before and after each variable
     620      ;
     621      ; data type: string
     622      ;
     623      ; example: plotvec = ",u,w,"
     624      ;
     625      ; default:   "plotvec" (no combining, the vector plot is created seperately)
     626      ;***************************************************
     627      if(.not. isvar("plotvec"))then
     628     
     629         
     630         plotvec = "plotvec"
     631         
    375632         
    376633      end if     
    377634      ;***************************************************
    378       ; OPTIONAL --plotvec-- VARIABLE FOR EXTRA VECTOR PLOT
    379       ;
    380       ; the vector plot can be combined with other variables; if desired inidcate the variables
    381       ; with one comma before and after it; e.g. plotvec=",u,v,w," (the vector plot
    382       ; will combined with the contour plots of u, v and w)
    383       ;
    384       ; data type: string
    385       ;
    386       ; default:   plotvec (no combining, the vector plot will be created seperately)
    387       ;***************************************************
    388       if(.not. isvar("plotvec"))then
    389      
    390          
    391          plotvec = "plotvec"
    392          
    393          
    394       end if     
    395       ;***************************************************
    396635      ; OPTIONAL --ref_mag-- REFERENCE VECTOR
    397636      ;
     
    399638      ;
    400639      ; data type: float
     640      ;
     641      ; example: ref_mag = 0.05
    401642      ;
    402643      ; default:   0.05
     
    420661
    421662      ;***************************************************
    422       ; REQUIRED --file_1-- INPUT FILE
    423       ;
    424       ; input file produced by PALM; please declare path and name
    425       ;
    426       ; data type: string
     663      ; REQUIRED --file_1-- 1ST INPUT FILE
     664      ;
     665      ; 1st input file produced by PALM; if there is more than one file for all timesteps of a job chain,
     666      ; declare path and file name without cycle number and ".nc" (first and last cycle numbers are given with the parameters 'start_f' and 'end_f');
     667      ; otherwise declare path and full file name
     668      ;
     669      ; data type: string
     670      ;
     671      ; example: file_1 = "/path_to_file/file_pr"      if you use more than one cycle number
     672      ;          file_1 = "/path_to_file/file_pr.4.nc" if you use one file (e.g. with cycle number 4)
    427673      ;
    428674      ; default:   "File in"
    429675      ;***************************************************
    430       if (.not. isvar("file_1"))then
    431      
    432      
    433          file_1 = "File in"
     676      if(.not. isvar("file_1"))then
     677
     678
     679         file_1 = "File in"
     680         
     681         
     682      end if
     683      ;***************************************************
     684      ; REQUIRED --start_f_1-- START OF CYCLIC NUMBER OF THE 1ST INPUT FILE
     685      ;
     686      ; if there is more than one file for all timesteps of a job chain, declare the first cyclicnumber which contains timesteps you'd like to use
     687      ;
     688      ; data type: integer
     689      ;
     690      ; example:   3 (if first file: example_ts.3.nc)
     691      ;
     692      ; default:   -1 
     693      ;***************************************************
     694      if (.not. isvar("start_f_1"))then 
     695     
     696         
     697         start_f_1 = -1
     698         
     699         
     700      end if
     701      ;***************************************************
     702      ; REQUIRED --end_f_1-- END OF CYCLIC NUMBER OF THE 1ST INPUT FILE
     703      ;
     704      ; if there is more than one file for all timesteps of a job chain, declare the last cyclicnumber which contains timesteps you'd like to use
     705      ;
     706      ; data type: integer
     707      ;
     708      ; example:   48 (if last file: example_ts.48.nc)
     709      ;
     710      ; default:   -1 
     711      ;***************************************************
     712      if (.not. isvar("end_f_1"))then 
     713     
     714         
     715         end_f_1 = -1
    434716         
    435717         
     
    470752      ; OPTIONAL --file_out-- OUTPUT FILE
    471753      ;
    472       ; please declare path and name
     754      ; declare path and name of the output file
    473755      ;
    474756      ; data type: string
     
    484766      end if     
    485767      ;***************************************************
    486       ; OPTIONAL --no_columns-- NUMBER OF ROWS
     768      ; OPTIONAL --no_columns-- NUMBER OF COLUMNS
    487769      ;
    488770      ; number of plots in one row
     
    500782      end if     
    501783      ;***************************************************
    502       ; OPTIONAL --no_lines-- NUMBER OF COLUMNS
     784      ; OPTIONAL --no_rows-- NUMBER OF ROWS
    503785      ;
    504786      ; number of plots in one column
     
    508790      ; default:   2
    509791      ;***************************************************
    510       if (.not. isvar("no_lines"))then 
     792      if (.not. isvar("no_rows"))then   
    511793     
    512794       
    513          no_lines = 2
     795         no_rows = 2
    514796         
    515797         
     
    518800      ; OPTIONAL --var-- OUTPUT VARIABLES
    519801      ;
    520       ; names of output variables; please inidcate them with one comma before and after each variable;
    521       ; e.g.: var=",u,v,w," (the windcomponents u, v and w will be plotted)
     802      ; names of output variables; inidcate them with one comma before and after each variable
    522803      ;
    523804      ; data type: string
    524       ;
    525       ; default:   "all" (all available variables will be plotted)
     805      ;
     806      ; example: var = ",u,v,w";
     807      ;
     808      ; default:   "all" (all available variables are plotted)
    526809      ;***************************************************
    527810      if (.not. isvar("var"))then
     
    535818      ; OPTIONAL ---no_files-- NUMBER OF INPUT FILES
    536819      ;
    537       ; up to 12 input files can be used; all input files must contain the same variables and dimensionalties;
    538       ; identical variables will be merged in one plot;
     820      ; up to 6 input files of different runs can be used; all input files must contain the same variables and dimensionalties;
     821      ; identical variables are merged in one plot;
    539822      ; combining or overlaying of any variables is not feasible;
    540       ; an extra legend (besides the time step legend) will be drawn where you can explicitily name the items
     823      ; an extra legend (besides the time step legend) is drawn where you can explicitily name the items
    541824      ;
    542825      ; data type: integer
     
    566849         
    567850         
    568       end if     
    569       ;***************************************************
    570       ; REQUIRED IF no_files > 1 --file_2-- 2nd INPUT FILE
    571       ;
    572       ; 2nd input file produced by PALM; please declare path and name
    573       ;
    574       ; data type: string
     851      end if
     852      ;***************************************************
     853      ; REQUIRED if no_files > 1 --file_2-- 2ND INPUT FILE
     854      ;
     855      ; 2nd input file produced by PALM; if there is more than one file for all timesteps of a job chain,
     856      ; declare path and file name without cycle number and ".nc" (first and last cycle numbers are given with the parameters 'start_f' and 'end_f');
     857      ; otherwise declare path and full file name
     858      ;
     859      ; data type: string
     860      ;
     861      ; example: file_2 = "/path_to_file/file_pr"      if you use more than one cycle number
     862      ;          file_2 = "/path_to_file/file_pr.4.nc" if you use one file (e.g. with cycle number 4)
    575863      ;
    576864      ; default:   "File in"
    577865      ;***************************************************
    578       if (.not. isvar("file_2"))then 
     866      if(.not. isvar("file_2"))then
     867
     868
     869         file_2 = "File in"
     870         
     871         
     872      end if
     873      ;***************************************************
     874      ; REQUIRED --start_f_2-- START OF CYCLIC NUMBER OF THE 2ND INPUT FILE
     875      ;
     876      ; if there is more than one file for all timesteps of a job chain, declare the first cyclicnumber which contains timesteps you'd like to use
     877      ;
     878      ; data type: integer
     879      ;
     880      ; example:   3 (if first file: example_ts.3.nc)
     881      ;
     882      ; default:   -1 
     883      ;***************************************************
     884      if (.not. isvar("start_f_2"))then 
    579885     
    580886         
    581          file_2 = "File in"
    582          
    583          
    584       end if     
     887         start_f_2 = -1
     888         
     889         
     890      end if
     891      ;***************************************************
     892      ; REQUIRED --end_f_2-- END OF CYCLIC NUMBER OF THE 2ND INPUT FILE
     893      ;
     894      ; if there is more than one file for all timesteps of a job chain, declare the last cyclicnumber which contains timesteps you'd like to use
     895      ;
     896      ; data type: integer
     897      ;
     898      ; example:   48 (if last file: example_ts.48.nc)
     899      ;
     900      ; default:   -1 
     901      ;***************************************************
     902      if (.not. isvar("end_f_2"))then 
     903     
     904         
     905         end_f_2 = -1
     906         
     907         
     908      end if             
    585909      ;***************************************************
    586910      ; OPTIONAL --name_legend_2-- 2nd LEGEND ITEM
     
    598922         
    599923         
    600       end if     
    601       ;***************************************************
    602       ; REQUIRED IF no_files > 2 --file_3-- 3rd INPUT FILE
    603       ;
    604       ; 3rd input file produced by PALM; please declare path and name
    605       ;
    606       ; data type: string
     924      end if
     925      ;***************************************************
     926      ; REQUIRED if no_files > 2 --file_3-- 3RD INPUT FILE
     927      ;
     928      ; 3rd input file produced by PALM; if there is more than one file for all timesteps of a job chain,
     929      ; declare path and file name without cycle number and ".nc" (first and last cycle numbers are given with the parameters 'start_f' and 'end_f');
     930      ; otherwise declare path and full file name
     931      ;
     932      ; data type: string
     933      ;
     934      ; example: file_3 = "/path_to_file/file_pr"      if you use more than one cycle number
     935      ;          file_3 = "/path_to_file/file_pr.4.nc" if you use one file (e.g. with cycle number 4)
    607936      ;
    608937      ; default:   "File in"
    609938      ;***************************************************
    610       if (.not. isvar("file_3"))then   
    611      
    612        
    613          file_3 = "File in"
    614          
    615          
    616       end if     
     939      if(.not. isvar("file_3"))then
     940
     941
     942         file_3 = "File in"
     943         
     944         
     945      end if
     946      ;***************************************************
     947      ; REQUIRED --start_f_3-- START OF CYCLIC NUMBER OF THE 3RD INPUT FILE
     948      ;
     949      ; if there is more than one file for all timesteps of a job chain, declare the first cyclicnumber which contains timesteps you'd like to use
     950      ;
     951      ; data type: integer
     952      ;
     953      ; example:   3 (if first file: example_ts.3.nc)
     954      ;
     955      ; default:   -1 
     956      ;***************************************************
     957      if (.not. isvar("start_f_3"))then 
     958     
     959         
     960         start_f_3 = -1
     961         
     962         
     963      end if
     964      ;***************************************************
     965      ; REQUIRED --end_f_3-- END OF CYCLIC NUMBER OF THE 3RD INPUT FILE
     966      ;
     967      ; if there is more than one file for all timesteps of a job chain, declare the last cyclicnumber which contains timesteps you'd like to use
     968      ;
     969      ; data type: integer
     970      ;
     971      ; example:   48 (if last file: example_ts.48.nc)
     972      ;
     973      ; default:   -1 
     974      ;***************************************************
     975      if (.not. isvar("end_f_3"))then 
     976     
     977         
     978         end_f_3 = -1
     979         
     980         
     981      end if   
    617982      ;***************************************************
    618983      ; OPTIONAL --name_legend_3-- 3rd LEGEND ITEM
     
    630995         
    631996         
    632       end if     
    633       ;***************************************************
    634       ; REQUIRED IF no_files > 3 --file_4-- 4th INPUT FILE
    635       ;
    636       ; 4th input file produced by PALM; please declare path and name
    637       ;
    638       ; data type: string
     997      end if
     998      ;***************************************************
     999      ; REQUIRED if no_files > 3 --file_4-- 4TH INPUT FILE
     1000      ;
     1001      ; 4th input file produced by PALM; if there is more than one file for all timesteps of a job chain,
     1002      ; declare path and file name without cycle number and ".nc" (first and last cycle numbers are given with the parameters 'start_f' and 'end_f');
     1003      ; otherwise declare path and full file name
     1004      ;
     1005      ; data type: string
     1006      ;
     1007      ; example: file_4 = "/path_to_file/file_pr"      if you use more than one cycle number
     1008      ;          file_4 = "/path_to_file/file_pr.4.nc" if you use one file (e.g. with cycle number 4)
    6391009      ;
    6401010      ; default:   "File in"
    6411011      ;***************************************************
    642       if (.not. isvar("file_4"))then
    643      
    644            
    645          file_4 = "File in"
    646          
    647          
    648       end if     
     1012      if(.not. isvar("file_4"))then
     1013
     1014
     1015         file_4 = "File in"
     1016         
     1017         
     1018      end if
     1019      ;***************************************************
     1020      ; REQUIRED --start_f_4-- START OF CYCLIC NUMBER OF THE 4TH INPUT FILE
     1021      ;
     1022      ; if there is more than one file for all timesteps of a job chain, declare the first cyclicnumber which contains timesteps you'd like to use
     1023      ;
     1024      ; data type: integer
     1025      ;
     1026      ; example:   3 (if first file: example_ts.3.nc)
     1027      ;
     1028      ; default:   -1 
     1029      ;***************************************************
     1030      if (.not. isvar("start_f_4"))then 
     1031     
     1032         
     1033         start_f_4 = -1
     1034         
     1035         
     1036      end if
     1037      ;***************************************************
     1038      ; REQUIRED --end_f_4-- END OF CYCLIC NUMBER OF THE 4TH INPUT FILE
     1039      ;
     1040      ; if there is more than one file for all timesteps of a job chain, declare the last cyclicnumber which contains timesteps you'd like to use
     1041      ;
     1042      ; data type: integer
     1043      ;
     1044      ; example:   48 (if last file: example_ts.48.nc)
     1045      ;
     1046      ; default:   -1 
     1047      ;***************************************************
     1048      if (.not. isvar("end_f_4"))then 
     1049     
     1050         
     1051         end_f_4 = -1
     1052         
     1053         
     1054      end if 
    6491055      ;***************************************************
    6501056      ; OPTIONAL --name_legend_4-- 4th LEGEND ITEM
     
    6621068         
    6631069         
    664       end if     
    665       ;***************************************************
    666       ; REQUIRED IF no_files > 4 --file_5-- 5th INPUT FILE
    667       ;
    668       ; 5th input file produced by PALM; please declare path and name
    669       ;
    670       ; data type: string
     1070      end if
     1071      ;***************************************************
     1072      ; REQUIRED if no_files > 4 --file_5-- 5TH INPUT FILE
     1073      ;
     1074      ; 5th input file produced by PALM; if there is more than one file for all timesteps of a job chain,
     1075      ; declare path and file name without cycle number and ".nc" (first and last cycle numbers are given with the parameters 'start_f' and 'end_f');
     1076      ; otherwise declare path and full file name
     1077      ;
     1078      ; data type: string
     1079      ;
     1080      ; example: file_5 = "/path_to_file/file_pr"      if you use more than one cycle number
     1081      ;          file_5 = "/path_to_file/file_pr.4.nc" if you use one file (e.g. with cycle number 4)
    6711082      ;
    6721083      ; default:   "File in"
    6731084      ;***************************************************
    674       if (.not. isvar("file_5"))then
    675      
    676            
    677          file_5 = "File in"
    678          
    679          
    680       end if     
     1085      if(.not. isvar("file_5"))then
     1086
     1087
     1088         file_5 = "File in"
     1089         
     1090         
     1091      end if
     1092      ;***************************************************
     1093      ; REQUIRED --start_f_5-- START OF CYCLIC NUMBER OF THE 5TH INPUT FILE
     1094      ;
     1095      ; if there is more than one file for all timesteps of a job chain, declare the first cyclicnumber which contains timesteps you'd like to use
     1096      ;
     1097      ; data type: integer
     1098      ;
     1099      ; example:   3 (if first file: example_ts.3.nc)
     1100      ;
     1101      ; default:   -1 
     1102      ;***************************************************
     1103      if (.not. isvar("start_f_5"))then 
     1104     
     1105         
     1106         start_f_5 = -1
     1107         
     1108         
     1109      end if
     1110      ;***************************************************
     1111      ; REQUIRED --end_f_5-- END OF CYCLIC NUMBER OF THE 5TH INPUT FILE
     1112      ;
     1113      ; if there is more than one file for all timesteps of a job chain, declare the last cyclicnumber which contains timesteps you'd like to use
     1114      ;
     1115      ; data type: integer
     1116      ;
     1117      ; example:   48 (if last file: example_ts.48.nc)
     1118      ;
     1119      ; default:   -1 
     1120      ;***************************************************
     1121      if (.not. isvar("end_f_5"))then 
     1122     
     1123         
     1124         end_f_5 = -1
     1125         
     1126         
     1127      end if     
    6811128      ;***************************************************
    6821129      ; OPTIONAL --name_legend_5-- 5th LEGEND ITEM
     
    6941141         
    6951142         
    696       end if     
    697       ;***************************************************
    698       ; REQUIRED IF no_files > 5 --file_6-- 6th INPUT FILE
    699       ;
    700       ; 6th input file produced by PALM; please declare path and name
    701       ;
    702       ; data type: string
     1143      end if
     1144      ;***************************************************
     1145      ; REQUIRED if no_files > 5 --file_6-- 6TH INPUT FILE
     1146      ;
     1147      ; 6th input file produced by PALM; if there is more than one file for all timesteps of a job chain,
     1148      ; declare path and file name without cycle number and ".nc" (first and last cycle numbers are given with the parameters 'start_f' and 'end_f');
     1149      ; otherwise declare path and full file name
     1150      ;
     1151      ; data type: string
     1152      ;
     1153      ; example: file_6 = "/path_to_file/file_pr"      if you use more than one cycle number
     1154      ;          file_6 = "/path_to_file/file_pr.4.nc" if you use one file (e.g. with cycle number 4)
    7031155      ;
    7041156      ; default:   "File in"
    7051157      ;***************************************************
    706       if (.not. isvar("file_6"))then   
    707      
    708          
    709          file_6 = "File in"
    710          
    711          
    712       end if     
    713       ;***************************************************
    714       ; OPTIONAL --name_legend_6-- 6th LEGEND ITEM
     1158      if(.not. isvar("file_6"))then
     1159
     1160
     1161         file_6 = "File in"
     1162         
     1163         
     1164      end if
     1165      ;***************************************************
     1166      ; REQUIRED --start_f_6-- START OF CYCLIC NUMBER OF THE 6TH INPUT FILE
     1167      ;
     1168      ; if there is more than one file for all timesteps of a job chain, declare the first cyclicnumber which contains timesteps you'd like to use
     1169      ;
     1170      ; data type: integer
     1171      ;
     1172      ; example:   3 (if first file: example_ts.3.nc)
     1173      ;
     1174      ; default:   -1 
     1175      ;***************************************************
     1176      if (.not. isvar("start_f_6"))then 
     1177     
     1178         
     1179         start_f_6 = -1
     1180         
     1181         
     1182      end if
     1183      ;***************************************************
     1184      ; REQUIRED --end_f_6-- END OF CYCLIC NUMBER OF THE 6TH INPUT FILE
     1185      ;
     1186      ; if there is more than one file for all timesteps of a job chain, declare the last cyclicnumber which contains timesteps you'd like to use
     1187      ;
     1188      ; data type: integer
     1189      ;
     1190      ; example:   48 (if last file: example_ts.48.nc)
     1191      ;
     1192      ; default:   -1 
     1193      ;***************************************************
     1194      if (.not. isvar("end_f_6"))then 
     1195     
     1196         
     1197         end_f_6 = -1
     1198         
     1199         
     1200      end if                     
     1201      ;***************************************************
     1202      ; OPTIONAL --name_legend_6-- 6TH LEGEND ITEM
    7151203      ;
    7161204      ; if no_files > 5: name of 6th item in extra legend
     
    7261214         
    7271215         
    728       end if     
    729       ;***************************************************
    730       ; REQUIRED IF no_files > 6 --file_7-- 7th INPUT FILE
    731       ;
    732       ; 7th input file produced by PALM; please declare path and name
    733       ;
    734       ; data type: string
    735       ;
    736       ; default:   "File in"
    737       ;***************************************************
    738       if (.not. isvar("file_7"))then
    739      
    740            
    741          file_7 = "File in"
    742          
    743          
    744       end if     
    745       ;***************************************************
    746       ; OPTIONAL --name_legend_7-- 7th LEGEND ITEM
    747       ;
    748       ; if no_files > 6: name of 7th item in extra legend
    749       ;
    750       ; data type: string
    751       ;
    752       ; default:   "File 7"
    753       ;***************************************************
    754       if (.not. isvar("name_legend_7"))then
    755      
    756            
    757          name_legend_7 = "File 7"
    758          
    759          
    760       end if     
    761       ;***************************************************
    762       ; REQUIRED IF no_files > 7 --file_8-- 8th INPUT FILE
    763       ;
    764       ; 8th input file produced by PALM; please declare path and name
    765       ;
    766       ; data type: string
    767       ;
    768       ; default:   "File in"
    769       ;***************************************************
    770       if (.not. isvar("file_8"))then
    771      
    772            
    773          file_8 = "File in"
    774          
    775          
    776       end if     
    777       ;***************************************************
    778       ; OPTIONAL --name_legend_8-- 8th LEGEND ITEM
    779       ;
    780       ; if no_files > 1: name of 8th item in extra legend
    781       ;
    782       ; data type: string
    783       ;
    784       ; default:   "File 8"
    785       ;***************************************************
    786       if (.not. isvar("name_legend_8"))then   
    787      
    788          
    789          name_legend_8 = "File 8"
    790          
    791          
    792       end if     
    793       ;***************************************************
    794       ; REQUIRED IF no_files > 8 --file_9-- 9th INPUT FILE
    795       ;
    796       ; 9th input file produced by PALM; please declare path and name
    797       ;
    798       ; data type: string
    799       ;
    800       ; default:   "File in"
    801       ;***************************************************
    802       if (.not. isvar("file_9"))then
    803      
    804            
    805          file_9 = "File in"
    806          
    807          
    808       end if     
    809       ;***************************************************
    810       ; OPTIONAL --name_legend_9-- 9th LEGEND ITEM
    811       ;
    812       ; if no_files > 8: name of 9th item in extra legend
    813       ;
    814       ; data type: string
    815       ;
    816       ; default:   "File 9"
    817       ;***************************************************
    818       if (.not. isvar("name_legend_9"))then
    819      
    820            
    821          name_legend_9 = "File 9"
    822          
    823          
    824       end if     
    825       ;***************************************************
    826       ; REQUIRED IF no_files > 9 --file_10-- 10th INPUT FILE
    827       ;
    828       ; 10th input file produced by PALM; please declare path and name
    829       ;
    830       ; data type: string
    831       ;
    832       ; default:   "File in"
    833       ;***************************************************
    834       if (.not. isvar("file_10"))then
    835      
    836            
    837          file_10 = "File in"
    838          
    839          
    840       end if     
    841       ;***************************************************
    842       ; OPTIONAL --name_legend_10-- 10th LEGEND ITEM
    843       ;
    844       ; if no_files > 9: name of 10th item in extra legend
    845       ;
    846       ; data type: string
    847       ;
    848       ; default:   "File 10"
    849       ;***************************************************
    850       if (.not. isvar("name_legend_10"))then 
    851      
    852          
    853          name_legend_10 = "File 10"
    854          
    855          
    856       end if     
    857       ;***************************************************
    858       ; REQUIRED IF no_files > 10 --file_11-- 11th INPUT FILE
    859       ;
    860       ; 11th input file produced by PALM; please declare path and name
    861       ;
    862       ; data type: string
    863       ;
    864       ; default:   "File in"
    865       ;***************************************************
    866       if (.not. isvar("file_11"))then
    867      
    868            
    869          file_11 = "File in"
    870          
    871          
    872       end if     
    873       ;***************************************************
    874       ; OPTIONAL --name_legend_11-- 11th LEGEND ITEM
    875       ;
    876       ; if no_files > 1: name of 11th item in extra legend
    877       ;
    878       ; data type: string
    879       ;
    880       ; default:   "File 11"
    881       ;***************************************************
    882       if (.not. isvar("name_legend_11"))then
    883      
    884            
    885          name_legend_11 = "File 11"
    886          
    887          
    888       end if     
    889       ;***************************************************
    890       ; REQUIRED IF no_files > 11 --file_12-- 12th INPUT FILE
    891       ;
    892       ; 12th input file produced by PALM; please declare path and name
    893       ;
    894       ; data type: string
    895       ;
    896       ; default:   "File in"
    897       ;***************************************************
    898       if (.not. isvar("file_12"))then
    899      
    900            
    901          file_12 = "File in"
    902          
    903          
    904       end if     
    905       ;***************************************************
    906       ; OPTIONAL --name_legend_12-- 12th LEGEND ITEM
    907       ;
    908       ; if no_files > 11: name of 12th item in extra legend
    909       ;
    910       ; data type: string
    911       ;
    912       ; default:   "File 12"
    913       ;***************************************************
    914       if (.not. isvar("name_legend_12"))then
    915      
    916            
    917          name_legend_12 = "File 12"
    918          
    919          
    920       end if     
     1216      end if             
    9211217      ;***************************************************
    9221218      ; OPTIONAL --start_time_step-- FIRST TIME STEP
    9231219      ;
    924       ; first time step of plot; you don't have to specify a precise time value from the input file,
    925       ; it will be rounded to the next existent time step
     1220      ; first time step of plot; there is no need to specify a precise time value from the input file,
     1221      ; 'start_time_step' is rounded to the next existent time step
    9261222      ;
    9271223      ; data type: double; dimension: [hour]
     1224      ;
     1225      ; example:   start_time_step = 0.5
    9281226      ;
    9291227      ; default:   -1.d (first available time step)
     
    9391237      ; OPTIONAL --end_time_step-- LAST TIME STEP
    9401238      ;
    941       ; last time step of plot; you don't have to specify a precise time value from the input file,
    942       ; it will be rounded to the next existent time step
     1239      ; last time step of plot; there is no need to specify a precise time value from the input file,
     1240      ; 'end_time_step' is rounded to the next existent time step
    9431241      ;
    9441242      ; data type: double; dimension: [hour]
    9451243      ;
     1244      ; example:   end_time_step = 2.0
     1245      ;
    9461246      ; default:   -1.d (last available time step)
    9471247      ;***************************************************
     
    9491249     
    9501250         
    951          end_time_step = -1.d     
    952          
    953          
    954       end if     
     1251         end_time_step = -1.d     
     1252         
     1253         
     1254      end if
     1255      ;***************************************************
     1256      ; OPTIONAL --start_x-- STARTVALUE IN X-DIRECTION 
     1257      ;
     1258      ; if 3D data is used for profiles you can choose a startvalue in x-direction for horizontal averaging
     1259      ;
     1260      ; data type: integer; dimension [grippoint]
     1261      ;
     1262      ; example: start_x = 0
     1263      ;
     1264      ; default:   0 (minimum of range)
     1265      ;***************************************************
     1266      if (.not. isvar("start_x"))then   
     1267     
     1268         
     1269         start_x = 0
     1270         
     1271         
     1272      end if
     1273      ;***************************************************
     1274      ; OPTIONAL --end_x-- ENDVALUE IN X-DIRECTION 
     1275      ;
     1276      ; if 3D data is used for profiles you can choose an endvalue in x-direction for horizontal averaging
     1277      ;
     1278      ; data type: integer; dimension [grippoint]
     1279      ;
     1280      ; example: end_x = 200
     1281      ;
     1282      ; default:   -1 (maximum of range excluding ghostpoint)
     1283      ;***************************************************
     1284      if (.not. isvar("end_x"))then   
     1285     
     1286         
     1287         end_x = -1
     1288         
     1289         
     1290      end if
     1291      ;***************************************************
     1292      ; OPTIONAL --start_y-- STARTVALUE IN Y-DIRECTION 
     1293      ;
     1294      ; if 3D data is used for profiles you can choose a startvalue in y-direction for horizontal averaging
     1295      ;
     1296      ; data type: integer; dimension [grippoint]
     1297      ;
     1298      ; example: start_y = 0
     1299      ;
     1300      ; default:   0 (minimum of range)
     1301      ;***************************************************
     1302      if (.not. isvar("start_y"))then   
     1303     
     1304         
     1305         start_y = 0
     1306         
     1307         
     1308      end if
     1309      ;***************************************************
     1310      ; OPTIONAL --end_y-- ENDVALUE IN Y-DIRECTION 
     1311      ;
     1312      ; if 3D data is used for profiles you can choose an endvalue in y-direction for horizontal averaging
     1313      ;
     1314      ; data type: integer; dimension [grippoint]
     1315      ;
     1316      ; example: end_y = 200
     1317      ;
     1318      ; default:   -1 (maximum of range excluding ghostpoint)
     1319      ;***************************************************
     1320      if (.not. isvar("end_y"))then   
     1321     
     1322         
     1323         end_y = -1
     1324         
     1325         
     1326      end if     
    9551327      ;***************************************************
    9561328      ; OPTIONAL --xs-- MINIMUM X-AXIS
     
    9591331      ;
    9601332      ; data type: float; dimension of variable
     1333      ;
     1334      ; example: xs = 300. (e.g. if var=",pt")
    9611335      ;
    9621336      ; default:   -1. (minimum of range)
     
    9761350      ; data type: float; dimension of variable
    9771351      ;
     1352      ; example: xe = 302. (e.g. if var=",pt")
     1353      ;
    9781354      ; default:   -1. (maximum of range)
    9791355      ;***************************************************
     
    9921368      ; data type: double; dimension: [meter]
    9931369      ;
     1370      ; example: min_z = 0.0
     1371      ;
    9941372      ; default:   -1.d (related to the minimum of z-axis in input file)
    9951373      ;***************************************************
     
    10081386      ; data type: double; dimension: [meter]
    10091387      ;
     1388      ; example: max_z = 1000.0
     1389      ;
    10101390      ; default:   -1.d (related to the maximum of z-axis in input file)
    10111391      ;***************************************************
     
    10181398      end if     
    10191399      ;***************************************************
    1020       ; OPTIONAL --logy-- LOGARITHMIC SCALE
    1021       ;
    1022       ; logarithmic scale for y-axis will be switched on [=1] or off [=0]
    1023       ;
    1024       ; data type: integer
     1400      ; OPTIONAL --log_z-- LOGARITHMIC SCALE OF Z-AXIS
     1401      ;
     1402      ; logarithmic scale for z-axis is switched on [=1] or off [=0]
     1403      ;
     1404      ; data type: integer
     1405      ;
     1406      ; example: log_z = 0
    10251407      ;
    10261408      ; default:   0
    10271409      ;***************************************************
    1028       if (.not. isvar("logy"))then
     1410      if (.not. isvar("log_z"))then
    10291411     
    10301412           
    1031          logy = 0
    1032          
    1033          
    1034       end if     
    1035       ;***************************************************
    1036       ; OPTIONAL --norm-- NORMALISE Y-AXIS
    1037       ;
    1038       ; value for normalising the y-axis
     1413         log_z = 0
     1414         
     1415         
     1416      end if     
     1417      ;***************************************************
     1418      ; OPTIONAL --norm_z-- NORMALISING Z-AXIS
     1419      ;
     1420      ; value for normalising the z-axis
    10391421      ;
    10401422      ; data type: float
    10411423      ;
    1042       ; default:   1 (no normalising)
    1043       ;***************************************************
    1044       if (.not. isvar("norm"))then
     1424      ; example: norm_z = 1.
     1425      ;
     1426      ; default:   1. (no normalising)
     1427      ;***************************************************
     1428      if (.not. isvar("norm_z"))then
    10451429     
    10461430           
    1047          norm = 1
     1431         norm_z = 1.
    10481432         
    10491433         
     
    10521436      ; OPTIONAL --over-- OVERLAYING
    10531437      ;
    1054       ; predefined overlaying of standard variables will be switched on [=1] or off [=0];
     1438      ; predefined overlaying of standard variables is switched on [=1] or off [=0];
    10551439      ; this feature cannot be used with statistic regions
    10561440      ;
    10571441      ; data type: integer
    10581442      ;
     1443      ; example: over = 0
     1444      ;
    10591445      ; default:   0
    10601446      ;***************************************************
     
    10721458      ;
    10731459      ; data type: integer
     1460      ;
     1461      ; example: combine = 0
    10741462      ;
    10751463      ; default:   0
     
    10891477      ; data type: integer
    10901478      ;
     1479      ; example: number_comb = 2
     1480      ;
    10911481      ; default:   -1
    10921482      ;***************************************************
     
    11011491      ; REQUIRED IF combine=1 --c_var-- VARIABLES FOR COMBINING
    11021492      ;
    1103       ; names of variables; please inidcate them with one comma before and after each variable;
    1104       ; e.g.: var=",u,v,w," (the windcomponents u, v and w will be combined)
    1105       ;
    1106       ; data type: string
     1493      ; names of variables; inidcate them with one comma before and after each variable;
     1494      ;
     1495      ; data type: string
     1496      ;
     1497      ; example: c_var = ",u,v,"
    11071498      ;
    11081499      ; default:   "c_variables"
     
    11221513      ; data type: integer
    11231514      ;
     1515      ; example: black = 0
     1516      ;
    11241517      ; default:   0
    11251518      ;***************************************************
     
    11381531      ; data type: integer
    11391532      ;
     1533      ; example: dash = 0
     1534      ;
    11401535      ; default:   0
    11411536      ;***************************************************
     
    11441539       
    11451540         dash = 0     
     1541         
     1542         
     1543      end if
     1544      ;***************************************************
     1545      ; OPTIONAL --font_size-- FONT SIZE OF STRINGS
     1546      ;
     1547      ; font size can be changed to fit to the current plot
     1548      ;
     1549      ; data type: float
     1550      ;
     1551      ; example: font_size = 0.02
     1552      ;
     1553      ; default:   0.02
     1554      ;***************************************************
     1555      if(.not. isvar("font_size"))then
     1556     
     1557         
     1558         font_size = 0.02
     1559         
     1560         
     1561      end if
     1562      ;***************************************************
     1563      ; OPTIONAL --font_size_legend-- FONT SIZE OF LEGEND STRINGS
     1564      ;
     1565      ; font size of the legend can be changed to fit to the current plot
     1566      ;
     1567      ; data type: float
     1568      ;
     1569      ; example: font_size_legend = 0.2
     1570      ;
     1571      ; default:   0.2
     1572      ;***************************************************
     1573      if(.not. isvar("font_size_legend"))then
     1574     
     1575         
     1576         font_size_legend = 0.2
    11461577         
    11471578         
     
    11601591      ; REQUIRED --file_1-- INPUT FILE
    11611592      ;
    1162       ; input file produced by PALM; please declare path and name
    1163       ;
    1164       ; data type: string
     1593      ; input file produced by PALM; if there is more than one file for all timesteps of a job chain,
     1594      ; declare path and file name without cycle number and ".nc" (first and last cycle numbers are given with the parameters 'start_f' and 'end_f');
     1595      ; otherwise declare path and full file name
     1596      ;
     1597      ; data type: string
     1598      ;
     1599      ; example: file_1 = "/path_to_file/file_sp"      if you use more than one cycle number
     1600      ;          file_1 = "/path_to_file/file_sp.4.nc" if you use one file (e.g. with cycle number 4)
    11651601      ;
    11661602      ; default:   "File in"
    11671603      ;***************************************************
    1168       if (.not. isvar("file_1"))then 
     1604      if(.not. isvar("file_1"))then
     1605
     1606
     1607         file_1 = "File in"
     1608         
     1609         
     1610      end if
     1611      ;***************************************************
     1612      ; REQUIRED if more than one file of a job chain --start_f-- START OF CYCLIC NUMBER
     1613      ;
     1614      ; if there is more than one file for all timesteps of a job chain, declare the first cyclicnumber which contains timesteps you'd like to use
     1615      ;
     1616      ; data type: integer
     1617      ;
     1618      ; example:   start_f = 3 (if first file: example_ts.3.nc)
     1619      ;
     1620      ; default:   -1 
     1621      ;***************************************************
     1622      if (.not. isvar("start_f"))then 
    11691623     
    11701624         
    1171          file_1 = "File in"
    1172          
    1173          
    1174       end if     
     1625         start_f = -1
     1626         
     1627         
     1628      end if
     1629      ;***************************************************
     1630      ; REQUIRED if more than one file of a job chain --end_f-- END OF CYCLIC NUMBER
     1631      ;
     1632      ; if there is more than one file for all timesteps of a job chain, declare the last cyclicnumber which contains timesteps you'd like to use
     1633      ;
     1634      ; data type: integer
     1635      ;
     1636      ; example:   end_f = 48 (if last file: example_ts.48.nc)
     1637      ;
     1638      ; default:   -1 
     1639      ;***************************************************
     1640      if (.not. isvar("end_f"))then 
     1641     
     1642         
     1643         end_f = -1
     1644         
     1645         
     1646      end if
    11751647      ;***************************************************
    11761648      ; OPTIONAL --format_out-- OUTPUT FORMAT
    11771649      ;
    1178       ; supported file formats: x11, pdf, ps, eps, epsi or ncgm
    1179       ;
    1180       ; data type: string
     1650      ; format of the output file; supported file formats: x11, pdf, ps, eps, epsi or ncgm
     1651      ;
     1652      ; data type: string
     1653      ;
     1654      ; example: format_out = "pdf"
    11811655      ;
    11821656      ; default:   "x11"
     
    11921666      ; OPTIONAL --file_out-- OUTPUT FILE
    11931667      ;
    1194       ; please declare path and name
    1195       ;
    1196       ; data type: string
     1668      ; path and name of the output file
     1669      ;
     1670      ; data type: string
     1671      ;
     1672      ; example: file_out = "/home/path_to_output/test_sp"
    11971673      ;
    11981674      ; default:   "~/test_sp"
     
    12061682      end if     
    12071683      ;***************************************************
    1208       ; OPTIONAL --no_columns-- NUMBER OF ROWS
     1684      ; OPTIONAL --no_columns-- NUMBER OF COLUMNS
    12091685      ;
    12101686      ; number of plots in one row
    12111687      ;
    12121688      ; data type: integer
     1689      ;
     1690      ; example: no_columns = 1
    12131691      ;
    12141692      ; default:   1
     
    12221700      end if     
    12231701      ;***************************************************
    1224       ; OPTIONAL --no_lines-- NUMBER OF COLUMNS
     1702      ; OPTIONAL --no_rows-- NUMBER OF ROWS
    12251703      ;
    12261704      ; number of plots in one column
    12271705      ;
    12281706      ; data type: integer
     1707      ;
     1708      ; example: no_rows = 2
    12291709      ;
    12301710      ; default:   2
    12311711      ;***************************************************
    1232       if (.not. isvar("no_lines"))then
    1233      
    1234          
    1235          no_lines = 2
     1712      if (.not. isvar("no_rows"))then
     1713     
     1714         
     1715         no_rows = 2
    12361716         
    12371717         
     
    12401720      ; OPTIONAL --var-- OUTPUT VARIABLES
    12411721      ;
    1242       ; names of output variables; please inidcate them with one comma before and after each variable;
    1243       ; e.g.: var=",u,v,w," (the windcomponents u, v and w will be plotted)
     1722      ; names of output variables; inidcate them with one comma before and after each variable
    12441723      ;
    12451724      ; data type: string
    1246       ;
    1247       ; default:   "all" (all available variables will be plotted)
     1725      ;
     1726      ; example: var = ",u,v,w"
     1727      ;
     1728      ; default:   "all" (all available variables are plotted)
    12481729      ;***************************************************
    12491730      if (.not. isvar("var"))then
     
    12571738      ; OPTIONAL --height_level-- HEIGHT LEVELS
    12581739      ;
    1259       ; defines if all heights or selected heights will be output; a 1D array with up to 10 elements
    1260       ; from 0 to 9 must be indicated, e.g.: height_level=(/0,3,7/) (level 0, 3 and 7 as in input file
    1261       ; will be output)
     1740      ; defines if all levels or selected levels are output; a 1D array with up to 100 elements
     1741      ; from 0 to 99 must be indicated
    12621742      ;
    12631743      ; data type: integer array
    12641744      ;
    1265       ; default:   (/-1/) (all levels will be output)
     1745      ; example: height_level=(/0,3,7/)
     1746      ;
     1747      ; default:   (/-1/) (all levels are output)
    12661748      ;***************************************************
    12671749      if (.not. isvar("height_level"))then
     
    12751757      ; OPTIONAL --sort-- TYPE OF SORTING
    12761758      ;
    1277       ; defines if all heights [="height"] or all timesteps [="time"] will be merged in one plot
    1278       ;
    1279       ; data type: string
     1759      ; defines if all heights [="height"] or all timesteps [="time"] are merged in one plot
     1760      ;
     1761      ; data type: string
     1762      ;
     1763      ; example: sort = "height"
    12801764      ;
    12811765      ; default:   "height"
     
    12911775      ; OPTIONAL --start_time_step-- FIRST TIME STEP
    12921776      ;
    1293       ; first time step of plot; you don't have to specify a precise time value from the input file,
    1294       ; it will be rounded to the next existent time step
     1777      ; first time step of plot; there is no need to specify a precise time value from the input file,
     1778      ; 'start_time_step' is rounded to the next existent time step
    12951779      ;
    12961780      ; data type: double; dimension: [hour]
    12971781      ;
     1782      ; example:   start_time_step = 0.5
     1783      ;
    12981784      ; default:   -1.d (first available time step)
    12991785      ;***************************************************
     
    13081794      ; OPTIONAL --end_time_step-- LAST TIME STEP
    13091795      ;
    1310       ; last time step of plot; you don't have to specify a precise time value from the input file,
    1311       ; it will be rounded to the next existent time step
     1796      ; last time step of plot; there is no need to specify a precise time value from the input file,
     1797      ; 'end_time_step' is rounded to the next existent time step
    13121798      ;
    13131799      ; data type: double; dimension: [hour]
    13141800      ;
     1801      ; example:   end_time_step = 2.0
     1802      ;
    13151803      ; default:   -1.d (last available time step)
    13161804      ;***************************************************
     
    13281816      ;
    13291817      ; data type: integer
     1818      ;
     1819      ; example: black = 0
    13301820      ;
    13311821      ; default:   0
     
    13451835      ; data type: integer
    13461836      ;
     1837      ; example: dash = 0
     1838      ;
    13471839      ; default:   0
    13481840      ;***************************************************
     
    13551847      end if
    13561848      ;***************************************************
    1357       ; OPTIONAL --logx-- LOGARITHMIC SCALE X-AXIS
    1358       ;
    1359       ; logarithmic scale for x-axis will be switched on [=1] or off [=0]
    1360       ;
    1361       ; data type: integer
     1849      ; OPTIONAL --log_x-- LOGARITHMIC SCALE X-AXIS
     1850      ;
     1851      ; logarithmic scale for x-axis is switched on [=1] or off [=0]
     1852      ;
     1853      ; data type: integer
     1854      ;
     1855      ; example: log_x = 1
    13621856      ;
    13631857      ; default:   1
    13641858      ;***************************************************
    1365       if (.not. isvar("logx"))then   
     1859      if (.not. isvar("log_x"))then   
    13661860     
    13671861         
    1368          logx = 1
     1862         log_x = 1
    13691863         
    13701864         
    13711865      end if     
    13721866      ;***************************************************
    1373       ; OPTIONAL --logy-- LOGARITHMIC SCALE Y-AXIS
    1374       ;
    1375       ; logarithmic scale for y-axis will be switched on [=1] or off [=0]
    1376       ;
    1377       ; data type: integer
     1867      ; OPTIONAL --log_y-- LOGARITHMIC SCALE Y-AXIS
     1868      ;
     1869      ; logarithmic scale for y-axis is switched on [=1] or off [=0]
     1870      ;
     1871      ; data type: integer
     1872      ;
     1873      ; example: log_y = 1
    13781874      ;
    13791875      ; default:   1
    13801876      ;***************************************************
    1381       if (.not. isvar("logy"))then   
     1877      if (.not. isvar("log_y"))then   
    13821878     
    13831879         
    1384          logy = 1
    1385          
    1386          
    1387       end if     
    1388       ;***************************************************
    1389       ; OPTIONAL --normx-- NORMALISE X-AXIS
     1880         log_y = 1
     1881         
     1882         
     1883      end if     
     1884      ;***************************************************
     1885      ; OPTIONAL --norm_x-- NORMALISING X-AXIS
    13901886      ;
    13911887      ; value for normalising the x-axis
     
    13931889      ; data type: float
    13941890      ;
     1891      ; example: norm_x = 1
     1892      ;
    13951893      ; default:   1. (no normalising)
    13961894      ;***************************************************
    1397       if (.not. isvar("normx"))then   
     1895      if (.not. isvar("norm_x"))then   
    13981896     
    13991897         
    1400          normx = 1.
    1401          
    1402          
    1403       end if     
    1404       ;***************************************************
    1405       ; OPTIONAL --normy-- NORMALISE Y-AXIS
     1898         norm_x = 1.
     1899         
     1900         
     1901      end if
     1902      ;***************************************************
     1903      ; OPTIONAL --norm_height-- NORMALISING X-AXIS WITH HEIGHT
     1904      ;
     1905      ; normalising the x-axis with the height of each spectra is switched on [=1] or off [=0]
     1906      ;
     1907      ; data type: integer
     1908      ;
     1909      ; example: norm_height = 0
     1910      ;
     1911      ; default:   0 (no normalising)
     1912      ;***************************************************
     1913      if (.not. isvar("norm_height"))then   
     1914     
     1915         
     1916         norm_height = 0
     1917         
     1918         
     1919      end if     
     1920      ;***************************************************
     1921      ; OPTIONAL --norm_y-- NORMALISE Y-AXIS
    14061922      ;
    14071923      ; value for normalising the y-axis
     
    14091925      ; data type: float
    14101926      ;
     1927      ; example: norm_y = 1.
     1928      ;
    14111929      ; default:   1. (no normalising)
    14121930      ;***************************************************
    1413       if (.not. isvar("normy"))then     
     1931      if (.not. isvar("norm_y"))then     
    14141932     
    14151933       
    1416          normy = 1.
    1417          
    1418          
    1419       end if     
     1934         norm_y = 1.
     1935         
     1936         
     1937      end if
     1938      ;***************************************************
     1939      ; OPTIONAL --unit_x--unit_y-- UNITS OF X-, Y-AXES
     1940      ;
     1941      ; unit of the x-axis can be set if norm_x is unequal 1, otherwise the unit is [1/m]; by default no unit is set for the y-axis
     1942      ;
     1943      ; data type: string,
     1944      ;
     1945      ; example: unit_x = "100/m"
     1946      ;
     1947      ; default:   " "
     1948      ;***************************************************
     1949      if(.not. isvar("unit_x"))then     
     1950         
     1951         unit_x = " "   
     1952         
     1953      end if
     1954      if(.not. isvar("unit_y"))then     
     1955         
     1956         unit_y = " "   
     1957         
     1958      end if
     1959      ;***************************************************
     1960      ; OPTIONAL --font_size-- FONT SIZE OF STRINGS
     1961      ;
     1962      ; font size can be changed to fit to the current plot
     1963      ;
     1964      ; data type: float
     1965      ;
     1966      ; example: font_size = 0.02
     1967      ;
     1968      ; default:   0.02
     1969      ;***************************************************
     1970      if(.not. isvar("font_size"))then
     1971     
     1972         
     1973         font_size = 0.02
     1974         
     1975         
     1976      end if
     1977      ;***************************************************
     1978      ; OPTIONAL --font_size_legend-- FONT SIZE OF LEGEND STRINGS
     1979      ;
     1980      ; font size of the legend can be changed to fit to the current plot
     1981      ;
     1982      ; data type: float
     1983      ;
     1984      ; example: font_size_legend = 0.2
     1985      ;
     1986      ; default:   0.2
     1987      ;***************************************************
     1988      if(.not. isvar("font_size_legend"))then
     1989     
     1990         
     1991         font_size_legend = 0.2
     1992         
     1993         
     1994      end if     
    14201995      ;***************************************************
    14211996     
     
    14312006      ; REQUIRED --file_1-- INPUT FILE
    14322007      ;
    1433       ; input file produced by PALM; please declare path and name
    1434       ;
    1435       ; data type: string
     2008      ; input file produced by PALM; if there is more than one file for all timesteps of a job chain,
     2009      ; declare path and file name without cycle number and ".nc" (first and last cycle numbers are given with the parameters 'start_f' and 'end_f');
     2010      ; otherwise declare path and full file name
     2011      ;
     2012      ; data type: string
     2013      ;
     2014      ; example: file_1 = "/path_to_file/file_ts"      if you use more than one cycle number
     2015      ;          file_1 = "/path_to_file/file_ts.4.nc" if you use one file (e.g. with cycle number 4)
    14362016      ;
    14372017      ; default:   "File in"
    14382018      ;***************************************************
    1439       if (.not. isvar("file_1"))then 
     2019      if(.not. isvar("file_1"))then
     2020
     2021
     2022         file_1 = "File in"
     2023         
     2024         
     2025      end if
     2026      ;***************************************************
     2027      ; REQUIRED if more than one file of a job chain --start_f-- START OF CYCLIC NUMBER
     2028      ;
     2029      ; if there is more than one file for all timesteps of a job chain, declare the first cyclicnumber which contains timesteps you'd like to use
     2030      ;
     2031      ; data type: integer
     2032      ;
     2033      ; example:   start_f = 3 (if first file: example_ts.3.nc)
     2034      ;
     2035      ; default:   -1 
     2036      ;***************************************************
     2037      if (.not. isvar("start_f"))then 
    14402038     
    14412039         
    1442          file_1 = "File in"
    1443          
    1444          
    1445       end if     
     2040         start_f = -1
     2041         
     2042         
     2043      end if
     2044      ;***************************************************
     2045      ; REQUIRED if more than one file of a job chain --end_f-- END OF CYCLIC NUMBER
     2046      ;
     2047      ; if there is more than one file for all timesteps of a job chain, declare the last cyclicnumber which contains timesteps you'd like to use
     2048      ;
     2049      ; data type: integer
     2050      ;
     2051      ; example:   end_f = 48 (if last file: example_ts.48.nc)
     2052      ;
     2053      ; default:   -1 
     2054      ;***************************************************
     2055      if (.not. isvar("end_f"))then 
     2056     
     2057         
     2058         end_f = -1
     2059         
     2060         
     2061      end if
    14462062      ;***************************************************
    14472063      ; OPTIONAL --format_out-- OUTPUT FORMAT
    14482064      ;
    1449       ; supported file formats: x11, pdf, ps, eps, epsi or ncgm
    1450       ;
    1451       ; data type: string
     2065      ; format of the output file; supported file formats: x11, pdf, ps, eps, epsi or ncgm
     2066      ;
     2067      ; data type: string
     2068      ;
     2069      ; example: format_out = "pdf"
    14522070      ;
    14532071      ; default:   "x11"
     
    14632081      ; OPTIONAL --file_out-- OUTPUT FILE
    14642082      ;
    1465       ; please declare path and name
    1466       ;
    1467       ; data type: string
     2083      ; path and name of the output file
     2084      ;
     2085      ; data type: string
     2086      ;
     2087      ; example: file_out = "/home/path_to_output/test_ts"
    14682088      ;
    14692089      ; default:   "~/test_ts"
     
    14772097      end if     
    14782098      ;***************************************************
    1479       ; OPTIONAL --no_columns-- NUMBER OF ROWS
     2099      ; OPTIONAL --no_columns-- NUMBER OF COLUMNS
    14802100      ;
    14812101      ; number of plots in one row
    14822102      ;
    14832103      ; data type: integer
     2104      ;
     2105      ; example: no_columns = 1
    14842106      ;
    14852107      ; default:   1
     
    14932115      end if     
    14942116      ;***************************************************
    1495       ; OPTIONAL --no_lines-- NUMBER OF COLUMNS
     2117      ; OPTIONAL --no_rows-- NUMBER OF ROWS
    14962118      ;
    14972119      ; number of plots in one column
    14982120      ;
    14992121      ; data type: integer
     2122      ;
     2123      ; example: no_rows = 2
    15002124      ;
    15012125      ; default:   2
    15022126      ;***************************************************
    1503       if (.not. isvar("no_lines"))then
    1504      
    1505          
    1506          no_lines = 2
     2127      if (.not. isvar("no_rows"))then
     2128     
     2129         
     2130         no_rows = 2
    15072131         
    15082132         
     
    15112135      ; OPTIONAL --var-- OUTPUT VARIABLES
    15122136      ;
    1513       ; names of output variables; please inidcate them with one comma before and after each variable;
    1514       ; e.g.: var=",u,v,w," (the windcomponents u, v and w will be plotted)
     2137      ; names of output variables; inidcate them with one comma before and after each variable
    15152138      ;
    15162139      ; data type: string
    1517       ;
    1518       ; default:   "all" (all available variables will be plotted)
     2140      ;
     2141      ; example: var = ",u,v,w"
     2142      ;
     2143      ; default:   "all" (all available variables are plotted)
    15192144      ;***************************************************
    15202145      if (.not. isvar("var"))then
     
    15282153      ; OPTIONAL --start_time_step-- FIRST TIME STEP
    15292154      ;
    1530       ; first time step of plot; you don't have to specify a precise time value from the input file,
    1531       ; it will be rounded to the next existent time step
     2155      ; first time step of plot; there is no need to specify a precise time value from the input file,
     2156      ; 'start_time_step' is rounded to the next existent time step
    15322157      ;
    15332158      ; data type: double; dimension: [hour]
    15342159      ;
     2160      ; example:   start_time_step = 0.5
     2161      ;
    15352162      ; default:   -1.d (first available time step)
    15362163      ;***************************************************
     
    15452172      ; OPTIONAL --end_time_step-- LAST TIME STEP
    15462173      ;
    1547       ; last time step of plot; you don't have to specify a precise time value from the input file,
    1548       ; it will be rounded to the next existent time step
     2174      ; last time step of plot; there is no need to specify a precise time value from the input file,
     2175      ; 'end_time_step' is rounded to the next existent time step
    15492176      ;
    15502177      ; data type: double; dimension: [hour]
    15512178      ;
     2179      ; example:   end_time_step = 2.0
     2180      ;
    15522181      ; default:   -1.d (last available time step)
    15532182      ;***************************************************
     
    15622191      ; OPTIONAL --over-- OVERLAYING
    15632192      ;
    1564       ; predefined overlaying of standard variables will be switched on [=1] or off [=0];
     2193      ; predefined overlaying of standard variables are switched on [=1] or off [=0];
    15652194      ; this feature cannot be used with statistic regions
    15662195      ;
    15672196      ; data type: integer
    15682197      ;
     2198      ; example: over = 0
     2199      ;
    15692200      ; default:   0
    15702201      ;***************************************************
     
    15752206         
    15762207         
    1577       end if
     2208      end if
     2209      ;***************************************************
     2210      ; OPTIONAL --font_size-- FONT SIZE OF STRINGS
     2211      ;
     2212      ; font size can be changed to fit to the current plot
     2213      ;
     2214      ; data type: float
     2215      ;
     2216      ; example: font_size = 0.07
     2217      ;
     2218      ; default:   0.07
     2219      ;***************************************************
     2220      if(.not. isvar("font_size"))then
     2221     
     2222         
     2223         font_size = 0.07
     2224         
     2225         
     2226      end if
     2227      ;***************************************************
     2228      ; OPTIONAL --norm_t-- NORMALISING X-AXIS OF TIMESERIES PLOTS (T-AXIS)
     2229      ;
     2230      ; normalising the x-axis with the given value; the unit can explicitly set; by default no normalising and unit is set to [s]
     2231      ;
     2232      ; data type: float
     2233      ;
     2234      ; example: norm_t = 3600
     2235      ;
     2236      ; default:   1.
     2237      ;***************************************************
     2238      if(.not. isvar("normt"))then
     2239     
     2240         
     2241         norm_t = 1.
     2242         
     2243         
     2244      end if
     2245      ;***************************************************
     2246      ; OPTIONAL --unit_t-- UNIT OF X-AXIS OF TIMESERIES PLOTS (T-AXIS)
     2247      ;
     2248      ; unit of the x-axis can be set; by default unit is set to [s]
     2249      ;
     2250      ; data type: string
     2251      ;
     2252      ; example: unit_t = "h"
     2253      ;
     2254      ; default:   " "
     2255      ;***************************************************
     2256      if(.not. isvar("unitt"))then
     2257     
     2258         
     2259         unit_t = " "
     2260         
     2261         
     2262      end if
    15782263      ;***************************************************
    15792264
Note: See TracChangeset for help on using the changeset viewer.