Home

Context Navigation

ncl_preferences.ncl @ 412

Last change on this file since 412 was 351, checked in by heinze, 15 years ago
New parameter time_stride is in profiles.ncl available.
File size: 65.9 KB

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

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

Download in other formats:

Original Format