Home

Context Navigation

ncl_preferences.ncl @ 218

Last change on this file since 218 was 218, checked in by letzel, 16 years ago
User can check user parameters and deduce further quantities in user_check_parameters Topography definition according to new user parameter topography_grid_convention (init_grid, modules, user_header, user_init_grid, user_parin) New subdirectory trunk/EXAMPLES with two examples highlighting the topography_grid_convention
File size: 65.0 KB

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

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

Download in other formats:

Original Format