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

Last change on this file since 494 was 418, checked in by heinze, 15 years ago

ncl_preferences is renamed .ncl.config.default

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