Home

Context Navigation

ncl_preferences.ncl @ 346

Last change on this file since 346 was 342, checked in by heinze, 16 years ago
Plots of profiles appear in the order specified by the variable var.
File size: 65.2 KB

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

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

Download in other formats:

Original Format