Home

Context Navigation

ncl_preferences.ncl @ 266

Last change on this file since 266 was 236, checked in by letzel, 16 years ago
NCL scripts in trunk/SCRIPTS/NCL updated
File size: 65.7 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	; REQUIRED --prof3d-- 3D OR PROFILE DATA
721	;
722	; must be switched on [=1] if 3D data are used or switched off [=0] if profile data are used for profiles
723	;
724	; data type: integer
725	;
726	; default: -1
727	;***************************************************
728	if (.not. isvar("prof3d"))then
729
730
731	prof3d = -1
732
733
734	end if
735	;***************************************************
736	; OPTIONAL --format_out-- OUTPUT FORMAT
737	;
738	; supported file formats: x11, pdf, ps, eps, epsi or ncgm
739	;
740	; data type: string
741	;
742	; default: "x11"
743	;***************************************************
744	if (.not. isvar("format_out"))then
745
746
747	format_out = "x11"
748
749
750	end if
751	;***************************************************
752	; OPTIONAL --file_out-- OUTPUT FILE
753	;
754	; declare path and name of the output file
755	;
756	; data type: string
757	;
758	; default: "~/test_pr"
759	;***************************************************
760	if (.not. isvar("file_out"))then
761
762
763	file_out = "~/test_pr"
764
765
766	end if
767	;***************************************************
768	; OPTIONAL --no_columns-- NUMBER OF COLUMNS
769	;
770	; number of plots in one row
771	;
772	; data type: integer
773	;
774	; default: 1
775	;***************************************************
776	if (.not. isvar("no_columns"))then
777
778
779	no_columns = 1
780
781
782	end if
783	;***************************************************
784	; OPTIONAL --no_rows-- NUMBER OF ROWS
785	;
786	; number of plots in one column
787	;
788	; data type: integer
789	;
790	; default: 2
791	;***************************************************
792	if (.not. isvar("no_rows"))then
793
794
795	no_rows = 2
796
797
798	end if
799	;***************************************************
800	; OPTIONAL --var-- OUTPUT VARIABLES
801	;
802	; names of output variables; inidcate them with one comma before and after each variable
803	;
804	; data type: string
805	;
806	; example: var = ",u,v,w";
807	;
808	; default: "all" (all available variables are plotted)
809	;***************************************************
810	if (.not. isvar("var"))then
811
812
813	var = "all"
814
815
816	end if
817	;***************************************************
818	; OPTIONAL --legend-- SWITCHES ON A LEGEND
819	;
820	; a legend is switched on [=1] or off [=0]
821	;
822	; data type: integer
823	;
824	; example: legend = 1
825	;
826	; default: 1
827	;***************************************************
828	if (.not. isvar("legend"))then
829
830
831	legend = 1
832
833
834	end if
835	;***************************************************
836	; OPTIONAL ---no_files-- NUMBER OF INPUT FILES
837	;
838	; up to 6 input files of different runs can be used; all input files must contain the same variables and dimensionalties;
839	; identical variables are merged in one plot;
840	; combining or overlaying of any variables is not feasible;
841	; an extra legend (besides the time step legend) is drawn where you can explicitily name the items
842	;
843	; data type: integer
844	;
845	; default: 1
846	;***************************************************
847	if (.not. isvar("no_files"))then
848
849
850	no_files = 1
851
852
853	end if
854	;***************************************************
855	; OPTIONAL --name_legend_1-- 1st LEGEND ITEM
856	;
857	; if no_files > 1: name of 1st item in extra legend
858	;
859	; data type: string
860	;
861	; default: "File 1"
862	;***************************************************
863	if (.not. isvar("name_legend_1"))then
864
865
866	name_legend_1 = "File 1"
867
868
869	end if
870	;***************************************************
871	; REQUIRED if no_files > 1 --file_2-- 2ND INPUT FILE
872	;
873	; 2nd input file produced by PALM; if there is more than one file for all timesteps of a job chain,
874	; 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');
875	; otherwise declare path and full file name
876	;
877	; data type: string
878	;
879	; example: file_2 = "/path_to_file/file_pr" if you use more than one cycle number
880	; file_2 = "/path_to_file/file_pr.4.nc" if you use one file (e.g. with cycle number 4)
881	;
882	; default: "File in"
883	;***************************************************
884	if(.not. isvar("file_2"))then
885
886
887	file_2 = "File in"
888
889
890	end if
891	;***************************************************
892	; REQUIRED --start_f_2-- START OF CYCLIC NUMBER OF THE 2ND INPUT FILE
893	;
894	; 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
895	;
896	; data type: integer
897	;
898	; example: 3 (if first file: example_ts.3.nc)
899	;
900	; default: -1
901	;***************************************************
902	if (.not. isvar("start_f_2"))then
903
904
905	start_f_2 = -1
906
907
908	end if
909	;***************************************************
910	; REQUIRED --end_f_2-- END OF CYCLIC NUMBER OF THE 2ND INPUT FILE
911	;
912	; 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
913	;
914	; data type: integer
915	;
916	; example: 48 (if last file: example_ts.48.nc)
917	;
918	; default: -1
919	;***************************************************
920	if (.not. isvar("end_f_2"))then
921
922
923	end_f_2 = -1
924
925
926	end if
927	;***************************************************
928	; OPTIONAL --name_legend_2-- 2nd LEGEND ITEM
929	;
930	; if no_files > 1: name of 2nd item in extra legend
931	;
932	; data type: string
933	;
934	; default: "File 2"
935	;***************************************************
936	if (.not. isvar("name_legend_2"))then
937
938
939	name_legend_2 = "File 2"
940
941
942	end if
943	;***************************************************
944	; REQUIRED if no_files > 2 --file_3-- 3RD INPUT FILE
945	;
946	; 3rd input file produced by PALM; if there is more than one file for all timesteps of a job chain,
947	; 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');
948	; otherwise declare path and full file name
949	;
950	; data type: string
951	;
952	; example: file_3 = "/path_to_file/file_pr" if you use more than one cycle number
953	; file_3 = "/path_to_file/file_pr.4.nc" if you use one file (e.g. with cycle number 4)
954	;
955	; default: "File in"
956	;***************************************************
957	if(.not. isvar("file_3"))then
958
959
960	file_3 = "File in"
961
962
963	end if
964	;***************************************************
965	; REQUIRED --start_f_3-- START OF CYCLIC NUMBER OF THE 3RD INPUT FILE
966	;
967	; 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
968	;
969	; data type: integer
970	;
971	; example: 3 (if first file: example_ts.3.nc)
972	;
973	; default: -1
974	;***************************************************
975	if (.not. isvar("start_f_3"))then
976
977
978	start_f_3 = -1
979
980
981	end if
982	;***************************************************
983	; REQUIRED --end_f_3-- END OF CYCLIC NUMBER OF THE 3RD INPUT FILE
984	;
985	; 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
986	;
987	; data type: integer
988	;
989	; example: 48 (if last file: example_ts.48.nc)
990	;
991	; default: -1
992	;***************************************************
993	if (.not. isvar("end_f_3"))then
994
995
996	end_f_3 = -1
997
998
999	end if
1000	;***************************************************
1001	; OPTIONAL --name_legend_3-- 3rd LEGEND ITEM
1002	;
1003	; if no_files > 2: name of 3rd item in extra legend
1004	;
1005	; data type: string
1006	;
1007	; default: "File 3"
1008	;***************************************************
1009	if (.not. isvar("name_legend_3"))then
1010
1011
1012	name_legend_3 = "File 3"
1013
1014
1015	end if
1016	;***************************************************
1017	; REQUIRED if no_files > 3 --file_4-- 4TH INPUT FILE
1018	;
1019	; 4th input file produced by PALM; if there is more than one file for all timesteps of a job chain,
1020	; 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');
1021	; otherwise declare path and full file name
1022	;
1023	; data type: string
1024	;
1025	; example: file_4 = "/path_to_file/file_pr" if you use more than one cycle number
1026	; file_4 = "/path_to_file/file_pr.4.nc" if you use one file (e.g. with cycle number 4)
1027	;
1028	; default: "File in"
1029	;***************************************************
1030	if(.not. isvar("file_4"))then
1031
1032
1033	file_4 = "File in"
1034
1035
1036	end if
1037	;***************************************************
1038	; REQUIRED --start_f_4-- START OF CYCLIC NUMBER OF THE 4TH INPUT FILE
1039	;
1040	; 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
1041	;
1042	; data type: integer
1043	;
1044	; example: 3 (if first file: example_ts.3.nc)
1045	;
1046	; default: -1
1047	;***************************************************
1048	if (.not. isvar("start_f_4"))then
1049
1050
1051	start_f_4 = -1
1052
1053
1054	end if
1055	;***************************************************
1056	; REQUIRED --end_f_4-- END OF CYCLIC NUMBER OF THE 4TH INPUT FILE
1057	;
1058	; 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
1059	;
1060	; data type: integer
1061	;
1062	; example: 48 (if last file: example_ts.48.nc)
1063	;
1064	; default: -1
1065	;***************************************************
1066	if (.not. isvar("end_f_4"))then
1067
1068
1069	end_f_4 = -1
1070
1071
1072	end if
1073	;***************************************************
1074	; OPTIONAL --name_legend_4-- 4th LEGEND ITEM
1075	;
1076	; if no_files > 3: name of 4th item in extra legend
1077	;
1078	; data type: string
1079	;
1080	; default: "File 4"
1081	;***************************************************
1082	if (.not. isvar("name_legend_4"))then
1083
1084
1085	name_legend_4 = "File 4"
1086
1087
1088	end if
1089	;***************************************************
1090	; REQUIRED if no_files > 4 --file_5-- 5TH INPUT FILE
1091	;
1092	; 5th input file produced by PALM; if there is more than one file for all timesteps of a job chain,
1093	; 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');
1094	; otherwise declare path and full file name
1095	;
1096	; data type: string
1097	;
1098	; example: file_5 = "/path_to_file/file_pr" if you use more than one cycle number
1099	; file_5 = "/path_to_file/file_pr.4.nc" if you use one file (e.g. with cycle number 4)
1100	;
1101	; default: "File in"
1102	;***************************************************
1103	if(.not. isvar("file_5"))then
1104
1105
1106	file_5 = "File in"
1107
1108
1109	end if
1110	;***************************************************
1111	; REQUIRED --start_f_5-- START OF CYCLIC NUMBER OF THE 5TH INPUT FILE
1112	;
1113	; 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
1114	;
1115	; data type: integer
1116	;
1117	; example: 3 (if first file: example_ts.3.nc)
1118	;
1119	; default: -1
1120	;***************************************************
1121	if (.not. isvar("start_f_5"))then
1122
1123
1124	start_f_5 = -1
1125
1126
1127	end if
1128	;***************************************************
1129	; REQUIRED --end_f_5-- END OF CYCLIC NUMBER OF THE 5TH INPUT FILE
1130	;
1131	; 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
1132	;
1133	; data type: integer
1134	;
1135	; example: 48 (if last file: example_ts.48.nc)
1136	;
1137	; default: -1
1138	;***************************************************
1139	if (.not. isvar("end_f_5"))then
1140
1141
1142	end_f_5 = -1
1143
1144
1145	end if
1146	;***************************************************
1147	; OPTIONAL --name_legend_5-- 5th LEGEND ITEM
1148	;
1149	; if no_files > 4: name of 5th item in extra legend
1150	;
1151	; data type: string
1152	;
1153	; default: "File 5"
1154	;***************************************************
1155	if (.not. isvar("name_legend_5"))then
1156
1157
1158	name_legend_5 = "File 5"
1159
1160
1161	end if
1162	;***************************************************
1163	; REQUIRED if no_files > 5 --file_6-- 6TH INPUT FILE
1164	;
1165	; 6th input file produced by PALM; if there is more than one file for all timesteps of a job chain,
1166	; 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');
1167	; otherwise declare path and full file name
1168	;
1169	; data type: string
1170	;
1171	; example: file_6 = "/path_to_file/file_pr" if you use more than one cycle number
1172	; file_6 = "/path_to_file/file_pr.4.nc" if you use one file (e.g. with cycle number 4)
1173	;
1174	; default: "File in"
1175	;***************************************************
1176	if(.not. isvar("file_6"))then
1177
1178
1179	file_6 = "File in"
1180
1181
1182	end if
1183	;***************************************************
1184	; REQUIRED --start_f_6-- START OF CYCLIC NUMBER OF THE 6TH INPUT FILE
1185	;
1186	; 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
1187	;
1188	; data type: integer
1189	;
1190	; example: 3 (if first file: example_ts.3.nc)
1191	;
1192	; default: -1
1193	;***************************************************
1194	if (.not. isvar("start_f_6"))then
1195
1196
1197	start_f_6 = -1
1198
1199
1200	end if
1201	;***************************************************
1202	; REQUIRED --end_f_6-- END OF CYCLIC NUMBER OF THE 6TH INPUT FILE
1203	;
1204	; 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
1205	;
1206	; data type: integer
1207	;
1208	; example: 48 (if last file: example_ts.48.nc)
1209	;
1210	; default: -1
1211	;***************************************************
1212	if (.not. isvar("end_f_6"))then
1213
1214
1215	end_f_6 = -1
1216
1217
1218	end if
1219	;***************************************************
1220	; OPTIONAL --name_legend_6-- 6TH LEGEND ITEM
1221	;
1222	; if no_files > 5: name of 6th item in extra legend
1223	;
1224	; data type: string
1225	;
1226	; default: "File 6"
1227	;***************************************************
1228	if (.not. isvar("name_legend_6"))then
1229
1230
1231	name_legend_6 = "File 6"
1232
1233
1234	end if
1235	;***************************************************
1236	; OPTIONAL --start_time_step-- FIRST TIME STEP
1237	;
1238	; first time step of plot; there is no need to specify a precise time value from the input file,
1239	; 'start_time_step' is rounded to the next existent time step
1240	;
1241	; data type: double; dimension: [hour]
1242	;
1243	; example: start_time_step = 0.5
1244	;
1245	; default: -1.d (first available time step)
1246	;***************************************************
1247	if (.not. isvar("start_time_step"))then
1248
1249
1250	start_time_step = -1.d
1251
1252
1253	end if
1254	;***************************************************
1255	; OPTIONAL --end_time_step-- LAST TIME STEP
1256	;
1257	; last time step of plot; there is no need to specify a precise time value from the input file,
1258	; 'end_time_step' is rounded to the next existent time step
1259	;
1260	; data type: double; dimension: [hour]
1261	;
1262	; example: end_time_step = 2.0
1263	;
1264	; default: -1.d (last available time step)
1265	;***************************************************
1266	if (.not. isvar("end_time_step"))then
1267
1268
1269	end_time_step = -1.d
1270
1271
1272	end if
1273	;***************************************************
1274	; OPTIONAL --start_x-- STARTVALUE IN X-DIRECTION
1275	;
1276	; if 3D data is used for profiles you can choose a startvalue in x-direction for horizontal averaging
1277	;
1278	; data type: integer; dimension [grippoint]
1279	;
1280	; example: start_x = 0
1281	;
1282	; default: 0 (minimum of range)
1283	;***************************************************
1284	if (.not. isvar("start_x"))then
1285
1286
1287	start_x = 0
1288
1289
1290	end if
1291	;***************************************************
1292	; OPTIONAL --end_x-- ENDVALUE IN X-DIRECTION
1293	;
1294	; if 3D data is used for profiles you can choose an endvalue in x-direction for horizontal averaging
1295	;
1296	; data type: integer; dimension [grippoint]
1297	;
1298	; example: end_x = 200
1299	;
1300	; default: -1 (maximum of range excluding ghostpoint)
1301	;***************************************************
1302	if (.not. isvar("end_x"))then
1303
1304
1305	end_x = -1
1306
1307
1308	end if
1309	;***************************************************
1310	; OPTIONAL --start_y-- STARTVALUE IN Y-DIRECTION
1311	;
1312	; if 3D data is used for profiles you can choose a startvalue in y-direction for horizontal averaging
1313	;
1314	; data type: integer; dimension [grippoint]
1315	;
1316	; example: start_y = 0
1317	;
1318	; default: 0 (minimum of range)
1319	;***************************************************
1320	if (.not. isvar("start_y"))then
1321
1322
1323	start_y = 0
1324
1325
1326	end if
1327	;***************************************************
1328	; OPTIONAL --end_y-- ENDVALUE IN Y-DIRECTION
1329	;
1330	; if 3D data is used for profiles you can choose an endvalue in y-direction for horizontal averaging
1331	;
1332	; data type: integer; dimension [grippoint]
1333	;
1334	; example: end_y = 200
1335	;
1336	; default: -1 (maximum of range excluding ghostpoint)
1337	;***************************************************
1338	if (.not. isvar("end_y"))then
1339
1340
1341	end_y = -1
1342
1343
1344	end if
1345	;***************************************************
1346	; OPTIONAL --xs-- MINIMUM X-AXIS
1347	;
1348	; value for minimum of x-axis
1349	;
1350	; data type: float; dimension of variable
1351	;
1352	; example: xs = 300. (e.g. if var=",pt")
1353	;
1354	; default: -1. (minimum of range)
1355	;***************************************************
1356	if (.not. isvar("xs"))then
1357
1358
1359	xs = -1.
1360
1361
1362	end if
1363	;***************************************************
1364	; OPTIONAL --xe-- MAXIMUM X-AXIS
1365	;
1366	; value for maximum of x-axis
1367	;
1368	; data type: float; dimension of variable
1369	;
1370	; example: xe = 302. (e.g. if var=",pt")
1371	;
1372	; default: -1. (maximum of range)
1373	;***************************************************
1374	if (.not. isvar("xe"))then
1375
1376
1377	xe = -1.
1378
1379
1380	end if
1381	;***************************************************
1382	; OPTIONAL --min_z-- MINIMUM Z-AXIS
1383	;
1384	; value for minimum height of z-axis
1385	;
1386	; data type: double; dimension: [meter]
1387	;
1388	; example: min_z = 0.0
1389	;
1390	; default: -1.d (related to the minimum of z-axis in input file)
1391	;***************************************************
1392	if (.not. isvar("min_z"))then
1393
1394
1395	min_z = -1.d
1396
1397
1398	end if
1399	;***************************************************
1400	; OPTIONAL --max_z-- MAXIMUM Z-AXIS
1401	;
1402	; value for maximum height of z-axis
1403	;
1404	; data type: double; dimension: [meter]
1405	;
1406	; example: max_z = 1000.0
1407	;
1408	; default: -1.d (related to the maximum of z-axis in input file)
1409	;***************************************************
1410	if (.not. isvar("max_z"))then
1411
1412
1413	max_z = -1.d
1414
1415
1416	end if
1417	;***************************************************
1418	; OPTIONAL --log_z-- LOGARITHMIC SCALE OF Z-AXIS
1419	;
1420	; logarithmic scale for z-axis is switched on [=1] or off [=0]
1421	;
1422	; data type: integer
1423	;
1424	; example: log_z = 0
1425	;
1426	; default: 0
1427	;***************************************************
1428	if (.not. isvar("log_z"))then
1429
1430
1431	log_z = 0
1432
1433
1434	end if
1435	;***************************************************
1436	; OPTIONAL --norm_z-- NORMALISING Z-AXIS
1437	;
1438	; value for normalising the z-axis
1439	;
1440	; data type: float
1441	;
1442	; example: norm_z = 1.
1443	;
1444	; default: 1. (no normalising)
1445	;***************************************************
1446	if (.not. isvar("norm_z"))then
1447
1448
1449	norm_z = 1.
1450
1451
1452	end if
1453	;***************************************************
1454	; OPTIONAL --over-- OVERLAYING
1455	;
1456	; predefined overlaying of standard variables is switched on [=1] or off [=0];
1457	; this feature cannot be used with statistic regions
1458	;
1459	; data type: integer
1460	;
1461	; example: over = 0
1462	;
1463	; default: 0
1464	;***************************************************
1465	if (.not. isvar("over"))then
1466
1467
1468	over = 0
1469
1470
1471	end if
1472	;***************************************************
1473	; OPTIONAL --combine-- COMBINING
1474	;
1475	; combining of two or three individual variables is possible and switched on [=1] or off[=0]
1476	;
1477	; data type: integer
1478	;
1479	; example: combine = 0
1480	;
1481	; default: 0
1482	;***************************************************
1483	if (.not. isvar("combine"))then
1484
1485
1486	combine = 0
1487
1488
1489	end if
1490	;***************************************************
1491	; REQUIRED IF combine=1 --number_comb-- NUMBER FOR COMBINING
1492	;
1493	; number of variables for combining; 2 or 3 is valid
1494	;
1495	; data type: integer
1496	;
1497	; example: number_comb = 2
1498	;
1499	; default: -1
1500	;***************************************************
1501	if (.not. isvar("number_comb"))then
1502
1503
1504	number_comb = -1
1505
1506
1507	end if
1508	;***************************************************
1509	; REQUIRED IF combine=1 --c_var-- VARIABLES FOR COMBINING
1510	;
1511	; names of variables; inidcate them with one comma before and after each variable;
1512	;
1513	; data type: string
1514	;
1515	; example: c_var = ",u,v,"
1516	;
1517	; default: "c_variables"
1518	;***************************************************
1519	if (.not. isvar("c_var"))then
1520
1521
1522	c_var = "c_variables"
1523
1524
1525	end if
1526	;***************************************************
1527	; OPTIONAL --black-- BLACK OR COLOR LINES
1528	;
1529	; colored [=0] or black and white [=1] plots
1530	;
1531	; data type: integer
1532	;
1533	; example: black = 0
1534	;
1535	; default: 0
1536	;***************************************************
1537	if (.not. isvar("black"))then
1538
1539
1540	black = 0
1541
1542
1543	end if
1544	;***************************************************
1545	; OPTIONAL --dash-- DASHED LINES
1546	;
1547	; use of different line patterns [=1] or continous lines [=0] for different time steps
1548	;
1549	; data type: integer
1550	;
1551	; example: dash = 0
1552	;
1553	; default: 0
1554	;***************************************************
1555	if (.not. isvar("dash"))then
1556
1557
1558	dash = 0
1559
1560
1561	end if
1562	;***************************************************
1563	; OPTIONAL --font_size-- FONT SIZE OF STRINGS
1564	;
1565	; font size can be changed to fit to the current plot
1566	;
1567	; data type: float
1568	;
1569	; example: font_size = 0.02
1570	;
1571	; default: 0.02
1572	;***************************************************
1573	if(.not. isvar("font_size"))then
1574
1575
1576	font_size = 0.02
1577
1578
1579	end if
1580	;***************************************************
1581	; OPTIONAL --font_size_legend-- FONT SIZE OF LEGEND STRINGS
1582	;
1583	; font size of the legend can be changed to fit to the current plot
1584	;
1585	; data type: float
1586	;
1587	; example: font_size_legend = 0.2
1588	;
1589	; default: 0.2
1590	;***************************************************
1591	if(.not. isvar("font_size_legend"))then
1592
1593
1594	font_size_legend = 0.2
1595
1596
1597	end if
1598	;***************************************************
1599
1600	end if
1601
1602	;******************************************************
1603	; parameters used by SPECTRA
1604	;******************************************************
1605
1606	if (script .EQ. "spectra") then
1607
1608	;***************************************************
1609	; REQUIRED --file_1-- INPUT FILE
1610	;
1611	; input file produced by PALM; if there is more than one file for all timesteps of a job chain,
1612	; 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');
1613	; otherwise declare path and full file name
1614	;
1615	; data type: string
1616	;
1617	; example: file_1 = "/path_to_file/file_sp" if you use more than one cycle number
1618	; file_1 = "/path_to_file/file_sp.4.nc" if you use one file (e.g. with cycle number 4)
1619	;
1620	; default: "File in"
1621	;***************************************************
1622	if(.not. isvar("file_1"))then
1623
1624
1625	file_1 = "File in"
1626
1627
1628	end if
1629	;***************************************************
1630	; REQUIRED if more than one file of a job chain --start_f-- START OF CYCLIC NUMBER
1631	;
1632	; 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
1633	;
1634	; data type: integer
1635	;
1636	; example: start_f = 3 (if first file: example_sp.3.nc) and 0 (if first file: example_sp.nc)
1637	;
1638	; default: -1
1639	;***************************************************
1640	if (.not. isvar("start_f"))then
1641
1642
1643	start_f = -1
1644
1645
1646	end if
1647	;***************************************************
1648	; REQUIRED if more than one file of a job chain --end_f-- END OF CYCLIC NUMBER
1649	;
1650	; 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
1651	;
1652	; data type: integer
1653	;
1654	; example: end_f = 48 (if last file: example_sp.48.nc) and 0 (if last file: example_sp.nc)
1655	;
1656	; default: -1
1657	;***************************************************
1658	if (.not. isvar("end_f"))then
1659
1660
1661	end_f = -1
1662
1663
1664	end if
1665	;***************************************************
1666	; OPTIONAL --format_out-- OUTPUT FORMAT
1667	;
1668	; format of the output file; supported file formats: x11, pdf, ps, eps, epsi or ncgm
1669	;
1670	; data type: string
1671	;
1672	; example: format_out = "pdf"
1673	;
1674	; default: "x11"
1675	;***************************************************
1676	if (.not. isvar("format_out"))then
1677
1678
1679	format_out = "x11"
1680
1681
1682	end if
1683	;***************************************************
1684	; OPTIONAL --file_out-- OUTPUT FILE
1685	;
1686	; path and name of the output file
1687	;
1688	; data type: string
1689	;
1690	; example: file_out = "/home/path_to_output/test_sp"
1691	;
1692	; default: "~/test_sp"
1693	;***************************************************
1694	if (.not. isvar("file_out"))then
1695
1696
1697	file_out = "~/test_sp"
1698
1699
1700	end if
1701	;***************************************************
1702	; OPTIONAL --no_columns-- NUMBER OF COLUMNS
1703	;
1704	; number of plots in one row
1705	;
1706	; data type: integer
1707	;
1708	; example: no_columns = 1
1709	;
1710	; default: 1
1711	;***************************************************
1712	if (.not. isvar("no_columns"))then
1713
1714
1715	no_columns = 1
1716
1717
1718	end if
1719	;***************************************************
1720	; OPTIONAL --no_rows-- NUMBER OF ROWS
1721	;
1722	; number of plots in one column
1723	;
1724	; data type: integer
1725	;
1726	; example: no_rows = 2
1727	;
1728	; default: 2
1729	;***************************************************
1730	if (.not. isvar("no_rows"))then
1731
1732
1733	no_rows = 2
1734
1735
1736	end if
1737	;***************************************************
1738	; OPTIONAL --var-- OUTPUT VARIABLES
1739	;
1740	; names of output variables; inidcate them with one comma before and after each variable
1741	;
1742	; data type: string
1743	;
1744	; example: var = ",u,v,w"
1745	;
1746	; default: "all" (all available variables are plotted)
1747	;***************************************************
1748	if (.not. isvar("var"))then
1749
1750
1751	var = "all"
1752
1753
1754	end if
1755	;***************************************************
1756	; OPTIONAL --height_level-- HEIGHT LEVELS
1757	;
1758	; defines if all levels or selected levels are output; a 1D array with up to 100 elements
1759	; from 0 to 99 must be indicated
1760	;
1761	; data type: integer array
1762	;
1763	; example: height_level=(/0,3,7/)
1764	;
1765	; default: (/-1/) (all levels are output)
1766	;***************************************************
1767	if (.not. isvar("height_level"))then
1768
1769
1770	height_level = (/-1/)
1771
1772
1773	end if
1774	;***************************************************
1775	; OPTIONAL --sort-- TYPE OF SORTING
1776	;
1777	; defines if all heights [="height"] or all timesteps [="time"] are merged in one plot
1778	;
1779	; data type: string
1780	;
1781	; example: sort = "height"
1782	;
1783	; default: "height"
1784	;***************************************************
1785	if (.not. isvar("sort"))then
1786
1787
1788	sort = "height"
1789
1790
1791	end if
1792	;***************************************************
1793	; OPTIONAL --start_time_step-- FIRST TIME STEP
1794	;
1795	; first time step of plot; there is no need to specify a precise time value from the input file,
1796	; 'start_time_step' is rounded to the next existent time step
1797	;
1798	; data type: double; dimension: [hour]
1799	;
1800	; example: start_time_step = 0.5
1801	;
1802	; default: -1.d (first available time step)
1803	;***************************************************
1804	if (.not. isvar("start_time_step"))then
1805
1806
1807	start_time_step = -1.d
1808
1809
1810	end if
1811	;***************************************************
1812	; OPTIONAL --end_time_step-- LAST TIME STEP
1813	;
1814	; last time step of plot; there is no need to specify a precise time value from the input file,
1815	; 'end_time_step' is rounded to the next existent time step
1816	;
1817	; data type: double; dimension: [hour]
1818	;
1819	; example: end_time_step = 2.0
1820	;
1821	; default: -1.d (last available time step)
1822	;***************************************************
1823	if (.not. isvar("end_time_step"))then
1824
1825
1826	end_time_step = -1.d
1827
1828
1829	end if
1830	;***************************************************
1831	; OPTIONAL --black-- BLACK OR COLOR LINES
1832	;
1833	; colored [=0] or black and white [=1] plots
1834	;
1835	; data type: integer
1836	;
1837	; example: black = 0
1838	;
1839	; default: 0
1840	;***************************************************
1841	if (.not. isvar("black"))then
1842
1843
1844	black = 0
1845
1846
1847	end if
1848	;***************************************************
1849	; OPTIONAL --dash-- DASHED LINES
1850	;
1851	; use of different line patterns [=1] or continous lines [=0]
1852	;
1853	; data type: integer
1854	;
1855	; example: dash = 0
1856	;
1857	; default: 0
1858	;***************************************************
1859	if (.not. isvar("dash"))then
1860
1861
1862	dash = 0
1863
1864
1865	end if
1866	;***************************************************
1867	; OPTIONAL --log_x-- LOGARITHMIC SCALE X-AXIS
1868	;
1869	; logarithmic scale for x-axis is switched on [=1] or off [=0]
1870	;
1871	; data type: integer
1872	;
1873	; example: log_x = 1
1874	;
1875	; default: 1
1876	;***************************************************
1877	if (.not. isvar("log_x"))then
1878
1879
1880	log_x = 1
1881
1882
1883	end if
1884	;***************************************************
1885	; OPTIONAL --log_y-- LOGARITHMIC SCALE Y-AXIS
1886	;
1887	; logarithmic scale for y-axis is switched on [=1] or off [=0]
1888	;
1889	; data type: integer
1890	;
1891	; example: log_y = 1
1892	;
1893	; default: 1
1894	;***************************************************
1895	if (.not. isvar("log_y"))then
1896
1897
1898	log_y = 1
1899
1900
1901	end if
1902	;***************************************************
1903	; OPTIONAL --norm_x-- NORMALISING X-AXIS
1904	;
1905	; value for normalising the x-axis
1906	;
1907	; data type: float
1908	;
1909	; example: norm_x = 1
1910	;
1911	; default: 1. (no normalising)
1912	;***************************************************
1913	if (.not. isvar("norm_x"))then
1914
1915
1916	norm_x = 1.
1917
1918
1919	end if
1920	;***************************************************
1921	; OPTIONAL --norm_height-- NORMALISING X-AXIS WITH HEIGHT
1922	;
1923	; normalising the x-axis with the height of each spectra is switched on [=1] or off [=0]
1924	;
1925	; data type: integer
1926	;
1927	; example: norm_height = 0
1928	;
1929	; default: 0 (no normalising)
1930	;***************************************************
1931	if (.not. isvar("norm_height"))then
1932
1933
1934	norm_height = 0
1935
1936
1937	end if
1938	;***************************************************
1939	; OPTIONAL --norm_y-- NORMALISE Y-AXIS
1940	;
1941	; value for normalising the y-axis
1942	;
1943	; data type: float
1944	;
1945	; example: norm_y = 1.
1946	;
1947	; default: 1. (no normalising)
1948	;***************************************************
1949	if (.not. isvar("norm_y"))then
1950
1951
1952	norm_y = 1.
1953
1954
1955	end if
1956	;***************************************************
1957	; OPTIONAL --unit_x--unit_y-- UNITS OF X-, Y-AXES
1958	;
1959	; 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
1960	;
1961	; data type: string,
1962	;
1963	; example: unit_x = "100/m"
1964	;
1965	; default: " "
1966	;***************************************************
1967	if(.not. isvar("unit_x"))then
1968
1969	unit_x = " "
1970
1971	end if
1972	if(.not. isvar("unit_y"))then
1973
1974	unit_y = " "
1975
1976	end if
1977	;***************************************************
1978	; OPTIONAL --font_size-- FONT SIZE OF STRINGS
1979	;
1980	; font size can be changed to fit to the current plot
1981	;
1982	; data type: float
1983	;
1984	; example: font_size = 0.02
1985	;
1986	; default: 0.02
1987	;***************************************************
1988	if(.not. isvar("font_size"))then
1989
1990
1991	font_size = 0.02
1992
1993
1994	end if
1995	;***************************************************
1996	; OPTIONAL --font_size_legend-- FONT SIZE OF LEGEND STRINGS
1997	;
1998	; font size of the legend can be changed to fit to the current plot
1999	;
2000	; data type: float
2001	;
2002	; example: font_size_legend = 0.2
2003	;
2004	; default: 0.2
2005	;***************************************************
2006	if(.not. isvar("font_size_legend"))then
2007
2008
2009	font_size_legend = 0.2
2010
2011
2012	end if
2013	;***************************************************
2014
2015	end if
2016
2017	;******************************************************
2018	; parameters used by TIMESERIES
2019	;******************************************************
2020
2021	if (script .EQ. "timeseries") then
2022
2023	;***************************************************
2024	; REQUIRED --file_1-- INPUT FILE
2025	;
2026	; input file produced by PALM; if there is more than one file for all timesteps of a job chain,
2027	; 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');
2028	; otherwise declare path and full file name
2029	;
2030	; data type: string
2031	;
2032	; example: file_1 = "/path_to_file/file_ts" if you use more than one cycle number
2033	; file_1 = "/path_to_file/file_ts.4.nc" if you use one file (e.g. with cycle number 4)
2034	;
2035	; default: "File in"
2036	;***************************************************
2037	if(.not. isvar("file_1"))then
2038
2039
2040	file_1 = "File in"
2041
2042
2043	end if
2044	;***************************************************
2045	; REQUIRED if more than one file of a job chain --start_f-- START OF CYCLIC NUMBER
2046	;
2047	; 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
2048	;
2049	; data type: integer
2050	;
2051	; example: start_f = 3 (if first file: example_ts.3.nc) and 0 (if first file: example_ts.nc)
2052	;
2053	; default: -1
2054	;***************************************************
2055	if (.not. isvar("start_f"))then
2056
2057
2058	start_f = -1
2059
2060
2061	end if
2062	;***************************************************
2063	; REQUIRED if more than one file of a job chain --end_f-- END OF CYCLIC NUMBER
2064	;
2065	; 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
2066	;
2067	; data type: integer
2068	;
2069	; example: end_f = 48 (if last file: example_ts.48.nc) and 0 (if last file: example_ts.nc)
2070	;
2071	; default: -1
2072	;***************************************************
2073	if (.not. isvar("end_f"))then
2074
2075
2076	end_f = -1
2077
2078
2079	end if
2080	;***************************************************
2081	; OPTIONAL --format_out-- OUTPUT FORMAT
2082	;
2083	; format of the output file; supported file formats: x11, pdf, ps, eps, epsi or ncgm
2084	;
2085	; data type: string
2086	;
2087	; example: format_out = "pdf"
2088	;
2089	; default: "x11"
2090	;***************************************************
2091	if (.not. isvar("format_out"))then
2092
2093
2094	format_out = "x11"
2095
2096
2097	end if
2098	;***************************************************
2099	; OPTIONAL --file_out-- OUTPUT FILE
2100	;
2101	; path and name of the output file
2102	;
2103	; data type: string
2104	;
2105	; example: file_out = "/home/path_to_output/test_ts"
2106	;
2107	; default: "~/test_ts"
2108	;***************************************************
2109	if (.not. isvar("file_out"))then
2110
2111
2112	file_out = "~/test_ts"
2113
2114
2115	end if
2116	;***************************************************
2117	; OPTIONAL --no_columns-- NUMBER OF COLUMNS
2118	;
2119	; number of plots in one row
2120	;
2121	; data type: integer
2122	;
2123	; example: no_columns = 1
2124	;
2125	; default: 1
2126	;***************************************************
2127	if (.not. isvar("no_columns"))then
2128
2129
2130	no_columns = 1
2131
2132
2133	end if
2134	;***************************************************
2135	; OPTIONAL --no_rows-- NUMBER OF ROWS
2136	;
2137	; number of plots in one column
2138	;
2139	; data type: integer
2140	;
2141	; example: no_rows = 2
2142	;
2143	; default: 2
2144	;***************************************************
2145	if (.not. isvar("no_rows"))then
2146
2147
2148	no_rows = 2
2149
2150
2151	end if
2152	;***************************************************
2153	; OPTIONAL --var-- OUTPUT VARIABLES
2154	;
2155	; names of output variables; inidcate them with one comma before and after each variable
2156	;
2157	; data type: string
2158	;
2159	; example: var = ",u,v,w"
2160	;
2161	; default: "all" (all available variables are plotted)
2162	;***************************************************
2163	if (.not. isvar("var"))then
2164
2165
2166	var = "all"
2167
2168
2169	end if
2170	;***************************************************
2171	; OPTIONAL --start_time_step-- FIRST TIME STEP
2172	;
2173	; first time step of plot; there is no need to specify a precise time value from the input file,
2174	; 'start_time_step' is rounded to the next existent time step
2175	;
2176	; data type: double; dimension: [hour]
2177	;
2178	; example: start_time_step = 0.5
2179	;
2180	; default: -1.d (first available time step)
2181	;***************************************************
2182	if (.not. isvar("start_time_step"))then
2183
2184
2185	start_time_step = -1.d
2186
2187
2188	end if
2189	;***************************************************
2190	; OPTIONAL --end_time_step-- LAST TIME STEP
2191	;
2192	; last time step of plot; there is no need to specify a precise time value from the input file,
2193	; 'end_time_step' is rounded to the next existent time step
2194	;
2195	; data type: double; dimension: [hour]
2196	;
2197	; example: end_time_step = 2.0
2198	;
2199	; default: -1.d (last available time step)
2200	;***************************************************
2201	if (.not. isvar("end_time_step"))then
2202
2203
2204	end_time_step = -1.d
2205
2206
2207	end if
2208	;***************************************************
2209	; OPTIONAL --over-- OVERLAYING
2210	;
2211	; predefined overlaying of standard variables are switched on [=1] or off [=0];
2212	; this feature cannot be used with statistic regions
2213	;
2214	; data type: integer
2215	;
2216	; example: over = 0
2217	;
2218	; default: 0
2219	;***************************************************
2220	if (.not. isvar("over"))then
2221
2222
2223	over = 0
2224
2225
2226	end if
2227	;***************************************************
2228	; OPTIONAL --font_size-- FONT SIZE OF STRINGS
2229	;
2230	; font size can be changed to fit to the current plot
2231	;
2232	; data type: float
2233	;
2234	; example: font_size = 0.07
2235	;
2236	; default: 0.07
2237	;***************************************************
2238	if(.not. isvar("font_size"))then
2239
2240
2241	font_size = 0.07
2242
2243
2244	end if
2245	;***************************************************
2246	; OPTIONAL --norm_t-- NORMALISING X-AXIS OF TIMESERIES PLOTS (T-AXIS)
2247	;
2248	; normalising the x-axis with the given value; the unit can explicitly set; by default no normalising and unit is set to [s]
2249	;
2250	; data type: float
2251	;
2252	; example: norm_t = 3600
2253	;
2254	; default: 1.
2255	;***************************************************
2256	if(.not. isvar("normt"))then
2257
2258
2259	norm_t = 1.
2260
2261
2262	end if
2263	;***************************************************
2264	; OPTIONAL --unit_t-- UNIT OF X-AXIS OF TIMESERIES PLOTS (T-AXIS)
2265	;
2266	; unit of the x-axis can be set; by default unit is set to [s]
2267	;
2268	; data type: string
2269	;
2270	; example: unit_t = "h"
2271	;
2272	; default: " "
2273	;***************************************************
2274	if(.not. isvar("unitt"))then
2275
2276
2277	unit_t = " "
2278
2279
2280	end if
2281	;***************************************************
2282
2283	end if
2284
2285	;*********************************************************
2286
2287	end

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

Download in other formats:

Original Format