Home

Context Navigation

.ncl.config.default @ 951

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

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

Download in other formats:

Original Format