Home

Context Navigation

.ncl.config.default @ 532

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

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

Download in other formats:

Original Format