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

Last change on this file since 927 was 782, checked in by heinze, 13 years ago

Bugfix concerning setting of font size in legends.

File size: 69.6 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      ;
825      ; default:   "all" (all available variables are plotted)
826      ;***************************************************
827      if (.not. isvar("var"))then
828     
829         
830         var = "all"
831         
832         
833      end if
834      ;***************************************************
835      ; OPTIONAL --start_time_step-- FIRST TIME STEP
836      ;
837      ; first time step of plot; there is no need to specify a precise
838      ; time value from the input file,
839      ; 'start_time_step' is rounded to the next existent time step
840      ;
841      ; data type: double; dimension: [hour]
842      ;
843      ; example:   start_time_step = 0.5d
844      ;
845      ; default:   -1.d (first available time step)
846      ;***************************************************
847      if (.not. isvar("start_time_step"))then   
848     
849       
850         start_time_step = -1.d
851         
852         
853      end if     
854      ;***************************************************
855      ; OPTIONAL --end_time_step-- LAST TIME STEP
856      ;
857      ; last time step of plot; there is no need to specify a precise
858      ; time value from the input file,
859      ; 'end_time_step' is rounded to the next existent time step
860      ;
861      ; data type: double; dimension: [hour]
862      ;
863      ; example:   end_time_step = 2.0d
864      ;
865      ; default:   -1.d (last available time step)
866      ;***************************************************
867      if (.not. isvar("end_time_step"))then
868     
869         
870         end_time_step = -1.d   
871         
872         
873      end if
874      ;***************************************************
875      ; OPTIONAL --time_stride-- TEMPORAL STRIDE
876      ;
877      ; temporal stride for the plots; if time_stride=1 all the profiles
878      ; between first_time_step and end_time_step are plotted
879      ; (including first_time_step and end_time_step),
880      ; if time_stride=3 every third time step between start_time_step and
881      ; end_time_step is plotted
882      ;
883      ; data type: integer
884      ;
885      ; example:   time_stride = 3
886      ;
887      ; default:   1
888      ;***************************************************
889      if (.not. isvar("time_stride"))then       
890     
891       
892         time_stride = 1
893         
894         
895      end if
896      ;***************************************************
897      ; OPTIONAL --xs-- MINIMUM X-AXIS
898      ;
899      ; value for minimum of x-axis
900      ;
901      ; data type: float; dimension of variable
902      ;
903      ; example: xs = 300. (e.g. if var=",pt")
904      ;
905      ; default:   -1. (minimum of range)
906      ;***************************************************
907      if (.not. isvar("xs"))then   
908     
909         
910         xs = -1.
911         
912         
913      end if     
914      ;***************************************************
915      ; OPTIONAL --xe-- MAXIMUM X-AXIS
916      ;
917      ; value for maximum of x-axis
918      ;
919      ; data type: float; dimension of variable
920      ;
921      ; example: xe = 302. (e.g. if var=",pt")
922      ;
923      ; default:   -1. (maximum of range)
924      ;***************************************************
925      if (.not. isvar("xe"))then 
926     
927         
928         xe = -1.
929         
930         
931      end if     
932      ;***************************************************
933      ; OPTIONAL --min_z-- MINIMUM Z-AXIS
934      ;
935      ; value for minimum height of z-axis
936      ;
937      ; data type: double; dimension: [meter]
938      ;
939      ; example: min_z = 0.0d
940      ;
941      ; default:   -1.d (related to the minimum of z-axis in input file)
942      ;***************************************************
943      if (.not. isvar("min_z"))then   
944     
945         
946         min_z = -1.d
947         
948         
949      end if     
950      ;***************************************************
951      ; OPTIONAL --max_z-- MAXIMUM Z-AXIS
952      ;
953      ; value for maximum height of z-axis
954      ;
955      ; data type: double; dimension: [meter]
956      ;
957      ; example: max_z = 1000.0d
958      ;
959      ; default:   -1.d (related to the maximum of z-axis in input file)
960      ;***************************************************
961      if (.not. isvar("max_z"))then   
962     
963         
964         max_z = -1.d
965         
966         
967      end if     
968      ;***************************************************
969      ; OPTIONAL --start_x-- STARTVALUE IN X-DIRECTION 
970      ;
971      ; if 3D data is used for profiles you can choose a startvalue
972      ; in x-direction for horizontal averaging
973      ;
974      ; data type: integer; dimension [grippoint]
975      ;
976      ; example: start_x = 0
977      ;
978      ; default:   0 (minimum of range)
979      ;***************************************************
980      if (.not. isvar("start_x"))then   
981     
982         
983         start_x = 0
984         
985         
986      end if
987      ;***************************************************
988      ; OPTIONAL --end_x-- ENDVALUE IN X-DIRECTION 
989      ;
990      ; if 3D data is used for profiles you can choose an endvalue
991      ; in x-direction for horizontal averaging
992      ;
993      ; data type: integer; dimension [grippoint]
994      ;
995      ; example: end_x = 200
996      ;
997      ; default:   -1 (maximum of range excluding ghostpoint)
998      ;***************************************************
999      if (.not. isvar("end_x"))then   
1000     
1001         
1002         end_x = -1
1003         
1004         
1005      end if
1006      ;***************************************************
1007      ; OPTIONAL --start_y-- STARTVALUE IN Y-DIRECTION 
1008      ;
1009      ; if 3D data is used for profiles you can choose a startvalue
1010      ; in y-direction for horizontal averaging
1011      ;
1012      ; data type: integer; dimension [grippoint]
1013      ;
1014      ; example: start_y = 0
1015      ;
1016      ; default:   0 (minimum of range)
1017      ;***************************************************
1018      if (.not. isvar("start_y"))then   
1019     
1020         
1021         start_y = 0
1022         
1023         
1024      end if
1025      ;***************************************************
1026      ; OPTIONAL --end_y-- ENDVALUE IN Y-DIRECTION 
1027      ;
1028      ; if 3D data is used for profiles you can choose an endvalue
1029      ; in y-direction for horizontal averaging
1030      ;
1031      ; data type: integer; dimension [grippoint]
1032      ;
1033      ; example: end_y = 200
1034      ;
1035      ; default:   -1 (maximum of range excluding ghostpoint)
1036      ;***************************************************
1037      if (.not. isvar("end_y"))then   
1038     
1039         
1040         end_y = -1
1041         
1042         
1043      end if     
1044      ;***************************************************
[218]1045      ; OPTIONAL --no_columns-- NUMBER OF COLUMNS
[190]1046      ;
1047      ; number of plots in one row
1048      ;
1049      ; data type: integer
1050      ;
1051      ; default:   1
1052      ;***************************************************
1053      if (.not. isvar("no_columns"))then         
1054     
1055     
1056         no_columns = 1
1057         
1058         
1059      end if     
1060      ;***************************************************
[218]1061      ; OPTIONAL --no_rows-- NUMBER OF ROWS
[190]1062      ;
1063      ; number of plots in one column
1064      ;
1065      ; data type: integer
1066      ;
1067      ; default:   2
1068      ;***************************************************
[218]1069      if (.not. isvar("no_rows"))then   
[190]1070     
1071       
[218]1072         no_rows = 2
[190]1073         
1074         
1075      end if     
1076      ;***************************************************
[566]1077      ; OPTIONAL --log_z-- LOGARITHMIC SCALE OF Z-AXIS
[190]1078      ;
[566]1079      ; logarithmic scale for z-axis is switched on [=1] or off [=0]
1080      ;
1081      ; data type: integer
1082      ;
1083      ; example: log_z = 0
1084      ;
1085      ; default:   0
1086      ;***************************************************
1087      if (.not. isvar("log_z"))then
1088     
1089           
1090         log_z = 0
1091         
1092         
1093      end if     
1094      ;***************************************************
1095      ; OPTIONAL --norm_z-- NORMALISING Z-AXIS
1096      ;
1097      ; value for normalising the z-axis
1098      ;
1099      ; data type: float
1100      ;
1101      ; example: norm_z = 1.
1102      ;
1103      ; default:   1. (no normalising)
1104      ;***************************************************
1105      if (.not. isvar("norm_z"))then
1106     
1107           
1108         norm_z = 1.
1109         
1110         
1111      end if     
1112      ;***************************************************
1113      ; OPTIONAL --over-- OVERLAYING
1114      ;
1115      ; predefined overlaying of standard variables is switched
1116      ; on [=1] or off [=0];
1117      ; this feature cannot be used with statistic regions
1118      ;
1119      ; data type: integer
1120      ;
1121      ; example: over = 0
1122      ;
1123      ; default:   0
1124      ;***************************************************
1125      if (.not. isvar("over"))then
1126     
1127         
1128         over = 0
1129         
1130         
1131      end if     
1132      ;***************************************************
1133      ; OPTIONAL --combine-- COMBINING
1134      ;
1135      ; combining of two or three individual variables is possible and
1136      ; switched on [=1] or off[=0]
1137      ;
1138      ; data type: integer
1139      ;
1140      ; example: combine = 0
1141      ;
1142      ; default:   0
1143      ;***************************************************
1144      if (.not. isvar("combine"))then 
1145     
1146         
1147         combine = 0
1148         
1149         
1150      end if     
1151      ;***************************************************
1152      ; REQUIRED IF combine=1 --number_comb-- NUMBER FOR COMBINING
1153      ;
1154      ; number of variables for combining; 2 or 3 is valid
1155      ;
1156      ; data type: integer
1157      ;
1158      ; example: number_comb = 2
1159      ;
1160      ; default:   -1
1161      ;***************************************************
1162      if (.not. isvar("number_comb"))then
1163     
1164           
1165         number_comb = -1
1166         
1167         
1168      end if     
1169      ;***************************************************
1170      ; REQUIRED IF combine=1 --c_var-- VARIABLES FOR COMBINING
1171      ;
1172      ; names of variables;
[418]1173      ; inidcate them with one comma before and after each variable
[526]1174      ; if the name contains double quotes (") use dq and the '+' operation to
1175      ; concatenate the strings
[190]1176      ;
[566]1177      ; data type: string
[190]1178      ;
[566]1179      ; example: c_var = ",u,v,w,w*pt*,w"+dq+"pt"+dq+",wpt,"
1180      ;
1181      ; default:   "c_variables"
[190]1182      ;***************************************************
[566]1183      if (.not. isvar("c_var"))then
[190]1184     
[566]1185           
1186         c_var = "c_variables"
[190]1187         
1188         
[566]1189      end if     
1190      ;***************************************************
1191      ; OPTIONAL --black-- BLACK OR COLOR LINES
1192      ;
1193      ; colored [=0] or black and white [=1] plots
1194      ;
1195      ; data type: integer
1196      ;
1197      ; example: black = 0
1198      ;
1199      ; default:   0
1200      ;***************************************************
1201      if (.not. isvar("black"))then
1202     
1203           
1204         black = 0
1205         
1206         
1207      end if     
1208      ;***************************************************
1209      ; OPTIONAL --dash-- DASHED LINES
1210      ;
1211      ; use of different line patterns [=1] or continous lines [=0]
1212      ; for different time steps
1213      ;
1214      ; data type: integer
1215      ;
1216      ; example: dash = 0
1217      ;
1218      ; default:   0
1219      ;***************************************************
1220      if (.not. isvar("dash"))then     
1221     
1222       
1223         dash = 0     
1224         
1225         
[219]1226      end if
1227      ;***************************************************
[566]1228      ; OPTIONAL --font_size-- FONT SIZE OF STRINGS
1229      ;
1230      ; font size can be changed to fit to the current plot
1231      ;
1232      ; data type: float
1233      ;
1234      ; example: font_size = 0.02
1235      ;
1236      ; default:   0.02
1237      ;***************************************************
1238      if(.not. isvar("font_size"))then
1239     
1240         
1241         font_size = 0.02
1242         
1243         
1244      end if
1245      ;***************************************************
[219]1246      ; OPTIONAL --legend-- SWITCHES ON A LEGEND
1247      ;
1248      ; a legend is switched on [=1] or off [=0]
1249      ;
1250      ; data type: integer
1251      ;
1252      ; example: legend = 1
1253      ;
1254      ; default:   1
1255      ;***************************************************
1256      if (.not. isvar("legend"))then
1257     
1258         
[236]1259         legend = 1
[219]1260         
1261         
[190]1262      end if   
1263      ;***************************************************
[566]1264      ; OPTIONAL --font_size_legend-- FONT SIZE OF LEGEND STRINGS
1265      ;
1266      ; font size of the legend can be changed to fit to the current plot
1267      ;
1268      ; data type: float
1269      ;
[782]1270      ; example: font_size_legend = 0.02
[566]1271      ;
[782]1272      ; default:   0.02
[566]1273      ;***************************************************
1274      if(.not. isvar("font_size_legend"))then
1275     
1276         
[782]1277         font_size_legend = 0.02
[566]1278         
1279         
1280      end if     
1281      ;***************************************************
[190]1282      ; OPTIONAL ---no_files-- NUMBER OF INPUT FILES
1283      ;
[418]1284      ; up to 6 input files of different runs can be used; all input files
1285      ; must contain the same variables and dimensionalties;
[218]1286      ; identical variables are merged in one plot;
[190]1287      ; combining or overlaying of any variables is not feasible;
[418]1288      ; an extra legend (besides the time step legend) is drawn where you
1289      ; can explicitily name the items
[190]1290      ;
1291      ; data type: integer
1292      ;
1293      ; default:   1
1294      ;***************************************************
1295      if (.not. isvar("no_files"))then
1296     
1297         
1298         no_files = 1
1299         
1300         
1301      end if     
1302      ;***************************************************
1303      ; OPTIONAL --name_legend_1-- 1st LEGEND ITEM
1304      ;
1305      ; if no_files > 1: name of 1st item in extra legend
1306      ;
1307      ; data type: string
1308      ;
1309      ; default:   "File 1"
1310      ;***************************************************
1311      if (.not. isvar("name_legend_1"))then
1312     
1313           
1314         name_legend_1 = "File 1"
1315         
1316         
[218]1317      end if
[190]1318      ;***************************************************
[218]1319      ; REQUIRED if no_files > 1 --file_2-- 2ND INPUT FILE
[190]1320      ;
[418]1321      ; 2nd input file produced by PALM; if there is more than one file
1322      ; for all timesteps of a job chain,
1323      ; declare path and file name without cycle number and ".nc" (first and
[513]1324      ; last cycle numbers are given with the parameters
1325      ; 'start_f' and 'end_f');
[218]1326      ; otherwise declare path and full file name
[190]1327      ;
1328      ; data type: string
1329      ;
[418]1330      ; example: file_2 = "/path_to_file/file_pr"     
1331      ;                    if you use more than one cycle number
1332      ;          file_2 = "/path_to_file/file_pr.4.nc"
1333      ;                    if you use one file (e.g. with cycle number 4)
[218]1334      ;
[190]1335      ; default:   "File in"
1336      ;***************************************************
[218]1337      if(.not. isvar("file_2"))then
1338
1339
1340         file_2 = "File in"
1341         
1342         
1343      end if
1344      ;***************************************************
[513]1345      ; OPTIONAL/REQUIRED --start_f_2-- if more than one file of a job chain
1346      ;                       START OF CYCLIC NUMBER OF THE 2ND INPUT FILE
[218]1347      ;
[418]1348      ; if there is more than one file for all timesteps of a job chain,
[513]1349      ; declare the first cyclicnumber which contains timesteps
1350      ; you'd like to use
[218]1351      ;
1352      ; data type: integer
1353      ;
1354      ; example:   3 (if first file: example_ts.3.nc)
1355      ;
1356      ; default:   -1 
1357      ;***************************************************
1358      if (.not. isvar("start_f_2"))then 
[190]1359     
1360         
[218]1361         start_f_2 = -1
[190]1362         
1363         
[218]1364      end if
[190]1365      ;***************************************************
[513]1366      ; REQUIRED --end_f_2-- if more than one file of a job chain
1367      ;                      END OF CYCLIC NUMBER OF THE 2ND INPUT FILE
[218]1368      ;
[418]1369      ; if there is more than one file for all timesteps of a job chain,
[513]1370      ; declare the last cyclicnumber which contains timesteps
1371      ; you'd like to use
[218]1372      ;
1373      ; data type: integer
1374      ;
1375      ; example:   48 (if last file: example_ts.48.nc)
1376      ;
1377      ; default:   -1 
1378      ;***************************************************
1379      if (.not. isvar("end_f_2"))then 
1380     
1381         
1382         end_f_2 = -1
1383         
1384         
1385      end if             
1386      ;***************************************************
[190]1387      ; OPTIONAL --name_legend_2-- 2nd LEGEND ITEM
1388      ;
1389      ; if no_files > 1: name of 2nd item in extra legend
1390      ;
1391      ; data type: string
1392      ;
1393      ; default:   "File 2"
1394      ;***************************************************
1395      if (.not. isvar("name_legend_2"))then   
1396     
1397         
1398         name_legend_2 = "File 2"
1399         
1400         
[218]1401      end if
[190]1402      ;***************************************************
[218]1403      ; REQUIRED if no_files > 2 --file_3-- 3RD INPUT FILE
[190]1404      ;
[418]1405      ; 3rd input file produced by PALM; if there is more than one file
1406      ; for all timesteps of a job chain,
1407      ; declare path and file name without cycle number and ".nc" (first and
[513]1408      ; last cycle numbers are given with the parameters
1409      ; 'start_f' and 'end_f');
[218]1410      ; otherwise declare path and full file name
[190]1411      ;
1412      ; data type: string
1413      ;
[418]1414      ; example: file_3 = "/path_to_file/file_pr"     
1415      ;                     if you use more than one cycle number
1416      ;          file_3 = "/path_to_file/file_pr.4.nc"
1417      ;                     if you use one file (e.g. with cycle number 4)
[218]1418      ;
[190]1419      ; default:   "File in"
1420      ;***************************************************
[218]1421      if(.not. isvar("file_3"))then
1422
1423
1424         file_3 = "File in"
[190]1425         
1426         
[218]1427      end if
[190]1428      ;***************************************************
[513]1429      ; OPTIONAL/REQUIRED --start_f_3-- if more than one file of a job chain
1430      ;                        START OF CYCLIC NUMBER OF THE 3RD INPUT FILE
[190]1431      ;
[418]1432      ; if there is more than one file for all timesteps of a job chain,
[513]1433      ; declare the first cyclicnumber which contains timesteps
1434      ; you'd like to use
[190]1435      ;
[218]1436      ; data type: integer
[190]1437      ;
[218]1438      ; example:   3 (if first file: example_ts.3.nc)
1439      ;
1440      ; default:   -1 
[190]1441      ;***************************************************
[218]1442      if (.not. isvar("start_f_3"))then 
[190]1443     
[218]1444         
1445         start_f_3 = -1
[190]1446         
1447         
[218]1448      end if
[190]1449      ;***************************************************
[513]1450      ; OPTIONAL/REQUIRED --end_f_3-- if more than one file of a job chain
1451      ;                    END OF CYCLIC NUMBER OF THE 3RD INPUT FILE
[190]1452      ;
[513]1453      ; if there is more than one file for all timesteps of a job chain,
1454      ; declare
[418]1455      ; the last cyclicnumber which contains timesteps you'd like to use
[190]1456      ;
[218]1457      ; data type: integer
[190]1458      ;
[218]1459      ; example:   48 (if last file: example_ts.48.nc)
1460      ;
1461      ; default:   -1 
[190]1462      ;***************************************************
[218]1463      if (.not. isvar("end_f_3"))then 
[190]1464     
[218]1465         
1466         end_f_3 = -1
[190]1467         
1468         
[218]1469      end if   
[190]1470      ;***************************************************
[218]1471      ; OPTIONAL --name_legend_3-- 3rd LEGEND ITEM
[190]1472      ;
[218]1473      ; if no_files > 2: name of 3rd item in extra legend
[190]1474      ;
1475      ; data type: string
1476      ;
[218]1477      ; default:   "File 3"
[190]1478      ;***************************************************
[218]1479      if (.not. isvar("name_legend_3"))then
[190]1480     
[218]1481           
1482         name_legend_3 = "File 3"
[190]1483         
1484         
[218]1485      end if
[190]1486      ;***************************************************
[218]1487      ; REQUIRED if no_files > 3 --file_4-- 4TH INPUT FILE
[190]1488      ;
[418]1489      ; 4th input file produced by PALM; if there is more than one file for
1490      ; all timesteps of a job chain,
1491      ; declare path and file name without cycle number and ".nc" (first and
[513]1492      ; last cycle numbers are given with the parameters
1493      ; 'start_f' and 'end_f');
[218]1494      ; otherwise declare path and full file name
[190]1495      ;
1496      ; data type: string
1497      ;
[418]1498      ; example: file_4 = "/path_to_file/file_pr"     
1499      ;                     if you use more than one cycle number
1500      ;          file_4 = "/path_to_file/file_pr.4.nc"
1501      ;                      if you use one file (e.g. with cycle number 4)
[218]1502      ;
[190]1503      ; default:   "File in"
1504      ;***************************************************
[218]1505      if(.not. isvar("file_4"))then
1506
1507
1508         file_4 = "File in"
[190]1509         
1510         
[218]1511      end if
[190]1512      ;***************************************************
[513]1513      ; OPTIONAL/REQUIRED --start_f_4-- if more than one file of a job chain
1514      ;                          START OF CYCLIC NUMBER OF THE 4TH INPUT FILE
[190]1515      ;
[418]1516      ; if there is more than one file for all timesteps of a job chain,
[513]1517      ; declare the first cyclicnumber which contains timesteps
1518      ; you'd like to use
[190]1519      ;
[218]1520      ; data type: integer
[190]1521      ;
[218]1522      ; example:   3 (if first file: example_ts.3.nc)
1523      ;
1524      ; default:   -1 
[190]1525      ;***************************************************
[218]1526      if (.not. isvar("start_f_4"))then 
[190]1527     
[218]1528         
1529         start_f_4 = -1
[190]1530         
1531         
[218]1532      end if
[190]1533      ;***************************************************
[513]1534      ; OPTIONAL/REQUIRED --end_f_4-- if more than one file of a job chain
1535      ;                   END OF CYCLIC NUMBER OF THE 4TH INPUT FILE
[190]1536      ;
[418]1537      ; if there is more than one file for all timesteps of a job chain,
[513]1538      ; declare the last cyclicnumber which contains timesteps
1539      ; you'd like to use
[190]1540      ;
[218]1541      ; data type: integer
[190]1542      ;
[218]1543      ; example:   48 (if last file: example_ts.48.nc)
1544      ;
1545      ; default:   -1 
[190]1546      ;***************************************************
[218]1547      if (.not. isvar("end_f_4"))then 
[190]1548     
[218]1549         
1550         end_f_4 = -1
[190]1551         
1552         
[218]1553      end if 
[190]1554      ;***************************************************
[218]1555      ; OPTIONAL --name_legend_4-- 4th LEGEND ITEM
[190]1556      ;
[218]1557      ; if no_files > 3: name of 4th item in extra legend
[190]1558      ;
1559      ; data type: string
1560      ;
[218]1561      ; default:   "File 4"
[190]1562      ;***************************************************
[218]1563      if (.not. isvar("name_legend_4"))then   
[190]1564     
[218]1565       
1566         name_legend_4 = "File 4"
[190]1567         
1568         
[218]1569      end if
[190]1570      ;***************************************************
[218]1571      ; REQUIRED if no_files > 4 --file_5-- 5TH INPUT FILE
[190]1572      ;
[418]1573      ; 5th input file produced by PALM; if there is more than one file for
1574      ; all timesteps of a job chain,
1575      ; declare path and file name without cycle number and ".nc" (first and
[513]1576      ; last cycle numbers are given with the parameters
1577      ; 'start_f' and 'end_f');
[218]1578      ; otherwise declare path and full file name
[190]1579      ;
1580      ; data type: string
1581      ;
[418]1582      ; example: file_5 = "/path_to_file/file_pr"     
1583      ;                     if you use more than one cycle number
1584      ;          file_5 = "/path_to_file/file_pr.4.nc"
1585      ;                     if you use one file (e.g. with cycle number 4)
[218]1586      ;
[190]1587      ; default:   "File in"
1588      ;***************************************************
[218]1589      if(.not. isvar("file_5"))then
1590
1591
1592         file_5 = "File in"
[190]1593         
1594         
[218]1595      end if
[190]1596      ;***************************************************
[513]1597      ; OPTIONAL/REQUIRED --start_f_5-- if more than one file of a job chain
1598      ;                     START OF CYCLIC NUMBER OF THE 5TH INPUT FILE
[190]1599      ;
[418]1600      ; if there is more than one file for all timesteps of a job chain,
[513]1601      ; declare the first cyclicnumber which contains timesteps
1602      ; you'd like to use
[190]1603      ;
[218]1604      ; data type: integer
[190]1605      ;
[218]1606      ; example:   3 (if first file: example_ts.3.nc)
1607      ;
1608      ; default:   -1 
[190]1609      ;***************************************************
[218]1610      if (.not. isvar("start_f_5"))then 
[190]1611     
[218]1612         
1613         start_f_5 = -1
[190]1614         
1615         
[218]1616      end if
[190]1617      ;***************************************************
[513]1618      ; OPTIONAL/REQUIRED --end_f_5-- if more than one file of a job chain
1619      ;END OF CYCLIC NUMBER OF THE 5TH INPUT FILE
[190]1620      ;
[418]1621      ; if there is more than one file for all timesteps of a job chain,
[513]1622      ; declare the last cyclicnumber which contains timesteps
1623      ; you'd like to use
[190]1624      ;
[218]1625      ; data type: integer
[190]1626      ;
[218]1627      ; example:   48 (if last file: example_ts.48.nc)
1628      ;
1629      ; default:   -1 
[190]1630      ;***************************************************
[218]1631      if (.not. isvar("end_f_5"))then 
[190]1632     
[218]1633         
1634         end_f_5 = -1
[190]1635         
1636         
[218]1637      end if     
[190]1638      ;***************************************************
[218]1639      ; OPTIONAL --name_legend_5-- 5th LEGEND ITEM
[190]1640      ;
[218]1641      ; if no_files > 4: name of 5th item in extra legend
[190]1642      ;
1643      ; data type: string
1644      ;
[218]1645      ; default:   "File 5"
[190]1646      ;***************************************************
[218]1647      if (.not. isvar("name_legend_5"))then
[190]1648     
[218]1649           
1650         name_legend_5 = "File 5"
[190]1651         
1652         
[218]1653      end if
[190]1654      ;***************************************************
[218]1655      ; REQUIRED if no_files > 5 --file_6-- 6TH INPUT FILE
[190]1656      ;
[418]1657      ; 6th input file produced by PALM; if there is more than one file for
1658      ; all timesteps of a job chain,
1659      ; declare path and file name without cycle number and ".nc" (first and
[513]1660      ; last cycle numbers are given with the parameters
1661      ; 'start_f' and 'end_f');
[218]1662      ; otherwise declare path and full file name
[190]1663      ;
1664      ; data type: string
1665      ;
[418]1666      ; example: file_6 = "/path_to_file/file_pr"     
1667      ;                     if you use more than one cycle number
1668      ;          file_6 = "/path_to_file/file_pr.4.nc"
1669      ;                     if you use one file (e.g. with cycle number 4)
[218]1670      ;
[190]1671      ; default:   "File in"
1672      ;***************************************************
[218]1673      if(.not. isvar("file_6"))then
1674
1675
1676         file_6 = "File in"
[190]1677         
1678         
[218]1679      end if
[190]1680      ;***************************************************
[513]1681      ; OPTIONAL/REQUIRED --start_f_6-- if more than one file of a job chain
1682      ;                     START OF CYCLIC NUMBER OF THE 6TH INPUT FILE
[190]1683      ;
[418]1684      ; if there is more than one file for all timesteps of a job chain,
[513]1685      ; declare the first cyclicnumber which contains timesteps
1686      ; you'd like to use
[190]1687      ;
[218]1688      ; data type: integer
[190]1689      ;
[218]1690      ; example:   3 (if first file: example_ts.3.nc)
1691      ;
1692      ; default:   -1 
[190]1693      ;***************************************************
[218]1694      if (.not. isvar("start_f_6"))then 
[190]1695     
[218]1696         
1697         start_f_6 = -1
[190]1698         
1699         
[218]1700      end if
[190]1701      ;***************************************************
[513]1702      ; OPTIONAL/REQUIRED --end_f_6-- if more than one file of a job chain
1703      ; END OF CYCLIC NUMBER OF THE 6TH INPUT FILE
[190]1704      ;
[418]1705      ; if there is more than one file for all timesteps of a job chain,
[513]1706      ; declare the last cyclicnumber which contains timesteps
1707      ; you'd like to use
[190]1708      ;
[218]1709      ; data type: integer
[190]1710      ;
[218]1711      ; example:   48 (if last file: example_ts.48.nc)
1712      ;
1713      ; default:   -1 
[190]1714      ;***************************************************
[218]1715      if (.not. isvar("end_f_6"))then 
[190]1716     
[218]1717         
1718         end_f_6 = -1
[190]1719         
1720         
[218]1721      end if                     
[190]1722      ;***************************************************
[218]1723      ; OPTIONAL --name_legend_6-- 6TH LEGEND ITEM
[190]1724      ;
[218]1725      ; if no_files > 5: name of 6th item in extra legend
[190]1726      ;
1727      ; data type: string
1728      ;
[218]1729      ; default:   "File 6"
[190]1730      ;***************************************************
[218]1731      if (.not. isvar("name_legend_6"))then
[190]1732     
[218]1733           
1734         name_legend_6 = "File 6"
[190]1735         
1736         
[218]1737      end if             
[190]1738      ;***************************************************
1739   end if
[194]1740
[190]1741   ;******************************************************
1742   ; parameters used by SPECTRA
1743   ;******************************************************
1744   
[194]1745   if (script .EQ. "spectra") then
1746
[190]1747      ;***************************************************
1748      ; REQUIRED --file_1-- INPUT FILE
1749      ;
[418]1750      ; input file produced by PALM; if there is more than one file for
1751      ; all timesteps of a job chain,
1752      ; declare path and file name without cycle number and ".nc" (first and
[513]1753      ; last cycle numbers are given with the parameters
1754      ; 'start_f' and 'end_f');
[218]1755      ; otherwise declare path and full file name
[190]1756      ;
1757      ; data type: string
1758      ;
[418]1759      ; example: file_1 = "/path_to_file/file_sp"     
1760      ;                     if you use more than one cycle number
1761      ;          file_1 = "/path_to_file/file_sp.4.nc"
1762      ;                     if you use one file (e.g. with cycle number 4)
[218]1763      ;
[190]1764      ; default:   "File in"
1765      ;***************************************************
[218]1766      if(.not. isvar("file_1"))then
1767
1768
1769         file_1 = "File in"
1770         
1771         
1772      end if
1773      ;***************************************************
[513]1774      ; OPTIONAL/REQUIRED if more than one file of a job chain
[418]1775      ;           --start_f-- START OF CYCLIC NUMBER
[218]1776      ;
[418]1777      ; if there is more than one file for all timesteps of a job chain,
[513]1778      ; declare the first cyclicnumber which contains timesteps
1779      ; you'd like to use
[218]1780      ;
1781      ; data type: integer
1782      ;
[418]1783      ; example:   start_f = 3 (if first file: example_sp.3.nc) and
1784      ;            start_f = 0 (if first file: example_sp.nc)
[218]1785      ;
1786      ; default:   -1 
1787      ;***************************************************
1788      if (.not. isvar("start_f"))then 
[190]1789     
1790         
[218]1791         start_f = -1
[190]1792         
1793         
[218]1794      end if
[190]1795      ;***************************************************
[513]1796      ; OPTIONAL/REQUIRED if more than one file of a job chain
[418]1797      ;          --end_f-- END OF CYCLIC NUMBER
[218]1798      ;
[418]1799      ; if there is more than one file for all timesteps of a job chain,
[513]1800      ; declare the last cyclicnumber which contains timesteps
1801      ; you'd like to use
[218]1802      ;
1803      ; data type: integer
1804      ;
[418]1805      ; example:   end_f = 48 (if last file: example_sp.48.nc) and
1806      ;            end_f =  0 (if last file: example_sp.nc)
[218]1807      ;
1808      ; default:   -1 
1809      ;***************************************************
1810      if (.not. isvar("end_f"))then 
1811     
1812         
1813         end_f = -1
1814         
1815         
1816      end if
1817      ;***************************************************
[190]1818      ; OPTIONAL --format_out-- OUTPUT FORMAT
1819      ;
[418]1820      ; format of the output file;
[532]1821      ; supported file formats: x11, pdf, ps, eps, epsi, ncgm or png
[190]1822      ;
1823      ; data type: string
1824      ;
[218]1825      ; example: format_out = "pdf"
1826      ;
[190]1827      ; default:   "x11"
1828      ;***************************************************
1829      if (.not. isvar("format_out"))then       
1830     
1831       
1832         format_out = "x11"
1833         
1834         
1835      end if     
1836      ;***************************************************
1837      ; OPTIONAL --file_out-- OUTPUT FILE
1838      ;
[218]1839      ; path and name of the output file
[190]1840      ;
1841      ; data type: string
1842      ;
[218]1843      ; example: file_out = "/home/path_to_output/test_sp"
1844      ;
[194]1845      ; default:   "~/test_sp"
[190]1846      ;***************************************************
1847      if (.not. isvar("file_out"))then 
1848     
1849       
1850         file_out = "~/test_sp"
1851         
1852         
[566]1853      end if   
[190]1854      ;***************************************************
[566]1855      ; OPTIONAL --var-- OUTPUT VARIABLES
[190]1856      ;
[566]1857      ; names of output variables;
1858      ; inidcate them with one comma before and after each variable
1859      ; if the name contains double quotes (") use dq and the '+' operation to
1860      ; concatenate the strings
[190]1861      ;
[566]1862      ; data type: string
[218]1863      ;
[566]1864      ; example: var = ",u,v,w,"
[190]1865      ;
[566]1866      ; default:   "all" (all available variables are plotted)
[190]1867      ;***************************************************
[566]1868      if (.not. isvar("var"))then
[190]1869     
[566]1870         
1871         var = "all"
[190]1872         
1873         
1874      end if     
1875      ;***************************************************
[566]1876      ; OPTIONAL --start_time_step-- FIRST TIME STEP
[190]1877      ;
[566]1878      ; first time step of plot; there is no need to specify a precise
1879      ; time value from the input file,
1880      ; 'start_time_step' is rounded to the next existent time step
[190]1881      ;
[566]1882      ; data type: double; dimension: [hour]
[190]1883      ;
[566]1884      ; example:   start_time_step = 0.5d
1885      ;
1886      ; default:   -1.d (first available time step)
[190]1887      ;***************************************************
[566]1888      if (.not. isvar("start_time_step"))then
[190]1889     
1890         
[566]1891         start_time_step = -1.d
[190]1892         
1893         
1894      end if     
1895      ;***************************************************
[566]1896      ; OPTIONAL --end_time_step-- LAST TIME STEP
[190]1897      ;
[566]1898      ; last time step of plot; there is no need to specify a precise
1899      ; time value from the input file,
1900      ; 'end_time_step' is rounded to the next existent time step
[190]1901      ;
[566]1902      ; data type: double; dimension: [hour]
[190]1903      ;
[566]1904      ; example:   end_time_step = 2.0d
1905      ;
1906      ; default:   -1.d (last available time step)
[190]1907      ;***************************************************
[566]1908      if (.not. isvar("end_time_step"))then
[190]1909     
1910         
[566]1911         end_time_step = -1.d
[190]1912         
1913         
[566]1914      end if   
[190]1915      ;***************************************************
1916      ; OPTIONAL --height_level-- HEIGHT LEVELS
1917      ;
[418]1918      ; defines if all levels or selected levels are output;
1919      ; a 1D array with up to 100 elements
[218]1920      ; from 0 to 99 must be indicated
[190]1921      ;
1922      ; data type: integer array
1923      ;
[218]1924      ; example: height_level=(/0,3,7/)
1925      ;
1926      ; default:   (/-1/) (all levels are output)
[190]1927      ;***************************************************
1928      if (.not. isvar("height_level"))then
1929     
1930         
1931         height_level = (/-1/)
1932         
1933         
1934      end if     
1935      ;***************************************************
1936      ; OPTIONAL --sort-- TYPE OF SORTING
1937      ;
[418]1938      ; defines if all heights [="height"] or all timesteps [="time"]
1939      ; are merged in one plot
[190]1940      ;
1941      ; data type: string
1942      ;
[218]1943      ; example: sort = "height"
1944      ;
[190]1945      ; default:   "height"
1946      ;***************************************************
1947      if (.not. isvar("sort"))then
1948     
1949         
1950         sort = "height"
1951         
1952         
1953      end if     
1954      ;***************************************************
[566]1955      ; OPTIONAL --no_columns-- NUMBER OF COLUMNS
[190]1956      ;
[566]1957      ; number of plots in one row
[190]1958      ;
[566]1959      ; data type: integer
1960      ;
1961      ; example: no_columns = 1
[190]1962      ;
[566]1963      ; default:   1
[190]1964      ;***************************************************
[566]1965      if (.not. isvar("no_columns"))then         
[190]1966     
[566]1967     
1968         no_columns = 1
[190]1969         
1970         
1971      end if     
1972      ;***************************************************
[566]1973      ; OPTIONAL --no_rows-- NUMBER OF ROWS
[190]1974      ;
[566]1975      ; number of plots in one column
[190]1976      ;
[566]1977      ; data type: integer
1978      ;
1979      ; example: no_rows = 2
[190]1980      ;
[566]1981      ; default:   2
[190]1982      ;***************************************************
[566]1983      if (.not. isvar("no_rows"))then
[190]1984     
1985         
[566]1986         no_rows = 2
[190]1987         
1988         
1989      end if     
1990      ;***************************************************
[218]1991      ; OPTIONAL --log_x-- LOGARITHMIC SCALE X-AXIS
[190]1992      ;
[218]1993      ; logarithmic scale for x-axis is switched on [=1] or off [=0]
[190]1994      ;
1995      ; data type: integer
1996      ;
[218]1997      ; example: log_x = 1
1998      ;
[190]1999      ; default:   1
2000      ;***************************************************
[218]2001      if (.not. isvar("log_x"))then   
[190]2002     
2003         
[218]2004         log_x = 1
[190]2005         
2006         
2007      end if     
2008      ;***************************************************
[218]2009      ; OPTIONAL --log_y-- LOGARITHMIC SCALE Y-AXIS
[190]2010      ;
[218]2011      ; logarithmic scale for y-axis is switched on [=1] or off [=0]
[190]2012      ;
2013      ; data type: integer
2014      ;
[218]2015      ; example: log_y = 1
2016      ;
[190]2017      ; default:   1
2018      ;***************************************************
[218]2019      if (.not. isvar("log_y"))then   
[190]2020     
2021         
[218]2022         log_y = 1
[190]2023         
2024         
2025      end if     
2026      ;***************************************************
[566]2027      ; OPTIONAL --norm_height-- NORMALISING X-AXIS WITH HEIGHT
[190]2028      ;
[566]2029      ; normalising the x-axis with the height of each spectra
2030      ; is switched on [=1] or off [=0]
[190]2031      ;
[566]2032      ; data type: integer
[190]2033      ;
[566]2034      ; example: norm_height = 0
[218]2035      ;
[566]2036      ; default:   0 (no normalising)
[190]2037      ;***************************************************
[566]2038      if (.not. isvar("norm_height"))then   
[190]2039     
2040         
[566]2041         norm_height = 0
[190]2042         
2043         
[566]2044      end if     
[218]2045      ;***************************************************
[566]2046      ; OPTIONAL --norm_x-- NORMALISING X-AXIS
[218]2047      ;
[566]2048      ; value for normalising the x-axis
[218]2049      ;
[566]2050      ; data type: float
[218]2051      ;
[566]2052      ; example: norm_x = 1.
[218]2053      ;
[566]2054      ; default:   1. (no normalising)
[218]2055      ;***************************************************
[566]2056      if (.not. isvar("norm_x"))then   
[218]2057     
2058         
[566]2059         norm_x = 1.
[218]2060         
2061         
[566]2062      end if
[190]2063      ;***************************************************
[218]2064      ; OPTIONAL --norm_y-- NORMALISE Y-AXIS
[190]2065      ;
2066      ; value for normalising the y-axis
2067      ;
2068      ; data type: float
2069      ;
[218]2070      ; example: norm_y = 1.
2071      ;
[190]2072      ; default:   1. (no normalising)
2073      ;***************************************************
[218]2074      if (.not. isvar("norm_y"))then     
[190]2075     
2076       
[218]2077         norm_y = 1.
[190]2078         
2079         
[218]2080      end if
[190]2081      ;***************************************************
[218]2082      ; OPTIONAL --unit_x--unit_y-- UNITS OF X-, Y-AXES
2083      ;
[513]2084      ; unit of the x-axis can be set if norm_x is unequal 1,
2085      ; otherwise the unit is [1/m]; by default no unit is set for the y-axis
[218]2086      ;
2087      ; data type: string,
2088      ;
2089      ; example: unit_x = "100/m"
2090      ;
2091      ; default:   " "
2092      ;***************************************************
2093      if(.not. isvar("unit_x"))then     
2094         
2095         unit_x = " "   
2096         
2097      end if
2098      if(.not. isvar("unit_y"))then     
2099         
2100         unit_y = " "   
2101         
2102      end if
2103      ;***************************************************
[566]2104      ; OPTIONAL --black-- BLACK OR COLOR LINES
2105      ;
2106      ; colored [=0] or black and white [=1] plots
2107      ;
2108      ; data type: integer
2109      ;
2110      ; example: black = 0
2111      ;
2112      ; default:   0
2113      ;***************************************************
2114      if (.not. isvar("black"))then
2115     
2116           
2117         black = 0
2118         
2119         
2120      end if     
2121      ;***************************************************
2122      ; OPTIONAL --dash-- DASHED LINES
2123      ;
2124      ; use of different line patterns [=1] or continous lines [=0]
2125      ;
2126      ; data type: integer
2127      ;
2128      ; example: dash = 0
2129      ;
2130      ; default:   0
2131      ;***************************************************
2132      if (.not. isvar("dash"))then   
2133     
2134         
2135         dash = 0
2136         
2137         
2138      end if
2139      ;***************************************************
[218]2140      ; OPTIONAL --font_size-- FONT SIZE OF STRINGS
2141      ;
2142      ; font size can be changed to fit to the current plot
2143      ;
2144      ; data type: float
2145      ;
2146      ; example: font_size = 0.02
2147      ;
2148      ; default:   0.02
2149      ;***************************************************
2150      if(.not. isvar("font_size"))then
2151     
2152         
2153         font_size = 0.02
2154         
2155         
2156      end if
2157      ;***************************************************
2158      ; OPTIONAL --font_size_legend-- FONT SIZE OF LEGEND STRINGS
2159      ;
2160      ; font size of the legend can be changed to fit to the current plot
2161      ;
2162      ; data type: float
2163      ;
[534]2164      ; example: font_size_legend = 0.02
[218]2165      ;
[534]2166      ; default:   0.02
[218]2167      ;***************************************************
2168      if(.not. isvar("font_size_legend"))then
2169     
2170         
[534]2171         font_size_legend = 0.02
[218]2172         
2173         
2174      end if     
2175      ;***************************************************
[190]2176     
2177   end if
[194]2178
[190]2179   ;******************************************************
2180   ; parameters used by TIMESERIES
2181   ;******************************************************
2182
[194]2183   if (script .EQ. "timeseries") then
2184
[190]2185      ;***************************************************
2186      ; REQUIRED --file_1-- INPUT FILE
2187      ;
[418]2188      ; input file produced by PALM; if there is more than one file for
2189      ; all timesteps of a job chain,
2190      ; declare path and file name without cycle number and ".nc" (first and
[513]2191      ; last cycle numbers are given with the parameters
2192      ; 'start_f' and 'end_f');
[218]2193      ; otherwise declare path and full file name
[190]2194      ;
2195      ; data type: string
2196      ;
[418]2197      ; example: file_1 = "/path_to_file/file_ts"     
2198      ;                      if you use more than one cycle number
2199      ;          file_1 = "/path_to_file/file_ts.4.nc"
2200      ;                      if you use one file (e.g. with cycle number 4)
[218]2201      ;
[190]2202      ; default:   "File in"
2203      ;***************************************************
[218]2204      if(.not. isvar("file_1"))then
2205
2206
[267]2207         file_1 = "File in"
[218]2208         
2209         
2210      end if
2211      ;***************************************************
[513]2212      ; OPTIONAL/REQUIRED if more than one file of a job chain
[418]2213      ;          --start_f-- START OF CYCLIC NUMBER
[218]2214      ;
[418]2215      ; if there is more than one file for all timesteps of a job chain,
[513]2216      ; declare the first cyclicnumber which contains timesteps
2217      ; you'd like to use
[218]2218      ;
2219      ; data type: integer
2220      ;
[418]2221      ; example:   start_f = 3 (if first file: example_ts.3.nc) and
2222      ;            start_f = 0 (if first file: example_ts.nc)
[218]2223      ;
2224      ; default:   -1 
2225      ;***************************************************
2226      if (.not. isvar("start_f"))then 
[190]2227     
2228         
[218]2229         start_f = -1
[190]2230         
2231         
[218]2232      end if
[190]2233      ;***************************************************
[513]2234      ; OPTIONAL/REQUIRED if more than one file of a job chain
[418]2235      ;           --end_f-- END OF CYCLIC NUMBER
[218]2236      ;
[418]2237      ; if there is more than one file for all timesteps of a job chain,
[513]2238      ; declare the last cyclicnumber which contains timesteps
2239      ; you'd like to use
[218]2240      ;
2241      ; data type: integer
2242      ;
[418]2243      ; example:   end_f = 48 (if last file: example_ts.48.nc) and
2244      ;            end_f = 0 (if last file: example_ts.nc)
[218]2245      ;
2246      ; default:   -1 
2247      ;***************************************************
2248      if (.not. isvar("end_f"))then 
2249     
2250         
2251         end_f = -1
2252         
2253         
2254      end if
2255      ;***************************************************
[190]2256      ; OPTIONAL --format_out-- OUTPUT FORMAT
2257      ;
[418]2258      ; format of the output file;
[532]2259      ; supported file formats: x11, pdf, ps, eps, epsi, ncgm or png
[190]2260      ;
2261      ; data type: string
2262      ;
[218]2263      ; example: format_out = "pdf"
2264      ;
[190]2265      ; default:   "x11"
2266      ;***************************************************
2267      if (.not. isvar("format_out"))then         
2268     
2269     
2270         format_out = "x11"
2271         
2272         
2273      end if     
2274      ;***************************************************
2275      ; OPTIONAL --file_out-- OUTPUT FILE
2276      ;
[218]2277      ; path and name of the output file
[190]2278      ;
2279      ; data type: string
2280      ;
[218]2281      ; example: file_out = "/home/path_to_output/test_ts"
2282      ;
[194]2283      ; default:   "~/test_ts"
[190]2284      ;***************************************************
2285      if (.not. isvar("file_out"))then 
2286     
2287       
2288         file_out = "~/test_ts"
2289         
2290         
[566]2291      end if
2292 ;***************************************************
[190]2293      ; OPTIONAL --var-- OUTPUT VARIABLES
2294      ;
[526]2295      ; names of output variables;
2296      ; inidcate them with one comma before and after each variable
2297      ; if the name contains double quotes (") use dq and the '+' operation to
2298      ; concatenate the strings
[190]2299      ;
2300      ; data type: string
[218]2301      ;
[526]2302      ; example: var = ",umax,vmax,wmax,w"+dq+"pt"+dq+"0,"
[190]2303      ;
[218]2304      ; default:   "all" (all available variables are plotted)
[190]2305      ;***************************************************
2306      if (.not. isvar("var"))then
2307     
2308         
2309         var = "all"
2310         
2311         
2312      end if     
2313      ;***************************************************
2314      ; OPTIONAL --start_time_step-- FIRST TIME STEP
2315      ;
[418]2316      ; first time step of plot; there is no need to specify a precise
2317      ; time value from the input file,
[218]2318      ; 'start_time_step' is rounded to the next existent time step
[190]2319      ;
2320      ; data type: double; dimension: [hour]
2321      ;
[418]2322      ; example:   start_time_step = 0.5d
[218]2323      ;
[190]2324      ; default:   -1.d (first available time step)
2325      ;***************************************************
2326      if (.not. isvar("start_time_step"))then
2327     
2328         
2329         start_time_step = -1.d
2330         
2331         
2332      end if     
2333      ;***************************************************
2334      ; OPTIONAL --end_time_step-- LAST TIME STEP
2335      ;
[418]2336      ; last time step of plot; there is no need to specify a precise
2337      ; time value from the input file,
[218]2338      ; 'end_time_step' is rounded to the next existent time step
[190]2339      ;
2340      ; data type: double; dimension: [hour]
2341      ;
[418]2342      ; example:   end_time_step = 2.0d
[218]2343      ;
[190]2344      ; default:   -1.d (last available time step)
2345      ;***************************************************
2346      if (.not. isvar("end_time_step"))then
2347     
2348         
2349         end_time_step = -1.d
2350         
[267]2351
[566]2352      end if             
[190]2353      ;***************************************************
2354      ; OPTIONAL --over-- OVERLAYING
2355      ;
[513]2356      ; predefined overlaying of standard variables are switched
2357      ; on [=1] or off [=0];
[190]2358      ; this feature cannot be used with statistic regions
2359      ;
2360      ; data type: integer
2361      ;
[218]2362      ; example: over = 0
2363      ;
[190]2364      ; default:   0
2365      ;***************************************************
2366      if (.not. isvar("over"))then
2367     
2368         
2369         over = 0
2370         
2371         
[218]2372      end if
2373      ;***************************************************
2374      ; OPTIONAL --norm_t-- NORMALISING X-AXIS OF TIMESERIES PLOTS (T-AXIS)
2375      ;
[513]2376      ; normalising the x-axis with the given value;
[566]2377      ; the unit can be set explicitly;
[418]2378      ; by default no normalising and unit is set to [s]
[218]2379      ;
2380      ; data type: float
2381      ;
[418]2382      ; example: norm_t = 3600.
[218]2383      ;
2384      ; default:   1.
2385      ;***************************************************
[513]2386      if(.not. isvar("norm_t"))then
[218]2387     
2388         
2389         norm_t = 1.
2390         
2391         
2392      end if
2393      ;***************************************************
2394      ; OPTIONAL --unit_t-- UNIT OF X-AXIS OF TIMESERIES PLOTS (T-AXIS)
2395      ;
2396      ; unit of the x-axis can be set; by default unit is set to [s]
2397      ;
2398      ; data type: string
2399      ;
2400      ; example: unit_t = "h"
2401      ;
2402      ; default:   " "
2403      ;***************************************************
[513]2404      if(.not. isvar("unit_t"))then
[218]2405     
2406         
2407         unit_t = " "
2408         
2409         
2410      end if
2411      ;***************************************************
[566]2412      ; OPTIONAL --no_columns-- NUMBER OF COLUMNS
2413      ;
2414      ; number of plots in one row
2415      ;
2416      ; data type: integer
2417      ;
2418      ; example: no_columns = 1
2419      ;
2420      ; default:   1
2421      ;***************************************************
2422      if (.not. isvar("no_columns"))then       
2423     
2424       
2425         no_columns = 1
2426         
2427         
2428      end if     
2429      ;***************************************************
2430      ; OPTIONAL --no_rows-- NUMBER OF ROWS
2431      ;
2432      ; number of plots in one column
2433      ;
2434      ; data type: integer
2435      ;
2436      ; example: no_rows = 2
2437      ;
2438      ; default:   2
2439      ;***************************************************
2440      if (.not. isvar("no_rows"))then
2441     
2442         
2443         no_rows = 2
2444         
2445         
2446      end if     
2447      ;***************************************************
2448      ; OPTIONAL --font_size-- FONT SIZE OF STRINGS
2449      ;
2450      ; font size can be changed to fit to the current plot
2451      ;
2452      ; data type: float
2453      ;
2454      ; example: font_size = 0.07
2455      ;
2456      ; default:   0.07
2457      ;***************************************************
2458      if(.not. isvar("font_size"))then
2459     
2460         
2461         font_size = 0.07
2462         
2463         
2464      end if
2465      ;***************************************************
[190]2466
2467   end if
2468         
2469;*********************************************************       
2470   
2471end
Note: See TracBrowser for help on using the repository browser.