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

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

ncl_preferences is renamed .ncl.config.default

File size: 68.1 KB
Line 
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
12begin
13
14;*********************************************************
15
16script=which_script
17   
18;*********************************************************
19   
20   ;******************************************************
21   ; parameters used by CROSS_SECTIONS
22   ;******************************************************
23   
24   if (script .EQ. "cross_section")then
25   
26      ;***************************************************
27      ; REQUIRED --file_1-- INPUT FILE
28      ;
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');
35      ; otherwise declare path and full file name
36      ;
37      ; data type: string
38      ;
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)
43      ;
44      ; default:   "File in"
45      ;***************************************************
46      if(.not. isvar("file_1"))then
47
48
49         file_1 = "File in"
50         
51         
52      end if
53      ;***************************************************
54      ; REQUIRED if more than one file of a job chain
55      ;          --start_f-- START OF CYCLIC NUMBER
56      ;
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
60      ;
61      ; data type: integer
62      ;
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)
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      ;
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
80      ;
81      ; data type: integer
82      ;
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)
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      ;
98      ; desired section has to be set to 1;
99      ; e.g.: xyc=1, xzc=0 and yzc=0 for xy-section
100      ;
101      ; data type: integer
102      ;
103      ; example: xyc = 1
104      ;          xzc = 0
105      ;          yzc = 0
106      ;
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      ;
127      ; format of the output file;
128      ; supported file formats: x11, pdf, ps, eps, epsi or ncgm
129      ;
130      ; data type: string
131      ;
132      ; example: format_out = "pdf"
133      ;
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      ;
146      ; path and name of the output file
147      ;
148      ; data type: string
149      ;
150      ; example: file_out = "/home/path_to_output/test_cs"
151      ;
152      ; default:   "~/test_cs"
153      ;***************************************************
154      if(.not. isvar("file_out"))then   
155     
156     
157         file_out = "~/test_cs"
158         
159         
160      end if     
161      ;***************************************************
162      ; OPTIONAL --no_columns-- NUMBER OF COLUMNS
163      ;
164      ; number of plots in one row
165      ;
166      ; data type: integer
167      ;
168      ; example: no_columns = 1
169      ;
170      ; default:   1
171      ;***************************************************
172      if(.not. isvar("no_columns"))then 
173     
174     
175         no_columns = 1
176         
177         
178      end if     
179      ;***************************************************
180      ; OPTIONAL --no_rows-- NUMBER OF ROWS
181      ;
182      ; number of plots in one column
183      ;
184      ; data type: integer
185      ;
186      ; example: no_rows = 2
187      ;
188      ; default:   2
189      ;***************************************************
190      if(.not. isvar("no_rows"))then     
191     
192     
193         no_rows = 2
194         
195         
196      end if     
197      ;***************************************************
198      ; OPTIONAL --sort-- TYPE OF SORTING
199      ;
200      ; sequence of plots; sorting either by time step [="time"] or
201      ; by layer/height [="layer"]
202      ;
203      ; data type: string
204      ;
205      ; example: sort = "layer"
206      ;
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      ;
219      ; names of output variables; inidcate them with one comma before and after
220      ; each variable
221      ;
222      ; data type: string
223      ;
224      ; example: var = ",u,v,w,"
225      ;
226      ; default:   "all" (all available variables are plotted)
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      ;
238      ; first time step of plot; there is no need to specify a precise time
239      ; value from the input file,
240      ; 'start_time_step' is rounded to the next existent time step
241      ;
242      ; data type: double; dimension: [hour]
243      ;
244      ; example:   start_time_step = 0.5d
245      ;
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      ;
258      ; last time step of plot; there is no need to specify a precise time
259      ; value from the input file,
260      ; 'end_time_step' is rounded to the next existent time step
261      ;
262      ; data type: double; dimension: [hour]
263      ;
264      ; example:   end_time_step = 2.0d
265      ;
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      ;
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
282      ;
283      ; data type: double; dimension: [meter]
284      ;
285      ; example: xs = 150.0d
286      ;
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      ;
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
303      ;
304      ; data type: double; dimension: [meter]
305      ;
306      ; example: xe = 300.0d
307      ;
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      ;
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
324      ;
325      ; data type: double; dimension: [meter]
326      ;
327      ; example: ys = 150.0d
328      ;
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      ;
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
345      ;
346      ; data type: double; dimension: [meter]
347      ;
348      ; example: ye = 300.0d
349      ;
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      ;
362      ; index for minimum of z-axis, selectable if there are NO preseted layers
363      ; for z (as for XZ or YZ 2D-DATA);
364      ; you cannot specify a meter value from the input file due to grid stretching 
365      ;
366      ; data type: integer
367      ;
368      ; example: zs = 0
369      ;
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      ;
382      ; index for maximum of z-axis, selectable if there are NO preseted layers
383      ; for z (as for XZ or YZ 2D-DATA);
384      ; you cannot specify a meter value from the input file due to grid stretching 
385      ;
386      ; data type: integer
387      ;
388      ; example: ze = 15
389      ;
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      ;
402      ; isoline plots [="Line"], contour plots [="Fill"] or both [="Both"]
403      ; can be choosen
404      ;
405      ; data type: string
406      ;
407      ; example: mode = "Fill"
408      ;
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      ;
421      ; area fill [="AreaFill"], raster fill [="RasterFill"] or
422      ; cell fill [="CellFill"] can be choosen
423      ;
424      ; data type: string
425      ;
426      ; example: fill_mode = "AreaFill"
427      ;
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      ;
440      ; aspect ratio of axis is kept [=1] or not[=0]
441      ;
442      ; data type: integer
443      ;
444      ; example: shape = 1
445      ;
446      ; default:   1
447      ;***************************************************
448      if(.not. isvar("shape"))then
449     
450         
451         shape = 1
452         
453         
454      end if
455      ;***************************************************
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      ;***************************************************
492      ; OPTIONAL --legend_label_stride-- REDUCTION OF THE NUMBER
493      ;                                  OF LABELS IN THE LEGEND
494      ;
495      ; if legend_label_stride is set to 1 every label is shown;
496      ; if set to 2 every second, ...
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      ;
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]
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      ;***************************************************
531      ; OPTIONAL if axes_explicit = 1 --major_ticks_x--
532      ;                              NUMBER OF MAJOR TICK MARKS AT X-AXIS
533      ;
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
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      ;***************************************************
553      ; OPTIONAL if axes_explicit = 1 --major_ticks_y--
554      ;                              NUMBER OF MAJOR TICK MARKS AT Y-AXIS
555      ;
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
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      ;***************************************************
575      ; OPTIONAL if axes_explicit = 1 --norm_x--norm_y--norm_z--
576      ;                              NORMALISE X-, Y- and Z-AXES
577      ;
578      ; normalising the axes with the given value; units can explicitly be set;
579      ; by default no normalising is used and units are [meter]
580      ;
581      ; data type: float
582      ;
583      ; default:   1.0
584      ;***************************************************
585      if(.not. isvar("norm_x"))then     
586         
587         norm_x = 1.0   
588         
589      end if
590      if(.not. isvar("norm_y"))then     
591         
592         norm_y = 1.0   
593         
594      end if
595      if(.not. isvar("norm_z"))then     
596         
597         norm_z = 1.0   
598         
599      end if
600      ;***************************************************
601      ; OPTIONAL if axes_explicit = 1 --unit_x--unit_y--unit_z--
602      ;                              UNITS OF X-, Y- and Z-AXES
603      ;
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]
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      ;***************************************************
629      ; OPTIONAL --vector-- VECTOR PLOT
630      ;
631      ; a vector plot for one defined vector is switched on [=1] or off [=0]
632      ;
633      ; data type: integer
634      ;
635      ; example: vector = 0
636      ;
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      ;
649      ; variables for the vector, vec1 hast to be set to component 1 and
650      ; vec2 to component 2;
651      ; inidcate them with one comma before and after each variable;
652      ;
653      ; data type: string
654      ;
655      ; example: vec1=",u,"
656      ;          vec2=",v,"
657      ;
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      ;
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
677      ;
678      ; data type: string
679      ;
680      ; example: plotvec = ",u,w,"
681      ;
682      ; default:   "plotvec" (no combining, the vector plot is created seperately)
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      ;
698      ; example: ref_mag = 0.05
699      ;
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      ;***************************************************
710
711   end if
712
713   ;******************************************************
714   ; parameters used by PROFILES
715   ;******************************************************
716   
717   if (script .EQ. "profiles")then
718
719      ;***************************************************
720      ; REQUIRED --file_1-- 1ST INPUT FILE
721      ;
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');
726      ; otherwise declare path and full file name
727      ;
728      ; data type: string
729      ;
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)
734      ;
735      ; default:   "File in"
736      ;***************************************************
737      if(.not. isvar("file_1"))then
738
739
740         file_1 = "File in"
741         
742         
743      end if
744      ;***************************************************
745      ; REQUIRED --start_f_1-- START OF CYCLIC NUMBER OF THE 1ST INPUT FILE
746      ;
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
749      ;
750      ; data type: integer
751      ;
752      ; example:   3 (if first file: example_pr.3.nc) and
753      ;            0 (if first file: example_pr.nc)
754      ;
755      ; default:   -1 
756      ;***************************************************
757      if (.not. isvar("start_f_1"))then 
758     
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      ;
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
769      ;
770      ; data type: integer
771      ;
772      ; example:   48 (if last file: example_pr.48.nc) and
773      ;             0 (if last file: example_pr.nc)
774      ;
775      ; default:   -1 
776      ;***************************************************
777      if (.not. isvar("end_f_1"))then 
778     
779         
780         end_f_1 = -1
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      ;
803      ; declare path and name of the output file
804      ;
805      ; data type: string
806      ;
807      ; default:   "~/test_pr"
808      ;***************************************************
809      if (.not. isvar("file_out"))then 
810     
811       
812         file_out = "~/test_pr"
813         
814         
815      end if     
816      ;***************************************************
817      ; OPTIONAL --no_columns-- NUMBER OF COLUMNS
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      ;***************************************************
833      ; OPTIONAL --no_rows-- NUMBER OF ROWS
834      ;
835      ; number of plots in one column
836      ;
837      ; data type: integer
838      ;
839      ; default:   2
840      ;***************************************************
841      if (.not. isvar("no_rows"))then   
842     
843       
844         no_rows = 2
845         
846         
847      end if     
848      ;***************************************************
849      ; OPTIONAL --var-- OUTPUT VARIABLES
850      ;
851      ; names of output variables;
852      ; inidcate them with one comma before and after each variable
853      ;
854      ; data type: string
855      ;
856      ; example: var = ",u,v,w,";
857      ;
858      ; default:   "all" (all available variables are plotted)
859      ;***************************************************
860      if (.not. isvar("var"))then
861     
862         
863         var = "all"
864         
865         
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         
881         legend = 1
882         
883         
884      end if   
885      ;***************************************************
886      ; OPTIONAL ---no_files-- NUMBER OF INPUT FILES
887      ;
888      ; up to 6 input files of different runs can be used; all input files
889      ; must contain the same variables and dimensionalties;
890      ; identical variables are merged in one plot;
891      ; combining or overlaying of any variables is not feasible;
892      ; an extra legend (besides the time step legend) is drawn where you
893      ; can explicitily name the items
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         
921      end if
922      ;***************************************************
923      ; REQUIRED if no_files > 1 --file_2-- 2ND INPUT FILE
924      ;
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');
929      ; otherwise declare path and full file name
930      ;
931      ; data type: string
932      ;
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)
937      ;
938      ; default:   "File in"
939      ;***************************************************
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      ;
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
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 
960     
961         
962         start_f_2 = -1
963         
964         
965      end if
966      ;***************************************************
967      ; REQUIRED --end_f_2-- END OF CYCLIC NUMBER OF THE 2ND INPUT FILE
968      ;
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
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      ;***************************************************
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         
1000      end if
1001      ;***************************************************
1002      ; REQUIRED if no_files > 2 --file_3-- 3RD INPUT FILE
1003      ;
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');
1008      ; otherwise declare path and full file name
1009      ;
1010      ; data type: string
1011      ;
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)
1016      ;
1017      ; default:   "File in"
1018      ;***************************************************
1019      if(.not. isvar("file_3"))then
1020
1021
1022         file_3 = "File in"
1023         
1024         
1025      end if
1026      ;***************************************************
1027      ; REQUIRED --start_f_3-- START OF CYCLIC NUMBER OF THE 3RD INPUT FILE
1028      ;
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
1031      ;
1032      ; data type: integer
1033      ;
1034      ; example:   3 (if first file: example_ts.3.nc)
1035      ;
1036      ; default:   -1 
1037      ;***************************************************
1038      if (.not. isvar("start_f_3"))then 
1039     
1040         
1041         start_f_3 = -1
1042         
1043         
1044      end if
1045      ;***************************************************
1046      ; REQUIRED --end_f_3-- END OF CYCLIC NUMBER OF THE 3RD INPUT FILE
1047      ;
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
1050      ;
1051      ; data type: integer
1052      ;
1053      ; example:   48 (if last file: example_ts.48.nc)
1054      ;
1055      ; default:   -1 
1056      ;***************************************************
1057      if (.not. isvar("end_f_3"))then 
1058     
1059         
1060         end_f_3 = -1
1061         
1062         
1063      end if   
1064      ;***************************************************
1065      ; OPTIONAL --name_legend_3-- 3rd LEGEND ITEM
1066      ;
1067      ; if no_files > 2: name of 3rd item in extra legend
1068      ;
1069      ; data type: string
1070      ;
1071      ; default:   "File 3"
1072      ;***************************************************
1073      if (.not. isvar("name_legend_3"))then
1074     
1075           
1076         name_legend_3 = "File 3"
1077         
1078         
1079      end if
1080      ;***************************************************
1081      ; REQUIRED if no_files > 3 --file_4-- 4TH INPUT FILE
1082      ;
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');
1087      ; otherwise declare path and full file name
1088      ;
1089      ; data type: string
1090      ;
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)
1095      ;
1096      ; default:   "File in"
1097      ;***************************************************
1098      if(.not. isvar("file_4"))then
1099
1100
1101         file_4 = "File in"
1102         
1103         
1104      end if
1105      ;***************************************************
1106      ; REQUIRED --start_f_4-- START OF CYCLIC NUMBER OF THE 4TH INPUT FILE
1107      ;
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
1110      ;
1111      ; data type: integer
1112      ;
1113      ; example:   3 (if first file: example_ts.3.nc)
1114      ;
1115      ; default:   -1 
1116      ;***************************************************
1117      if (.not. isvar("start_f_4"))then 
1118     
1119         
1120         start_f_4 = -1
1121         
1122         
1123      end if
1124      ;***************************************************
1125      ; REQUIRED --end_f_4-- END OF CYCLIC NUMBER OF THE 4TH INPUT FILE
1126      ;
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
1129      ;
1130      ; data type: integer
1131      ;
1132      ; example:   48 (if last file: example_ts.48.nc)
1133      ;
1134      ; default:   -1 
1135      ;***************************************************
1136      if (.not. isvar("end_f_4"))then 
1137     
1138         
1139         end_f_4 = -1
1140         
1141         
1142      end if 
1143      ;***************************************************
1144      ; OPTIONAL --name_legend_4-- 4th LEGEND ITEM
1145      ;
1146      ; if no_files > 3: name of 4th item in extra legend
1147      ;
1148      ; data type: string
1149      ;
1150      ; default:   "File 4"
1151      ;***************************************************
1152      if (.not. isvar("name_legend_4"))then   
1153     
1154       
1155         name_legend_4 = "File 4"
1156         
1157         
1158      end if
1159      ;***************************************************
1160      ; REQUIRED if no_files > 4 --file_5-- 5TH INPUT FILE
1161      ;
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');
1166      ; otherwise declare path and full file name
1167      ;
1168      ; data type: string
1169      ;
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)
1174      ;
1175      ; default:   "File in"
1176      ;***************************************************
1177      if(.not. isvar("file_5"))then
1178
1179
1180         file_5 = "File in"
1181         
1182         
1183      end if
1184      ;***************************************************
1185      ; REQUIRED --start_f_5-- START OF CYCLIC NUMBER OF THE 5TH INPUT FILE
1186      ;
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
1189      ;
1190      ; data type: integer
1191      ;
1192      ; example:   3 (if first file: example_ts.3.nc)
1193      ;
1194      ; default:   -1 
1195      ;***************************************************
1196      if (.not. isvar("start_f_5"))then 
1197     
1198         
1199         start_f_5 = -1
1200         
1201         
1202      end if
1203      ;***************************************************
1204      ; REQUIRED --end_f_5-- END OF CYCLIC NUMBER OF THE 5TH INPUT FILE
1205      ;
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
1208      ;
1209      ; data type: integer
1210      ;
1211      ; example:   48 (if last file: example_ts.48.nc)
1212      ;
1213      ; default:   -1 
1214      ;***************************************************
1215      if (.not. isvar("end_f_5"))then 
1216     
1217         
1218         end_f_5 = -1
1219         
1220         
1221      end if     
1222      ;***************************************************
1223      ; OPTIONAL --name_legend_5-- 5th LEGEND ITEM
1224      ;
1225      ; if no_files > 4: name of 5th item in extra legend
1226      ;
1227      ; data type: string
1228      ;
1229      ; default:   "File 5"
1230      ;***************************************************
1231      if (.not. isvar("name_legend_5"))then
1232     
1233           
1234         name_legend_5 = "File 5"
1235         
1236         
1237      end if
1238      ;***************************************************
1239      ; REQUIRED if no_files > 5 --file_6-- 6TH INPUT FILE
1240      ;
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');
1245      ; otherwise declare path and full file name
1246      ;
1247      ; data type: string
1248      ;
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)
1253      ;
1254      ; default:   "File in"
1255      ;***************************************************
1256      if(.not. isvar("file_6"))then
1257
1258
1259         file_6 = "File in"
1260         
1261         
1262      end if
1263      ;***************************************************
1264      ; REQUIRED --start_f_6-- START OF CYCLIC NUMBER OF THE 6TH INPUT FILE
1265      ;
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
1268      ;
1269      ; data type: integer
1270      ;
1271      ; example:   3 (if first file: example_ts.3.nc)
1272      ;
1273      ; default:   -1 
1274      ;***************************************************
1275      if (.not. isvar("start_f_6"))then 
1276     
1277         
1278         start_f_6 = -1
1279         
1280         
1281      end if
1282      ;***************************************************
1283      ; REQUIRED --end_f_6-- END OF CYCLIC NUMBER OF THE 6TH INPUT FILE
1284      ;
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
1287      ;
1288      ; data type: integer
1289      ;
1290      ; example:   48 (if last file: example_ts.48.nc)
1291      ;
1292      ; default:   -1 
1293      ;***************************************************
1294      if (.not. isvar("end_f_6"))then 
1295     
1296         
1297         end_f_6 = -1
1298         
1299         
1300      end if                     
1301      ;***************************************************
1302      ; OPTIONAL --name_legend_6-- 6TH LEGEND ITEM
1303      ;
1304      ; if no_files > 5: name of 6th item in extra legend
1305      ;
1306      ; data type: string
1307      ;
1308      ; default:   "File 6"
1309      ;***************************************************
1310      if (.not. isvar("name_legend_6"))then
1311     
1312           
1313         name_legend_6 = "File 6"
1314         
1315         
1316      end if             
1317      ;***************************************************
1318      ; OPTIONAL --start_time_step-- FIRST TIME STEP
1319      ;
1320      ; first time step of plot; there is no need to specify a precise
1321      ; time value from the input file,
1322      ; 'start_time_step' is rounded to the next existent time step
1323      ;
1324      ; data type: double; dimension: [hour]
1325      ;
1326      ; example:   start_time_step = 0.5d
1327      ;
1328      ; default:   -1.d (first available time step)
1329      ;***************************************************
1330      if (.not. isvar("start_time_step"))then   
1331     
1332       
1333         start_time_step = -1.d
1334         
1335         
1336      end if     
1337      ;***************************************************
1338      ; OPTIONAL --end_time_step-- LAST TIME STEP
1339      ;
1340      ; last time step of plot; there is no need to specify a precise
1341      ; time value from the input file,
1342      ; 'end_time_step' is rounded to the next existent time step
1343      ;
1344      ; data type: double; dimension: [hour]
1345      ;
1346      ; example:   end_time_step = 2.0d
1347      ;
1348      ; default:   -1.d (last available time step)
1349      ;***************************************************
1350      if (.not. isvar("end_time_step"))then
1351     
1352         
1353         end_time_step = -1.d   
1354         
1355         
1356      end if
1357      ;***************************************************
1358      ; OPTIONAL --time_stride-- TEMPORAL STRIDE
1359      ;
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
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      ;***************************************************
1380      ; OPTIONAL --start_x-- STARTVALUE IN X-DIRECTION 
1381      ;
1382      ; if 3D data is used for profiles you can choose a startvalue
1383      ; in x-direction for horizontal averaging
1384      ;
1385      ; data type: integer; dimension [grippoint]
1386      ;
1387      ; example: start_x = 0
1388      ;
1389      ; default:   0 (minimum of range)
1390      ;***************************************************
1391      if (.not. isvar("start_x"))then   
1392     
1393         
1394         start_x = 0
1395         
1396         
1397      end if
1398      ;***************************************************
1399      ; OPTIONAL --end_x-- ENDVALUE IN X-DIRECTION 
1400      ;
1401      ; if 3D data is used for profiles you can choose an endvalue
1402      ; in x-direction for horizontal averaging
1403      ;
1404      ; data type: integer; dimension [grippoint]
1405      ;
1406      ; example: end_x = 200
1407      ;
1408      ; default:   -1 (maximum of range excluding ghostpoint)
1409      ;***************************************************
1410      if (.not. isvar("end_x"))then   
1411     
1412         
1413         end_x = -1
1414         
1415         
1416      end if
1417      ;***************************************************
1418      ; OPTIONAL --start_y-- STARTVALUE IN Y-DIRECTION 
1419      ;
1420      ; if 3D data is used for profiles you can choose a startvalue
1421      ; in y-direction for horizontal averaging
1422      ;
1423      ; data type: integer; dimension [grippoint]
1424      ;
1425      ; example: start_y = 0
1426      ;
1427      ; default:   0 (minimum of range)
1428      ;***************************************************
1429      if (.not. isvar("start_y"))then   
1430     
1431         
1432         start_y = 0
1433         
1434         
1435      end if
1436      ;***************************************************
1437      ; OPTIONAL --end_y-- ENDVALUE IN Y-DIRECTION 
1438      ;
1439      ; if 3D data is used for profiles you can choose an endvalue
1440      ; in y-direction for horizontal averaging
1441      ;
1442      ; data type: integer; dimension [grippoint]
1443      ;
1444      ; example: end_y = 200
1445      ;
1446      ; default:   -1 (maximum of range excluding ghostpoint)
1447      ;***************************************************
1448      if (.not. isvar("end_y"))then   
1449     
1450         
1451         end_y = -1
1452         
1453         
1454      end if     
1455      ;***************************************************
1456      ; OPTIONAL --xs-- MINIMUM X-AXIS
1457      ;
1458      ; value for minimum of x-axis
1459      ;
1460      ; data type: float; dimension of variable
1461      ;
1462      ; example: xs = 300. (e.g. if var=",pt")
1463      ;
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      ;
1480      ; example: xe = 302. (e.g. if var=",pt")
1481      ;
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      ;
1498      ; example: min_z = 0.0d
1499      ;
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      ;
1516      ; example: max_z = 1000.0d
1517      ;
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      ;***************************************************
1528      ; OPTIONAL --log_z-- LOGARITHMIC SCALE OF Z-AXIS
1529      ;
1530      ; logarithmic scale for z-axis is switched on [=1] or off [=0]
1531      ;
1532      ; data type: integer
1533      ;
1534      ; example: log_z = 0
1535      ;
1536      ; default:   0
1537      ;***************************************************
1538      if (.not. isvar("log_z"))then
1539     
1540           
1541         log_z = 0
1542         
1543         
1544      end if     
1545      ;***************************************************
1546      ; OPTIONAL --norm_z-- NORMALISING Z-AXIS
1547      ;
1548      ; value for normalising the z-axis
1549      ;
1550      ; data type: float
1551      ;
1552      ; example: norm_z = 1.
1553      ;
1554      ; default:   1. (no normalising)
1555      ;***************************************************
1556      if (.not. isvar("norm_z"))then
1557     
1558           
1559         norm_z = 1.
1560         
1561         
1562      end if     
1563      ;***************************************************
1564      ; OPTIONAL --over-- OVERLAYING
1565      ;
1566      ; predefined overlaying of standard variables is switched on [=1] or off [=0];
1567      ; this feature cannot be used with statistic regions
1568      ;
1569      ; data type: integer
1570      ;
1571      ; example: over = 0
1572      ;
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      ;
1585      ; combining of two or three individual variables is possible and
1586      ; switched on [=1] or off[=0]
1587      ;
1588      ; data type: integer
1589      ;
1590      ; example: combine = 0
1591      ;
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      ;
1608      ; example: number_comb = 2
1609      ;
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      ;
1622      ; names of variables; inidcate them with one comma before and
1623      ; after each variable;
1624      ;
1625      ; data type: string
1626      ;
1627      ; example: c_var = ",u,v,"
1628      ;
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      ;
1645      ; example: black = 0
1646      ;
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      ;
1659      ; use of different line patterns [=1] or continous lines [=0]
1660      ; for different time steps
1661      ;
1662      ; data type: integer
1663      ;
1664      ; example: dash = 0
1665      ;
1666      ; default:   0
1667      ;***************************************************
1668      if (.not. isvar("dash"))then     
1669     
1670       
1671         dash = 0     
1672         
1673         
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         
1710      end if     
1711      ;***************************************************
1712       
1713   end if
1714
1715   ;******************************************************
1716   ; parameters used by SPECTRA
1717   ;******************************************************
1718   
1719   if (script .EQ. "spectra") then
1720
1721      ;***************************************************
1722      ; REQUIRED --file_1-- INPUT FILE
1723      ;
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');
1728      ; otherwise declare path and full file name
1729      ;
1730      ; data type: string
1731      ;
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)
1736      ;
1737      ; default:   "File in"
1738      ;***************************************************
1739      if(.not. isvar("file_1"))then
1740
1741
1742         file_1 = "File in"
1743         
1744         
1745      end if
1746      ;***************************************************
1747      ; REQUIRED if more than one file of a job chain
1748      ;           --start_f-- START OF CYCLIC NUMBER
1749      ;
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
1752      ;
1753      ; data type: integer
1754      ;
1755      ; example:   start_f = 3 (if first file: example_sp.3.nc) and
1756      ;            start_f = 0 (if first file: example_sp.nc)
1757      ;
1758      ; default:   -1 
1759      ;***************************************************
1760      if (.not. isvar("start_f"))then 
1761     
1762         
1763         start_f = -1
1764         
1765         
1766      end if
1767      ;***************************************************
1768      ; REQUIRED if more than one file of a job chain
1769      ;          --end_f-- END OF CYCLIC NUMBER
1770      ;
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
1773      ;
1774      ; data type: integer
1775      ;
1776      ; example:   end_f = 48 (if last file: example_sp.48.nc) and
1777      ;            end_f =  0 (if last file: example_sp.nc)
1778      ;
1779      ; default:   -1 
1780      ;***************************************************
1781      if (.not. isvar("end_f"))then 
1782     
1783         
1784         end_f = -1
1785         
1786         
1787      end if
1788      ;***************************************************
1789      ; OPTIONAL --format_out-- OUTPUT FORMAT
1790      ;
1791      ; format of the output file;
1792      ; supported file formats: x11, pdf, ps, eps, epsi or ncgm
1793      ;
1794      ; data type: string
1795      ;
1796      ; example: format_out = "pdf"
1797      ;
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      ;
1810      ; path and name of the output file
1811      ;
1812      ; data type: string
1813      ;
1814      ; example: file_out = "/home/path_to_output/test_sp"
1815      ;
1816      ; default:   "~/test_sp"
1817      ;***************************************************
1818      if (.not. isvar("file_out"))then 
1819     
1820       
1821         file_out = "~/test_sp"
1822         
1823         
1824      end if     
1825      ;***************************************************
1826      ; OPTIONAL --no_columns-- NUMBER OF COLUMNS
1827      ;
1828      ; number of plots in one row
1829      ;
1830      ; data type: integer
1831      ;
1832      ; example: no_columns = 1
1833      ;
1834      ; default:   1
1835      ;***************************************************
1836      if (.not. isvar("no_columns"))then         
1837     
1838     
1839         no_columns = 1
1840         
1841         
1842      end if     
1843      ;***************************************************
1844      ; OPTIONAL --no_rows-- NUMBER OF ROWS
1845      ;
1846      ; number of plots in one column
1847      ;
1848      ; data type: integer
1849      ;
1850      ; example: no_rows = 2
1851      ;
1852      ; default:   2
1853      ;***************************************************
1854      if (.not. isvar("no_rows"))then
1855     
1856         
1857         no_rows = 2
1858         
1859         
1860      end if     
1861      ;***************************************************
1862      ; OPTIONAL --var-- OUTPUT VARIABLES
1863      ;
1864      ; names of output variables; inidcate them with one comma before
1865      ; and after each variable
1866      ;
1867      ; data type: string
1868      ;
1869      ; example: var = ",u,v,w,"
1870      ;
1871      ; default:   "all" (all available variables are plotted)
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      ;
1883      ; defines if all levels or selected levels are output;
1884      ; a 1D array with up to 100 elements
1885      ; from 0 to 99 must be indicated
1886      ;
1887      ; data type: integer array
1888      ;
1889      ; example: height_level=(/0,3,7/)
1890      ;
1891      ; default:   (/-1/) (all levels are output)
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      ;
1903      ; defines if all heights [="height"] or all timesteps [="time"]
1904      ; are merged in one plot
1905      ;
1906      ; data type: string
1907      ;
1908      ; example: sort = "height"
1909      ;
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      ;
1922      ; first time step of plot; there is no need to specify a precise
1923      ; time value from the input file,
1924      ; 'start_time_step' is rounded to the next existent time step
1925      ;
1926      ; data type: double; dimension: [hour]
1927      ;
1928      ; example:   start_time_step = 0.5d
1929      ;
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      ;
1942      ; last time step of plot; there is no need to specify a precise
1943      ; time value from the input file,
1944      ; 'end_time_step' is rounded to the next existent time step
1945      ;
1946      ; data type: double; dimension: [hour]
1947      ;
1948      ; example:   end_time_step = 2.0d
1949      ;
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      ;
1966      ; example: black = 0
1967      ;
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      ;
1984      ; example: dash = 0
1985      ;
1986      ; default:   0
1987      ;***************************************************
1988      if (.not. isvar("dash"))then   
1989     
1990         
1991         dash = 0
1992         
1993         
1994      end if
1995      ;***************************************************
1996      ; OPTIONAL --log_x-- LOGARITHMIC SCALE X-AXIS
1997      ;
1998      ; logarithmic scale for x-axis is switched on [=1] or off [=0]
1999      ;
2000      ; data type: integer
2001      ;
2002      ; example: log_x = 1
2003      ;
2004      ; default:   1
2005      ;***************************************************
2006      if (.not. isvar("log_x"))then   
2007     
2008         
2009         log_x = 1
2010         
2011         
2012      end if     
2013      ;***************************************************
2014      ; OPTIONAL --log_y-- LOGARITHMIC SCALE Y-AXIS
2015      ;
2016      ; logarithmic scale for y-axis is switched on [=1] or off [=0]
2017      ;
2018      ; data type: integer
2019      ;
2020      ; example: log_y = 1
2021      ;
2022      ; default:   1
2023      ;***************************************************
2024      if (.not. isvar("log_y"))then   
2025     
2026         
2027         log_y = 1
2028         
2029         
2030      end if     
2031      ;***************************************************
2032      ; OPTIONAL --norm_x-- NORMALISING X-AXIS
2033      ;
2034      ; value for normalising the x-axis
2035      ;
2036      ; data type: float
2037      ;
2038      ; example: norm_x = 1
2039      ;
2040      ; default:   1. (no normalising)
2041      ;***************************************************
2042      if (.not. isvar("norm_x"))then   
2043     
2044         
2045         norm_x = 1.
2046         
2047         
2048      end if
2049      ;***************************************************
2050      ; OPTIONAL --norm_height-- NORMALISING X-AXIS WITH HEIGHT
2051      ;
2052      ; normalising the x-axis with the height of each spectra
2053      ; is switched on [=1] or off [=0]
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         
2067      end if     
2068      ;***************************************************
2069      ; OPTIONAL --norm_y-- NORMALISE Y-AXIS
2070      ;
2071      ; value for normalising the y-axis
2072      ;
2073      ; data type: float
2074      ;
2075      ; example: norm_y = 1.
2076      ;
2077      ; default:   1. (no normalising)
2078      ;***************************************************
2079      if (.not. isvar("norm_y"))then     
2080     
2081       
2082         norm_y = 1.
2083         
2084         
2085      end if
2086      ;***************************************************
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      ;***************************************************
2144     
2145   end if
2146
2147   ;******************************************************
2148   ; parameters used by TIMESERIES
2149   ;******************************************************
2150
2151   if (script .EQ. "timeseries") then
2152
2153      ;***************************************************
2154      ; REQUIRED --file_1-- INPUT FILE
2155      ;
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');
2160      ; otherwise declare path and full file name
2161      ;
2162      ; data type: string
2163      ;
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)
2168      ;
2169      ; default:   "File in"
2170      ;***************************************************
2171      if(.not. isvar("file_1"))then
2172
2173
2174         file_1 = "File in"
2175         
2176         
2177      end if
2178      ;***************************************************
2179      ; REQUIRED if more than one file of a job chain
2180      ;          --start_f-- START OF CYCLIC NUMBER
2181      ;
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
2184      ;
2185      ; data type: integer
2186      ;
2187      ; example:   start_f = 3 (if first file: example_ts.3.nc) and
2188      ;            start_f = 0 (if first file: example_ts.nc)
2189      ;
2190      ; default:   -1 
2191      ;***************************************************
2192      if (.not. isvar("start_f"))then 
2193     
2194         
2195         start_f = -1
2196         
2197         
2198      end if
2199      ;***************************************************
2200      ; REQUIRED if more than one file of a job chain
2201      ;           --end_f-- END OF CYCLIC NUMBER
2202      ;
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
2205      ;
2206      ; data type: integer
2207      ;
2208      ; example:   end_f = 48 (if last file: example_ts.48.nc) and
2209      ;            end_f = 0 (if last file: example_ts.nc)
2210      ;
2211      ; default:   -1 
2212      ;***************************************************
2213      if (.not. isvar("end_f"))then 
2214     
2215         
2216         end_f = -1
2217         
2218         
2219      end if
2220      ;***************************************************
2221      ; OPTIONAL --format_out-- OUTPUT FORMAT
2222      ;
2223      ; format of the output file;
2224      ; supported file formats: x11, pdf, ps, eps, epsi or ncgm
2225      ;
2226      ; data type: string
2227      ;
2228      ; example: format_out = "pdf"
2229      ;
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      ;
2242      ; path and name of the output file
2243      ;
2244      ; data type: string
2245      ;
2246      ; example: file_out = "/home/path_to_output/test_ts"
2247      ;
2248      ; default:   "~/test_ts"
2249      ;***************************************************
2250      if (.not. isvar("file_out"))then 
2251     
2252       
2253         file_out = "~/test_ts"
2254         
2255         
2256      end if     
2257      ;***************************************************
2258      ; OPTIONAL --no_columns-- NUMBER OF COLUMNS
2259      ;
2260      ; number of plots in one row
2261      ;
2262      ; data type: integer
2263      ;
2264      ; example: no_columns = 1
2265      ;
2266      ; default:   1
2267      ;***************************************************
2268      if (.not. isvar("no_columns"))then       
2269     
2270       
2271         no_columns = 1
2272         
2273         
2274      end if     
2275      ;***************************************************
2276      ; OPTIONAL --no_rows-- NUMBER OF ROWS
2277      ;
2278      ; number of plots in one column
2279      ;
2280      ; data type: integer
2281      ;
2282      ; example: no_rows = 2
2283      ;
2284      ; default:   2
2285      ;***************************************************
2286      if (.not. isvar("no_rows"))then
2287     
2288         
2289         no_rows = 1
2290         
2291         
2292      end if     
2293      ;***************************************************
2294      ; OPTIONAL --var-- OUTPUT VARIABLES
2295      ;
2296      ; names of output variables; inidcate them with one comma before
2297      ; and after each variable
2298      ;
2299      ; data type: string
2300      ;
2301      ; example: var = ",u,v,w,"
2302      ;
2303      ; default:   "all" (all available variables are plotted)
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      ;
2315      ; first time step of plot; there is no need to specify a precise
2316      ; time value from the input file,
2317      ; 'start_time_step' is rounded to the next existent time step
2318      ;
2319      ; data type: double; dimension: [hour]
2320      ;
2321      ; example:   start_time_step = 0.5d
2322      ;
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      ;
2335      ; last time step of plot; there is no need to specify a precise
2336      ; time value from the input file,
2337      ; 'end_time_step' is rounded to the next existent time step
2338      ;
2339      ; data type: double; dimension: [hour]
2340      ;
2341      ; example:   end_time_step = 2.0d
2342      ;
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         
2350
2351      end if     
2352      ;***************************************************
2353      ; OPTIONAL --over-- OVERLAYING
2354      ;
2355      ; predefined overlaying of standard variables are switched on [=1] or off [=0];
2356      ; this feature cannot be used with statistic regions
2357      ;
2358      ; data type: integer
2359      ;
2360      ; example: over = 0
2361      ;
2362      ; default:   0
2363      ;***************************************************
2364      if (.not. isvar("over"))then
2365     
2366         
2367         over = 0
2368         
2369         
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         
2388      end if
2389      ;***************************************************
2390      ; OPTIONAL --norm_t-- NORMALISING X-AXIS OF TIMESERIES PLOTS (T-AXIS)
2391      ;
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]
2394      ;
2395      ; data type: float
2396      ;
2397      ; example: norm_t = 3600.
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      ;***************************************************
2427
2428   end if
2429         
2430;*********************************************************       
2431   
2432end
Note: See TracBrowser for help on using the repository browser.