Home

Context Navigation

ncl_preferences.ncl @ 218

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

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

Download in other formats:

Original Format