Home

Context Navigation

ncl_preferences.ncl @ 350

Last change on this file since 350 was 342, checked in by heinze, 15 years ago
Plots of profiles appear in the order specified by the variable var.
File size: 65.2 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 --start_x-- STARTVALUE IN X-DIRECTION
1259	;
1260	; if 3D data is used for profiles you can choose a startvalue in x-direction for horizontal averaging
1261	;
1262	; data type: integer; dimension [grippoint]
1263	;
1264	; example: start_x = 0
1265	;
1266	; default: 0 (minimum of range)
1267	;***************************************************
1268	if (.not. isvar("start_x"))then
1269
1270
1271	start_x = 0
1272
1273
1274	end if
1275	;***************************************************
1276	; OPTIONAL --end_x-- ENDVALUE IN X-DIRECTION
1277	;
1278	; if 3D data is used for profiles you can choose an endvalue in x-direction for horizontal averaging
1279	;
1280	; data type: integer; dimension [grippoint]
1281	;
1282	; example: end_x = 200
1283	;
1284	; default: -1 (maximum of range excluding ghostpoint)
1285	;***************************************************
1286	if (.not. isvar("end_x"))then
1287
1288
1289	end_x = -1
1290
1291
1292	end if
1293	;***************************************************
1294	; OPTIONAL --start_y-- STARTVALUE IN Y-DIRECTION
1295	;
1296	; if 3D data is used for profiles you can choose a startvalue in y-direction for horizontal averaging
1297	;
1298	; data type: integer; dimension [grippoint]
1299	;
1300	; example: start_y = 0
1301	;
1302	; default: 0 (minimum of range)
1303	;***************************************************
1304	if (.not. isvar("start_y"))then
1305
1306
1307	start_y = 0
1308
1309
1310	end if
1311	;***************************************************
1312	; OPTIONAL --end_y-- ENDVALUE IN Y-DIRECTION
1313	;
1314	; if 3D data is used for profiles you can choose an endvalue in y-direction for horizontal averaging
1315	;
1316	; data type: integer; dimension [grippoint]
1317	;
1318	; example: end_y = 200
1319	;
1320	; default: -1 (maximum of range excluding ghostpoint)
1321	;***************************************************
1322	if (.not. isvar("end_y"))then
1323
1324
1325	end_y = -1
1326
1327
1328	end if
1329	;***************************************************
1330	; OPTIONAL --xs-- MINIMUM X-AXIS
1331	;
1332	; value for minimum of x-axis
1333	;
1334	; data type: float; dimension of variable
1335	;
1336	; example: xs = 300. (e.g. if var=",pt")
1337	;
1338	; default: -1. (minimum of range)
1339	;***************************************************
1340	if (.not. isvar("xs"))then
1341
1342
1343	xs = -1.
1344
1345
1346	end if
1347	;***************************************************
1348	; OPTIONAL --xe-- MAXIMUM X-AXIS
1349	;
1350	; value for maximum of x-axis
1351	;
1352	; data type: float; dimension of variable
1353	;
1354	; example: xe = 302. (e.g. if var=",pt")
1355	;
1356	; default: -1. (maximum of range)
1357	;***************************************************
1358	if (.not. isvar("xe"))then
1359
1360
1361	xe = -1.
1362
1363
1364	end if
1365	;***************************************************
1366	; OPTIONAL --min_z-- MINIMUM Z-AXIS
1367	;
1368	; value for minimum height of z-axis
1369	;
1370	; data type: double; dimension: [meter]
1371	;
1372	; example: min_z = 0.0
1373	;
1374	; default: -1.d (related to the minimum of z-axis in input file)
1375	;***************************************************
1376	if (.not. isvar("min_z"))then
1377
1378
1379	min_z = -1.d
1380
1381
1382	end if
1383	;***************************************************
1384	; OPTIONAL --max_z-- MAXIMUM Z-AXIS
1385	;
1386	; value for maximum height of z-axis
1387	;
1388	; data type: double; dimension: [meter]
1389	;
1390	; example: max_z = 1000.0
1391	;
1392	; default: -1.d (related to the maximum of z-axis in input file)
1393	;***************************************************
1394	if (.not. isvar("max_z"))then
1395
1396
1397	max_z = -1.d
1398
1399
1400	end if
1401	;***************************************************
1402	; OPTIONAL --log_z-- LOGARITHMIC SCALE OF Z-AXIS
1403	;
1404	; logarithmic scale for z-axis is switched on [=1] or off [=0]
1405	;
1406	; data type: integer
1407	;
1408	; example: log_z = 0
1409	;
1410	; default: 0
1411	;***************************************************
1412	if (.not. isvar("log_z"))then
1413
1414
1415	log_z = 0
1416
1417
1418	end if
1419	;***************************************************
1420	; OPTIONAL --norm_z-- NORMALISING Z-AXIS
1421	;
1422	; value for normalising the z-axis
1423	;
1424	; data type: float
1425	;
1426	; example: norm_z = 1.
1427	;
1428	; default: 1. (no normalising)
1429	;***************************************************
1430	if (.not. isvar("norm_z"))then
1431
1432
1433	norm_z = 1.
1434
1435
1436	end if
1437	;***************************************************
1438	; OPTIONAL --over-- OVERLAYING
1439	;
1440	; predefined overlaying of standard variables is switched on [=1] or off [=0];
1441	; this feature cannot be used with statistic regions
1442	;
1443	; data type: integer
1444	;
1445	; example: over = 0
1446	;
1447	; default: 0
1448	;***************************************************
1449	if (.not. isvar("over"))then
1450
1451
1452	over = 0
1453
1454
1455	end if
1456	;***************************************************
1457	; OPTIONAL --combine-- COMBINING
1458	;
1459	; combining of two or three individual variables is possible and switched on [=1] or off[=0]
1460	;
1461	; data type: integer
1462	;
1463	; example: combine = 0
1464	;
1465	; default: 0
1466	;***************************************************
1467	if (.not. isvar("combine"))then
1468
1469
1470	combine = 0
1471
1472
1473	end if
1474	;***************************************************
1475	; REQUIRED IF combine=1 --number_comb-- NUMBER FOR COMBINING
1476	;
1477	; number of variables for combining; 2 or 3 is valid
1478	;
1479	; data type: integer
1480	;
1481	; example: number_comb = 2
1482	;
1483	; default: -1
1484	;***************************************************
1485	if (.not. isvar("number_comb"))then
1486
1487
1488	number_comb = -1
1489
1490
1491	end if
1492	;***************************************************
1493	; REQUIRED IF combine=1 --c_var-- VARIABLES FOR COMBINING
1494	;
1495	; names of variables; inidcate them with one comma before and after each variable;
1496	;
1497	; data type: string
1498	;
1499	; example: c_var = ",u,v,"
1500	;
1501	; default: "c_variables"
1502	;***************************************************
1503	if (.not. isvar("c_var"))then
1504
1505
1506	c_var = "c_variables"
1507
1508
1509	end if
1510	;***************************************************
1511	; OPTIONAL --black-- BLACK OR COLOR LINES
1512	;
1513	; colored [=0] or black and white [=1] plots
1514	;
1515	; data type: integer
1516	;
1517	; example: black = 0
1518	;
1519	; default: 0
1520	;***************************************************
1521	if (.not. isvar("black"))then
1522
1523
1524	black = 0
1525
1526
1527	end if
1528	;***************************************************
1529	; OPTIONAL --dash-- DASHED LINES
1530	;
1531	; use of different line patterns [=1] or continous lines [=0] for different time steps
1532	;
1533	; data type: integer
1534	;
1535	; example: dash = 0
1536	;
1537	; default: 0
1538	;***************************************************
1539	if (.not. isvar("dash"))then
1540
1541
1542	dash = 0
1543
1544
1545	end if
1546	;***************************************************
1547	; OPTIONAL --font_size-- FONT SIZE OF STRINGS
1548	;
1549	; font size can be changed to fit to the current plot
1550	;
1551	; data type: float
1552	;
1553	; example: font_size = 0.02
1554	;
1555	; default: 0.02
1556	;***************************************************
1557	if(.not. isvar("font_size"))then
1558
1559
1560	font_size = 0.02
1561
1562
1563	end if
1564	;***************************************************
1565	; OPTIONAL --font_size_legend-- FONT SIZE OF LEGEND STRINGS
1566	;
1567	; font size of the legend can be changed to fit to the current plot
1568	;
1569	; data type: float
1570	;
1571	; example: font_size_legend = 0.2
1572	;
1573	; default: 0.2
1574	;***************************************************
1575	if(.not. isvar("font_size_legend"))then
1576
1577
1578	font_size_legend = 0.2
1579
1580
1581	end if
1582	;***************************************************
1583
1584	end if
1585
1586	;******************************************************
1587	; parameters used by SPECTRA
1588	;******************************************************
1589
1590	if (script .EQ. "spectra") then
1591
1592	;***************************************************
1593	; REQUIRED --file_1-- INPUT FILE
1594	;
1595	; input file produced by PALM; if there is more than one file for all timesteps of a job chain,
1596	; 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');
1597	; otherwise declare path and full file name
1598	;
1599	; data type: string
1600	;
1601	; example: file_1 = "/path_to_file/file_sp" if you use more than one cycle number
1602	; file_1 = "/path_to_file/file_sp.4.nc" if you use one file (e.g. with cycle number 4)
1603	;
1604	; default: "File in"
1605	;***************************************************
1606	if(.not. isvar("file_1"))then
1607
1608
1609	file_1 = "File in"
1610
1611
1612	end if
1613	;***************************************************
1614	; REQUIRED if more than one file of a job chain --start_f-- START OF CYCLIC NUMBER
1615	;
1616	; 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
1617	;
1618	; data type: integer
1619	;
1620	; example: start_f = 3 (if first file: example_sp.3.nc) and 0 (if first file: example_sp.nc)
1621	;
1622	; default: -1
1623	;***************************************************
1624	if (.not. isvar("start_f"))then
1625
1626
1627	start_f = -1
1628
1629
1630	end if
1631	;***************************************************
1632	; REQUIRED if more than one file of a job chain --end_f-- END OF CYCLIC NUMBER
1633	;
1634	; 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
1635	;
1636	; data type: integer
1637	;
1638	; example: end_f = 48 (if last file: example_sp.48.nc) and 0 (if last file: example_sp.nc)
1639	;
1640	; default: -1
1641	;***************************************************
1642	if (.not. isvar("end_f"))then
1643
1644
1645	end_f = -1
1646
1647
1648	end if
1649	;***************************************************
1650	; OPTIONAL --format_out-- OUTPUT FORMAT
1651	;
1652	; format of the output file; supported file formats: x11, pdf, ps, eps, epsi or ncgm
1653	;
1654	; data type: string
1655	;
1656	; example: format_out = "pdf"
1657	;
1658	; default: "x11"
1659	;***************************************************
1660	if (.not. isvar("format_out"))then
1661
1662
1663	format_out = "x11"
1664
1665
1666	end if
1667	;***************************************************
1668	; OPTIONAL --file_out-- OUTPUT FILE
1669	;
1670	; path and name of the output file
1671	;
1672	; data type: string
1673	;
1674	; example: file_out = "/home/path_to_output/test_sp"
1675	;
1676	; default: "~/test_sp"
1677	;***************************************************
1678	if (.not. isvar("file_out"))then
1679
1680
1681	file_out = "~/test_sp"
1682
1683
1684	end if
1685	;***************************************************
1686	; OPTIONAL --no_columns-- NUMBER OF COLUMNS
1687	;
1688	; number of plots in one row
1689	;
1690	; data type: integer
1691	;
1692	; example: no_columns = 1
1693	;
1694	; default: 1
1695	;***************************************************
1696	if (.not. isvar("no_columns"))then
1697
1698
1699	no_columns = 1
1700
1701
1702	end if
1703	;***************************************************
1704	; OPTIONAL --no_rows-- NUMBER OF ROWS
1705	;
1706	; number of plots in one column
1707	;
1708	; data type: integer
1709	;
1710	; example: no_rows = 2
1711	;
1712	; default: 2
1713	;***************************************************
1714	if (.not. isvar("no_rows"))then
1715
1716
1717	no_rows = 2
1718
1719
1720	end if
1721	;***************************************************
1722	; OPTIONAL --var-- OUTPUT VARIABLES
1723	;
1724	; names of output variables; inidcate them with one comma before and after each variable
1725	;
1726	; data type: string
1727	;
1728	; example: var = ",u,v,w,"
1729	;
1730	; default: "all" (all available variables are plotted)
1731	;***************************************************
1732	if (.not. isvar("var"))then
1733
1734
1735	var = "all"
1736
1737
1738	end if
1739	;***************************************************
1740	; OPTIONAL --height_level-- HEIGHT LEVELS
1741	;
1742	; defines if all levels or selected levels are output; a 1D array with up to 100 elements
1743	; from 0 to 99 must be indicated
1744	;
1745	; data type: integer array
1746	;
1747	; example: height_level=(/0,3,7/)
1748	;
1749	; default: (/-1/) (all levels are output)
1750	;***************************************************
1751	if (.not. isvar("height_level"))then
1752
1753
1754	height_level = (/-1/)
1755
1756
1757	end if
1758	;***************************************************
1759	; OPTIONAL --sort-- TYPE OF SORTING
1760	;
1761	; defines if all heights [="height"] or all timesteps [="time"] are merged in one plot
1762	;
1763	; data type: string
1764	;
1765	; example: sort = "height"
1766	;
1767	; default: "height"
1768	;***************************************************
1769	if (.not. isvar("sort"))then
1770
1771
1772	sort = "height"
1773
1774
1775	end if
1776	;***************************************************
1777	; OPTIONAL --start_time_step-- FIRST TIME STEP
1778	;
1779	; first time step of plot; there is no need to specify a precise time value from the input file,
1780	; 'start_time_step' is rounded to the next existent time step
1781	;
1782	; data type: double; dimension: [hour]
1783	;
1784	; example: start_time_step = 0.5
1785	;
1786	; default: -1.d (first available time step)
1787	;***************************************************
1788	if (.not. isvar("start_time_step"))then
1789
1790
1791	start_time_step = -1.d
1792
1793
1794	end if
1795	;***************************************************
1796	; OPTIONAL --end_time_step-- LAST TIME STEP
1797	;
1798	; last time step of plot; there is no need to specify a precise time value from the input file,
1799	; 'end_time_step' is rounded to the next existent time step
1800	;
1801	; data type: double; dimension: [hour]
1802	;
1803	; example: end_time_step = 2.0
1804	;
1805	; default: -1.d (last available time step)
1806	;***************************************************
1807	if (.not. isvar("end_time_step"))then
1808
1809
1810	end_time_step = -1.d
1811
1812
1813	end if
1814	;***************************************************
1815	; OPTIONAL --black-- BLACK OR COLOR LINES
1816	;
1817	; colored [=0] or black and white [=1] plots
1818	;
1819	; data type: integer
1820	;
1821	; example: black = 0
1822	;
1823	; default: 0
1824	;***************************************************
1825	if (.not. isvar("black"))then
1826
1827
1828	black = 0
1829
1830
1831	end if
1832	;***************************************************
1833	; OPTIONAL --dash-- DASHED LINES
1834	;
1835	; use of different line patterns [=1] or continous lines [=0]
1836	;
1837	; data type: integer
1838	;
1839	; example: dash = 0
1840	;
1841	; default: 0
1842	;***************************************************
1843	if (.not. isvar("dash"))then
1844
1845
1846	dash = 0
1847
1848
1849	end if
1850	;***************************************************
1851	; OPTIONAL --log_x-- LOGARITHMIC SCALE X-AXIS
1852	;
1853	; logarithmic scale for x-axis is switched on [=1] or off [=0]
1854	;
1855	; data type: integer
1856	;
1857	; example: log_x = 1
1858	;
1859	; default: 1
1860	;***************************************************
1861	if (.not. isvar("log_x"))then
1862
1863
1864	log_x = 1
1865
1866
1867	end if
1868	;***************************************************
1869	; OPTIONAL --log_y-- LOGARITHMIC SCALE Y-AXIS
1870	;
1871	; logarithmic scale for y-axis is switched on [=1] or off [=0]
1872	;
1873	; data type: integer
1874	;
1875	; example: log_y = 1
1876	;
1877	; default: 1
1878	;***************************************************
1879	if (.not. isvar("log_y"))then
1880
1881
1882	log_y = 1
1883
1884
1885	end if
1886	;***************************************************
1887	; OPTIONAL --norm_x-- NORMALISING X-AXIS
1888	;
1889	; value for normalising the x-axis
1890	;
1891	; data type: float
1892	;
1893	; example: norm_x = 1
1894	;
1895	; default: 1. (no normalising)
1896	;***************************************************
1897	if (.not. isvar("norm_x"))then
1898
1899
1900	norm_x = 1.
1901
1902
1903	end if
1904	;***************************************************
1905	; OPTIONAL --norm_height-- NORMALISING X-AXIS WITH HEIGHT
1906	;
1907	; normalising the x-axis with the height of each spectra is switched on [=1] or off [=0]
1908	;
1909	; data type: integer
1910	;
1911	; example: norm_height = 0
1912	;
1913	; default: 0 (no normalising)
1914	;***************************************************
1915	if (.not. isvar("norm_height"))then
1916
1917
1918	norm_height = 0
1919
1920
1921	end if
1922	;***************************************************
1923	; OPTIONAL --norm_y-- NORMALISE Y-AXIS
1924	;
1925	; value for normalising the y-axis
1926	;
1927	; data type: float
1928	;
1929	; example: norm_y = 1.
1930	;
1931	; default: 1. (no normalising)
1932	;***************************************************
1933	if (.not. isvar("norm_y"))then
1934
1935
1936	norm_y = 1.
1937
1938
1939	end if
1940	;***************************************************
1941	; OPTIONAL --unit_x--unit_y-- UNITS OF X-, Y-AXES
1942	;
1943	; 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
1944	;
1945	; data type: string,
1946	;
1947	; example: unit_x = "100/m"
1948	;
1949	; default: " "
1950	;***************************************************
1951	if(.not. isvar("unit_x"))then
1952
1953	unit_x = " "
1954
1955	end if
1956	if(.not. isvar("unit_y"))then
1957
1958	unit_y = " "
1959
1960	end if
1961	;***************************************************
1962	; OPTIONAL --font_size-- FONT SIZE OF STRINGS
1963	;
1964	; font size can be changed to fit to the current plot
1965	;
1966	; data type: float
1967	;
1968	; example: font_size = 0.02
1969	;
1970	; default: 0.02
1971	;***************************************************
1972	if(.not. isvar("font_size"))then
1973
1974
1975	font_size = 0.02
1976
1977
1978	end if
1979	;***************************************************
1980	; OPTIONAL --font_size_legend-- FONT SIZE OF LEGEND STRINGS
1981	;
1982	; font size of the legend can be changed to fit to the current plot
1983	;
1984	; data type: float
1985	;
1986	; example: font_size_legend = 0.2
1987	;
1988	; default: 0.2
1989	;***************************************************
1990	if(.not. isvar("font_size_legend"))then
1991
1992
1993	font_size_legend = 0.2
1994
1995
1996	end if
1997	;***************************************************
1998
1999	end if
2000
2001	;******************************************************
2002	; parameters used by TIMESERIES
2003	;******************************************************
2004
2005	if (script .EQ. "timeseries") then
2006
2007	;***************************************************
2008	; REQUIRED --file_1-- INPUT FILE
2009	;
2010	; input file produced by PALM; if there is more than one file for all timesteps of a job chain,
2011	; 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');
2012	; otherwise declare path and full file name
2013	;
2014	; data type: string
2015	;
2016	; example: file_1 = "/path_to_file/file_ts" if you use more than one cycle number
2017	; file_1 = "/path_to_file/file_ts.4.nc" if you use one file (e.g. with cycle number 4)
2018	;
2019	; default: "File in"
2020	;***************************************************
2021	if(.not. isvar("file_1"))then
2022
2023
2024	file_1 = "File in"
2025
2026
2027	end if
2028	;***************************************************
2029	; REQUIRED if more than one file of a job chain --start_f-- START OF CYCLIC NUMBER
2030	;
2031	; 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
2032	;
2033	; data type: integer
2034	;
2035	; example: start_f = 3 (if first file: example_ts.3.nc) and 0 (if first file: example_ts.nc)
2036	;
2037	; default: -1
2038	;***************************************************
2039	if (.not. isvar("start_f"))then
2040
2041
2042	start_f = -1
2043
2044
2045	end if
2046	;***************************************************
2047	; REQUIRED if more than one file of a job chain --end_f-- END OF CYCLIC NUMBER
2048	;
2049	; 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
2050	;
2051	; data type: integer
2052	;
2053	; example: end_f = 48 (if last file: example_ts.48.nc) and 0 (if last file: example_ts.nc)
2054	;
2055	; default: -1
2056	;***************************************************
2057	if (.not. isvar("end_f"))then
2058
2059
2060	end_f = -1
2061
2062
2063	end if
2064	;***************************************************
2065	; OPTIONAL --format_out-- OUTPUT FORMAT
2066	;
2067	; format of the output file; supported file formats: x11, pdf, ps, eps, epsi or ncgm
2068	;
2069	; data type: string
2070	;
2071	; example: format_out = "pdf"
2072	;
2073	; default: "x11"
2074	;***************************************************
2075	if (.not. isvar("format_out"))then
2076
2077
2078	format_out = "x11"
2079
2080
2081	end if
2082	;***************************************************
2083	; OPTIONAL --file_out-- OUTPUT FILE
2084	;
2085	; path and name of the output file
2086	;
2087	; data type: string
2088	;
2089	; example: file_out = "/home/path_to_output/test_ts"
2090	;
2091	; default: "~/test_ts"
2092	;***************************************************
2093	if (.not. isvar("file_out"))then
2094
2095
2096	file_out = "~/test_ts"
2097
2098
2099	end if
2100	;***************************************************
2101	; OPTIONAL --no_columns-- NUMBER OF COLUMNS
2102	;
2103	; number of plots in one row
2104	;
2105	; data type: integer
2106	;
2107	; example: no_columns = 1
2108	;
2109	; default: 1
2110	;***************************************************
2111	if (.not. isvar("no_columns"))then
2112
2113
2114	no_columns = 1
2115
2116
2117	end if
2118	;***************************************************
2119	; OPTIONAL --no_rows-- NUMBER OF ROWS
2120	;
2121	; number of plots in one column
2122	;
2123	; data type: integer
2124	;
2125	; example: no_rows = 2
2126	;
2127	; default: 2
2128	;***************************************************
2129	if (.not. isvar("no_rows"))then
2130
2131
2132	no_rows = 1
2133
2134
2135	end if
2136	;***************************************************
2137	; OPTIONAL --var-- OUTPUT VARIABLES
2138	;
2139	; names of output variables; inidcate them with one comma before and after each variable
2140	;
2141	; data type: string
2142	;
2143	; example: var = ",u,v,w,"
2144	;
2145	; default: "all" (all available variables are plotted)
2146	;***************************************************
2147	if (.not. isvar("var"))then
2148
2149
2150	var = "all"
2151
2152
2153	end if
2154	;***************************************************
2155	; OPTIONAL --start_time_step-- FIRST TIME STEP
2156	;
2157	; first time step of plot; there is no need to specify a precise time value from the input file,
2158	; 'start_time_step' is rounded to the next existent time step
2159	;
2160	; data type: double; dimension: [hour]
2161	;
2162	; example: start_time_step = 0.5
2163	;
2164	; default: -1.d (first available time step)
2165	;***************************************************
2166	if (.not. isvar("start_time_step"))then
2167
2168
2169	start_time_step = -1.d
2170
2171
2172	end if
2173	;***************************************************
2174	; OPTIONAL --end_time_step-- LAST TIME STEP
2175	;
2176	; last time step of plot; there is no need to specify a precise time value from the input file,
2177	; 'end_time_step' is rounded to the next existent time step
2178	;
2179	; data type: double; dimension: [hour]
2180	;
2181	; example: end_time_step = 2.0
2182	;
2183	; default: -1.d (last available time step)
2184	;***************************************************
2185	if (.not. isvar("end_time_step"))then
2186
2187
2188	end_time_step = -1.d
2189
2190
2191	end if
2192	;***************************************************
2193	; OPTIONAL --over-- OVERLAYING
2194	;
2195	; predefined overlaying of standard variables are switched on [=1] or off [=0];
2196	; this feature cannot be used with statistic regions
2197	;
2198	; data type: integer
2199	;
2200	; example: over = 0
2201	;
2202	; default: 0
2203	;***************************************************
2204	if (.not. isvar("over"))then
2205
2206
2207	over = 0
2208
2209
2210	end if
2211	;***************************************************
2212	; OPTIONAL --font_size-- FONT SIZE OF STRINGS
2213	;
2214	; font size can be changed to fit to the current plot
2215	;
2216	; data type: float
2217	;
2218	; example: font_size = 0.07
2219	;
2220	; default: 0.07
2221	;***************************************************
2222	if(.not. isvar("font_size"))then
2223
2224
2225	font_size = 0.07
2226
2227
2228	end if
2229	;***************************************************
2230	; OPTIONAL --norm_t-- NORMALISING X-AXIS OF TIMESERIES PLOTS (T-AXIS)
2231	;
2232	; normalising the x-axis with the given value; the unit can explicitly set; by default no normalising and unit is set to [s]
2233	;
2234	; data type: float
2235	;
2236	; example: norm_t = 3600
2237	;
2238	; default: 1.
2239	;***************************************************
2240	if(.not. isvar("normt"))then
2241
2242
2243	norm_t = 1.
2244
2245
2246	end if
2247	;***************************************************
2248	; OPTIONAL --unit_t-- UNIT OF X-AXIS OF TIMESERIES PLOTS (T-AXIS)
2249	;
2250	; unit of the x-axis can be set; by default unit is set to [s]
2251	;
2252	; data type: string
2253	;
2254	; example: unit_t = "h"
2255	;
2256	; default: " "
2257	;***************************************************
2258	if(.not. isvar("unitt"))then
2259
2260
2261	unit_t = " "
2262
2263
2264	end if
2265	;***************************************************
2266
2267	end if
2268
2269	;*********************************************************
2270
2271	end

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

Download in other formats:

Original Format