Home

Context Navigation

.ncl.config.default @ 515

Last change on this file since 515 was 513, checked in by heinze, 15 years ago
Using the NCL scripts by means of the shell script palmplot
File size: 69.4 KB

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

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

Download in other formats:

Original Format