source: palm/trunk/SCRIPTS/NCL/.ncl.config.default @ 1836

Last change on this file since 1836 was 958, checked in by hoffmann, 12 years ago

minor changes in palmplot pr

File size: 69.9 KB
RevLine 
[418]1;*********************************************************
2; Copying .hluresfile in the home directory
3;*********************************************************
4
5command = "ksh 'if [ ! -f $HOME/.hluresfile ]; then; "+\
6          "cp $PALM_BIN/NCL/.hluresfile $HOME; fi'"
7system(command)
8   
9;*********************************************************
10
11
[190]12begin
13
14;*********************************************************
[194]15
16script=which_script
[526]17
18   ;******************************************************
19   ; BE AWARE: NCl uses double quotes (") to denote strings
20   ;       to retrieve the double quote character - use dq:
21
22    dq=str_get_dq()
23
24   ;******************************************************
[190]25   
26;*********************************************************
[194]27   
[190]28   ;******************************************************
29   ; parameters used by CROSS_SECTIONS
30   ;******************************************************
31   
[194]32   if (script .EQ. "cross_section")then
[190]33   
34      ;***************************************************
35      ; REQUIRED --file_1-- INPUT FILE
36      ;
[418]37      ; input file produced by PALM; if there is more than one file for
38      ; all timesteps
39      ; of a job chain,
40      ; declare path and file name without cycle number and ".nc" (first and
41      ; last cycle numbers are given with the parameters 'start_f' and
42      ; 'end_f');
[218]43      ; otherwise declare path and full file name
[190]44      ;
45      ; data type: string
46      ;
[418]47      ; example: file_1 = "/path_to_file/file_xy_av"      if you use more than
48      ;                                                   one cycle number
49      ;          file_1 = "/path_to_file/file_xy_av.4.nc" if you use one file
50      ;                                              (e.g. with cycle number 4)
[218]51      ;
[190]52      ; default:   "File in"
53      ;***************************************************
54      if(.not. isvar("file_1"))then
[194]55
56
[267]57         file_1 = "File in"
[190]58         
[194]59         
[190]60      end if
61      ;***************************************************
[513]62      ; OPTIONAL/REQUIRED if more than one file of a job chain
[418]63      ;          --start_f-- START OF CYCLIC NUMBER
[190]64      ;
[418]65      ; if there is more than one file for all timesteps of a job chain,
66      ; declare the first cyclicnumber which contains timesteps
67      ; you'd like to use
[218]68      ;
69      ; data type: integer
70      ;
[418]71      ; example:   start_f = 3 (if first file: example_xy_av.3.nc) and
72      ;            start_f = 0 (if first file: example_xy_av.nc)
[218]73      ;
74      ; default:   -1 
75      ;***************************************************
76      if (.not. isvar("start_f"))then 
77     
78         
79         start_f = -1
80         
81         
82      end if
83      ;***************************************************
[513]84      ; OPTIONAL/REQUIRED if more than one file of a job chain
85      ;          --end_f-- END OF CYCLIC NUMBER
[218]86      ;
[513]87      ; if there is more than one file for all timesteps of a job chain,
88      ; declare the last cyclicnumber which contains timesteps you'd like
89      ; to use
[218]90      ;
91      ; data type: integer
92      ;
[418]93      ; example:   end_f = 48 (if last file: example_xy_av.48.nc) and
94      ;            end_f = 0 (if last file: example_xy_av.nc)
[218]95      ;
96      ; default:   -1 
97      ;***************************************************
98      if (.not. isvar("end_f"))then 
99     
100         
101         end_f = -1
102         
103         
104      end if
105      ;***************************************************
[190]106      ; OPTIONAL --format_out-- OUTPUT FORMAT
107      ;
[418]108      ; format of the output file;
[532]109      ; supported file formats: x11, pdf, ps, eps, epsi ncgm or png
[190]110      ;
111      ; data type: string
112      ;
[218]113      ; example: format_out = "pdf"
114      ;
[190]115      ; default:   "x11"
116      ;***************************************************
117      if(.not. isvar("format_out"))then 
118     
119     
120         format_out = "x11"
121         
122         
123      end if
124      ;***************************************************
125      ; OPTIONAL --file_out-- OUTPUT FILE
126      ;
[218]127      ; path and name of the output file
[190]128      ;
129      ; data type: string
130      ;
[218]131      ; example: file_out = "/home/path_to_output/test_cs"
132      ;
[194]133      ; default:   "~/test_cs"
[190]134      ;***************************************************
135      if(.not. isvar("file_out"))then   
136     
137     
138         file_out = "~/test_cs"
139         
140         
141      end if     
142      ;***************************************************
143      ; OPTIONAL --var-- OUTPUT VARIABLES
144      ;
[526]145      ; names of output variables;
146      ; inidcate them with one comma before and after each variable;
147      ; if the name contains double quotes (") use dq and the '+' operation to
148      ; concatenate the strings
[190]149      ;
150      ; data type: string
[218]151      ;
[342]152      ; example: var = ",u,v,w,"
[190]153      ;
[218]154      ; default:   "all" (all available variables are plotted)
[190]155      ;***************************************************
156      if(.not. isvar("var"))then
157     
158         
159         var = "all"
160         
161         
162      end if     
163      ;***************************************************
164      ; OPTIONAL --start_time_step-- FIRST TIME STEP
165      ;
[418]166      ; first time step of plot; there is no need to specify a precise time
167      ; value from the input file,
[218]168      ; 'start_time_step' is rounded to the next existent time step
[190]169      ;
170      ; data type: double; dimension: [hour]
171      ;
[418]172      ; example:   start_time_step = 0.5d
[218]173      ;
[190]174      ; default:   -1.d (first available time step)
175      ;***************************************************
176      if(.not. isvar("start_time_step"))then
177     
178         
179         start_time_step = -1.d
180         
181         
182      end if     
183      ;***************************************************
184      ; OPTIONAL --end_time_step-- LAST TIME STEP
185      ;
[418]186      ; last time step of plot; there is no need to specify a precise time
187      ; value from the input file,
[218]188      ; 'end_time_step' is rounded to the next existent time step
[190]189      ;
190      ; data type: double; dimension: [hour]
191      ;
[418]192      ; example:   end_time_step = 2.0d
[218]193      ;
[190]194      ; default:   -1.d (last available time step)
195      ;***************************************************
196      if(.not. isvar("end_time_step"))then       
197     
198     
199         end_time_step = -1.d
200         
201         
202      end if     
203      ;***************************************************
204      ; OPTIONAL --xs-- MINIMUM X-AXIS
205      ;
[513]206      ; value for minimum of x-axis, selectable if there are
[566]207      ; NO preset layers
[418]208      ; for x (as for XY or XZ 2D-DATA);
209      ; there is no need to specify a precise meter value from the input file,
210      ; it is rounded to the next existent value
[190]211      ;
212      ; data type: double; dimension: [meter]
213      ;
[418]214      ; example: xs = 150.0d
[218]215      ;
[190]216      ; default:   -1.d (minimum of range)
217      ;***************************************************
218      if(.not. isvar("xs"))then     
219     
220     
221         xs = -1.d
222         
223         
224      end if     
225      ;***************************************************
226      ; OPTIONAL --xe-- MAXIMUM X-AXIS
227      ;
[513]228      ; value for maximum of x-axis, selectable if there are
[566]229      ; NO preset layers
[418]230      ; for x (as for XY or XZ 2D-DATA);
231      ; there is no need to specify a precise meter value from the input file,
232      ; it is rounded to the next existent value
[190]233      ;
234      ; data type: double; dimension: [meter]
235      ;
[418]236      ; example: xe = 300.0d
[218]237      ;
[190]238      ; default:   -1.d (maximum of range)
239      ;***************************************************
240      if(.not. isvar("xe"))then     
241     
242     
243         xe = -1.d
244         
245         
246      end if     
247      ;***************************************************
248      ; OPTIONAL --ys-- MINIMUM Y-AXIS
249      ;
[513]250      ; value for minimum of y-axis, selectable if there are
[566]251      ; NO preset layers
[418]252      ; for y (as for XY or YZ 2D-DATA);
253      ; there is no need to specify a precise meter value from the input file,
254      ; it is rounded to the next existent value
[190]255      ;
256      ; data type: double; dimension: [meter]
257      ;
[418]258      ; example: ys = 150.0d
[218]259      ;
[190]260      ; default:   -1.d (minimum of range)
261      ;***************************************************
262      if(.not. isvar("ys"))then     
263     
264     
265         ys = -1.d
266         
267         
268      end if     
269      ;***************************************************
270      ; OPTIONAL --ye-- MAXIMUM Y-AXIS
271      ;
[566]272      ; value for maximum of y-axis, selectable if there are NO preset layers
[418]273      ; for y (as for XY or YZ 2D-DATA);
274      ; there is no need to specify a precise meter value from the input file,
275      ; it is rounded to the next existent value
[190]276      ;
277      ; data type: double; dimension: [meter]
278      ;
[418]279      ; example: ye = 300.0d
[218]280      ;
[190]281      ; default:   -1.d (maximum of range)
282      ;***************************************************
283      if(.not. isvar("ye"))then   
284     
285         
286         ye = -1.d
287         
288         
289      end if     
290      ;***************************************************
291      ; OPTIONAL --zs-- MINIMUM Z-AXIS
292      ;
[513]293      ; index for minimum of z-axis, selectable if there are
[566]294      ; NO preset layers
[418]295      ; for z (as for XZ or YZ 2D-DATA);
[513]296      ; you cannot specify a meter value from the input file due to
297      ; grid stretching 
[190]298      ;
299      ; data type: integer
300      ;
[218]301      ; example: zs = 0
302      ;
[190]303      ; default:   -1 (minimum of range)
304      ;***************************************************
305      if(.not. isvar("zs"))then     
306     
307     
308         zs = -1
309         
310         
311      end if     
312      ;***************************************************
313      ; OPTIONAL --ze-- MAXIMUM Z-AXIS
314      ;
[513]315      ; index for maximum of z-axis, selectable if there are
[566]316      ; NO preset layers
[418]317      ; for z (as for XZ or YZ 2D-DATA);
[513]318      ; you cannot specify a meter value from the input file due to
319      ; grid stretching 
[190]320      ;
321      ; data type: integer
322      ;
[218]323      ; example: ze = 15
324      ;
[190]325      ; default:   -1 (maximum of range)
326      ;***************************************************
327      if(.not. isvar("ze"))then
328     
329           
330         ze = -1
331         
332         
333      end if     
334      ;***************************************************
[566]335      ; OPTIONAL --no_columns-- NUMBER OF COLUMNS
336      ;
337      ; number of plots in one row
338      ;
339      ; data type: integer
340      ;
341      ; example: no_columns = 1
342      ;
343      ; default:   1
344      ;***************************************************
345      if(.not. isvar("no_columns"))then 
346     
347     
348         no_columns = 1
349         
350         
351      end if     
352      ;***************************************************
353      ; OPTIONAL --no_rows-- NUMBER OF ROWS
354      ;
355      ; number of plots in one column
356      ;
357      ; data type: integer
358      ;
359      ; example: no_rows = 2
360      ;
361      ; default:   2
362      ;***************************************************
363      if(.not. isvar("no_rows"))then     
364     
365     
366         no_rows = 2
367         
368         
369      end if     
370      ;***************************************************
371      ; OPTIONAL --sort-- TYPE OF SORTING
372      ;
373      ; sequence of plots; sorting either by time step [="time"] or
374      ; by layer/height [="layer"]
375      ;
376      ; data type: string
377      ;
378      ; example: sort = "layer"
379      ;
380      ; default:   "layer"
381      ;***************************************************
382      if(.not. isvar("sort"))then       
383     
384     
385         sort = "layer"
386         
387         
388      end if     
389      ;***************************************************
[190]390      ; OPTIONAL --mode-- TYPE OF PLOT
391      ;
[418]392      ; isoline plots [="Line"], contour plots [="Fill"] or both [="Both"]
393      ; can be choosen
[190]394      ;
395      ; data type: string
396      ;
[218]397      ; example: mode = "Fill"
398      ;
[190]399      ; default:   "Fill"
400      ;***************************************************
401      if(.not. isvar("mode"))then       
402     
403     
404         mode = "Fill"
405         
406         
407      end if     
408      ;***************************************************
409      ; OPTIONAL --fill_mode-- TYPE OF FILLING
410      ;
[418]411      ; area fill [="AreaFill"], raster fill [="RasterFill"] or
412      ; cell fill [="CellFill"] can be choosen
[190]413      ;
414      ; data type: string
415      ;
[218]416      ; example: fill_mode = "AreaFill"
417      ;
[517]418      ; default:    "RasterFill"
[190]419      ;***************************************************
420      if(.not. isvar("fill_mode"))then   
421     
422     
[517]423         fill_mode = "RasterFill"
[190]424         
425         
426      end if     
427      ;***************************************************
428      ; OPTIONAL --shape-- ASPECT RATIO
429      ;
[218]430      ; aspect ratio of axis is kept [=1] or not[=0]
[190]431      ;
432      ; data type: integer
433      ;
[218]434      ; example: shape = 1
435      ;
[190]436      ; default:   1
437      ;***************************************************
438      if(.not. isvar("shape"))then
439     
440         
441         shape = 1
442         
443         
[218]444      end if
[190]445      ;***************************************************
[218]446      ; OPTIONAL --font_size-- FONT SIZE OF STRINGS
447      ;
448      ; font size can be changed to fit to the current plot
449      ;
450      ; data type: float
451      ;
452      ; example: font_size = 0.04
453      ;
454      ; default:   0.04
455      ;***************************************************
456      if(.not. isvar("font_size"))then
457     
458         
459         font_size = 0.04
460         
461         
462      end if
463      ;***************************************************
464      ; OPTIONAL --font_size_legend-- FONT SIZE OF LEGEND STRINGS
465      ;
466      ; font size of the legend can be changed to fit to the current plot
467      ;
468      ; data type: float
469      ;
[517]470      ; example: font_size_legend = 0.04
[218]471      ;
[517]472      ; default:   0.04
[218]473      ;***************************************************
474      if(.not. isvar("font_size_legend"))then
475     
476         
[517]477         font_size_legend = 0.04
[218]478         
479         
480      end if   
481      ;***************************************************
[418]482      ; OPTIONAL --legend_label_stride-- REDUCTION OF THE NUMBER
483      ;                                  OF LABELS IN THE LEGEND
[218]484      ;
[418]485      ; if legend_label_stride is set to 1 every label is shown;
486      ; if set to 2 every second, ...
[218]487      ;
488      ; data type: integer
489      ;
490      ; example: legend_label_stride = 2
491      ;
492      ; default:   2
493      ;***************************************************
494      if(.not. isvar("legend_label_reduction"))then
495     
496         
497         legend_label_stride = 2
498         
499         
500      end if
501      ;***************************************************
[513]502      ; OPTIONAL --axes_explicit--
503      ;          EXPLICIT SETTINGS OF X- AND Y-AXES ARE SWITCHED ON
[218]504      ;
[418]505      ; to use several parameters to set up the axes, 'axes_explicit'
506      ; must be switched on [=1]; if switched off [=0] units of axes are [m]
[218]507      ;
508      ; data type: integer
509      ;
510      ; example: axes_explicit = 0
511      ;
512      ; default:   0
513      ;***************************************************
514      if(.not. isvar("axes_explicit"))then
515     
516         
517         axes_explicit = 0
518         
519         
520      end if
521      ;***************************************************
[418]522      ; OPTIONAL if axes_explicit = 1 --major_ticks_x--
523      ;                              NUMBER OF MAJOR TICK MARKS AT X-AXIS
[218]524      ;
[418]525      ; major_ticks_x can be set to any integer, best segmentation of x-axis
526      ; can be achieved, if
527      ; (number of gripdpoints in x-direction)/(major_ticks_x-1)=INTEGER,
528      ; e.g. 128/8 = 16 -> major_ticks_x should be set to 9
[218]529      ;
530      ; data type: integer
531      ;
532      ; example: major_ticks_x = 6
533      ;
534      ; default:   6
535      ;***************************************************
536      if(.not. isvar("major_ticks_x"))then
537     
538         
539         major_ticks_x = 6
540         
541         
542      end if
543      ;***************************************************
[418]544      ; OPTIONAL if axes_explicit = 1 --major_ticks_y--
545      ;                              NUMBER OF MAJOR TICK MARKS AT Y-AXIS
[218]546      ;
[418]547      ; major_ticks_y can be set to any integer, best segmentation of y-axis
548      ; can be achieved, if
549      ; (number of gripdpoints in y-direction)/(major_ticks_y-1)=INTEGER,
550      ; e.g. 128/8 = 16 -> major_ticks_y should be set to 9
[218]551      ;
552      ; data type: integer
553      ;
554      ; example: major_ticks_y = 6
555      ;
556      ; default:   6
557      ;***************************************************
558      if(.not. isvar("major_ticks_y"))then
559     
560         
561         major_ticks_y = 6
562         
563         
564      end if
565      ;***************************************************
[418]566      ; OPTIONAL if axes_explicit = 1 --norm_x--norm_y--norm_z--
567      ;                              NORMALISE X-, Y- and Z-AXES
[218]568      ;
[513]569      ; normalising the axes with the given value;
570      ; units can explicitly be set;
[418]571      ; by default no normalising is used and units are [meter]
[218]572      ;
573      ; data type: float
574      ;
[418]575      ; default:   1.0
[218]576      ;***************************************************
577      if(.not. isvar("norm_x"))then     
578         
[418]579         norm_x = 1.0   
[218]580         
581      end if
582      if(.not. isvar("norm_y"))then     
583         
[418]584         norm_y = 1.0   
[218]585         
586      end if
587      if(.not. isvar("norm_z"))then     
588         
[418]589         norm_z = 1.0   
[218]590         
591      end if
592      ;***************************************************
[418]593      ; OPTIONAL if axes_explicit = 1 --unit_x--unit_y--unit_z--
594      ;                              UNITS OF X-, Y- and Z-AXES
[218]595      ;
[418]596      ; units of the axes can be set; by default no units are given,
597      ; if norm_x, norm_y or norm_z are unequal 1., otherwise units are [meter]
[218]598      ;
599      ; data type: string,
600      ;
601      ; example: unit_x = "m/s"
602      ;
603      ; default:   " "
604      ;***************************************************
605      if(.not. isvar("unit_x"))then     
606         
607         unit_x = " "   
608         
609      end if
610      if(.not. isvar("unit_y"))then     
611         
612         unit_y = " "   
613         
614      end if
615      if(.not. isvar("unit_z"))then     
616         
617         unit_z = " "   
618         
619      end if
620      ;***************************************************
[190]621      ; OPTIONAL --vector-- VECTOR PLOT
622      ;
[218]623      ; a vector plot for one defined vector is switched on [=1] or off [=0]
[190]624      ;
625      ; data type: integer
626      ;
[218]627      ; example: vector = 0
628      ;
[190]629      ; default:   0
630      ;***************************************************
631      if(.not. isvar("vector"))then
632     
633         
634         vector = 0
635         
636         
637      end if     
638      ;***************************************************
639      ; REQUIRED IF vector=1 --vec1, vec2-- VECTOR FOR VECTOR PLOT
640      ;
[418]641      ; variables for the vector, vec1 hast to be set to component 1 and
642      ; vec2 to component 2;
[218]643      ; inidcate them with one comma before and after each variable;
[190]644      ;
645      ; data type: string
646      ;
[218]647      ; example: vec1=",u,"
648      ;          vec2=",v,"
649      ;
[190]650      ; default:   vec1 = "component1", vec2 = "component2"
651      ;***************************************************
652      if(.not. isvar("vec1"))then     
653         
654         vec1 = "component1"
655         
656      end if
657      if (.not. isvar("vec2"))then
658     
659         vec2 = "component2"     
660         
661      end if     
662      ;***************************************************
663      ; OPTIONAL --plotvec-- VARIABLE FOR EXTRA VECTOR PLOT
664      ;
[418]665      ; the vector plot can be combined with other plots of variables;
666      ; if this is desired inidcate the variables,
667      ; where a vector plot shall overlay; inidcate them with one comma
668      ; before and after each variable
[190]669      ;
670      ; data type: string
671      ;
[218]672      ; example: plotvec = ",u,w,"
673      ;
[513]674      ; default:   "plotvec" (no combining,
675      ;                       the vector plot is created seperately)
[190]676      ;***************************************************
677      if(.not. isvar("plotvec"))then
678     
679         
680         plotvec = "plotvec"
681         
682         
683      end if     
684      ;***************************************************
685      ; OPTIONAL --ref_mag-- REFERENCE VECTOR
686      ;
687      ; length of the reference vector in a vector plot
688      ;
689      ; data type: float
690      ;
[218]691      ; example: ref_mag = 0.05
692      ;
[190]693      ; default:   0.05
694      ;***************************************************
695      if(.not. isvar("ref_mag"))then
696     
697         
698         ref_mag = 0.05
699         
700         
701      end if             
702      ;***************************************************
[194]703
[190]704   end if
705
706   ;******************************************************
707   ; parameters used by PROFILES
708   ;******************************************************
709   
[194]710   if (script .EQ. "profiles")then
711
[190]712      ;***************************************************
[218]713      ; REQUIRED --file_1-- 1ST INPUT FILE
[190]714      ;
[418]715      ; 1st input file produced by PALM; if there is more than one file for
716      ; all timesteps of a job chain, declare path and file name without
717      ; cycle number and ".nc" (first and last cycle numbers are given with
718      ; the parameters 'start_f' and 'end_f');
[218]719      ; otherwise declare path and full file name
[190]720      ;
721      ; data type: string
722      ;
[418]723      ; example: file_1 = "/path_to_file/file_pr"     
724      ;                     if you use more than one cycle number
725      ;          file_1 = "/path_to_file/file_pr.4.nc"
726      ;                     if you use one file (e.g. with cycle number 4)
[218]727      ;
[190]728      ; default:   "File in"
729      ;***************************************************
[218]730      if(.not. isvar("file_1"))then
731
732
[267]733         file_1 = "File in"
[218]734         
735         
736      end if
737      ;***************************************************
[566]738      ; OPTIONAL/REQUIRED --start_f-- if more than one file of a job chain
[513]739      ;                    START OF CYCLIC NUMBER OF THE 1ST INPUT FILE
[218]740      ;
[418]741      ; if there is more than one file for all timesteps of a job chain,
[513]742      ; declare the first cyclicnumber which contains timesteps
743      ; you'd like to use
[218]744      ;
745      ; data type: integer
746      ;
[418]747      ; example:   3 (if first file: example_pr.3.nc) and
748      ;            0 (if first file: example_pr.nc)
[218]749      ;
750      ; default:   -1 
751      ;***************************************************
[566]752      if (.not. isvar("start_f"))then 
[190]753     
[218]754         
[566]755         start_f = -1
[218]756         
757         
758      end if
759      ;***************************************************
[566]760      ; OPTIONAL/REQUIRED --end_f--  if more than one file of a job chain
[513]761      ;            END OF CYCLIC NUMBER OF THE 1ST INPUT FILE
[218]762      ;
[418]763      ; if there is more than one file for all timesteps of a job chain,
[513]764      ; declare the last cyclicnumber which contains timesteps
765      ; you'd like to use
[218]766      ;
767      ; data type: integer
768      ;
[418]769      ; example:   48 (if last file: example_pr.48.nc) and
770      ;             0 (if last file: example_pr.nc)
[218]771      ;
772      ; default:   -1 
773      ;***************************************************
[566]774      if (.not. isvar("end_f"))then 
[190]775     
[218]776         
[566]777         end_f = -1
[190]778         
779         
780      end if
781      ;***************************************************
782      ; OPTIONAL --format_out-- OUTPUT FORMAT
783      ;
[532]784      ; supported file formats: x11, pdf, ps, eps, epsi, ncgm or png
[190]785      ;
786      ; data type: string
787      ;
788      ; default:   "x11"
789      ;***************************************************
790      if (.not. isvar("format_out"))then         
791     
792     
793         format_out = "x11"
794         
795         
796      end if     
797      ;***************************************************
798      ; OPTIONAL --file_out-- OUTPUT FILE
799      ;
[218]800      ; declare path and name of the output file
[190]801      ;
802      ; data type: string
803      ;
[194]804      ; default:   "~/test_pr"
[190]805      ;***************************************************
806      if (.not. isvar("file_out"))then 
807     
808       
809         file_out = "~/test_pr"
810         
811         
812      end if     
813      ;***************************************************
[566]814      ; OPTIONAL --var-- OUTPUT VARIABLES
815      ;
816      ; names of output variables;
817      ; inidcate them with one comma before and after each variable
818      ; if the name contains double quotes (") use dq and the '+' operation to
819      ; concatenate the strings
820      ;
821      ; data type: string
822      ;
823      ; example: var = ",u,v,w,w*pt*,w"+dq+"pt"+dq+",wpt,"
824      ;
[958]825      ; default:   ",cross_profiles," (profiles are plotted together according
826      ;                                to the d3par parameter cross_profiles)
[566]827      ;***************************************************
828      if (.not. isvar("var"))then
829     
830         
[951]831         var = ",cross_profiles,"
[566]832         
833         
834      end if
835      ;***************************************************
836      ; OPTIONAL --start_time_step-- FIRST TIME STEP
837      ;
838      ; first time step of plot; there is no need to specify a precise
839      ; time value from the input file,
840      ; 'start_time_step' is rounded to the next existent time step
841      ;
842      ; data type: double; dimension: [hour]
843      ;
844      ; example:   start_time_step = 0.5d
845      ;
846      ; default:   -1.d (first available time step)
847      ;***************************************************
848      if (.not. isvar("start_time_step"))then   
849     
850       
851         start_time_step = -1.d
852         
853         
854      end if     
855      ;***************************************************
856      ; OPTIONAL --end_time_step-- LAST TIME STEP
857      ;
858      ; last time step of plot; there is no need to specify a precise
859      ; time value from the input file,
860      ; 'end_time_step' is rounded to the next existent time step
861      ;
862      ; data type: double; dimension: [hour]
863      ;
864      ; example:   end_time_step = 2.0d
865      ;
866      ; default:   -1.d (last available time step)
867      ;***************************************************
868      if (.not. isvar("end_time_step"))then
869     
870         
871         end_time_step = -1.d   
872         
873         
874      end if
875      ;***************************************************
876      ; OPTIONAL --time_stride-- TEMPORAL STRIDE
877      ;
878      ; temporal stride for the plots; if time_stride=1 all the profiles
879      ; between first_time_step and end_time_step are plotted
880      ; (including first_time_step and end_time_step),
881      ; if time_stride=3 every third time step between start_time_step and
882      ; end_time_step is plotted
883      ;
884      ; data type: integer
885      ;
886      ; example:   time_stride = 3
887      ;
888      ; default:   1
889      ;***************************************************
890      if (.not. isvar("time_stride"))then       
891     
892       
893         time_stride = 1
894         
895         
896      end if
897      ;***************************************************
898      ; OPTIONAL --xs-- MINIMUM X-AXIS
899      ;
900      ; value for minimum of x-axis
901      ;
902      ; data type: float; dimension of variable
903      ;
904      ; example: xs = 300. (e.g. if var=",pt")
905      ;
906      ; default:   -1. (minimum of range)
907      ;***************************************************
908      if (.not. isvar("xs"))then   
909     
910         
911         xs = -1.
912         
913         
914      end if     
915      ;***************************************************
916      ; OPTIONAL --xe-- MAXIMUM X-AXIS
917      ;
918      ; value for maximum of x-axis
919      ;
920      ; data type: float; dimension of variable
921      ;
922      ; example: xe = 302. (e.g. if var=",pt")
923      ;
924      ; default:   -1. (maximum of range)
925      ;***************************************************
926      if (.not. isvar("xe"))then 
927     
928         
929         xe = -1.
930         
931         
932      end if     
933      ;***************************************************
934      ; OPTIONAL --min_z-- MINIMUM Z-AXIS
935      ;
936      ; value for minimum height of z-axis
937      ;
938      ; data type: double; dimension: [meter]
939      ;
940      ; example: min_z = 0.0d
941      ;
942      ; default:   -1.d (related to the minimum of z-axis in input file)
943      ;***************************************************
944      if (.not. isvar("min_z"))then   
945     
946         
947         min_z = -1.d
948         
949         
950      end if     
951      ;***************************************************
952      ; OPTIONAL --max_z-- MAXIMUM Z-AXIS
953      ;
954      ; value for maximum height of z-axis
955      ;
956      ; data type: double; dimension: [meter]
957      ;
958      ; example: max_z = 1000.0d
959      ;
960      ; default:   -1.d (related to the maximum of z-axis in input file)
961      ;***************************************************
962      if (.not. isvar("max_z"))then   
963     
964         
965         max_z = -1.d
966         
967         
968      end if     
969      ;***************************************************
970      ; OPTIONAL --start_x-- STARTVALUE IN X-DIRECTION 
971      ;
972      ; if 3D data is used for profiles you can choose a startvalue
973      ; in x-direction for horizontal averaging
974      ;
975      ; data type: integer; dimension [grippoint]
976      ;
977      ; example: start_x = 0
978      ;
979      ; default:   0 (minimum of range)
980      ;***************************************************
981      if (.not. isvar("start_x"))then   
982     
983         
984         start_x = 0
985         
986         
987      end if
988      ;***************************************************
989      ; OPTIONAL --end_x-- ENDVALUE IN X-DIRECTION 
990      ;
991      ; if 3D data is used for profiles you can choose an endvalue
992      ; in x-direction for horizontal averaging
993      ;
994      ; data type: integer; dimension [grippoint]
995      ;
996      ; example: end_x = 200
997      ;
998      ; default:   -1 (maximum of range excluding ghostpoint)
999      ;***************************************************
1000      if (.not. isvar("end_x"))then   
1001     
1002         
1003         end_x = -1
1004         
1005         
1006      end if
1007      ;***************************************************
1008      ; OPTIONAL --start_y-- STARTVALUE IN Y-DIRECTION 
1009      ;
1010      ; if 3D data is used for profiles you can choose a startvalue
1011      ; in y-direction for horizontal averaging
1012      ;
1013      ; data type: integer; dimension [grippoint]
1014      ;
1015      ; example: start_y = 0
1016      ;
1017      ; default:   0 (minimum of range)
1018      ;***************************************************
1019      if (.not. isvar("start_y"))then   
1020     
1021         
1022         start_y = 0
1023         
1024         
1025      end if
1026      ;***************************************************
1027      ; OPTIONAL --end_y-- ENDVALUE IN Y-DIRECTION 
1028      ;
1029      ; if 3D data is used for profiles you can choose an endvalue
1030      ; in y-direction for horizontal averaging
1031      ;
1032      ; data type: integer; dimension [grippoint]
1033      ;
1034      ; example: end_y = 200
1035      ;
1036      ; default:   -1 (maximum of range excluding ghostpoint)
1037      ;***************************************************
1038      if (.not. isvar("end_y"))then   
1039     
1040         
1041         end_y = -1
1042         
1043         
1044      end if     
1045      ;***************************************************
[218]1046      ; OPTIONAL --no_columns-- NUMBER OF COLUMNS
[190]1047      ;
1048      ; number of plots in one row
1049      ;
1050      ; data type: integer
1051      ;
[958]1052      ; default:   -1      (profiles are plotted according
1053      ;                     to the d3par parameter profile_columns)
[190]1054      ;***************************************************
1055      if (.not. isvar("no_columns"))then         
1056     
1057     
[951]1058         no_columns = -1
[190]1059         
1060         
1061      end if     
1062      ;***************************************************
[218]1063      ; OPTIONAL --no_rows-- NUMBER OF ROWS
[190]1064      ;
1065      ; number of plots in one column
1066      ;
1067      ; data type: integer
1068      ;
[958]1069      ; default:   -1      (profiles are plotted according
1070      ;                     to the d3par parameter profile_rows)
[190]1071      ;***************************************************
[218]1072      if (.not. isvar("no_rows"))then   
[190]1073     
1074       
[951]1075         no_rows = -1
[190]1076         
1077         
1078      end if     
1079      ;***************************************************
[566]1080      ; OPTIONAL --log_z-- LOGARITHMIC SCALE OF Z-AXIS
[190]1081      ;
[566]1082      ; logarithmic scale for z-axis is switched on [=1] or off [=0]
1083      ;
1084      ; data type: integer
1085      ;
1086      ; example: log_z = 0
1087      ;
1088      ; default:   0
1089      ;***************************************************
1090      if (.not. isvar("log_z"))then
1091     
1092           
1093         log_z = 0
1094         
1095         
1096      end if     
1097      ;***************************************************
1098      ; OPTIONAL --norm_z-- NORMALISING Z-AXIS
1099      ;
1100      ; value for normalising the z-axis
1101      ;
1102      ; data type: float
1103      ;
1104      ; example: norm_z = 1.
1105      ;
1106      ; default:   1. (no normalising)
1107      ;***************************************************
1108      if (.not. isvar("norm_z"))then
1109     
1110           
1111         norm_z = 1.
1112         
1113         
1114      end if     
1115      ;***************************************************
1116      ; OPTIONAL --over-- OVERLAYING
1117      ;
1118      ; predefined overlaying of standard variables is switched
1119      ; on [=1] or off [=0];
1120      ; this feature cannot be used with statistic regions
1121      ;
1122      ; data type: integer
1123      ;
1124      ; example: over = 0
1125      ;
1126      ; default:   0
1127      ;***************************************************
1128      if (.not. isvar("over"))then
1129     
1130         
1131         over = 0
1132         
1133         
1134      end if     
1135      ;***************************************************
1136      ; OPTIONAL --combine-- COMBINING
1137      ;
1138      ; combining of two or three individual variables is possible and
1139      ; switched on [=1] or off[=0]
1140      ;
1141      ; data type: integer
1142      ;
1143      ; example: combine = 0
1144      ;
1145      ; default:   0
1146      ;***************************************************
1147      if (.not. isvar("combine"))then 
1148     
1149         
1150         combine = 0
1151         
1152         
1153      end if     
1154      ;***************************************************
1155      ; REQUIRED IF combine=1 --number_comb-- NUMBER FOR COMBINING
1156      ;
1157      ; number of variables for combining; 2 or 3 is valid
1158      ;
1159      ; data type: integer
1160      ;
1161      ; example: number_comb = 2
1162      ;
1163      ; default:   -1
1164      ;***************************************************
1165      if (.not. isvar("number_comb"))then
1166     
1167           
1168         number_comb = -1
1169         
1170         
1171      end if     
1172      ;***************************************************
1173      ; REQUIRED IF combine=1 --c_var-- VARIABLES FOR COMBINING
1174      ;
1175      ; names of variables;
[418]1176      ; inidcate them with one comma before and after each variable
[526]1177      ; if the name contains double quotes (") use dq and the '+' operation to
1178      ; concatenate the strings
[190]1179      ;
[566]1180      ; data type: string
[190]1181      ;
[566]1182      ; example: c_var = ",u,v,w,w*pt*,w"+dq+"pt"+dq+",wpt,"
1183      ;
1184      ; default:   "c_variables"
[190]1185      ;***************************************************
[566]1186      if (.not. isvar("c_var"))then
[190]1187     
[566]1188           
1189         c_var = "c_variables"
[190]1190         
1191         
[566]1192      end if     
1193      ;***************************************************
1194      ; OPTIONAL --black-- BLACK OR COLOR LINES
1195      ;
1196      ; colored [=0] or black and white [=1] plots
1197      ;
1198      ; data type: integer
1199      ;
1200      ; example: black = 0
1201      ;
1202      ; default:   0
1203      ;***************************************************
1204      if (.not. isvar("black"))then
1205     
1206           
1207         black = 0
1208         
1209         
1210      end if     
1211      ;***************************************************
1212      ; OPTIONAL --dash-- DASHED LINES
1213      ;
1214      ; use of different line patterns [=1] or continous lines [=0]
1215      ; for different time steps
1216      ;
1217      ; data type: integer
1218      ;
1219      ; example: dash = 0
1220      ;
1221      ; default:   0
1222      ;***************************************************
1223      if (.not. isvar("dash"))then     
1224     
1225       
1226         dash = 0     
1227         
1228         
[219]1229      end if
1230      ;***************************************************
[566]1231      ; OPTIONAL --font_size-- FONT SIZE OF STRINGS
1232      ;
1233      ; font size can be changed to fit to the current plot
1234      ;
1235      ; data type: float
1236      ;
1237      ; example: font_size = 0.02
1238      ;
1239      ; default:   0.02
1240      ;***************************************************
1241      if(.not. isvar("font_size"))then
1242     
1243         
1244         font_size = 0.02
1245         
1246         
1247      end if
1248      ;***************************************************
[219]1249      ; OPTIONAL --legend-- SWITCHES ON A LEGEND
1250      ;
1251      ; a legend is switched on [=1] or off [=0]
1252      ;
1253      ; data type: integer
1254      ;
1255      ; example: legend = 1
1256      ;
1257      ; default:   1
1258      ;***************************************************
1259      if (.not. isvar("legend"))then
1260     
1261         
[236]1262         legend = 1
[219]1263         
1264         
[190]1265      end if   
1266      ;***************************************************
[566]1267      ; OPTIONAL --font_size_legend-- FONT SIZE OF LEGEND STRINGS
1268      ;
1269      ; font size of the legend can be changed to fit to the current plot
1270      ;
1271      ; data type: float
1272      ;
[782]1273      ; example: font_size_legend = 0.02
[566]1274      ;
[782]1275      ; default:   0.02
[566]1276      ;***************************************************
1277      if(.not. isvar("font_size_legend"))then
1278     
1279         
[782]1280         font_size_legend = 0.02
[566]1281         
1282         
1283      end if     
1284      ;***************************************************
[190]1285      ; OPTIONAL ---no_files-- NUMBER OF INPUT FILES
1286      ;
[418]1287      ; up to 6 input files of different runs can be used; all input files
1288      ; must contain the same variables and dimensionalties;
[218]1289      ; identical variables are merged in one plot;
[190]1290      ; combining or overlaying of any variables is not feasible;
[418]1291      ; an extra legend (besides the time step legend) is drawn where you
1292      ; can explicitily name the items
[190]1293      ;
1294      ; data type: integer
1295      ;
1296      ; default:   1
1297      ;***************************************************
1298      if (.not. isvar("no_files"))then
1299     
1300         
1301         no_files = 1
1302         
1303         
1304      end if     
1305      ;***************************************************
1306      ; OPTIONAL --name_legend_1-- 1st LEGEND ITEM
1307      ;
1308      ; if no_files > 1: name of 1st item in extra legend
1309      ;
1310      ; data type: string
1311      ;
1312      ; default:   "File 1"
1313      ;***************************************************
1314      if (.not. isvar("name_legend_1"))then
1315     
1316           
1317         name_legend_1 = "File 1"
1318         
1319         
[218]1320      end if
[190]1321      ;***************************************************
[218]1322      ; REQUIRED if no_files > 1 --file_2-- 2ND INPUT FILE
[190]1323      ;
[418]1324      ; 2nd input file produced by PALM; if there is more than one file
1325      ; for all timesteps of a job chain,
1326      ; declare path and file name without cycle number and ".nc" (first and
[513]1327      ; last cycle numbers are given with the parameters
1328      ; 'start_f' and 'end_f');
[218]1329      ; otherwise declare path and full file name
[190]1330      ;
1331      ; data type: string
1332      ;
[418]1333      ; example: file_2 = "/path_to_file/file_pr"     
1334      ;                    if you use more than one cycle number
1335      ;          file_2 = "/path_to_file/file_pr.4.nc"
1336      ;                    if you use one file (e.g. with cycle number 4)
[218]1337      ;
[190]1338      ; default:   "File in"
1339      ;***************************************************
[218]1340      if(.not. isvar("file_2"))then
1341
1342
1343         file_2 = "File in"
1344         
1345         
1346      end if
1347      ;***************************************************
[513]1348      ; OPTIONAL/REQUIRED --start_f_2-- if more than one file of a job chain
1349      ;                       START OF CYCLIC NUMBER OF THE 2ND INPUT FILE
[218]1350      ;
[418]1351      ; if there is more than one file for all timesteps of a job chain,
[513]1352      ; declare the first cyclicnumber which contains timesteps
1353      ; you'd like to use
[218]1354      ;
1355      ; data type: integer
1356      ;
1357      ; example:   3 (if first file: example_ts.3.nc)
1358      ;
1359      ; default:   -1 
1360      ;***************************************************
1361      if (.not. isvar("start_f_2"))then 
[190]1362     
1363         
[218]1364         start_f_2 = -1
[190]1365         
1366         
[218]1367      end if
[190]1368      ;***************************************************
[513]1369      ; REQUIRED --end_f_2-- if more than one file of a job chain
1370      ;                      END OF CYCLIC NUMBER OF THE 2ND INPUT FILE
[218]1371      ;
[418]1372      ; if there is more than one file for all timesteps of a job chain,
[513]1373      ; declare the last cyclicnumber which contains timesteps
1374      ; you'd like to use
[218]1375      ;
1376      ; data type: integer
1377      ;
1378      ; example:   48 (if last file: example_ts.48.nc)
1379      ;
1380      ; default:   -1 
1381      ;***************************************************
1382      if (.not. isvar("end_f_2"))then 
1383     
1384         
1385         end_f_2 = -1
1386         
1387         
1388      end if             
1389      ;***************************************************
[190]1390      ; OPTIONAL --name_legend_2-- 2nd LEGEND ITEM
1391      ;
1392      ; if no_files > 1: name of 2nd item in extra legend
1393      ;
1394      ; data type: string
1395      ;
1396      ; default:   "File 2"
1397      ;***************************************************
1398      if (.not. isvar("name_legend_2"))then   
1399     
1400         
1401         name_legend_2 = "File 2"
1402         
1403         
[218]1404      end if
[190]1405      ;***************************************************
[218]1406      ; REQUIRED if no_files > 2 --file_3-- 3RD INPUT FILE
[190]1407      ;
[418]1408      ; 3rd input file produced by PALM; if there is more than one file
1409      ; for all timesteps of a job chain,
1410      ; declare path and file name without cycle number and ".nc" (first and
[513]1411      ; last cycle numbers are given with the parameters
1412      ; 'start_f' and 'end_f');
[218]1413      ; otherwise declare path and full file name
[190]1414      ;
1415      ; data type: string
1416      ;
[418]1417      ; example: file_3 = "/path_to_file/file_pr"     
1418      ;                     if you use more than one cycle number
1419      ;          file_3 = "/path_to_file/file_pr.4.nc"
1420      ;                     if you use one file (e.g. with cycle number 4)
[218]1421      ;
[190]1422      ; default:   "File in"
1423      ;***************************************************
[218]1424      if(.not. isvar("file_3"))then
1425
1426
1427         file_3 = "File in"
[190]1428         
1429         
[218]1430      end if
[190]1431      ;***************************************************
[513]1432      ; OPTIONAL/REQUIRED --start_f_3-- if more than one file of a job chain
1433      ;                        START OF CYCLIC NUMBER OF THE 3RD INPUT FILE
[190]1434      ;
[418]1435      ; if there is more than one file for all timesteps of a job chain,
[513]1436      ; declare the first cyclicnumber which contains timesteps
1437      ; you'd like to use
[190]1438      ;
[218]1439      ; data type: integer
[190]1440      ;
[218]1441      ; example:   3 (if first file: example_ts.3.nc)
1442      ;
1443      ; default:   -1 
[190]1444      ;***************************************************
[218]1445      if (.not. isvar("start_f_3"))then 
[190]1446     
[218]1447         
1448         start_f_3 = -1
[190]1449         
1450         
[218]1451      end if
[190]1452      ;***************************************************
[513]1453      ; OPTIONAL/REQUIRED --end_f_3-- if more than one file of a job chain
1454      ;                    END OF CYCLIC NUMBER OF THE 3RD INPUT FILE
[190]1455      ;
[513]1456      ; if there is more than one file for all timesteps of a job chain,
1457      ; declare
[418]1458      ; the last cyclicnumber which contains timesteps you'd like to use
[190]1459      ;
[218]1460      ; data type: integer
[190]1461      ;
[218]1462      ; example:   48 (if last file: example_ts.48.nc)
1463      ;
1464      ; default:   -1 
[190]1465      ;***************************************************
[218]1466      if (.not. isvar("end_f_3"))then 
[190]1467     
[218]1468         
1469         end_f_3 = -1
[190]1470         
1471         
[218]1472      end if   
[190]1473      ;***************************************************
[218]1474      ; OPTIONAL --name_legend_3-- 3rd LEGEND ITEM
[190]1475      ;
[218]1476      ; if no_files > 2: name of 3rd item in extra legend
[190]1477      ;
1478      ; data type: string
1479      ;
[218]1480      ; default:   "File 3"
[190]1481      ;***************************************************
[218]1482      if (.not. isvar("name_legend_3"))then
[190]1483     
[218]1484           
1485         name_legend_3 = "File 3"
[190]1486         
1487         
[218]1488      end if
[190]1489      ;***************************************************
[218]1490      ; REQUIRED if no_files > 3 --file_4-- 4TH INPUT FILE
[190]1491      ;
[418]1492      ; 4th input file produced by PALM; if there is more than one file for
1493      ; all timesteps of a job chain,
1494      ; declare path and file name without cycle number and ".nc" (first and
[513]1495      ; last cycle numbers are given with the parameters
1496      ; 'start_f' and 'end_f');
[218]1497      ; otherwise declare path and full file name
[190]1498      ;
1499      ; data type: string
1500      ;
[418]1501      ; example: file_4 = "/path_to_file/file_pr"     
1502      ;                     if you use more than one cycle number
1503      ;          file_4 = "/path_to_file/file_pr.4.nc"
1504      ;                      if you use one file (e.g. with cycle number 4)
[218]1505      ;
[190]1506      ; default:   "File in"
1507      ;***************************************************
[218]1508      if(.not. isvar("file_4"))then
1509
1510
1511         file_4 = "File in"
[190]1512         
1513         
[218]1514      end if
[190]1515      ;***************************************************
[513]1516      ; OPTIONAL/REQUIRED --start_f_4-- if more than one file of a job chain
1517      ;                          START OF CYCLIC NUMBER OF THE 4TH INPUT FILE
[190]1518      ;
[418]1519      ; if there is more than one file for all timesteps of a job chain,
[513]1520      ; declare the first cyclicnumber which contains timesteps
1521      ; you'd like to use
[190]1522      ;
[218]1523      ; data type: integer
[190]1524      ;
[218]1525      ; example:   3 (if first file: example_ts.3.nc)
1526      ;
1527      ; default:   -1 
[190]1528      ;***************************************************
[218]1529      if (.not. isvar("start_f_4"))then 
[190]1530     
[218]1531         
1532         start_f_4 = -1
[190]1533         
1534         
[218]1535      end if
[190]1536      ;***************************************************
[513]1537      ; OPTIONAL/REQUIRED --end_f_4-- if more than one file of a job chain
1538      ;                   END OF CYCLIC NUMBER OF THE 4TH INPUT FILE
[190]1539      ;
[418]1540      ; if there is more than one file for all timesteps of a job chain,
[513]1541      ; declare the last cyclicnumber which contains timesteps
1542      ; you'd like to use
[190]1543      ;
[218]1544      ; data type: integer
[190]1545      ;
[218]1546      ; example:   48 (if last file: example_ts.48.nc)
1547      ;
1548      ; default:   -1 
[190]1549      ;***************************************************
[218]1550      if (.not. isvar("end_f_4"))then 
[190]1551     
[218]1552         
1553         end_f_4 = -1
[190]1554         
1555         
[218]1556      end if 
[190]1557      ;***************************************************
[218]1558      ; OPTIONAL --name_legend_4-- 4th LEGEND ITEM
[190]1559      ;
[218]1560      ; if no_files > 3: name of 4th item in extra legend
[190]1561      ;
1562      ; data type: string
1563      ;
[218]1564      ; default:   "File 4"
[190]1565      ;***************************************************
[218]1566      if (.not. isvar("name_legend_4"))then   
[190]1567     
[218]1568       
1569         name_legend_4 = "File 4"
[190]1570         
1571         
[218]1572      end if
[190]1573      ;***************************************************
[218]1574      ; REQUIRED if no_files > 4 --file_5-- 5TH INPUT FILE
[190]1575      ;
[418]1576      ; 5th input file produced by PALM; if there is more than one file for
1577      ; all timesteps of a job chain,
1578      ; declare path and file name without cycle number and ".nc" (first and
[513]1579      ; last cycle numbers are given with the parameters
1580      ; 'start_f' and 'end_f');
[218]1581      ; otherwise declare path and full file name
[190]1582      ;
1583      ; data type: string
1584      ;
[418]1585      ; example: file_5 = "/path_to_file/file_pr"     
1586      ;                     if you use more than one cycle number
1587      ;          file_5 = "/path_to_file/file_pr.4.nc"
1588      ;                     if you use one file (e.g. with cycle number 4)
[218]1589      ;
[190]1590      ; default:   "File in"
1591      ;***************************************************
[218]1592      if(.not. isvar("file_5"))then
1593
1594
1595         file_5 = "File in"
[190]1596         
1597         
[218]1598      end if
[190]1599      ;***************************************************
[513]1600      ; OPTIONAL/REQUIRED --start_f_5-- if more than one file of a job chain
1601      ;                     START OF CYCLIC NUMBER OF THE 5TH INPUT FILE
[190]1602      ;
[418]1603      ; if there is more than one file for all timesteps of a job chain,
[513]1604      ; declare the first cyclicnumber which contains timesteps
1605      ; you'd like to use
[190]1606      ;
[218]1607      ; data type: integer
[190]1608      ;
[218]1609      ; example:   3 (if first file: example_ts.3.nc)
1610      ;
1611      ; default:   -1 
[190]1612      ;***************************************************
[218]1613      if (.not. isvar("start_f_5"))then 
[190]1614     
[218]1615         
1616         start_f_5 = -1
[190]1617         
1618         
[218]1619      end if
[190]1620      ;***************************************************
[513]1621      ; OPTIONAL/REQUIRED --end_f_5-- if more than one file of a job chain
1622      ;END OF CYCLIC NUMBER OF THE 5TH INPUT FILE
[190]1623      ;
[418]1624      ; if there is more than one file for all timesteps of a job chain,
[513]1625      ; declare the last cyclicnumber which contains timesteps
1626      ; you'd like to use
[190]1627      ;
[218]1628      ; data type: integer
[190]1629      ;
[218]1630      ; example:   48 (if last file: example_ts.48.nc)
1631      ;
1632      ; default:   -1 
[190]1633      ;***************************************************
[218]1634      if (.not. isvar("end_f_5"))then 
[190]1635     
[218]1636         
1637         end_f_5 = -1
[190]1638         
1639         
[218]1640      end if     
[190]1641      ;***************************************************
[218]1642      ; OPTIONAL --name_legend_5-- 5th LEGEND ITEM
[190]1643      ;
[218]1644      ; if no_files > 4: name of 5th item in extra legend
[190]1645      ;
1646      ; data type: string
1647      ;
[218]1648      ; default:   "File 5"
[190]1649      ;***************************************************
[218]1650      if (.not. isvar("name_legend_5"))then
[190]1651     
[218]1652           
1653         name_legend_5 = "File 5"
[190]1654         
1655         
[218]1656      end if
[190]1657      ;***************************************************
[218]1658      ; REQUIRED if no_files > 5 --file_6-- 6TH INPUT FILE
[190]1659      ;
[418]1660      ; 6th input file produced by PALM; if there is more than one file for
1661      ; all timesteps of a job chain,
1662      ; declare path and file name without cycle number and ".nc" (first and
[513]1663      ; last cycle numbers are given with the parameters
1664      ; 'start_f' and 'end_f');
[218]1665      ; otherwise declare path and full file name
[190]1666      ;
1667      ; data type: string
1668      ;
[418]1669      ; example: file_6 = "/path_to_file/file_pr"     
1670      ;                     if you use more than one cycle number
1671      ;          file_6 = "/path_to_file/file_pr.4.nc"
1672      ;                     if you use one file (e.g. with cycle number 4)
[218]1673      ;
[190]1674      ; default:   "File in"
1675      ;***************************************************
[218]1676      if(.not. isvar("file_6"))then
1677
1678
1679         file_6 = "File in"
[190]1680         
1681         
[218]1682      end if
[190]1683      ;***************************************************
[513]1684      ; OPTIONAL/REQUIRED --start_f_6-- if more than one file of a job chain
1685      ;                     START OF CYCLIC NUMBER OF THE 6TH INPUT FILE
[190]1686      ;
[418]1687      ; if there is more than one file for all timesteps of a job chain,
[513]1688      ; declare the first cyclicnumber which contains timesteps
1689      ; you'd like to use
[190]1690      ;
[218]1691      ; data type: integer
[190]1692      ;
[218]1693      ; example:   3 (if first file: example_ts.3.nc)
1694      ;
1695      ; default:   -1 
[190]1696      ;***************************************************
[218]1697      if (.not. isvar("start_f_6"))then 
[190]1698     
[218]1699         
1700         start_f_6 = -1
[190]1701         
1702         
[218]1703      end if
[190]1704      ;***************************************************
[513]1705      ; OPTIONAL/REQUIRED --end_f_6-- if more than one file of a job chain
1706      ; END OF CYCLIC NUMBER OF THE 6TH INPUT FILE
[190]1707      ;
[418]1708      ; if there is more than one file for all timesteps of a job chain,
[513]1709      ; declare the last cyclicnumber which contains timesteps
1710      ; you'd like to use
[190]1711      ;
[218]1712      ; data type: integer
[190]1713      ;
[218]1714      ; example:   48 (if last file: example_ts.48.nc)
1715      ;
1716      ; default:   -1 
[190]1717      ;***************************************************
[218]1718      if (.not. isvar("end_f_6"))then 
[190]1719     
[218]1720         
1721         end_f_6 = -1
[190]1722         
1723         
[218]1724      end if                     
[190]1725      ;***************************************************
[218]1726      ; OPTIONAL --name_legend_6-- 6TH LEGEND ITEM
[190]1727      ;
[218]1728      ; if no_files > 5: name of 6th item in extra legend
[190]1729      ;
1730      ; data type: string
1731      ;
[218]1732      ; default:   "File 6"
[190]1733      ;***************************************************
[218]1734      if (.not. isvar("name_legend_6"))then
[190]1735     
[218]1736           
1737         name_legend_6 = "File 6"
[190]1738         
1739         
[218]1740      end if             
[190]1741      ;***************************************************
1742   end if
[194]1743
[190]1744   ;******************************************************
1745   ; parameters used by SPECTRA
1746   ;******************************************************
1747   
[194]1748   if (script .EQ. "spectra") then
1749
[190]1750      ;***************************************************
1751      ; REQUIRED --file_1-- INPUT FILE
1752      ;
[418]1753      ; input file produced by PALM; if there is more than one file for
1754      ; all timesteps of a job chain,
1755      ; declare path and file name without cycle number and ".nc" (first and
[513]1756      ; last cycle numbers are given with the parameters
1757      ; 'start_f' and 'end_f');
[218]1758      ; otherwise declare path and full file name
[190]1759      ;
1760      ; data type: string
1761      ;
[418]1762      ; example: file_1 = "/path_to_file/file_sp"     
1763      ;                     if you use more than one cycle number
1764      ;          file_1 = "/path_to_file/file_sp.4.nc"
1765      ;                     if you use one file (e.g. with cycle number 4)
[218]1766      ;
[190]1767      ; default:   "File in"
1768      ;***************************************************
[218]1769      if(.not. isvar("file_1"))then
1770
1771
1772         file_1 = "File in"
1773         
1774         
1775      end if
1776      ;***************************************************
[513]1777      ; OPTIONAL/REQUIRED if more than one file of a job chain
[418]1778      ;           --start_f-- START OF CYCLIC NUMBER
[218]1779      ;
[418]1780      ; if there is more than one file for all timesteps of a job chain,
[513]1781      ; declare the first cyclicnumber which contains timesteps
1782      ; you'd like to use
[218]1783      ;
1784      ; data type: integer
1785      ;
[418]1786      ; example:   start_f = 3 (if first file: example_sp.3.nc) and
1787      ;            start_f = 0 (if first file: example_sp.nc)
[218]1788      ;
1789      ; default:   -1 
1790      ;***************************************************
1791      if (.not. isvar("start_f"))then 
[190]1792     
1793         
[218]1794         start_f = -1
[190]1795         
1796         
[218]1797      end if
[190]1798      ;***************************************************
[513]1799      ; OPTIONAL/REQUIRED if more than one file of a job chain
[418]1800      ;          --end_f-- END OF CYCLIC NUMBER
[218]1801      ;
[418]1802      ; if there is more than one file for all timesteps of a job chain,
[513]1803      ; declare the last cyclicnumber which contains timesteps
1804      ; you'd like to use
[218]1805      ;
1806      ; data type: integer
1807      ;
[418]1808      ; example:   end_f = 48 (if last file: example_sp.48.nc) and
1809      ;            end_f =  0 (if last file: example_sp.nc)
[218]1810      ;
1811      ; default:   -1 
1812      ;***************************************************
1813      if (.not. isvar("end_f"))then 
1814     
1815         
1816         end_f = -1
1817         
1818         
1819      end if
1820      ;***************************************************
[190]1821      ; OPTIONAL --format_out-- OUTPUT FORMAT
1822      ;
[418]1823      ; format of the output file;
[532]1824      ; supported file formats: x11, pdf, ps, eps, epsi, ncgm or png
[190]1825      ;
1826      ; data type: string
1827      ;
[218]1828      ; example: format_out = "pdf"
1829      ;
[190]1830      ; default:   "x11"
1831      ;***************************************************
1832      if (.not. isvar("format_out"))then       
1833     
1834       
1835         format_out = "x11"
1836         
1837         
1838      end if     
1839      ;***************************************************
1840      ; OPTIONAL --file_out-- OUTPUT FILE
1841      ;
[218]1842      ; path and name of the output file
[190]1843      ;
1844      ; data type: string
1845      ;
[218]1846      ; example: file_out = "/home/path_to_output/test_sp"
1847      ;
[194]1848      ; default:   "~/test_sp"
[190]1849      ;***************************************************
1850      if (.not. isvar("file_out"))then 
1851     
1852       
1853         file_out = "~/test_sp"
1854         
1855         
[566]1856      end if   
[190]1857      ;***************************************************
[566]1858      ; OPTIONAL --var-- OUTPUT VARIABLES
[190]1859      ;
[566]1860      ; names of output variables;
1861      ; inidcate them with one comma before and after each variable
1862      ; if the name contains double quotes (") use dq and the '+' operation to
1863      ; concatenate the strings
[190]1864      ;
[566]1865      ; data type: string
[218]1866      ;
[566]1867      ; example: var = ",u,v,w,"
[190]1868      ;
[566]1869      ; default:   "all" (all available variables are plotted)
[190]1870      ;***************************************************
[566]1871      if (.not. isvar("var"))then
[190]1872     
[566]1873         
1874         var = "all"
[190]1875         
1876         
1877      end if     
1878      ;***************************************************
[566]1879      ; OPTIONAL --start_time_step-- FIRST TIME STEP
[190]1880      ;
[566]1881      ; first time step of plot; there is no need to specify a precise
1882      ; time value from the input file,
1883      ; 'start_time_step' is rounded to the next existent time step
[190]1884      ;
[566]1885      ; data type: double; dimension: [hour]
[190]1886      ;
[566]1887      ; example:   start_time_step = 0.5d
1888      ;
1889      ; default:   -1.d (first available time step)
[190]1890      ;***************************************************
[566]1891      if (.not. isvar("start_time_step"))then
[190]1892     
1893         
[566]1894         start_time_step = -1.d
[190]1895         
1896         
1897      end if     
1898      ;***************************************************
[566]1899      ; OPTIONAL --end_time_step-- LAST TIME STEP
[190]1900      ;
[566]1901      ; last time step of plot; there is no need to specify a precise
1902      ; time value from the input file,
1903      ; 'end_time_step' is rounded to the next existent time step
[190]1904      ;
[566]1905      ; data type: double; dimension: [hour]
[190]1906      ;
[566]1907      ; example:   end_time_step = 2.0d
1908      ;
1909      ; default:   -1.d (last available time step)
[190]1910      ;***************************************************
[566]1911      if (.not. isvar("end_time_step"))then
[190]1912     
1913         
[566]1914         end_time_step = -1.d
[190]1915         
1916         
[566]1917      end if   
[190]1918      ;***************************************************
1919      ; OPTIONAL --height_level-- HEIGHT LEVELS
1920      ;
[418]1921      ; defines if all levels or selected levels are output;
1922      ; a 1D array with up to 100 elements
[218]1923      ; from 0 to 99 must be indicated
[190]1924      ;
1925      ; data type: integer array
1926      ;
[218]1927      ; example: height_level=(/0,3,7/)
1928      ;
1929      ; default:   (/-1/) (all levels are output)
[190]1930      ;***************************************************
1931      if (.not. isvar("height_level"))then
1932     
1933         
1934         height_level = (/-1/)
1935         
1936         
1937      end if     
1938      ;***************************************************
1939      ; OPTIONAL --sort-- TYPE OF SORTING
1940      ;
[418]1941      ; defines if all heights [="height"] or all timesteps [="time"]
1942      ; are merged in one plot
[190]1943      ;
1944      ; data type: string
1945      ;
[218]1946      ; example: sort = "height"
1947      ;
[190]1948      ; default:   "height"
1949      ;***************************************************
1950      if (.not. isvar("sort"))then
1951     
1952         
1953         sort = "height"
1954         
1955         
1956      end if     
1957      ;***************************************************
[566]1958      ; OPTIONAL --no_columns-- NUMBER OF COLUMNS
[190]1959      ;
[566]1960      ; number of plots in one row
[190]1961      ;
[566]1962      ; data type: integer
1963      ;
1964      ; example: no_columns = 1
[190]1965      ;
[566]1966      ; default:   1
[190]1967      ;***************************************************
[566]1968      if (.not. isvar("no_columns"))then         
[190]1969     
[566]1970     
1971         no_columns = 1
[190]1972         
1973         
1974      end if     
1975      ;***************************************************
[566]1976      ; OPTIONAL --no_rows-- NUMBER OF ROWS
[190]1977      ;
[566]1978      ; number of plots in one column
[190]1979      ;
[566]1980      ; data type: integer
1981      ;
1982      ; example: no_rows = 2
[190]1983      ;
[566]1984      ; default:   2
[190]1985      ;***************************************************
[566]1986      if (.not. isvar("no_rows"))then
[190]1987     
1988         
[566]1989         no_rows = 2
[190]1990         
1991         
1992      end if     
1993      ;***************************************************
[218]1994      ; OPTIONAL --log_x-- LOGARITHMIC SCALE X-AXIS
[190]1995      ;
[218]1996      ; logarithmic scale for x-axis is switched on [=1] or off [=0]
[190]1997      ;
1998      ; data type: integer
1999      ;
[218]2000      ; example: log_x = 1
2001      ;
[190]2002      ; default:   1
2003      ;***************************************************
[218]2004      if (.not. isvar("log_x"))then   
[190]2005     
2006         
[218]2007         log_x = 1
[190]2008         
2009         
2010      end if     
2011      ;***************************************************
[218]2012      ; OPTIONAL --log_y-- LOGARITHMIC SCALE Y-AXIS
[190]2013      ;
[218]2014      ; logarithmic scale for y-axis is switched on [=1] or off [=0]
[190]2015      ;
2016      ; data type: integer
2017      ;
[218]2018      ; example: log_y = 1
2019      ;
[190]2020      ; default:   1
2021      ;***************************************************
[218]2022      if (.not. isvar("log_y"))then   
[190]2023     
2024         
[218]2025         log_y = 1
[190]2026         
2027         
2028      end if     
2029      ;***************************************************
[566]2030      ; OPTIONAL --norm_height-- NORMALISING X-AXIS WITH HEIGHT
[190]2031      ;
[566]2032      ; normalising the x-axis with the height of each spectra
2033      ; is switched on [=1] or off [=0]
[190]2034      ;
[566]2035      ; data type: integer
[190]2036      ;
[566]2037      ; example: norm_height = 0
[218]2038      ;
[566]2039      ; default:   0 (no normalising)
[190]2040      ;***************************************************
[566]2041      if (.not. isvar("norm_height"))then   
[190]2042     
2043         
[566]2044         norm_height = 0
[190]2045         
2046         
[566]2047      end if     
[218]2048      ;***************************************************
[566]2049      ; OPTIONAL --norm_x-- NORMALISING X-AXIS
[218]2050      ;
[566]2051      ; value for normalising the x-axis
[218]2052      ;
[566]2053      ; data type: float
[218]2054      ;
[566]2055      ; example: norm_x = 1.
[218]2056      ;
[566]2057      ; default:   1. (no normalising)
[218]2058      ;***************************************************
[566]2059      if (.not. isvar("norm_x"))then   
[218]2060     
2061         
[566]2062         norm_x = 1.
[218]2063         
2064         
[566]2065      end if
[190]2066      ;***************************************************
[218]2067      ; OPTIONAL --norm_y-- NORMALISE Y-AXIS
[190]2068      ;
2069      ; value for normalising the y-axis
2070      ;
2071      ; data type: float
2072      ;
[218]2073      ; example: norm_y = 1.
2074      ;
[190]2075      ; default:   1. (no normalising)
2076      ;***************************************************
[218]2077      if (.not. isvar("norm_y"))then     
[190]2078     
2079       
[218]2080         norm_y = 1.
[190]2081         
2082         
[218]2083      end if
[190]2084      ;***************************************************
[218]2085      ; OPTIONAL --unit_x--unit_y-- UNITS OF X-, Y-AXES
2086      ;
[513]2087      ; unit of the x-axis can be set if norm_x is unequal 1,
2088      ; otherwise the unit is [1/m]; by default no unit is set for the y-axis
[218]2089      ;
2090      ; data type: string,
2091      ;
2092      ; example: unit_x = "100/m"
2093      ;
2094      ; default:   " "
2095      ;***************************************************
2096      if(.not. isvar("unit_x"))then     
2097         
2098         unit_x = " "   
2099         
2100      end if
2101      if(.not. isvar("unit_y"))then     
2102         
2103         unit_y = " "   
2104         
2105      end if
2106      ;***************************************************
[566]2107      ; OPTIONAL --black-- BLACK OR COLOR LINES
2108      ;
2109      ; colored [=0] or black and white [=1] plots
2110      ;
2111      ; data type: integer
2112      ;
2113      ; example: black = 0
2114      ;
2115      ; default:   0
2116      ;***************************************************
2117      if (.not. isvar("black"))then
2118     
2119           
2120         black = 0
2121         
2122         
2123      end if     
2124      ;***************************************************
2125      ; OPTIONAL --dash-- DASHED LINES
2126      ;
2127      ; use of different line patterns [=1] or continous lines [=0]
2128      ;
2129      ; data type: integer
2130      ;
2131      ; example: dash = 0
2132      ;
2133      ; default:   0
2134      ;***************************************************
2135      if (.not. isvar("dash"))then   
2136     
2137         
2138         dash = 0
2139         
2140         
2141      end if
2142      ;***************************************************
[218]2143      ; OPTIONAL --font_size-- FONT SIZE OF STRINGS
2144      ;
2145      ; font size can be changed to fit to the current plot
2146      ;
2147      ; data type: float
2148      ;
2149      ; example: font_size = 0.02
2150      ;
2151      ; default:   0.02
2152      ;***************************************************
2153      if(.not. isvar("font_size"))then
2154     
2155         
2156         font_size = 0.02
2157         
2158         
2159      end if
2160      ;***************************************************
2161      ; OPTIONAL --font_size_legend-- FONT SIZE OF LEGEND STRINGS
2162      ;
2163      ; font size of the legend can be changed to fit to the current plot
2164      ;
2165      ; data type: float
2166      ;
[534]2167      ; example: font_size_legend = 0.02
[218]2168      ;
[534]2169      ; default:   0.02
[218]2170      ;***************************************************
2171      if(.not. isvar("font_size_legend"))then
2172     
2173         
[534]2174         font_size_legend = 0.02
[218]2175         
2176         
2177      end if     
2178      ;***************************************************
[190]2179     
2180   end if
[194]2181
[190]2182   ;******************************************************
2183   ; parameters used by TIMESERIES
2184   ;******************************************************
2185
[194]2186   if (script .EQ. "timeseries") then
2187
[190]2188      ;***************************************************
2189      ; REQUIRED --file_1-- INPUT FILE
2190      ;
[418]2191      ; input file produced by PALM; if there is more than one file for
2192      ; all timesteps of a job chain,
2193      ; declare path and file name without cycle number and ".nc" (first and
[513]2194      ; last cycle numbers are given with the parameters
2195      ; 'start_f' and 'end_f');
[218]2196      ; otherwise declare path and full file name
[190]2197      ;
2198      ; data type: string
2199      ;
[418]2200      ; example: file_1 = "/path_to_file/file_ts"     
2201      ;                      if you use more than one cycle number
2202      ;          file_1 = "/path_to_file/file_ts.4.nc"
2203      ;                      if you use one file (e.g. with cycle number 4)
[218]2204      ;
[190]2205      ; default:   "File in"
2206      ;***************************************************
[218]2207      if(.not. isvar("file_1"))then
2208
2209
[267]2210         file_1 = "File in"
[218]2211         
2212         
2213      end if
2214      ;***************************************************
[513]2215      ; OPTIONAL/REQUIRED if more than one file of a job chain
[418]2216      ;          --start_f-- START OF CYCLIC NUMBER
[218]2217      ;
[418]2218      ; if there is more than one file for all timesteps of a job chain,
[513]2219      ; declare the first cyclicnumber which contains timesteps
2220      ; you'd like to use
[218]2221      ;
2222      ; data type: integer
2223      ;
[418]2224      ; example:   start_f = 3 (if first file: example_ts.3.nc) and
2225      ;            start_f = 0 (if first file: example_ts.nc)
[218]2226      ;
2227      ; default:   -1 
2228      ;***************************************************
2229      if (.not. isvar("start_f"))then 
[190]2230     
2231         
[218]2232         start_f = -1
[190]2233         
2234         
[218]2235      end if
[190]2236      ;***************************************************
[513]2237      ; OPTIONAL/REQUIRED if more than one file of a job chain
[418]2238      ;           --end_f-- END OF CYCLIC NUMBER
[218]2239      ;
[418]2240      ; if there is more than one file for all timesteps of a job chain,
[513]2241      ; declare the last cyclicnumber which contains timesteps
2242      ; you'd like to use
[218]2243      ;
2244      ; data type: integer
2245      ;
[418]2246      ; example:   end_f = 48 (if last file: example_ts.48.nc) and
2247      ;            end_f = 0 (if last file: example_ts.nc)
[218]2248      ;
2249      ; default:   -1 
2250      ;***************************************************
2251      if (.not. isvar("end_f"))then 
2252     
2253         
2254         end_f = -1
2255         
2256         
2257      end if
2258      ;***************************************************
[190]2259      ; OPTIONAL --format_out-- OUTPUT FORMAT
2260      ;
[418]2261      ; format of the output file;
[532]2262      ; supported file formats: x11, pdf, ps, eps, epsi, ncgm or png
[190]2263      ;
2264      ; data type: string
2265      ;
[218]2266      ; example: format_out = "pdf"
2267      ;
[190]2268      ; default:   "x11"
2269      ;***************************************************
2270      if (.not. isvar("format_out"))then         
2271     
2272     
2273         format_out = "x11"
2274         
2275         
2276      end if     
2277      ;***************************************************
2278      ; OPTIONAL --file_out-- OUTPUT FILE
2279      ;
[218]2280      ; path and name of the output file
[190]2281      ;
2282      ; data type: string
2283      ;
[218]2284      ; example: file_out = "/home/path_to_output/test_ts"
2285      ;
[194]2286      ; default:   "~/test_ts"
[190]2287      ;***************************************************
2288      if (.not. isvar("file_out"))then 
2289     
2290       
2291         file_out = "~/test_ts"
2292         
2293         
[566]2294      end if
2295 ;***************************************************
[190]2296      ; OPTIONAL --var-- OUTPUT VARIABLES
2297      ;
[526]2298      ; names of output variables;
2299      ; inidcate them with one comma before and after each variable
2300      ; if the name contains double quotes (") use dq and the '+' operation to
2301      ; concatenate the strings
[190]2302      ;
2303      ; data type: string
[218]2304      ;
[526]2305      ; example: var = ",umax,vmax,wmax,w"+dq+"pt"+dq+"0,"
[190]2306      ;
[218]2307      ; default:   "all" (all available variables are plotted)
[190]2308      ;***************************************************
2309      if (.not. isvar("var"))then
2310     
2311         
2312         var = "all"
2313         
2314         
2315      end if     
2316      ;***************************************************
2317      ; OPTIONAL --start_time_step-- FIRST TIME STEP
2318      ;
[418]2319      ; first time step of plot; there is no need to specify a precise
2320      ; time value from the input file,
[218]2321      ; 'start_time_step' is rounded to the next existent time step
[190]2322      ;
2323      ; data type: double; dimension: [hour]
2324      ;
[418]2325      ; example:   start_time_step = 0.5d
[218]2326      ;
[190]2327      ; default:   -1.d (first available time step)
2328      ;***************************************************
2329      if (.not. isvar("start_time_step"))then
2330     
2331         
2332         start_time_step = -1.d
2333         
2334         
2335      end if     
2336      ;***************************************************
2337      ; OPTIONAL --end_time_step-- LAST TIME STEP
2338      ;
[418]2339      ; last time step of plot; there is no need to specify a precise
2340      ; time value from the input file,
[218]2341      ; 'end_time_step' is rounded to the next existent time step
[190]2342      ;
2343      ; data type: double; dimension: [hour]
2344      ;
[418]2345      ; example:   end_time_step = 2.0d
[218]2346      ;
[190]2347      ; default:   -1.d (last available time step)
2348      ;***************************************************
2349      if (.not. isvar("end_time_step"))then
2350     
2351         
2352         end_time_step = -1.d
2353         
[267]2354
[566]2355      end if             
[190]2356      ;***************************************************
2357      ; OPTIONAL --over-- OVERLAYING
2358      ;
[513]2359      ; predefined overlaying of standard variables are switched
2360      ; on [=1] or off [=0];
[190]2361      ; this feature cannot be used with statistic regions
2362      ;
2363      ; data type: integer
2364      ;
[218]2365      ; example: over = 0
2366      ;
[190]2367      ; default:   0
2368      ;***************************************************
2369      if (.not. isvar("over"))then
2370     
2371         
2372         over = 0
2373         
2374         
[218]2375      end if
2376      ;***************************************************
2377      ; OPTIONAL --norm_t-- NORMALISING X-AXIS OF TIMESERIES PLOTS (T-AXIS)
2378      ;
[513]2379      ; normalising the x-axis with the given value;
[566]2380      ; the unit can be set explicitly;
[418]2381      ; by default no normalising and unit is set to [s]
[218]2382      ;
2383      ; data type: float
2384      ;
[418]2385      ; example: norm_t = 3600.
[218]2386      ;
2387      ; default:   1.
2388      ;***************************************************
[513]2389      if(.not. isvar("norm_t"))then
[218]2390     
2391         
2392         norm_t = 1.
2393         
2394         
2395      end if
2396      ;***************************************************
2397      ; OPTIONAL --unit_t-- UNIT OF X-AXIS OF TIMESERIES PLOTS (T-AXIS)
2398      ;
2399      ; unit of the x-axis can be set; by default unit is set to [s]
2400      ;
2401      ; data type: string
2402      ;
2403      ; example: unit_t = "h"
2404      ;
2405      ; default:   " "
2406      ;***************************************************
[513]2407      if(.not. isvar("unit_t"))then
[218]2408     
2409         
2410         unit_t = " "
2411         
2412         
2413      end if
2414      ;***************************************************
[566]2415      ; OPTIONAL --no_columns-- NUMBER OF COLUMNS
2416      ;
2417      ; number of plots in one row
2418      ;
2419      ; data type: integer
2420      ;
2421      ; example: no_columns = 1
2422      ;
2423      ; default:   1
2424      ;***************************************************
2425      if (.not. isvar("no_columns"))then       
2426     
2427       
2428         no_columns = 1
2429         
2430         
2431      end if     
2432      ;***************************************************
2433      ; OPTIONAL --no_rows-- NUMBER OF ROWS
2434      ;
2435      ; number of plots in one column
2436      ;
2437      ; data type: integer
2438      ;
2439      ; example: no_rows = 2
2440      ;
2441      ; default:   2
2442      ;***************************************************
2443      if (.not. isvar("no_rows"))then
2444     
2445         
2446         no_rows = 2
2447         
2448         
2449      end if     
2450      ;***************************************************
2451      ; OPTIONAL --font_size-- FONT SIZE OF STRINGS
2452      ;
2453      ; font size can be changed to fit to the current plot
2454      ;
2455      ; data type: float
2456      ;
2457      ; example: font_size = 0.07
2458      ;
2459      ; default:   0.07
2460      ;***************************************************
2461      if(.not. isvar("font_size"))then
2462     
2463         
2464         font_size = 0.07
2465         
2466         
2467      end if
2468      ;***************************************************
[190]2469
2470   end if
2471         
2472;*********************************************************       
2473   
2474end
Note: See TracBrowser for help on using the repository browser.