Home

Context Navigation

.ncl.config.default @ 431

Last change on this file since 431 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
5	command = "ksh 'if [ ! -f $HOME/.hluresfile ]; then; "+\
6	"cp $PALM_BIN/NCL/.hluresfile $HOME; fi'"
7	system(command)
8
9	;*********************************************************
10
11
12	begin
13
14	;*********************************************************
15
16	script=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
2432	end

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format