Home

Context Navigation

.ncl.config.default @ 3248

Last change on this file since 3248 was 958, checked in by hoffmann, 12 years ago
minor changes in palmplot pr
File size: 69.9 KB

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

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

Download in other formats:

Original Format