Home

Context Navigation

ncl_preferences.ncl @ 329

Last change on this file since 329 was 267, checked in by weinreis, 16 years ago
NCL scripts profiles cross_sections ncl_preference and UserGuide? modified
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