Home

Context Navigation

ncl_preferences.ncl @ 392

Last change on this file since 392 was 351, checked in by heinze, 15 years ago
New parameter time_stride is in profiles.ncl available.
File size: 65.9 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	;***************************************************
[351]	1258	; OPTIONAL --time_stride-- TEMPORAL STRIDE
	1259	;
	1260	; temporal stride for the plots; if time_stride=1 all the profiles between first_time_step
	1261	; and end_time_step are plotted (including first_time_step and end_time_step),
	1262	; if time_stride=3 every third time step between start_time_step and end_time_step is plotted
	1263	;
	1264	; data type: integer
	1265	;
	1266	; example: time_stride = 3
	1267	;
	1268	; default: 1
	1269	;***************************************************
	1270	if (.not. isvar("time_stride"))then
	1271
	1272
	1273	time_stride = 1
	1274
	1275
	1276	end if
	1277	;***************************************************
[218]	1278	; OPTIONAL --start_x-- STARTVALUE IN X-DIRECTION
[190]	1279	;
[218]	1280	; if 3D data is used for profiles you can choose a startvalue in x-direction for horizontal averaging
[190]	1281	;
[218]	1282	; data type: integer; dimension [grippoint]
[190]	1283	;
[218]	1284	; example: start_x = 0
	1285	;
	1286	; default: 0 (minimum of range)
[190]	1287	;***************************************************
[218]	1288	if (.not. isvar("start_x"))then
[190]	1289
[218]	1290
	1291	start_x = 0
[190]	1292
	1293
[218]	1294	end if
[190]	1295	;***************************************************
[218]	1296	; OPTIONAL --end_x-- ENDVALUE IN X-DIRECTION
[190]	1297	;
[218]	1298	; if 3D data is used for profiles you can choose an endvalue in x-direction for horizontal averaging
[190]	1299	;
[218]	1300	; data type: integer; dimension [grippoint]
[190]	1301	;
[218]	1302	; example: end_x = 200
	1303	;
	1304	; default: -1 (maximum of range excluding ghostpoint)
[190]	1305	;***************************************************
[218]	1306	if (.not. isvar("end_x"))then
[190]	1307
[218]	1308
	1309	end_x = -1
[190]	1310
	1311
[218]	1312	end if
[190]	1313	;***************************************************
[218]	1314	; OPTIONAL --start_y-- STARTVALUE IN Y-DIRECTION
[190]	1315	;
[218]	1316	; if 3D data is used for profiles you can choose a startvalue in y-direction for horizontal averaging
[190]	1317	;
[218]	1318	; data type: integer; dimension [grippoint]
[190]	1319	;
[218]	1320	; example: start_y = 0
	1321	;
	1322	; default: 0 (minimum of range)
[190]	1323	;***************************************************
[218]	1324	if (.not. isvar("start_y"))then
[190]	1325
[218]	1326
	1327	start_y = 0
[190]	1328
	1329
[218]	1330	end if
[190]	1331	;***************************************************
[218]	1332	; OPTIONAL --end_y-- ENDVALUE IN Y-DIRECTION
[190]	1333	;
[218]	1334	; if 3D data is used for profiles you can choose an endvalue in y-direction for horizontal averaging
[190]	1335	;
[218]	1336	; data type: integer; dimension [grippoint]
[190]	1337	;
[218]	1338	; example: end_y = 200
	1339	;
	1340	; default: -1 (maximum of range excluding ghostpoint)
[190]	1341	;***************************************************
[218]	1342	if (.not. isvar("end_y"))then
[190]	1343
[218]	1344
	1345	end_y = -1
[190]	1346
	1347
[218]	1348	end if
[190]	1349	;***************************************************
	1350	; OPTIONAL --xs-- MINIMUM X-AXIS
	1351	;
	1352	; value for minimum of x-axis
	1353	;
	1354	; data type: float; dimension of variable
	1355	;
[218]	1356	; example: xs = 300. (e.g. if var=",pt")
	1357	;
[190]	1358	; default: -1. (minimum of range)
	1359	;***************************************************
	1360	if (.not. isvar("xs"))then
	1361
	1362
	1363	xs = -1.
	1364
	1365
	1366	end if
	1367	;***************************************************
	1368	; OPTIONAL --xe-- MAXIMUM X-AXIS
	1369	;
	1370	; value for maximum of x-axis
	1371	;
	1372	; data type: float; dimension of variable
	1373	;
[218]	1374	; example: xe = 302. (e.g. if var=",pt")
	1375	;
[190]	1376	; default: -1. (maximum of range)
	1377	;***************************************************
	1378	if (.not. isvar("xe"))then
	1379
	1380
	1381	xe = -1.
	1382
	1383
	1384	end if
	1385	;***************************************************
	1386	; OPTIONAL --min_z-- MINIMUM Z-AXIS
	1387	;
	1388	; value for minimum height of z-axis
	1389	;
	1390	; data type: double; dimension: [meter]
	1391	;
[218]	1392	; example: min_z = 0.0
	1393	;
[190]	1394	; default: -1.d (related to the minimum of z-axis in input file)
	1395	;***************************************************
	1396	if (.not. isvar("min_z"))then
	1397
	1398
	1399	min_z = -1.d
	1400
	1401
	1402	end if
	1403	;***************************************************
	1404	; OPTIONAL --max_z-- MAXIMUM Z-AXIS
	1405	;
	1406	; value for maximum height of z-axis
	1407	;
	1408	; data type: double; dimension: [meter]
	1409	;
[218]	1410	; example: max_z = 1000.0
	1411	;
[190]	1412	; default: -1.d (related to the maximum of z-axis in input file)
	1413	;***************************************************
	1414	if (.not. isvar("max_z"))then
	1415
	1416
	1417	max_z = -1.d
	1418
	1419
	1420	end if
	1421	;***************************************************
[218]	1422	; OPTIONAL --log_z-- LOGARITHMIC SCALE OF Z-AXIS
[190]	1423	;
[218]	1424	; logarithmic scale for z-axis is switched on [=1] or off [=0]
[190]	1425	;
	1426	; data type: integer
	1427	;
[218]	1428	; example: log_z = 0
	1429	;
[190]	1430	; default: 0
	1431	;***************************************************
[218]	1432	if (.not. isvar("log_z"))then
[190]	1433
	1434
[218]	1435	log_z = 0
[190]	1436
	1437
	1438	end if
	1439	;***************************************************
[218]	1440	; OPTIONAL --norm_z-- NORMALISING Z-AXIS
[190]	1441	;
[218]	1442	; value for normalising the z-axis
[190]	1443	;
	1444	; data type: float
	1445	;
[218]	1446	; example: norm_z = 1.
	1447	;
	1448	; default: 1. (no normalising)
[190]	1449	;***************************************************
[218]	1450	if (.not. isvar("norm_z"))then
[190]	1451
	1452
[218]	1453	norm_z = 1.
[190]	1454
	1455
	1456	end if
	1457	;***************************************************
	1458	; OPTIONAL --over-- OVERLAYING
	1459	;
[218]	1460	; predefined overlaying of standard variables is switched on [=1] or off [=0];
[190]	1461	; this feature cannot be used with statistic regions
	1462	;
	1463	; data type: integer
	1464	;
[218]	1465	; example: over = 0
	1466	;
[190]	1467	; default: 0
	1468	;***************************************************
	1469	if (.not. isvar("over"))then
	1470
	1471
	1472	over = 0
	1473
	1474
	1475	end if
	1476	;***************************************************
	1477	; OPTIONAL --combine-- COMBINING
	1478	;
	1479	; combining of two or three individual variables is possible and switched on [=1] or off[=0]
	1480	;
	1481	; data type: integer
	1482	;
[218]	1483	; example: combine = 0
	1484	;
[190]	1485	; default: 0
	1486	;***************************************************
	1487	if (.not. isvar("combine"))then
	1488
	1489
	1490	combine = 0
	1491
	1492
	1493	end if
	1494	;***************************************************
	1495	; REQUIRED IF combine=1 --number_comb-- NUMBER FOR COMBINING
	1496	;
	1497	; number of variables for combining; 2 or 3 is valid
	1498	;
	1499	; data type: integer
	1500	;
[218]	1501	; example: number_comb = 2
	1502	;
[190]	1503	; default: -1
	1504	;***************************************************
	1505	if (.not. isvar("number_comb"))then
	1506
	1507
	1508	number_comb = -1
	1509
	1510
	1511	end if
	1512	;***************************************************
	1513	; REQUIRED IF combine=1 --c_var-- VARIABLES FOR COMBINING
	1514	;
[218]	1515	; names of variables; inidcate them with one comma before and after each variable;
[190]	1516	;
	1517	; data type: string
	1518	;
[218]	1519	; example: c_var = ",u,v,"
	1520	;
[190]	1521	; default: "c_variables"
	1522	;***************************************************
	1523	if (.not. isvar("c_var"))then
	1524
	1525
	1526	c_var = "c_variables"
	1527
	1528
	1529	end if
	1530	;***************************************************
	1531	; OPTIONAL --black-- BLACK OR COLOR LINES
	1532	;
	1533	; colored [=0] or black and white [=1] plots
	1534	;
	1535	; data type: integer
	1536	;
[218]	1537	; example: black = 0
	1538	;
[190]	1539	; default: 0
	1540	;***************************************************
	1541	if (.not. isvar("black"))then
	1542
	1543
	1544	black = 0
	1545
	1546
	1547	end if
	1548	;***************************************************
	1549	; OPTIONAL --dash-- DASHED LINES
	1550	;
	1551	; use of different line patterns [=1] or continous lines [=0] for different time steps
	1552	;
	1553	; data type: integer
	1554	;
[218]	1555	; example: dash = 0
	1556	;
[190]	1557	; default: 0
	1558	;***************************************************
	1559	if (.not. isvar("dash"))then
	1560
	1561
	1562	dash = 0
	1563
	1564
[218]	1565	end if
	1566	;***************************************************
	1567	; OPTIONAL --font_size-- FONT SIZE OF STRINGS
	1568	;
	1569	; font size can be changed to fit to the current plot
	1570	;
	1571	; data type: float
	1572	;
	1573	; example: font_size = 0.02
	1574	;
	1575	; default: 0.02
	1576	;***************************************************
	1577	if(.not. isvar("font_size"))then
	1578
	1579
	1580	font_size = 0.02
	1581
	1582
	1583	end if
	1584	;***************************************************
	1585	; OPTIONAL --font_size_legend-- FONT SIZE OF LEGEND STRINGS
	1586	;
	1587	; font size of the legend can be changed to fit to the current plot
	1588	;
	1589	; data type: float
	1590	;
	1591	; example: font_size_legend = 0.2
	1592	;
	1593	; default: 0.2
	1594	;***************************************************
	1595	if(.not. isvar("font_size_legend"))then
	1596
	1597
	1598	font_size_legend = 0.2
	1599
	1600
[190]	1601	end if
	1602	;***************************************************
	1603
	1604	end if
[194]	1605
[190]	1606	;******************************************************
	1607	; parameters used by SPECTRA
	1608	;******************************************************
	1609
[194]	1610	if (script .EQ. "spectra") then
	1611
[190]	1612	;***************************************************
	1613	; REQUIRED --file_1-- INPUT FILE
	1614	;
[218]	1615	; input file produced by PALM; if there is more than one file for all timesteps of a job chain,
	1616	; 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');
	1617	; otherwise declare path and full file name
[190]	1618	;
	1619	; data type: string
	1620	;
[218]	1621	; example: file_1 = "/path_to_file/file_sp" if you use more than one cycle number
	1622	; file_1 = "/path_to_file/file_sp.4.nc" if you use one file (e.g. with cycle number 4)
	1623	;
[190]	1624	; default: "File in"
	1625	;***************************************************
[218]	1626	if(.not. isvar("file_1"))then
	1627
	1628
	1629	file_1 = "File in"
	1630
	1631
	1632	end if
	1633	;***************************************************
	1634	; REQUIRED if more than one file of a job chain --start_f-- START OF CYCLIC NUMBER
	1635	;
	1636	; 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
	1637	;
	1638	; data type: integer
	1639	;
[236]	1640	; example: start_f = 3 (if first file: example_sp.3.nc) and 0 (if first file: example_sp.nc)
[218]	1641	;
	1642	; default: -1
	1643	;***************************************************
	1644	if (.not. isvar("start_f"))then
[190]	1645
	1646
[218]	1647	start_f = -1
[190]	1648
	1649
[218]	1650	end if
[190]	1651	;***************************************************
[218]	1652	; REQUIRED if more than one file of a job chain --end_f-- END OF CYCLIC NUMBER
	1653	;
	1654	; 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
	1655	;
	1656	; data type: integer
	1657	;
[236]	1658	; example: end_f = 48 (if last file: example_sp.48.nc) and 0 (if last file: example_sp.nc)
[218]	1659	;
	1660	; default: -1
	1661	;***************************************************
	1662	if (.not. isvar("end_f"))then
	1663
	1664
	1665	end_f = -1
	1666
	1667
	1668	end if
	1669	;***************************************************
[190]	1670	; OPTIONAL --format_out-- OUTPUT FORMAT
	1671	;
[218]	1672	; format of the output file; supported file formats: x11, pdf, ps, eps, epsi or ncgm
[190]	1673	;
	1674	; data type: string
	1675	;
[218]	1676	; example: format_out = "pdf"
	1677	;
[190]	1678	; default: "x11"
	1679	;***************************************************
	1680	if (.not. isvar("format_out"))then
	1681
	1682
	1683	format_out = "x11"
	1684
	1685
	1686	end if
	1687	;***************************************************
	1688	; OPTIONAL --file_out-- OUTPUT FILE
	1689	;
[218]	1690	; path and name of the output file
[190]	1691	;
	1692	; data type: string
	1693	;
[218]	1694	; example: file_out = "/home/path_to_output/test_sp"
	1695	;
[194]	1696	; default: "~/test_sp"
[190]	1697	;***************************************************
	1698	if (.not. isvar("file_out"))then
	1699
	1700
	1701	file_out = "~/test_sp"
	1702
	1703
	1704	end if
	1705	;***************************************************
[218]	1706	; OPTIONAL --no_columns-- NUMBER OF COLUMNS
[190]	1707	;
	1708	; number of plots in one row
	1709	;
	1710	; data type: integer
[218]	1711	;
	1712	; example: no_columns = 1
[190]	1713	;
	1714	; default: 1
	1715	;***************************************************
	1716	if (.not. isvar("no_columns"))then
	1717
	1718
	1719	no_columns = 1
	1720
	1721
	1722	end if
	1723	;***************************************************
[218]	1724	; OPTIONAL --no_rows-- NUMBER OF ROWS
[190]	1725	;
	1726	; number of plots in one column
	1727	;
	1728	; data type: integer
[218]	1729	;
	1730	; example: no_rows = 2
[190]	1731	;
	1732	; default: 2
	1733	;***************************************************
[218]	1734	if (.not. isvar("no_rows"))then
[190]	1735
	1736
[218]	1737	no_rows = 2
[190]	1738
	1739
	1740	end if
	1741	;***************************************************
	1742	; OPTIONAL --var-- OUTPUT VARIABLES
	1743	;
[218]	1744	; names of output variables; inidcate them with one comma before and after each variable
[190]	1745	;
	1746	; data type: string
[218]	1747	;
[342]	1748	; example: var = ",u,v,w,"
[190]	1749	;
[218]	1750	; default: "all" (all available variables are plotted)
[190]	1751	;***************************************************
	1752	if (.not. isvar("var"))then
	1753
	1754
	1755	var = "all"
	1756
	1757
	1758	end if
	1759	;***************************************************
	1760	; OPTIONAL --height_level-- HEIGHT LEVELS
	1761	;
[218]	1762	; defines if all levels or selected levels are output; a 1D array with up to 100 elements
	1763	; from 0 to 99 must be indicated
[190]	1764	;
	1765	; data type: integer array
	1766	;
[218]	1767	; example: height_level=(/0,3,7/)
	1768	;
	1769	; default: (/-1/) (all levels are output)
[190]	1770	;***************************************************
	1771	if (.not. isvar("height_level"))then
	1772
	1773
	1774	height_level = (/-1/)
	1775
	1776
	1777	end if
	1778	;***************************************************
	1779	; OPTIONAL --sort-- TYPE OF SORTING
	1780	;
[218]	1781	; defines if all heights [="height"] or all timesteps [="time"] are merged in one plot
[190]	1782	;
	1783	; data type: string
	1784	;
[218]	1785	; example: sort = "height"
	1786	;
[190]	1787	; default: "height"
	1788	;***************************************************
	1789	if (.not. isvar("sort"))then
	1790
	1791
	1792	sort = "height"
	1793
	1794
	1795	end if
	1796	;***************************************************
	1797	; OPTIONAL --start_time_step-- FIRST TIME STEP
	1798	;
[218]	1799	; first time step of plot; there is no need to specify a precise time value from the input file,
	1800	; 'start_time_step' is rounded to the next existent time step
[190]	1801	;
	1802	; data type: double; dimension: [hour]
	1803	;
[218]	1804	; example: start_time_step = 0.5
	1805	;
[190]	1806	; default: -1.d (first available time step)
	1807	;***************************************************
	1808	if (.not. isvar("start_time_step"))then
	1809
	1810
	1811	start_time_step = -1.d
	1812
	1813
	1814	end if
	1815	;***************************************************
	1816	; OPTIONAL --end_time_step-- LAST TIME STEP
	1817	;
[218]	1818	; last time step of plot; there is no need to specify a precise time value from the input file,
	1819	; 'end_time_step' is rounded to the next existent time step
[190]	1820	;
	1821	; data type: double; dimension: [hour]
	1822	;
[218]	1823	; example: end_time_step = 2.0
	1824	;
[190]	1825	; default: -1.d (last available time step)
	1826	;***************************************************
	1827	if (.not. isvar("end_time_step"))then
	1828
	1829
	1830	end_time_step = -1.d
	1831
	1832
	1833	end if
	1834	;***************************************************
	1835	; OPTIONAL --black-- BLACK OR COLOR LINES
	1836	;
	1837	; colored [=0] or black and white [=1] plots
	1838	;
	1839	; data type: integer
	1840	;
[218]	1841	; example: black = 0
	1842	;
[190]	1843	; default: 0
	1844	;***************************************************
	1845	if (.not. isvar("black"))then
	1846
	1847
	1848	black = 0
	1849
	1850
	1851	end if
	1852	;***************************************************
	1853	; OPTIONAL --dash-- DASHED LINES
	1854	;
	1855	; use of different line patterns [=1] or continous lines [=0]
	1856	;
	1857	; data type: integer
	1858	;
[218]	1859	; example: dash = 0
	1860	;
[190]	1861	; default: 0
	1862	;***************************************************
	1863	if (.not. isvar("dash"))then
	1864
	1865
	1866	dash = 0
	1867
	1868
	1869	end if
	1870	;***************************************************
[218]	1871	; OPTIONAL --log_x-- LOGARITHMIC SCALE X-AXIS
[190]	1872	;
[218]	1873	; logarithmic scale for x-axis is switched on [=1] or off [=0]
[190]	1874	;
	1875	; data type: integer
	1876	;
[218]	1877	; example: log_x = 1
	1878	;
[190]	1879	; default: 1
	1880	;***************************************************
[218]	1881	if (.not. isvar("log_x"))then
[190]	1882
	1883
[218]	1884	log_x = 1
[190]	1885
	1886
	1887	end if
	1888	;***************************************************
[218]	1889	; OPTIONAL --log_y-- LOGARITHMIC SCALE Y-AXIS
[190]	1890	;
[218]	1891	; logarithmic scale for y-axis is switched on [=1] or off [=0]
[190]	1892	;
	1893	; data type: integer
	1894	;
[218]	1895	; example: log_y = 1
	1896	;
[190]	1897	; default: 1
	1898	;***************************************************
[218]	1899	if (.not. isvar("log_y"))then
[190]	1900
	1901
[218]	1902	log_y = 1
[190]	1903
	1904
	1905	end if
	1906	;***************************************************
[218]	1907	; OPTIONAL --norm_x-- NORMALISING X-AXIS
[190]	1908	;
	1909	; value for normalising the x-axis
	1910	;
	1911	; data type: float
	1912	;
[218]	1913	; example: norm_x = 1
	1914	;
[190]	1915	; default: 1. (no normalising)
	1916	;***************************************************
[218]	1917	if (.not. isvar("norm_x"))then
[190]	1918
	1919
[218]	1920	norm_x = 1.
[190]	1921
	1922
[218]	1923	end if
	1924	;***************************************************
	1925	; OPTIONAL --norm_height-- NORMALISING X-AXIS WITH HEIGHT
	1926	;
	1927	; normalising the x-axis with the height of each spectra is switched on [=1] or off [=0]
	1928	;
	1929	; data type: integer
	1930	;
	1931	; example: norm_height = 0
	1932	;
	1933	; default: 0 (no normalising)
	1934	;***************************************************
	1935	if (.not. isvar("norm_height"))then
	1936
	1937
	1938	norm_height = 0
	1939
	1940
[190]	1941	end if
	1942	;***************************************************
[218]	1943	; OPTIONAL --norm_y-- NORMALISE Y-AXIS
[190]	1944	;
	1945	; value for normalising the y-axis
	1946	;
	1947	; data type: float
	1948	;
[218]	1949	; example: norm_y = 1.
	1950	;
[190]	1951	; default: 1. (no normalising)
	1952	;***************************************************
[218]	1953	if (.not. isvar("norm_y"))then
[190]	1954
	1955
[218]	1956	norm_y = 1.
[190]	1957
	1958
[218]	1959	end if
[190]	1960	;***************************************************
[218]	1961	; OPTIONAL --unit_x--unit_y-- UNITS OF X-, Y-AXES
	1962	;
	1963	; 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
	1964	;
	1965	; data type: string,
	1966	;
	1967	; example: unit_x = "100/m"
	1968	;
	1969	; default: " "
	1970	;***************************************************
	1971	if(.not. isvar("unit_x"))then
	1972
	1973	unit_x = " "
	1974
	1975	end if
	1976	if(.not. isvar("unit_y"))then
	1977
	1978	unit_y = " "
	1979
	1980	end if
	1981	;***************************************************
	1982	; OPTIONAL --font_size-- FONT SIZE OF STRINGS
	1983	;
	1984	; font size can be changed to fit to the current plot
	1985	;
	1986	; data type: float
	1987	;
	1988	; example: font_size = 0.02
	1989	;
	1990	; default: 0.02
	1991	;***************************************************
	1992	if(.not. isvar("font_size"))then
	1993
	1994
	1995	font_size = 0.02
	1996
	1997
	1998	end if
	1999	;***************************************************
	2000	; OPTIONAL --font_size_legend-- FONT SIZE OF LEGEND STRINGS
	2001	;
	2002	; font size of the legend can be changed to fit to the current plot
	2003	;
	2004	; data type: float
	2005	;
	2006	; example: font_size_legend = 0.2
	2007	;
	2008	; default: 0.2
	2009	;***************************************************
	2010	if(.not. isvar("font_size_legend"))then
	2011
	2012
	2013	font_size_legend = 0.2
	2014
	2015
	2016	end if
	2017	;***************************************************
[190]	2018
	2019	end if
[194]	2020
[190]	2021	;******************************************************
	2022	; parameters used by TIMESERIES
	2023	;******************************************************
	2024
[194]	2025	if (script .EQ. "timeseries") then
	2026
[190]	2027	;***************************************************
	2028	; REQUIRED --file_1-- INPUT FILE
	2029	;
[218]	2030	; input file produced by PALM; if there is more than one file for all timesteps of a job chain,
	2031	; 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');
	2032	; otherwise declare path and full file name
[190]	2033	;
	2034	; data type: string
	2035	;
[218]	2036	; example: file_1 = "/path_to_file/file_ts" if you use more than one cycle number
	2037	; file_1 = "/path_to_file/file_ts.4.nc" if you use one file (e.g. with cycle number 4)
	2038	;
[190]	2039	; default: "File in"
	2040	;***************************************************
[218]	2041	if(.not. isvar("file_1"))then
	2042
	2043
[267]	2044	file_1 = "File in"
[218]	2045
	2046
	2047	end if
	2048	;***************************************************
	2049	; REQUIRED if more than one file of a job chain --start_f-- START OF CYCLIC NUMBER
	2050	;
	2051	; 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
	2052	;
	2053	; data type: integer
	2054	;
[236]	2055	; example: start_f = 3 (if first file: example_ts.3.nc) and 0 (if first file: example_ts.nc)
[218]	2056	;
	2057	; default: -1
	2058	;***************************************************
	2059	if (.not. isvar("start_f"))then
[190]	2060
	2061
[218]	2062	start_f = -1
[190]	2063
	2064
[218]	2065	end if
[190]	2066	;***************************************************
[218]	2067	; REQUIRED if more than one file of a job chain --end_f-- END OF CYCLIC NUMBER
	2068	;
	2069	; 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
	2070	;
	2071	; data type: integer
	2072	;
[236]	2073	; example: end_f = 48 (if last file: example_ts.48.nc) and 0 (if last file: example_ts.nc)
[218]	2074	;
	2075	; default: -1
	2076	;***************************************************
	2077	if (.not. isvar("end_f"))then
	2078
	2079
	2080	end_f = -1
	2081
	2082
	2083	end if
	2084	;***************************************************
[190]	2085	; OPTIONAL --format_out-- OUTPUT FORMAT
	2086	;
[218]	2087	; format of the output file; supported file formats: x11, pdf, ps, eps, epsi or ncgm
[190]	2088	;
	2089	; data type: string
	2090	;
[218]	2091	; example: format_out = "pdf"
	2092	;
[190]	2093	; default: "x11"
	2094	;***************************************************
	2095	if (.not. isvar("format_out"))then
	2096
	2097
	2098	format_out = "x11"
	2099
	2100
	2101	end if
	2102	;***************************************************
	2103	; OPTIONAL --file_out-- OUTPUT FILE
	2104	;
[218]	2105	; path and name of the output file
[190]	2106	;
	2107	; data type: string
	2108	;
[218]	2109	; example: file_out = "/home/path_to_output/test_ts"
	2110	;
[194]	2111	; default: "~/test_ts"
[190]	2112	;***************************************************
	2113	if (.not. isvar("file_out"))then
	2114
	2115
	2116	file_out = "~/test_ts"
	2117
	2118
	2119	end if
	2120	;***************************************************
[218]	2121	; OPTIONAL --no_columns-- NUMBER OF COLUMNS
[190]	2122	;
	2123	; number of plots in one row
	2124	;
	2125	; data type: integer
[218]	2126	;
	2127	; example: no_columns = 1
[190]	2128	;
	2129	; default: 1
	2130	;***************************************************
	2131	if (.not. isvar("no_columns"))then
	2132
	2133
	2134	no_columns = 1
	2135
	2136
	2137	end if
	2138	;***************************************************
[218]	2139	; OPTIONAL --no_rows-- NUMBER OF ROWS
[190]	2140	;
	2141	; number of plots in one column
	2142	;
	2143	; data type: integer
[218]	2144	;
	2145	; example: no_rows = 2
[190]	2146	;
	2147	; default: 2
	2148	;***************************************************
[218]	2149	if (.not. isvar("no_rows"))then
[190]	2150
	2151
[267]	2152	no_rows = 1
[190]	2153
	2154
	2155	end if
	2156	;***************************************************
	2157	; OPTIONAL --var-- OUTPUT VARIABLES
	2158	;
[218]	2159	; names of output variables; inidcate them with one comma before and after each variable
[190]	2160	;
	2161	; data type: string
[218]	2162	;
[342]	2163	; example: var = ",u,v,w,"
[190]	2164	;
[218]	2165	; default: "all" (all available variables are plotted)
[190]	2166	;***************************************************
	2167	if (.not. isvar("var"))then
	2168
	2169
	2170	var = "all"
	2171
	2172
	2173	end if
	2174	;***************************************************
	2175	; OPTIONAL --start_time_step-- FIRST TIME STEP
	2176	;
[218]	2177	; first time step of plot; there is no need to specify a precise time value from the input file,
	2178	; 'start_time_step' is rounded to the next existent time step
[190]	2179	;
	2180	; data type: double; dimension: [hour]
	2181	;
[218]	2182	; example: start_time_step = 0.5
	2183	;
[190]	2184	; default: -1.d (first available time step)
	2185	;***************************************************
	2186	if (.not. isvar("start_time_step"))then
	2187
	2188
	2189	start_time_step = -1.d
	2190
	2191
	2192	end if
	2193	;***************************************************
	2194	; OPTIONAL --end_time_step-- LAST TIME STEP
	2195	;
[218]	2196	; last time step of plot; there is no need to specify a precise time value from the input file,
	2197	; 'end_time_step' is rounded to the next existent time step
[190]	2198	;
	2199	; data type: double; dimension: [hour]
	2200	;
[218]	2201	; example: end_time_step = 2.0
	2202	;
[190]	2203	; default: -1.d (last available time step)
	2204	;***************************************************
	2205	if (.not. isvar("end_time_step"))then
	2206
	2207
	2208	end_time_step = -1.d
	2209
[267]	2210
[190]	2211	end if
	2212	;***************************************************
	2213	; OPTIONAL --over-- OVERLAYING
	2214	;
[218]	2215	; predefined overlaying of standard variables are switched on [=1] or off [=0];
[190]	2216	; this feature cannot be used with statistic regions
	2217	;
	2218	; data type: integer
	2219	;
[218]	2220	; example: over = 0
	2221	;
[190]	2222	; default: 0
	2223	;***************************************************
	2224	if (.not. isvar("over"))then
	2225
	2226
	2227	over = 0
	2228
	2229
[218]	2230	end if
	2231	;***************************************************
	2232	; OPTIONAL --font_size-- FONT SIZE OF STRINGS
	2233	;
	2234	; font size can be changed to fit to the current plot
	2235	;
	2236	; data type: float
	2237	;
	2238	; example: font_size = 0.07
	2239	;
	2240	; default: 0.07
	2241	;***************************************************
	2242	if(.not. isvar("font_size"))then
	2243
	2244
	2245	font_size = 0.07
	2246
	2247
[190]	2248	end if
	2249	;***************************************************
[218]	2250	; OPTIONAL --norm_t-- NORMALISING X-AXIS OF TIMESERIES PLOTS (T-AXIS)
	2251	;
	2252	; normalising the x-axis with the given value; the unit can explicitly set; by default no normalising and unit is set to [s]
	2253	;
	2254	; data type: float
	2255	;
	2256	; example: norm_t = 3600
	2257	;
	2258	; default: 1.
	2259	;***************************************************
	2260	if(.not. isvar("normt"))then
	2261
	2262
	2263	norm_t = 1.
	2264
	2265
	2266	end if
	2267	;***************************************************
	2268	; OPTIONAL --unit_t-- UNIT OF X-AXIS OF TIMESERIES PLOTS (T-AXIS)
	2269	;
	2270	; unit of the x-axis can be set; by default unit is set to [s]
	2271	;
	2272	; data type: string
	2273	;
	2274	; example: unit_t = "h"
	2275	;
	2276	; default: " "
	2277	;***************************************************
	2278	if(.not. isvar("unitt"))then
	2279
	2280
	2281	unit_t = " "
	2282
	2283
	2284	end if
	2285	;***************************************************
[190]	2286
	2287	end if
	2288
	2289	;*********************************************************
	2290
	2291	end

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

Download in other formats:

Original Format