Home

Context Navigation

netcdf.f90 @ 77

Last change on this file since 77 was 77, checked in by raasch, 18 years ago

New:
---

particle reflection from vertical walls implemented, particle SGS model adjusted to walls

Wall functions for vertical walls now include diabatic conditions. New subroutines wall_fluxes, wall_fluxes_e. New 4D-array rif_wall.

new d3par-parameter netcdf_64bit_3d to switch on 64bit offset only for 3D files

new d3par-parameter dt_max to define the maximum value for the allowed timestep

new inipar-parameter loop_optimization to control the loop optimization method

new inipar-parameter pt_refrence. If given, this value is used as the reference that in buoyancy terms (otherwise, the instantaneous horizontally averaged temperature is used).

new user interface user_advec_particles

new initializing action "by_user" calls user_init_3d_model and allows the initial setting of all 3d arrays

topography height informations are stored on arrays zu_s_inner and zw_w_inner and output to the 2d/3d NetCDF files

samples added to the user interface which show how to add user-define time series quantities.

calculation/output of precipitation amount, precipitation rate and z0 (by setting "pra*", "prr*", "z0*" with data_output). The time interval on which the precipitation amount is defined is set by new d3par-parameter precipitation_amount_interval

unit 9 opened for debug output (file DEBUG_<pe#>)

Makefile, advec_particles, average_3d_data, buoyancy, calc_precipitation, check_open, check_parameters, data_output_2d, diffusion_e, diffusion_u, diffusion_v, diffusion_w, diffusivities, header, impact_of_latent_heat, init_particles, init_3d_model, modules, netcdf, parin, production_e, read_var_list, read_3d_binary, sum_up_3d_data, user_interface, write_var_list, write_3d_binary

New: wall_fluxes

Changed:

General revision of non-cyclic horizontal boundary conditions: radiation boundary conditions are now used instead of Neumann conditions at the outflow (calculation needs velocity values for t-dt, which are stored on new arrays u_m_l, u_m_r, etc.), calculation of mean outflow is not needed any more, volume flow control is added for the outflow boundary (currently only for the north boundary!!), additional gridpoints along x and y (uxrp, vynp) are not needed any more, routine "boundary_conds" now operates on timelevel t+dt and is not split in two parts (main, uvw_outflow) any more, Neumann boundary conditions at inflow/outflow in case of non-cyclic boundary conditions for all 2d-arrays that are handled by exchange_horiz_2d

The FFT-method for solving the Poisson-equation is now working with Neumann boundary conditions both at the bottom and the top. This requires adjustments of the tridiagonal coefficients and subtracting the horizontally averaged mean from the vertical velocity field.

+age_m in particle_type

Particles-package is now part of the default code ("-p particles" is not needed any more).

Move call of user_actions( 'after_integration' ) below increment of times
and counters. user_actions is now called for each statistic region and has as an argument the number of the respective region (sr)

d3par-parameter data_output_ts removed. Timeseries output for "profil" removed. Timeseries are now switched on by dt_dots. Timeseries data is collected in flow_statistics.

Initial velocities at nzb+1 are regarded for volume flow control in case they have been set zero before (to avoid small timesteps); see new internal parameters u/v_nzb_p1_for_vfc.

q is not allowed to become negative (prognostic_equations).

poisfft_init is only called if fft-solver is switched on (init_pegrid).

d3par-parameter moisture renamed to humidity.

Subversion global revision number is read from mrun and added to the run description header and to the run control (_rc) file.

vtk directives removed from main program.

The uitility routine interpret_config reads PALM environment variables from NAMELIST instead using the system call GETENV.

advec_u_pw, advec_u_up, advec_v_pw, advec_v_up, asselin_filter, check_parameters, coriolis, data_output_dvrp, data_output_ptseries, data_output_ts, data_output_2d, data_output_3d, diffusion_u, diffusion_v, exchange_horiz, exchange_horiz_2d, flow_statistics, header, init_grid, init_particles, init_pegrid, init_rankine, init_pt_anomaly, init_1d_model, init_3d_model, modules, palm, package_parin, parin, poisfft, poismg, prandtl_fluxes, pres, production_e, prognostic_equations, read_var_list, read_3d_binary, sor, swap_timelevel, time_integration, write_var_list, write_3d_binary

Errors:

Bugfix: preset of tendencies te_em, te_um, te_vm in init_1d_model

Bugfix in sample for reading user defined data from restart file (user_init)

Bugfix in setting diffusivities for cases with the outflow damping layer extending over more than one subdomain (init_3d_model)

Check for possible negative humidities in the initial humidity profile.

in Makefile, default suffixes removed from the suffix list to avoid calling of m2c in
# case of .mod files

Makefile
check_parameters, init_1d_model, init_3d_model, user_interface

Property svn:keywords set to Id

File size: 142.5 KB

Rev	Line
[1]	1	#if defined( __ibmy_special )
	2	@PROCESS NOOPTimize
	3	#endif
	4	SUBROUTINE define_netcdf_header( callmode, extend, av )
	5
	6	!------------------------------------------------------------------------------!
	7	! Current revisions:
	8	! ------------------
[77]	9	!
[1]	10	!
	11	! Former revisions:
	12	! -----------------
[3]	13	! $Id: netcdf.f90 77 2007-03-29 04:26:56Z raasch $
[77]	14	!
	15	! 48 2007-03-06 12:28:36Z raasch
	16	! Output topography height information (zu_s_inner, zw_s_inner) to 2d-xy and 3d
	17	! datasets
	18	!
[3]	19	! RCS Log replace by Id keyword, revision history cleaned up
	20	!
[1]	21	! Revision 1.12 2006/09/26 19:35:16 raasch
	22	! Bugfix yv coordinates for yz cross sections
	23	!
	24	! Revision 1.1 2005/05/18 15:37:16 raasch
	25	! Initial revision
	26	!
	27	!
	28	! Description:
	29	! ------------
	30	! In case of extend = .FALSE.:
	31	! Define all necessary dimensions, axes and variables for the different
	32	! NetCDF datasets. This subroutine is called from check_open after a new
	33	! dataset is created. It leaves the open NetCDF files ready to write.
	34	!
	35	! In case of extend = .TRUE.:
	36	! Find out if dimensions and variables of an existing file match the values
	37	! of the actual run. If so, get all necessary informations (ids, etc.) from
	38	! this file.
	39	!
	40	! Parameter av can assume values 0 (non-averaged data) and 1 (time averaged
	41	! data)
	42	!------------------------------------------------------------------------------!
	43	#if defined( __netcdf )
	44
	45	USE arrays_3d
	46	USE constants
	47	USE control_parameters
	48	USE grid_variables
	49	USE indices
	50	USE netcdf_control
	51	USE pegrid
	52	USE particle_attributes
	53	USE profil_parameter
	54	USE spectrum
	55	USE statistics
	56
	57
	58	IMPLICIT NONE
	59
	60	CHARACTER (LEN=2) :: suffix
	61	CHARACTER (LEN=2), INTENT (IN) :: callmode
	62	CHARACTER (LEN=3) :: suffix1
	63	CHARACTER (LEN=4) :: grid_x, grid_y, grid_z
	64	CHARACTER (LEN=6) :: mode
	65	CHARACTER (LEN=10) :: netcdf_var_name, netcdf_var_name_base, &
	66	precision, var
	67	CHARACTER (LEN=80) :: time_average_text
	68	CHARACTER (LEN=2000) :: var_list, var_list_old
	69
	70	INTEGER :: av, i, id_x, id_y, id_z, j, ns, ns_old, nz_old
	71
	72	INTEGER, DIMENSION(1) :: id_dim_time_old, id_dim_x_yz_old, &
	73	id_dim_y_xz_old, id_dim_zu_sp_old, &
	74	id_dim_zu_xy_old, id_dim_zu_3d_old
	75
	76	LOGICAL :: found
	77
	78	LOGICAL, INTENT (INOUT) :: extend
	79
	80	LOGICAL, SAVE :: init_netcdf = .FALSE.
	81
	82	REAL, DIMENSION(1) :: last_time_coordinate
	83
	84	REAL, DIMENSION(:), ALLOCATABLE :: netcdf_data
	85
	86
	87	!
	88	!-- Initializing actions (return to calling routine check_parameters afterwards)
	89	IF ( .NOT. init_netcdf ) THEN
	90	!
	91	!-- Check and set accuracy for NetCDF output. First set default value
	92	nc_precision = NF90_REAL4
	93
	94	i = 1
	95	DO WHILE ( netcdf_precision(i) /= ' ' )
	96	j = INDEX( netcdf_precision(i), '_' )
	97	IF ( j == 0 ) THEN
	98	IF ( myid == 0 ) THEN
	99	PRINT*, '+++ define_netcdf_header: netcdf_precision must ', &
	100	'contain a "_" netcdf_precision(', i, ')="', &
	101	TRIM( netcdf_precision(i) ),'"'
	102	ENDIF
	103	CALL local_stop
	104	ENDIF
	105
	106	var = netcdf_precision(i)(1:j-1)
	107	precision = netcdf_precision(i)(j+1:)
	108
	109	IF ( precision == 'NF90_REAL4' ) THEN
	110	j = NF90_REAL4
	111	ELSEIF ( precision == 'NF90_REAL8' ) THEN
	112	j = NF90_REAL8
	113	ELSE
	114	IF ( myid == 0 ) THEN
	115	PRINT*, '+++ define_netcdf_header: illegal netcdf precision: ',&
	116	'netcdf_precision(', i, ')="', &
	117	TRIM( netcdf_precision(i) ),'"'
	118	ENDIF
	119	CALL local_stop
	120	ENDIF
	121
	122	SELECT CASE ( var )
	123	CASE ( 'xy' )
	124	nc_precision(1) = j
	125	CASE ( 'xz' )
	126	nc_precision(2) = j
	127	CASE ( 'yz' )
	128	nc_precision(3) = j
	129	CASE ( '2d' )
	130	nc_precision(1:3) = j
	131	CASE ( '3d' )
	132	nc_precision(4) = j
	133	CASE ( 'pr' )
	134	nc_precision(5) = j
	135	CASE ( 'ts' )
	136	nc_precision(6) = j
	137	CASE ( 'sp' )
	138	nc_precision(7) = j
	139	CASE ( 'prt' )
	140	nc_precision(8) = j
	141	CASE ( 'all' )
	142	nc_precision = j
	143
	144	CASE DEFAULT
	145	IF ( myid == 0 ) THEN
	146	PRINT*, '+++ define_netcdf_header: unknown variable in ', &
	147	'inipar assignment: netcdf_precision(', i, ')="',&
	148	TRIM( netcdf_precision(i) ),'"'
	149	ENDIF
	150	CALL local_stop
	151
	152	END SELECT
	153
	154	i = i + 1
	155	IF ( i > 10 ) EXIT
	156	ENDDO
	157
	158	init_netcdf = .TRUE.
	159
	160	RETURN
	161
	162	ENDIF
	163
	164	!
	165	!-- Determine the mode to be processed
	166	IF ( extend ) THEN
	167	mode = callmode // '_ext'
	168	ELSE
	169	mode = callmode // '_new'
	170	ENDIF
	171
	172	!
	173	!-- Select the mode to be processed. Possibilities are xy, xz, yz, pr and ts.
	174	SELECT CASE ( mode )
	175
	176	CASE ( '3d_new' )
	177
	178	!
	179	!-- Define some global attributes of the dataset
	180	nc_stat = NF90_PUT_ATT( id_set_3d(av), NF90_GLOBAL, 'Conventions', &
	181	'COARDS' )
	182	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 62 )
	183
	184	IF ( av == 0 ) THEN
	185	time_average_text = ' '
	186	ELSE
	187	WRITE (time_average_text, '('', '',F7.1,'' s average'')') &
	188	averaging_interval
	189	ENDIF
	190	nc_stat = NF90_PUT_ATT( id_set_3d(av), NF90_GLOBAL, 'title', &
	191	TRIM( run_description_header ) // &
	192	TRIM( time_average_text ) )
	193	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 63 )
	194	IF ( av == 1 ) THEN
	195	WRITE ( time_average_text,'(F7.1,'' s avg'')' ) averaging_interval
	196	nc_stat = NF90_PUT_ATT( id_set_3d(av), NF90_GLOBAL, 'time_avg', &
	197	TRIM( time_average_text ) )
	198	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 63 )
	199	ENDIF
	200
	201	!
	202	!-- Define time coordinate for volume data (unlimited dimension)
	203	nc_stat = NF90_DEF_DIM( id_set_3d(av), 'time', NF90_UNLIMITED, &
	204	id_dim_time_3d(av) )
	205	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 64 )
	206
	207	nc_stat = NF90_DEF_VAR( id_set_3d(av), 'time', NF90_DOUBLE, &
	208	id_dim_time_3d(av), id_var_time_3d(av) )
	209	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 65 )
	210
	211	nc_stat = NF90_PUT_ATT( id_set_3d(av), id_var_time_3d(av), 'units', &
	212	'seconds')
	213	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 66 )
	214
	215	!
	216	!-- Define spatial dimensions and coordinates:
	217	!-- Define vertical coordinate grid (zu grid)
	218	nc_stat = NF90_DEF_DIM( id_set_3d(av), 'zu_3d', nz_do3d-nzb+1, &
	219	id_dim_zu_3d(av) )
	220	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 67 )
	221
	222	nc_stat = NF90_DEF_VAR( id_set_3d(av), 'zu_3d', NF90_DOUBLE, &
	223	id_dim_zu_3d(av), id_var_zu_3d(av) )
	224	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 68 )
	225
	226	nc_stat = NF90_PUT_ATT( id_set_3d(av), id_var_zu_3d(av), 'units', &
	227	'meters' )
	228	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 69 )
	229
	230	!
	231	!-- Define vertical coordinate grid (zw grid)
	232	nc_stat = NF90_DEF_DIM( id_set_3d(av), 'zw_3d', nz_do3d-nzb+1, &
	233	id_dim_zw_3d(av) )
	234	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 70 )
	235
	236	nc_stat = NF90_DEF_VAR( id_set_3d(av), 'zw_3d', NF90_DOUBLE, &
	237	id_dim_zw_3d(av), id_var_zw_3d(av) )
	238	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 71 )
	239
	240	nc_stat = NF90_PUT_ATT( id_set_3d(av), id_var_zw_3d(av), 'units', &
	241	'meters' )
	242	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 72 )
	243
	244	!
	245	!-- Define x-axis (for scalar position)
	246	nc_stat = NF90_DEF_DIM( id_set_3d(av), 'x', nx+2, id_dim_x_3d(av) )
	247	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 73 )
	248
	249	nc_stat = NF90_DEF_VAR( id_set_3d(av), 'x', NF90_DOUBLE, &
	250	id_dim_x_3d(av), id_var_x_3d(av) )
	251	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 74 )
	252
	253	nc_stat = NF90_PUT_ATT( id_set_3d(av), id_var_x_3d(av), 'units', &
	254	'meters' )
	255	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 75 )
	256
	257	!
	258	!-- Define x-axis (for u position)
	259	nc_stat = NF90_DEF_DIM( id_set_3d(av), 'xu', nx+2, id_dim_xu_3d(av) )
	260	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 358 )
	261
	262	nc_stat = NF90_DEF_VAR( id_set_3d(av), 'xu', NF90_DOUBLE, &
	263	id_dim_xu_3d(av), id_var_xu_3d(av) )
	264	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 359 )
	265
	266	nc_stat = NF90_PUT_ATT( id_set_3d(av), id_var_xu_3d(av), 'units', &
	267	'meters' )
	268	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 360 )
	269
	270	!
	271	!-- Define y-axis (for scalar position)
	272	nc_stat = NF90_DEF_DIM( id_set_3d(av), 'y', ny+2, id_dim_y_3d(av) )
	273	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 76 )
	274
	275	nc_stat = NF90_DEF_VAR( id_set_3d(av), 'y', NF90_DOUBLE, &
	276	id_dim_y_3d(av), id_var_y_3d(av) )
	277	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 77 )
	278
	279	nc_stat = NF90_PUT_ATT( id_set_3d(av), id_var_y_3d(av), 'units', &
	280	'meters' )
	281	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 78 )
	282
	283	!
	284	!-- Define y-axis (for v position)
	285	nc_stat = NF90_DEF_DIM( id_set_3d(av), 'yv', ny+2, id_dim_yv_3d(av) )
	286	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 361 )
	287
	288	nc_stat = NF90_DEF_VAR( id_set_3d(av), 'yv', NF90_DOUBLE, &
	289	id_dim_yv_3d(av), id_var_yv_3d(av) )
	290	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 362 )
	291
	292	nc_stat = NF90_PUT_ATT( id_set_3d(av), id_var_yv_3d(av), 'units', &
	293	'meters' )
	294	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 363 )
	295
	296	!
[48]	297	!-- In case of non-flat topography define 2d-arrays containing the height
	298	!-- informations
	299	IF ( TRIM( topography ) /= 'flat' ) THEN
	300	!
	301	!-- Define zusi = zu(nzb_s_inner)
	302	nc_stat = NF90_DEF_VAR( id_set_3d(av), 'zusi', NF90_DOUBLE, &
	303	(/ id_dim_x_3d(av), id_dim_y_3d(av) /), &
	304	id_var_zusi_3d(av) )
	305	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 413 )
	306
	307	nc_stat = NF90_PUT_ATT( id_set_3d(av), id_var_zusi_3d(av), &
	308	'units', 'meters' )
	309	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 414 )
	310
	311	nc_stat = NF90_PUT_ATT( id_set_3d(av), id_var_zusi_3d(av), &
	312	'long_name', 'zu(nzb_s_inner)' )
	313	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 415 )
	314
	315	!
	316	!-- Define zwwi = zw(nzb_w_inner)
	317	nc_stat = NF90_DEF_VAR( id_set_3d(av), 'zwwi', NF90_DOUBLE, &
	318	(/ id_dim_x_3d(av), id_dim_y_3d(av) /), &
	319	id_var_zwwi_3d(av) )
	320	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 416 )
	321
	322	nc_stat = NF90_PUT_ATT( id_set_3d(av), id_var_zwwi_3d(av), &
	323	'units', 'meters' )
	324	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 417 )
	325
	326	nc_stat = NF90_PUT_ATT( id_set_3d(av), id_var_zwwi_3d(av), &
	327	'long_name', 'zw(nzb_w_inner)' )
	328	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 418 )
	329
	330	ENDIF
	331
	332
	333	!
[1]	334	!-- Define the variables
	335	var_list = ';'
	336	i = 1
	337
	338	DO WHILE ( do3d(av,i)(1:1) /= ' ' )
	339
	340	!
	341	!-- Check for the grid
	342	found = .TRUE.
	343	SELECT CASE ( do3d(av,i) )
	344	!
	345	!-- Most variables are defined on the scalar grid
	346	CASE ( 'e', 'p', 'pc', 'pr', 'pt', 'q', 'ql', 'ql_c', 'ql_v', &
	347	'ql_vp', 'qv', 's', 'vpt' )
	348
	349	grid_x = 'x'
	350	grid_y = 'y'
	351	grid_z = 'zu'
	352	!
	353	!-- u grid
	354	CASE ( 'u' )
	355
	356	grid_x = 'xu'
	357	grid_y = 'y'
	358	grid_z = 'zu'
	359	!
	360	!-- v grid
	361	CASE ( 'v' )
	362
	363	grid_x = 'x'
	364	grid_y = 'yv'
	365	grid_z = 'zu'
	366	!
	367	!-- w grid
	368	CASE ( 'w' )
	369
	370	grid_x = 'x'
	371	grid_y = 'y'
	372	grid_z = 'zw'
	373
	374	CASE DEFAULT
	375	!
	376	!-- Check for user-defined quantities
	377	CALL user_define_netcdf_grid( do3d(av,i), found, grid_x, &
	378	grid_y, grid_z )
	379
	380	END SELECT
	381
	382	!
	383	!-- Select the respective dimension ids
	384	IF ( grid_x == 'x' ) THEN
	385	id_x = id_dim_x_3d(av)
	386	ELSEIF ( grid_x == 'xu' ) THEN
	387	id_x = id_dim_xu_3d(av)
	388	ENDIF
	389
	390	IF ( grid_y == 'y' ) THEN
	391	id_y = id_dim_y_3d(av)
	392	ELSEIF ( grid_y == 'yv' ) THEN
	393	id_y = id_dim_yv_3d(av)
	394	ENDIF
	395
	396	IF ( grid_z == 'zu' ) THEN
	397	id_z = id_dim_zu_3d(av)
	398	ELSEIF ( grid_z == 'zw' ) THEN
	399	id_z = id_dim_zw_3d(av)
	400	ENDIF
	401
	402	!
	403	!-- Define the grid
	404	nc_stat = NF90_DEF_VAR( id_set_3d(av), do3d(av,i), &
	405	nc_precision(4), &
	406	(/ id_x, id_y, id_z, id_dim_time_3d(av) /), &
	407	id_var_do3d(av,i) )
	408
	409	IF ( .NOT. found ) THEN
	410	PRINT*, '+++ define_netcdf_header: no grid defined for', &
	411	' variable ', do3d(av,i)
	412	ENDIF
	413
	414	var_list = TRIM( var_list ) // TRIM( do3d(av,i) ) // ';'
	415
	416	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 79 )
	417	!
	418	!-- Store the 'real' name of the variable (with *, for example)
	419	!-- in the long_name attribute. This is evaluated by Ferret,
	420	!-- for example.
	421	nc_stat = NF90_PUT_ATT( id_set_3d(av), id_var_do3d(av,i), &
	422	'long_name', do3d(av,i) )
	423	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 80 )
	424	!
	425	!-- Define the variable's unit
	426	nc_stat = NF90_PUT_ATT( id_set_3d(av), id_var_do3d(av,i), &
	427	'units', TRIM( do3d_unit(av,i) ) )
	428	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 357 )
	429
	430	i = i + 1
	431
	432	ENDDO
	433
	434	!
	435	!-- No arrays to output
	436	IF ( i == 1 ) RETURN
	437
	438	!
	439	!-- Write the list of variables as global attribute (this is used by
	440	!-- restart runs and by combine_plot_fields)
	441	nc_stat = NF90_PUT_ATT( id_set_3d(av), NF90_GLOBAL, 'VAR_LIST', &
	442	var_list )
	443	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 81 )
	444
	445	!
	446	!-- Leave NetCDF define mode
	447	nc_stat = NF90_ENDDEF( id_set_3d(av) )
	448	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 82 )
	449
	450	!
	451	!-- Write data for x and xu axis (shifted by -dx/2)
	452	ALLOCATE( netcdf_data(0:nx+1) )
	453
	454	DO i = 0, nx+1
	455	netcdf_data(i) = i * dx
	456	ENDDO
	457
	458	nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_x_3d(av), netcdf_data, &
	459	start = (/ 1 /), count = (/ nx+2 /) )
	460	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 83 )
	461
	462	DO i = 0, nx+1
	463	netcdf_data(i) = ( i - 0.5 ) * dx
	464	ENDDO
	465
	466	nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_xu_3d(av), &
	467	netcdf_data, start = (/ 1 /), &
	468	count = (/ nx+2 /) )
	469	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 385 )
	470
	471	DEALLOCATE( netcdf_data )
	472
	473	!
	474	!-- Write data for y and yv axis (shifted by -dy/2)
	475	ALLOCATE( netcdf_data(0:ny+1) )
	476
	477	DO i = 0, ny+1
	478	netcdf_data(i) = i * dy
	479	ENDDO
	480
	481	nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_y_3d(av), netcdf_data, &
	482	start = (/ 1 /), count = (/ ny+2 /))
	483	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 84 )
	484
	485	DO i = 0, ny+1
	486	netcdf_data(i) = ( i - 0.5 ) * dy
	487	ENDDO
	488
	489	nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_yv_3d(av), &
	490	netcdf_data, start = (/ 1 /), &
	491	count = (/ ny+2 /))
	492	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 387 )
	493
	494	DEALLOCATE( netcdf_data )
	495
	496	!
	497	!-- Write zu and zw data (vertical axes)
	498	nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zu_3d(av), &
	499	zu(nzb:nz_do3d), start = (/ 1 /), &
	500	count = (/ nz_do3d-nzb+1 /) )
	501	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 85 )
	502
	503	nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zw_3d(av), &
	504	zw(nzb:nz_do3d), start = (/ 1 /), &
	505	count = (/ nz_do3d-nzb+1 /) )
	506	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 86 )
	507
	508
[48]	509	!
	510	!-- In case of non-flat topography write height information
	511	IF ( TRIM( topography ) /= 'flat' ) THEN
	512
	513	nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zusi_3d(av), &
	514	zu_s_inner(0:nx+1,0:ny+1), &
	515	start = (/ 1, 1 /), &
	516	count = (/ nx+2, ny+2 /) )
	517	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 419 )
	518
	519	nc_stat = NF90_PUT_VAR( id_set_3d(av), id_var_zwwi_3d(av), &
	520	zw_w_inner(0:nx+1,0:ny+1), &
	521	start = (/ 1, 1 /), &
	522	count = (/ nx+2, ny+2 /) )
	523	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 420 )
	524
	525	ENDIF
	526
[1]	527	CASE ( '3d_ext' )
	528
	529	!
	530	!-- Get the list of variables and compare with the actual run.
	531	!-- First var_list_old has to be reset, since GET_ATT does not assign
	532	!-- trailing blanks.
	533	var_list_old = ' '
	534	nc_stat = NF90_GET_ATT( id_set_3d(av), NF90_GLOBAL, 'VAR_LIST', &
	535	var_list_old )
	536	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 87 )
	537
	538	var_list = ';'
	539	i = 1
	540	DO WHILE ( do3d(av,i)(1:1) /= ' ' )
	541	var_list = TRIM(var_list) // TRIM( do3d(av,i) ) // ';'
	542	i = i + 1
	543	ENDDO
	544
	545	IF ( av == 0 ) THEN
	546	var = '(3d)'
	547	ELSE
	548	var = '(3d_av)'
	549	ENDIF
	550
	551	IF ( TRIM( var_list ) /= TRIM( var_list_old ) ) THEN
	552	PRINT*, '+++ WARNING: NetCDF file for volume data ' // &
	553	TRIM( var ) // ' from previuos run found,'
	554	PRINT*, ' but this file cannot be extended due to' // &
	555	' variable mismatch.'
	556	PRINT*, ' New file is created instead.'
	557	extend = .FALSE.
	558	RETURN
	559	ENDIF
	560
	561	!
	562	!-- Get and compare the number of vertical gridpoints
	563	nc_stat = NF90_INQ_VARID( id_set_3d(av), 'zu_3d', id_var_zu_3d(av) )
	564	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 88 )
	565
	566	nc_stat = NF90_INQUIRE_VARIABLE( id_set_3d(av), id_var_zu_3d(av), &
	567	dimids = id_dim_zu_3d_old )
	568	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 89 )
	569	id_dim_zu_3d(av) = id_dim_zu_3d_old(1)
	570
	571	nc_stat = NF90_INQUIRE_DIMENSION( id_set_3d(av), id_dim_zu_3d(av), &
	572	len = nz_old )
	573	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 90 )
	574
	575	IF ( nz_do3d-nzb+1 /= nz_old ) THEN
	576	PRINT*, '+++ WARNING: NetCDF file for volume data ' // &
	577	TRIM( var ) // ' from previuos run found,'
	578	PRINT*, ' but this file cannot be extended due to' // &
	579	' mismatch in number of'
	580	PRINT*, ' vertical grid points (nz_do3d).'
	581	PRINT*, ' New file is created instead.'
	582	extend = .FALSE.
	583	RETURN
	584	ENDIF
	585
	586	!
	587	!-- Get the id of the time coordinate (unlimited coordinate) and its
	588	!-- last index on the file. The next time level is pl3d..count+1.
	589	!-- The current time must be larger than the last output time
	590	!-- on the file.
	591	nc_stat = NF90_INQ_VARID( id_set_3d(av), 'time', id_var_time_3d(av) )
	592	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 91 )
	593
	594	nc_stat = NF90_INQUIRE_VARIABLE( id_set_3d(av), id_var_time_3d(av), &
	595	dimids = id_dim_time_old )
	596	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 92 )
	597	id_dim_time_3d(av) = id_dim_time_old(1)
	598
	599	nc_stat = NF90_INQUIRE_DIMENSION( id_set_3d(av), id_dim_time_3d(av), &
	600	len = do3d_time_count(av) )
	601	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 93 )
	602
	603	nc_stat = NF90_GET_VAR( id_set_3d(av), id_var_time_3d(av), &
	604	last_time_coordinate, &
	605	start = (/ do3d_time_count(av) /), &
	606	count = (/ 1 /) )
	607	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 94 )
	608
	609	IF ( last_time_coordinate(1) >= simulated_time ) THEN
	610	PRINT*, '+++ WARNING: NetCDF file for volume data ' // &
	611	TRIM( var ) // ' from previuos run found,'
	612	PRINT*, ' but this file cannot be extended becaus' // &
	613	'e the current output time'
	614	PRINT*, ' is less or equal than the last output t' // &
	615	'ime on this file.'
	616	PRINT*, ' New file is created instead.'
	617	do3d_time_count(av) = 0
	618	extend = .FALSE.
	619	RETURN
	620	ENDIF
	621
	622	!
	623	!-- Dataset seems to be extendable.
	624	!-- Now get the variable ids.
	625	i = 1
	626	DO WHILE ( do3d(av,i)(1:1) /= ' ' )
	627	nc_stat = NF90_INQ_VARID( id_set_3d(av), TRIM( do3d(av,i) ), &
	628	id_var_do3d(av,i) )
	629	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 95 )
	630	i = i + 1
	631	ENDDO
	632
	633	!
	634	!-- Change the titel attribute on file
	635	nc_stat = NF90_PUT_ATT( id_set_3d(av), NF90_GLOBAL, 'title', &
	636	TRIM( run_description_header ) )
	637	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 96 )
	638	PRINT, '** NetCDF file for volume data ' // TRIM( var ) // &
	639	' from previous run found.'
	640	PRINT*, ' This file will be extended.'
	641
	642
	643	CASE ( 'xy_new' )
	644
	645	!
	646	!-- Define some global attributes of the dataset
	647	nc_stat = NF90_PUT_ATT( id_set_xy(av), NF90_GLOBAL, 'Conventions', &
	648	'COARDS' )
	649	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 97 )
	650
	651	IF ( av == 0 ) THEN
	652	time_average_text = ' '
	653	ELSE
	654	WRITE (time_average_text, '('', '',F7.1,'' s average'')') &
	655	averaging_interval
	656	ENDIF
	657	nc_stat = NF90_PUT_ATT( id_set_xy(av), NF90_GLOBAL, 'title', &
	658	TRIM( run_description_header ) // &
	659	TRIM( time_average_text ) )
	660	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 98 )
	661	IF ( av == 1 ) THEN
	662	WRITE ( time_average_text,'(F7.1,'' s avg'')' ) averaging_interval
	663	nc_stat = NF90_PUT_ATT( id_set_xy(av), NF90_GLOBAL, 'time_avg', &
	664	TRIM( time_average_text ) )
	665	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 98 )
	666	ENDIF
	667
	668	!
	669	!-- Define time coordinate for xy sections (unlimited dimension)
	670	nc_stat = NF90_DEF_DIM( id_set_xy(av), 'time', NF90_UNLIMITED, &
	671	id_dim_time_xy(av) )
	672	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 99 )
	673
	674	nc_stat = NF90_DEF_VAR( id_set_xy(av), 'time', NF90_DOUBLE, &
	675	id_dim_time_xy(av), id_var_time_xy(av) )
	676	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 100 )
	677
	678	nc_stat = NF90_PUT_ATT( id_set_xy(av), id_var_time_xy(av), 'units', &
	679	'seconds')
	680	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 101 )
	681
	682	!
	683	!-- Define the spatial dimensions and coordinates for xy-sections.
	684	!-- First, determine the number of horizontal sections to be written.
	685	IF ( section(1,1) == -9999 ) THEN
	686	RETURN
	687	ELSE
	688	ns = 1
	689	DO WHILE ( section(ns,1) /= -9999 .AND. ns <= 100 )
	690	ns = ns + 1
	691	ENDDO
	692	ns = ns - 1
	693	ENDIF
	694
	695	!
	696	!-- Define vertical coordinate grid (zu grid)
	697	nc_stat = NF90_DEF_DIM( id_set_xy(av), 'zu_xy', ns, id_dim_zu_xy(av) )
	698	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 102 )
	699
	700	nc_stat = NF90_DEF_VAR( id_set_xy(av), 'zu_xy', NF90_DOUBLE, &
	701	id_dim_zu_xy(av), id_var_zu_xy(av) )
	702	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 103 )
	703
	704	nc_stat = NF90_PUT_ATT( id_set_xy(av), id_var_zu_xy(av), 'units', &
	705	'meters' )
	706	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 104 )
	707
	708	!
	709	!-- Define vertical coordinate grid (zw grid)
	710	nc_stat = NF90_DEF_DIM( id_set_xy(av), 'zw_xy', ns, id_dim_zw_xy(av) )
	711	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 105 )
	712
	713	nc_stat = NF90_DEF_VAR( id_set_xy(av), 'zw_xy', NF90_DOUBLE, &
	714	id_dim_zw_xy(av), id_var_zw_xy(av) )
	715	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 106 )
	716
	717	nc_stat = NF90_PUT_ATT( id_set_xy(av), id_var_zw_xy(av), 'units', &
	718	'meters' )
	719	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 107 )
	720
	721	!
	722	!-- Define a pseudo vertical coordinate grid for the surface variables
	723	!-- u* and t* to store their height level
	724	nc_stat = NF90_DEF_DIM( id_set_xy(av), 'zu1_xy', 1, &
	725	id_dim_zu1_xy(av) )
	726	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 108 )
	727
	728	nc_stat = NF90_DEF_VAR( id_set_xy(av), 'zu1_xy', NF90_DOUBLE, &
	729	id_dim_zu1_xy(av), id_var_zu1_xy(av) )
	730	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 109 )
	731
	732	nc_stat = NF90_PUT_ATT( id_set_xy(av), id_var_zu1_xy(av), 'units', &
	733	'meters' )
	734	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 110 )
	735
	736	!
	737	!-- Define a variable to store the layer indices of the horizontal cross
	738	!-- sections, too
	739	nc_stat = NF90_DEF_VAR( id_set_xy(av), 'ind_z_xy', NF90_DOUBLE, &
	740	id_dim_zu_xy(av), id_var_ind_z_xy(av) )
	741	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 111 )
	742
	743	nc_stat = NF90_PUT_ATT( id_set_xy(av), id_var_ind_z_xy(av), 'units', &
	744	'gridpoints')
	745	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 112 )
	746
	747	!
	748	!-- Define x-axis (for scalar position)
	749	nc_stat = NF90_DEF_DIM( id_set_xy(av), 'x', nx+2, id_dim_x_xy(av) )
	750	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 113 )
	751
	752	nc_stat = NF90_DEF_VAR( id_set_xy(av), 'x', NF90_DOUBLE, &
	753	id_dim_x_xy(av), id_var_x_xy(av) )
	754	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 114 )
	755
	756	nc_stat = NF90_PUT_ATT( id_set_xy(av), id_var_x_xy(av), 'units', &
	757	'meters' )
	758	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 115 )
	759
	760	!
	761	!-- Define x-axis (for u position)
	762	nc_stat = NF90_DEF_DIM( id_set_xy(av), 'xu', nx+2, id_dim_xu_xy(av) )
	763	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 388 )
	764
	765	nc_stat = NF90_DEF_VAR( id_set_xy(av), 'xu', NF90_DOUBLE, &
	766	id_dim_xu_xy(av), id_var_xu_xy(av) )
	767	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 389 )
	768
	769	nc_stat = NF90_PUT_ATT( id_set_xy(av), id_var_xu_xy(av), 'units', &
	770	'meters' )
	771	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 390 )
	772
	773	!
	774	!-- Define y-axis (for scalar position)
	775	nc_stat = NF90_DEF_DIM( id_set_xy(av), 'y', ny+2, id_dim_y_xy(av) )
	776	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 116 )
	777
	778	nc_stat = NF90_DEF_VAR( id_set_xy(av), 'y', NF90_DOUBLE, &
	779	id_dim_y_xy(av), id_var_y_xy(av) )
	780	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 117 )
	781
	782	nc_stat = NF90_PUT_ATT( id_set_xy(av), id_var_y_xy(av), 'units', &
	783	'meters' )
	784	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 118 )
	785
	786	!
	787	!-- Define y-axis (for scalar position)
	788	nc_stat = NF90_DEF_DIM( id_set_xy(av), 'yv', ny+2, id_dim_yv_xy(av) )
	789	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 364 )
	790
	791	nc_stat = NF90_DEF_VAR( id_set_xy(av), 'yv', NF90_DOUBLE, &
	792	id_dim_yv_xy(av), id_var_yv_xy(av) )
	793	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 365 )
	794
	795	nc_stat = NF90_PUT_ATT( id_set_xy(av), id_var_yv_xy(av), 'units', &
	796	'meters' )
	797	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 366 )
	798
	799	!
[48]	800	!-- In case of non-flat topography define 2d-arrays containing the height
	801	!-- informations
	802	IF ( TRIM( topography ) /= 'flat' ) THEN
	803	!
	804	!-- Define zusi = zu(nzb_s_inner)
	805	nc_stat = NF90_DEF_VAR( id_set_xy(av), 'zusi', NF90_DOUBLE, &
	806	(/ id_dim_x_xy(av), id_dim_y_xy(av) /), &
	807	id_var_zusi_xy(av) )
	808	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 421 )
	809
	810	nc_stat = NF90_PUT_ATT( id_set_xy(av), id_var_zusi_xy(av), &
	811	'units', 'meters' )
	812	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 422 )
	813
	814	nc_stat = NF90_PUT_ATT( id_set_xy(av), id_var_zusi_xy(av), &
	815	'long_name', 'zu(nzb_s_inner)' )
	816	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 423 )
	817
	818	!
	819	!-- Define zwwi = zw(nzb_w_inner)
	820	nc_stat = NF90_DEF_VAR( id_set_xy(av), 'zwwi', NF90_DOUBLE, &
	821	(/ id_dim_x_xy(av), id_dim_y_xy(av) /), &
	822	id_var_zwwi_xy(av) )
	823	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 424 )
	824
	825	nc_stat = NF90_PUT_ATT( id_set_xy(av), id_var_zwwi_xy(av), &
	826	'units', 'meters' )
	827	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 425 )
	828
	829	nc_stat = NF90_PUT_ATT( id_set_xy(av), id_var_zwwi_xy(av), &
	830	'long_name', 'zw(nzb_w_inner)' )
	831	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 426 )
	832
	833	ENDIF
	834
	835
	836	!
[1]	837	!-- Define the variables
	838	var_list = ';'
	839	i = 1
	840
	841	DO WHILE ( do2d(av,i)(1:1) /= ' ' )
	842
	843	IF ( INDEX( do2d(av,i), 'xy' ) /= 0 ) THEN
	844	!
	845	!-- If there is a star in the variable name (u* or t*), it is a
	846	!-- surface variable. Define it with id_dim_zu1_xy.
	847	IF ( INDEX( do2d(av,i), '*' ) /= 0 ) THEN
	848	!
	849	!-- First, remove those characters not allowed by NetCDF
	850	netcdf_var_name = do2d(av,i)
	851	CALL clean_netcdf_varname( netcdf_var_name )
	852
	853	nc_stat = NF90_DEF_VAR( id_set_xy(av), netcdf_var_name, &
	854	nc_precision(1), &
	855	(/ id_dim_x_xy(av), id_dim_y_xy(av),&
	856	id_dim_zu1_xy(av), &
	857	id_dim_time_xy(av) /), &
	858	id_var_do2d(av,i) )
	859
	860	var_list = TRIM(var_list) // TRIM(netcdf_var_name) // ';'
	861
	862	ELSE
	863
	864	!
	865	!-- Check for the grid
	866	found = .TRUE.
	867	SELECT CASE ( do2d(av,i) )
	868	!
	869	!-- Most variables are defined on the zu grid
	870	CASE ( 'e_xy', 'p_xy', 'pc_xy', 'pr_xy', 'pt_xy', 'q_xy',&
	871	'ql_xy', 'ql_c_xy', 'ql_v_xy', 'ql_vp_xy', &
	872	'qv_xy', 's_xy', 'vpt_xy' )
	873
	874	grid_x = 'x'
	875	grid_y = 'y'
	876	grid_z = 'zu'
	877	!
	878	!-- u grid
	879	CASE ( 'u_xy' )
	880
	881	grid_x = 'xu'
	882	grid_y = 'y'
	883	grid_z = 'zu'
	884	!
	885	!-- v grid
	886	CASE ( 'v_xy' )
	887
	888	grid_x = 'x'
	889	grid_y = 'yv'
	890	grid_z = 'zu'
	891	!
	892	!-- w grid
	893	CASE ( 'w_xy' )
	894
	895	grid_x = 'x'
	896	grid_y = 'y'
	897	grid_z = 'zw'
	898
	899	CASE DEFAULT
	900	!
	901	!-- Check for user-defined quantities
	902	CALL user_define_netcdf_grid( do2d(av,i), found, &
	903	grid_x, grid_y, grid_z )
	904
	905	END SELECT
	906
	907	!
	908	!-- Select the respective dimension ids
	909	IF ( grid_x == 'x' ) THEN
	910	id_x = id_dim_x_xy(av)
	911	ELSEIF ( grid_x == 'xu' ) THEN
	912	id_x = id_dim_xu_xy(av)
	913	ENDIF
	914
	915	IF ( grid_y == 'y' ) THEN
	916	id_y = id_dim_y_xy(av)
	917	ELSEIF ( grid_y == 'yv' ) THEN
	918	id_y = id_dim_yv_xy(av)
	919	ENDIF
	920
	921	IF ( grid_z == 'zu' ) THEN
	922	id_z = id_dim_zu_xy(av)
	923	ELSEIF ( grid_z == 'zw' ) THEN
	924	id_z = id_dim_zw_xy(av)
	925	ENDIF
	926
	927	!
	928	!-- Define the grid
	929	nc_stat = NF90_DEF_VAR( id_set_xy(av), do2d(av,i), &
	930	nc_precision(1), &
	931	(/ id_x, id_y, id_z, id_dim_time_xy(av) /), &
	932	id_var_do2d(av,i) )
	933
	934	IF ( .NOT. found ) THEN
	935	PRINT*, '+++ define_netcdf_header: no grid defined ', &
	936	'for variable ', do2d(av,i)
	937	ENDIF
	938
	939	var_list = TRIM( var_list ) // TRIM( do2d(av,i) ) // ';'
	940
	941	ENDIF
	942
	943	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 119 )
	944	!
	945	!-- Store the 'real' name of the variable (with *, for example)
	946	!-- in the long_name attribute. This is evaluated by Ferret,
	947	!-- for example.
	948	nc_stat = NF90_PUT_ATT( id_set_xy(av), id_var_do2d(av,i), &
	949	'long_name', do2d(av,i) )
	950	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 120 )
	951	!
	952	!-- Define the variable's unit
	953	nc_stat = NF90_PUT_ATT( id_set_xy(av), id_var_do2d(av,i), &
	954	'units', TRIM( do2d_unit(av,i) ) )
	955	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 354 )
	956	ENDIF
	957
	958	i = i + 1
	959
	960	ENDDO
	961
	962	!
	963	!-- No arrays to output. Close the netcdf file and return.
	964	IF ( i == 1 ) RETURN
	965
	966	!
	967	!-- Write the list of variables as global attribute (this is used by
	968	!-- restart runs and by combine_plot_fields)
	969	nc_stat = NF90_PUT_ATT( id_set_xy(av), NF90_GLOBAL, 'VAR_LIST', &
	970	var_list )
	971	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 121 )
	972
	973	!
	974	!-- Leave NetCDF define mode
	975	nc_stat = NF90_ENDDEF( id_set_xy(av) )
	976	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 122 )
	977
	978	!
	979	!-- Write axis data: z_xy, x, y
	980	ALLOCATE( netcdf_data(1:ns) )
	981
	982	!
	983	!-- Write zu data
	984	DO i = 1, ns
	985	IF( section(i,1) == -1 ) THEN
	986	netcdf_data(i) = -1.0 ! section averaged along z
	987	ELSE
	988	netcdf_data(i) = zu( section(i,1) )
	989	ENDIF
	990	ENDDO
	991	nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zu_xy(av), &
	992	netcdf_data, start = (/ 1 /), &
	993	count = (/ ns /) )
	994	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 123 )
	995
	996	!
	997	!-- Write zw data
	998	DO i = 1, ns
	999	IF( section(i,1) == -1 ) THEN
	1000	netcdf_data(i) = -1.0 ! section averaged along z
	1001	ELSE
	1002	netcdf_data(i) = zw( section(i,1) )
	1003	ENDIF
	1004	ENDDO
	1005	nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zw_xy(av), &
	1006	netcdf_data, start = (/ 1 /), &
	1007	count = (/ ns /) )
	1008	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 124 )
	1009
	1010	!
	1011	!-- Write gridpoint number data
	1012	netcdf_data(1:ns) = section(1:ns,1)
	1013	nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_ind_z_xy(av), &
	1014	netcdf_data, start = (/ 1 /), &
	1015	count = (/ ns /) )
	1016	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 125 )
	1017
	1018	DEALLOCATE( netcdf_data )
	1019
	1020	!
	1021	!-- Write the cross section height u, t
	1022	nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zu1_xy(av), &
	1023	(/ zu(nzb+1) /), start = (/ 1 /), &
	1024	count = (/ 1 /) )
	1025	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 126 )
	1026
	1027	!
	1028	!-- Write data for x and xu axis (shifted by -dx/2)
	1029	ALLOCATE( netcdf_data(0:nx+1) )
	1030
	1031	DO i = 0, nx+1
	1032	netcdf_data(i) = i * dx
	1033	ENDDO
	1034
	1035	nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_x_xy(av), netcdf_data, &
	1036	start = (/ 1 /), count = (/ nx+2 /) )
	1037	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 127 )
	1038
	1039	DO i = 0, nx+1
	1040	netcdf_data(i) = ( i - 0.5 ) * dx
	1041	ENDDO
	1042
	1043	nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_xu_xy(av), &
	1044	netcdf_data, start = (/ 1 /), &
	1045	count = (/ nx+2 /) )
	1046	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 367 )
	1047
	1048	DEALLOCATE( netcdf_data )
	1049
	1050	!
	1051	!-- Write data for y and yv axis (shifted by -dy/2)
	1052	ALLOCATE( netcdf_data(0:ny+1) )
	1053
	1054	DO i = 0, ny+1
	1055	netcdf_data(i) = i * dy
	1056	ENDDO
	1057
	1058	nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_y_xy(av), netcdf_data, &
	1059	start = (/ 1 /), count = (/ ny+2 /))
	1060	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 128 )
	1061
	1062	DO i = 0, ny+1
	1063	netcdf_data(i) = ( i - 0.5 ) * dy
	1064	ENDDO
	1065
	1066	nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_yv_xy(av), &
	1067	netcdf_data, start = (/ 1 /), &
	1068	count = (/ ny+2 /))
	1069	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 368 )
	1070
	1071	DEALLOCATE( netcdf_data )
	1072
[48]	1073	!
	1074	!-- In case of non-flat topography write height information
	1075	IF ( TRIM( topography ) /= 'flat' ) THEN
[1]	1076
[48]	1077	nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zusi_xy(av), &
	1078	zu_s_inner(0:nx+1,0:ny+1), &
	1079	start = (/ 1, 1 /), &
	1080	count = (/ nx+2, ny+2 /) )
	1081	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 427 )
	1082
	1083	nc_stat = NF90_PUT_VAR( id_set_xy(av), id_var_zwwi_xy(av), &
	1084	zw_w_inner(0:nx+1,0:ny+1), &
	1085	start = (/ 1, 1 /), &
	1086	count = (/ nx+2, ny+2 /) )
	1087	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 428 )
	1088
	1089	ENDIF
	1090
	1091
[1]	1092	CASE ( 'xy_ext' )
	1093
	1094	!
	1095	!-- Get the list of variables and compare with the actual run.
	1096	!-- First var_list_old has to be reset, since GET_ATT does not assign
	1097	!-- trailing blanks.
	1098	var_list_old = ' '
	1099	nc_stat = NF90_GET_ATT( id_set_xy(av), NF90_GLOBAL, 'VAR_LIST', &
	1100	var_list_old )
	1101	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 129 )
	1102
	1103	var_list = ';'
	1104	i = 1
	1105	DO WHILE ( do2d(av,i)(1:1) /= ' ' )
	1106	IF ( INDEX( do2d(av,i), 'xy' ) /= 0 ) THEN
	1107	netcdf_var_name = do2d(av,i)
	1108	CALL clean_netcdf_varname( netcdf_var_name )
	1109	var_list = TRIM(var_list) // TRIM(netcdf_var_name) // ';'
	1110	ENDIF
	1111	i = i + 1
	1112	ENDDO
	1113
	1114	IF ( av == 0 ) THEN
	1115	var = '(xy)'
	1116	ELSE
	1117	var = '(xy_av)'
	1118	ENDIF
	1119
	1120	IF ( TRIM( var_list ) /= TRIM( var_list_old ) ) THEN
	1121	PRINT*, '+++ WARNING: NetCDF file for cross-sections ' // &
	1122	TRIM( var ) // ' from previuos run found,'
	1123	PRINT*, ' but this file cannot be extended due to' // &
	1124	' variable mismatch.'
	1125	PRINT*, ' New file is created instead.'
	1126	extend = .FALSE.
	1127	RETURN
	1128	ENDIF
	1129
	1130	!
	1131	!-- Calculate the number of current sections
	1132	ns = 1
	1133	DO WHILE ( section(ns,1) /= -9999 .AND. ns <= 100 )
	1134	ns = ns + 1
	1135	ENDDO
	1136	ns = ns - 1
	1137
	1138	!
	1139	!-- Get and compare the number of horizontal cross sections
	1140	nc_stat = NF90_INQ_VARID( id_set_xy(av), 'zu_xy', id_var_zu_xy(av) )
	1141	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 130 )
	1142
	1143	nc_stat = NF90_INQUIRE_VARIABLE( id_set_xy(av), id_var_zu_xy(av), &
	1144	dimids = id_dim_zu_xy_old )
	1145	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 131 )
	1146	id_dim_zu_xy(av) = id_dim_zu_xy_old(1)
	1147
	1148	nc_stat = NF90_INQUIRE_DIMENSION( id_set_xy(av), id_dim_zu_xy(av), &
	1149	len = ns_old )
	1150	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 132 )
	1151
	1152	IF ( ns /= ns_old ) THEN
	1153	PRINT*, '+++ WARNING: NetCDF file for cross-sections ' // &
	1154	TRIM( var ) // ' from previuos run found,'
	1155	PRINT*, ' but this file cannot be extended due to' // &
	1156	' mismatch in number of'
	1157	PRINT*, ' cross sections.'
	1158	PRINT*, ' New file is created instead.'
	1159	extend = .FALSE.
	1160	RETURN
	1161	ENDIF
	1162
	1163	!
	1164	!-- Get and compare the heights of the cross sections
	1165	ALLOCATE( netcdf_data(1:ns_old) )
	1166
	1167	nc_stat = NF90_GET_VAR( id_set_xy(av), id_var_zu_xy(av), netcdf_data )
	1168	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 133 )
	1169
	1170	DO i = 1, ns
	1171	IF ( section(i,1) /= -1 ) THEN
	1172	IF ( zu(section(i,1)) /= netcdf_data(i) ) THEN
	1173	PRINT*, '+++ WARNING: NetCDF file for cross-sections ' // &
	1174	TRIM( var ) // ' from previuos run found,'
	1175	PRINT*, ' but this file cannot be extended' // &
	1176	' due to mismatch in cross'
	1177	PRINT*, ' section levels.'
	1178	PRINT*, ' New file is created instead.'
	1179	extend = .FALSE.
	1180	RETURN
	1181	ENDIF
	1182	ELSE
	1183	IF ( -1.0 /= netcdf_data(i) ) THEN
	1184	PRINT*, '+++ WARNING: NetCDF file for cross-sections ' // &
	1185	TRIM( var ) // ' from previuos run found,'
	1186	PRINT*, ' but this file cannot be extended' // &
	1187	' due to mismatch in cross'
	1188	PRINT*, ' section levels.'
	1189	PRINT*, ' New file is created instead.'
	1190	extend = .FALSE.
	1191	RETURN
	1192	ENDIF
	1193	ENDIF
	1194	ENDDO
	1195
	1196	DEALLOCATE( netcdf_data )
	1197
	1198	!
	1199	!-- Get the id of the time coordinate (unlimited coordinate) and its
	1200	!-- last index on the file. The next time level is do2d..count+1.
	1201	!-- The current time must be larger than the last output time
	1202	!-- on the file.
	1203	nc_stat = NF90_INQ_VARID( id_set_xy(av), 'time', id_var_time_xy(av) )
	1204	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 134 )
	1205
	1206	nc_stat = NF90_INQUIRE_VARIABLE( id_set_xy(av), id_var_time_xy(av), &
	1207	dimids = id_dim_time_old )
	1208	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 135 )
	1209	id_dim_time_xy(av) = id_dim_time_old(1)
	1210
	1211	nc_stat = NF90_INQUIRE_DIMENSION( id_set_xy(av), id_dim_time_xy(av), &
	1212	len = do2d_xy_time_count(av) )
	1213	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 136 )
	1214
	1215	nc_stat = NF90_GET_VAR( id_set_xy(av), id_var_time_xy(av), &
	1216	last_time_coordinate, &
	1217	start = (/ do2d_xy_time_count(av) /), &
	1218	count = (/ 1 /) )
	1219	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 137 )
	1220
	1221	IF ( last_time_coordinate(1) >= simulated_time ) THEN
	1222	PRINT*, '+++ WARNING: NetCDF file for cross sections ' // &
	1223	TRIM( var ) // ' from previuos run found,'
	1224	PRINT*, ' but this file cannot be extended becaus' // &
	1225	'e the current output time'
	1226	PRINT*, ' is less or equal than the last output t' // &
	1227	'ime on this file.'
	1228	PRINT*, ' New file is created instead.'
	1229	do2d_xy_time_count(av) = 0
	1230	extend = .FALSE.
	1231	RETURN
	1232	ENDIF
	1233
	1234	!
	1235	!-- Dataset seems to be extendable.
	1236	!-- Now get the variable ids.
	1237	i = 1
	1238	DO WHILE ( do2d(av,i)(1:1) /= ' ' )
	1239	IF ( INDEX( do2d(av,i), 'xy' ) /= 0 ) THEN
	1240	netcdf_var_name = do2d(av,i)
	1241	CALL clean_netcdf_varname( netcdf_var_name )
	1242	nc_stat = NF90_INQ_VARID( id_set_xy(av), netcdf_var_name, &
	1243	id_var_do2d(av,i) )
	1244	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 138 )
	1245	ENDIF
	1246	i = i + 1
	1247	ENDDO
	1248
	1249	!
	1250	!-- Change the titel attribute on file
	1251	nc_stat = NF90_PUT_ATT( id_set_xy(av), NF90_GLOBAL, 'title', &
	1252	TRIM( run_description_header ) )
	1253	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 139 )
	1254	PRINT, '** NetCDF file for cross-sections ' // TRIM( var ) // &
	1255	' from previous run found.'
	1256	PRINT*, ' This file will be extended.'
	1257
	1258
	1259	CASE ( 'xz_new' )
	1260
	1261	!
	1262	!-- Define some global attributes of the dataset
	1263	nc_stat = NF90_PUT_ATT( id_set_xz(av), NF90_GLOBAL, 'Conventions', &
	1264	'COARDS' )
	1265	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 140 )
	1266
	1267	IF ( av == 0 ) THEN
	1268	time_average_text = ' '
	1269	ELSE
	1270	WRITE (time_average_text, '('', '',F7.1,'' s average'')') &
	1271	averaging_interval
	1272	ENDIF
	1273	nc_stat = NF90_PUT_ATT( id_set_xz(av), NF90_GLOBAL, 'title', &
	1274	TRIM( run_description_header ) // &
	1275	TRIM( time_average_text ) )
	1276	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 141 )
	1277	IF ( av == 1 ) THEN
	1278	WRITE ( time_average_text,'(F7.1,'' s avg'')' ) averaging_interval
	1279	nc_stat = NF90_PUT_ATT( id_set_xz(av), NF90_GLOBAL, 'time_avg', &
	1280	TRIM( time_average_text ) )
	1281	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 141 )
	1282	ENDIF
	1283
	1284	!
	1285	!-- Define time coordinate for xz sections (unlimited dimension)
	1286	nc_stat = NF90_DEF_DIM( id_set_xz(av), 'time', NF90_UNLIMITED, &
	1287	id_dim_time_xz(av) )
	1288	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 142 )
	1289
	1290	nc_stat = NF90_DEF_VAR( id_set_xz(av), 'time', NF90_DOUBLE, &
	1291	id_dim_time_xz(av), id_var_time_xz(av) )
	1292	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 143 )
	1293
	1294	nc_stat = NF90_PUT_ATT( id_set_xz(av), id_var_time_xz(av), 'units', &
	1295	'seconds')
	1296	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 144 )
	1297
	1298	!
	1299	!-- Define the spatial dimensions and coordinates for xz-sections.
	1300	!-- First, determine the number of vertical sections to be written.
	1301	IF ( section(1,2) == -9999 ) THEN
	1302	RETURN
	1303	ELSE
	1304	ns = 1
	1305	DO WHILE ( section(ns,2) /= -9999 .AND. ns <= 100 )
	1306	ns = ns + 1
	1307	ENDDO
	1308	ns = ns - 1
	1309	ENDIF
	1310
	1311	!
	1312	!-- Define y-axis (for scalar position)
	1313	nc_stat = NF90_DEF_DIM( id_set_xz(av), 'y_xz', ns, id_dim_y_xz(av) )
	1314	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 145 )
	1315
	1316	nc_stat = NF90_DEF_VAR( id_set_xz(av), 'y_xz', NF90_DOUBLE, &
	1317	id_dim_y_xz(av), id_var_y_xz(av) )
	1318	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 146 )
	1319
	1320	nc_stat = NF90_PUT_ATT( id_set_xz(av), id_var_y_xz(av), 'units', &
	1321	'meters' )
	1322	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 147 )
	1323
	1324	!
	1325	!-- Define y-axis (for v position)
	1326	nc_stat = NF90_DEF_DIM( id_set_xz(av), 'yv_xz', ns, id_dim_yv_xz(av) )
	1327	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 369 )
	1328
	1329	nc_stat = NF90_DEF_VAR( id_set_xz(av), 'yv_xz', NF90_DOUBLE, &
	1330	id_dim_yv_xz(av), id_var_yv_xz(av) )
	1331	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 370 )
	1332
	1333	nc_stat = NF90_PUT_ATT( id_set_xz(av), id_var_yv_xz(av), 'units', &
	1334	'meters' )
	1335	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 371 )
	1336
	1337	!
	1338	!-- Define a variable to store the layer indices of the vertical cross
	1339	!-- sections
	1340	nc_stat = NF90_DEF_VAR( id_set_xz(av), 'ind_y_xz', NF90_DOUBLE, &
	1341	id_dim_y_xz(av), id_var_ind_y_xz(av) )
	1342	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 148 )
	1343
	1344	nc_stat = NF90_PUT_ATT( id_set_xz(av), id_var_ind_y_xz(av), 'units', &
	1345	'gridpoints')
	1346	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 149 )
	1347
	1348	!
	1349	!-- Define x-axis (for scalar position)
	1350	nc_stat = NF90_DEF_DIM( id_set_xz(av), 'x', nx+2, id_dim_x_xz(av) )
	1351	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 150 )
	1352
	1353	nc_stat = NF90_DEF_VAR( id_set_xz(av), 'x', NF90_DOUBLE, &
	1354	id_dim_x_xz(av), id_var_x_xz(av) )
	1355	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 151 )
	1356
	1357	nc_stat = NF90_PUT_ATT( id_set_xz(av), id_var_x_xz(av), 'units', &
	1358	'meters' )
	1359	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 152 )
	1360
	1361	!
	1362	!-- Define x-axis (for u position)
	1363	nc_stat = NF90_DEF_DIM( id_set_xz(av), 'xu', nx+2, id_dim_xu_xz(av) )
	1364	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 372 )
	1365
	1366	nc_stat = NF90_DEF_VAR( id_set_xz(av), 'xu', NF90_DOUBLE, &
	1367	id_dim_xu_xz(av), id_var_xu_xz(av) )
	1368	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 373 )
	1369
	1370	nc_stat = NF90_PUT_ATT( id_set_xz(av), id_var_xu_xz(av), 'units', &
	1371	'meters' )
	1372	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 374 )
	1373
	1374	!
	1375	!-- Define the two z-axes (zu and zw)
	1376	nc_stat = NF90_DEF_DIM( id_set_xz(av), 'zu', nz+2, id_dim_zu_xz(av) )
	1377	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 153 )
	1378
	1379	nc_stat = NF90_DEF_VAR( id_set_xz(av), 'zu', NF90_DOUBLE, &
	1380	id_dim_zu_xz(av), id_var_zu_xz(av) )
	1381	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 154 )
	1382
	1383	nc_stat = NF90_PUT_ATT( id_set_xz(av), id_var_zu_xz(av), 'units', &
	1384	'meters' )
	1385	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 155 )
	1386
	1387	nc_stat = NF90_DEF_DIM( id_set_xz(av), 'zw', nz+2, id_dim_zw_xz(av) )
	1388	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 156 )
	1389
	1390	nc_stat = NF90_DEF_VAR( id_set_xz(av), 'zw', NF90_DOUBLE, &
	1391	id_dim_zw_xz(av), id_var_zw_xz(av) )
	1392	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 157 )
	1393
	1394	nc_stat = NF90_PUT_ATT( id_set_xz(av), id_var_zw_xz(av), 'units', &
	1395	'meters' )
	1396	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 158 )
	1397
	1398
	1399	!
	1400	!-- Define the variables
	1401	var_list = ';'
	1402	i = 1
	1403
	1404	DO WHILE ( do2d(av,i)(1:1) /= ' ' )
	1405
	1406	IF ( INDEX( do2d(av,i), 'xz' ) /= 0 ) THEN
	1407
	1408	!
	1409	!-- Check for the grid
	1410	found = .TRUE.
	1411	SELECT CASE ( do2d(av,i) )
	1412	!
	1413	!-- Most variables are defined on the zu grid
	1414	CASE ( 'e_xz', 'p_xz', 'pc_xz', 'pr_xz', 'pt_xz', 'q_xz', &
	1415	'ql_xz', 'ql_c_xz', 'ql_v_xz', 'ql_vp_xz', 'qv_xz', &
	1416	's_xz', 'vpt_xz' )
	1417
	1418	grid_x = 'x'
	1419	grid_y = 'y'
	1420	grid_z = 'zu'
	1421	!
	1422	!-- u grid
	1423	CASE ( 'u_xz' )
	1424
	1425	grid_x = 'xu'
	1426	grid_y = 'y'
	1427	grid_z = 'zu'
	1428	!
	1429	!-- v grid
	1430	CASE ( 'v_xz' )
	1431
	1432	grid_x = 'x'
	1433	grid_y = 'yv'
	1434	grid_z = 'zu'
	1435	!
	1436	!-- w grid
	1437	CASE ( 'w_xz' )
	1438
	1439	grid_x = 'x'
	1440	grid_y = 'y'
	1441	grid_z = 'zw'
	1442
	1443	CASE DEFAULT
	1444	!
	1445	!-- Check for user-defined quantities
	1446	CALL user_define_netcdf_grid( do2d(av,i), found, &
	1447	grid_x, grid_y, grid_z )
	1448
	1449	END SELECT
	1450
	1451	!
	1452	!-- Select the respective dimension ids
	1453	IF ( grid_x == 'x' ) THEN
	1454	id_x = id_dim_x_xz(av)
	1455	ELSEIF ( grid_x == 'xu' ) THEN
	1456	id_x = id_dim_xu_xz(av)
	1457	ENDIF
	1458
	1459	IF ( grid_y == 'y' ) THEN
	1460	id_y = id_dim_y_xz(av)
	1461	ELSEIF ( grid_y == 'yv' ) THEN
	1462	id_y = id_dim_yv_xz(av)
	1463	ENDIF
	1464
	1465	IF ( grid_z == 'zu' ) THEN
	1466	id_z = id_dim_zu_xz(av)
	1467	ELSEIF ( grid_z == 'zw' ) THEN
	1468	id_z = id_dim_zw_xz(av)
	1469	ENDIF
	1470
	1471	!
	1472	!-- Define the grid
	1473	nc_stat = NF90_DEF_VAR( id_set_xz(av), do2d(av,i), &
	1474	nc_precision(2), &
	1475	(/ id_x, id_y, id_z, id_dim_time_xz(av) /), &
	1476	id_var_do2d(av,i) )
	1477
	1478	IF ( .NOT. found ) THEN
	1479	PRINT*, '+++ define_netcdf_header: no grid defined for', &
	1480	' variable ', do2d(av,i)
	1481	ENDIF
	1482
	1483	var_list = TRIM( var_list ) // TRIM( do2d(av,i) ) // ';'
	1484
	1485	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 159 )
	1486	!
	1487	!-- Store the 'real' name of the variable (with *, for example)
	1488	!-- in the long_name attribute. This is evaluated by Ferret,
	1489	!-- for example.
	1490	nc_stat = NF90_PUT_ATT( id_set_xz(av), id_var_do2d(av,i), &
	1491	'long_name', do2d(av,i) )
	1492	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 160 )
	1493	!
	1494	!-- Define the variable's unit
	1495	nc_stat = NF90_PUT_ATT( id_set_xz(av), id_var_do2d(av,i), &
	1496	'units', TRIM( do2d_unit(av,i) ) )
	1497	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 355 )
	1498	ENDIF
	1499
	1500	i = i + 1
	1501
	1502	ENDDO
	1503
	1504	!
	1505	!-- No arrays to output. Close the netcdf file and return.
	1506	IF ( i == 1 ) RETURN
	1507
	1508	!
	1509	!-- Write the list of variables as global attribute (this is used by
	1510	!-- restart runs and by combine_plot_fields)
	1511	nc_stat = NF90_PUT_ATT( id_set_xz(av), NF90_GLOBAL, 'VAR_LIST', &
	1512	var_list )
	1513	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 161 )
	1514
	1515	!
	1516	!-- Leave NetCDF define mode
	1517	nc_stat = NF90_ENDDEF( id_set_xz(av) )
	1518	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 162 )
	1519
	1520	!
	1521	!-- Write axis data: y_xz, x, zu, zw
	1522	ALLOCATE( netcdf_data(1:ns) )
	1523
	1524	!
	1525	!-- Write y_xz data
	1526	DO i = 1, ns
	1527	IF( section(i,2) == -1 ) THEN
	1528	netcdf_data(i) = -1.0 ! section averaged along y
	1529	ELSE
	1530	netcdf_data(i) = section(i,2) * dy
	1531	ENDIF
	1532	ENDDO
	1533	nc_stat = NF90_PUT_VAR( id_set_xz(av), id_var_y_xz(av), netcdf_data, &
	1534	start = (/ 1 /), count = (/ ns /) )
	1535	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 163 )
	1536
	1537	!
	1538	!-- Write yv_xz data
	1539	DO i = 1, ns
	1540	IF( section(i,2) == -1 ) THEN
	1541	netcdf_data(i) = -1.0 ! section averaged along y
	1542	ELSE
	1543	netcdf_data(i) = ( section(i,2) - 0.5 ) * dy
	1544	ENDIF
	1545	ENDDO
	1546	nc_stat = NF90_PUT_VAR( id_set_xz(av), id_var_yv_xz(av), &
	1547	netcdf_data, start = (/ 1 /), &
	1548	count = (/ ns /) )
	1549	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 375 )
	1550
	1551	!
	1552	!-- Write gridpoint number data
	1553	netcdf_data(1:ns) = section(1:ns,2)
	1554	nc_stat = NF90_PUT_VAR( id_set_xz(av), id_var_ind_y_xz(av), &
	1555	netcdf_data, start = (/ 1 /), &
	1556	count = (/ ns /) )
	1557	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 164 )
	1558
	1559	DEALLOCATE( netcdf_data )
	1560
	1561	!
	1562	!-- Write data for x and xu axis (shifted by -dx/2)
	1563	ALLOCATE( netcdf_data(0:nx+1) )
	1564
	1565	DO i = 0, nx+1
	1566	netcdf_data(i) = i * dx
	1567	ENDDO
	1568
	1569	nc_stat = NF90_PUT_VAR( id_set_xz(av), id_var_x_xz(av), netcdf_data, &
	1570	start = (/ 1 /), count = (/ nx+2 /) )
	1571	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 165 )
	1572
	1573	DO i = 0, nx+1
	1574	netcdf_data(i) = ( i - 0.5 ) * dx
	1575	ENDDO
	1576
	1577	nc_stat = NF90_PUT_VAR( id_set_xz(av), id_var_xu_xz(av), &
	1578	netcdf_data, start = (/ 1 /), &
	1579	count = (/ nx+2 /) )
	1580	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 376 )
	1581
	1582	DEALLOCATE( netcdf_data )
	1583
	1584	!
	1585	!-- Write zu and zw data (vertical axes)
	1586	ALLOCATE( netcdf_data(0:nz+1) )
	1587
	1588	netcdf_data(0:nz+1) = zu(nzb:nzt+1)
	1589	nc_stat = NF90_PUT_VAR( id_set_xz(av), id_var_zu_xz(av), &
	1590	netcdf_data, start = (/ 1 /), &
	1591	count = (/ nz+2 /) )
	1592	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 166 )
	1593
	1594	netcdf_data(0:nz+1) = zw(nzb:nzt+1)
	1595	nc_stat = NF90_PUT_VAR( id_set_xz(av), id_var_zw_xz(av), &
	1596	netcdf_data, start = (/ 1 /), &
	1597	count = (/ nz+2 /) )
	1598	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 167 )
	1599
	1600	DEALLOCATE( netcdf_data )
	1601
	1602
	1603	CASE ( 'xz_ext' )
	1604
	1605	!
	1606	!-- Get the list of variables and compare with the actual run.
	1607	!-- First var_list_old has to be reset, since GET_ATT does not assign
	1608	!-- trailing blanks.
	1609	var_list_old = ' '
	1610	nc_stat = NF90_GET_ATT( id_set_xz(av), NF90_GLOBAL, 'VAR_LIST', &
	1611	var_list_old )
	1612	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 168 )
	1613
	1614	var_list = ';'
	1615	i = 1
	1616	DO WHILE ( do2d(av,i)(1:1) /= ' ' )
	1617	IF ( INDEX( do2d(av,i), 'xz' ) /= 0 ) THEN
	1618	netcdf_var_name = do2d(av,i)
	1619	CALL clean_netcdf_varname( netcdf_var_name )
	1620	var_list = TRIM(var_list) // TRIM(netcdf_var_name) // ';'
	1621	ENDIF
	1622	i = i + 1
	1623	ENDDO
	1624
	1625	IF ( av == 0 ) THEN
	1626	var = '(xz)'
	1627	ELSE
	1628	var = '(xz_av)'
	1629	ENDIF
	1630
	1631	IF ( TRIM( var_list ) /= TRIM( var_list_old ) ) THEN
	1632	PRINT*, '+++ WARNING: NetCDF file for cross-sections ' // &
	1633	TRIM( var ) // ' from previuos run found,'
	1634	PRINT*, ' but this file cannot be extended due to' // &
	1635	' variable mismatch.'
	1636	PRINT*, ' New file is created instead.'
	1637	extend = .FALSE.
	1638	RETURN
	1639	ENDIF
	1640
	1641	!
	1642	!-- Calculate the number of current sections
	1643	ns = 1
	1644	DO WHILE ( section(ns,2) /= -9999 .AND. ns <= 100 )
	1645	ns = ns + 1
	1646	ENDDO
	1647	ns = ns - 1
	1648
	1649	!
	1650	!-- Get and compare the number of vertical cross sections
	1651	nc_stat = NF90_INQ_VARID( id_set_xz(av), 'y_xz', id_var_y_xz(av) )
	1652	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 169 )
	1653
	1654	nc_stat = NF90_INQUIRE_VARIABLE( id_set_xz(av), id_var_y_xz(av), &
	1655	dimids = id_dim_y_xz_old )
	1656	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 170 )
	1657	id_dim_y_xz(av) = id_dim_y_xz_old(1)
	1658
	1659	nc_stat = NF90_INQUIRE_DIMENSION( id_set_xz(av), id_dim_y_xz(av), &
	1660	len = ns_old )
	1661	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 171 )
	1662
	1663	IF ( ns /= ns_old ) THEN
	1664	PRINT*, '+++ WARNING: NetCDF file for cross-sections ' // &
	1665	TRIM( var ) // ' from previuos run found,'
	1666	PRINT*, ' but this file cannot be extended due to' // &
	1667	' mismatch in number of'
	1668	PRINT*, ' cross sections.'
	1669	PRINT*, ' New file is created instead.'
	1670	extend = .FALSE.
	1671	RETURN
	1672	ENDIF
	1673
	1674	!
	1675	!-- Get and compare the heights of the cross sections
	1676	ALLOCATE( netcdf_data(1:ns_old) )
	1677
	1678	nc_stat = NF90_GET_VAR( id_set_xz(av), id_var_y_xz(av), netcdf_data )
	1679	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 172 )
	1680
	1681	DO i = 1, ns
	1682	IF ( section(i,2) /= -1 ) THEN
	1683	IF ( ( section(i,2) * dy ) /= netcdf_data(i) ) THEN
	1684	PRINT*, '+++ WARNING: NetCDF file for cross-sections ' // &
	1685	TRIM( var ) // ' from previuos run found,'
	1686	PRINT*, ' but this file cannot be extended' // &
	1687	' due to mismatch in cross'
	1688	PRINT*, ' section indices.'
	1689	PRINT*, ' New file is created instead.'
	1690	extend = .FALSE.
	1691	RETURN
	1692	ENDIF
	1693	ELSE
	1694	IF ( -1.0 /= netcdf_data(i) ) THEN
	1695	PRINT*, '+++ WARNING: NetCDF file for cross-sections ' // &
	1696	TRIM( var ) // ' from previuos run found,'
	1697	PRINT*, ' but this file cannot be extended' // &
	1698	' due to mismatch in cross'
	1699	PRINT*, ' section indices.'
	1700	PRINT*, ' New file is created instead.'
	1701	extend = .FALSE.
	1702	RETURN
	1703	ENDIF
	1704	ENDIF
	1705	ENDDO
	1706
	1707	DEALLOCATE( netcdf_data )
	1708
	1709	!
	1710	!-- Get the id of the time coordinate (unlimited coordinate) and its
	1711	!-- last index on the file. The next time level is do2d..count+1.
	1712	!-- The current time must be larger than the last output time
	1713	!-- on the file.
	1714	nc_stat = NF90_INQ_VARID( id_set_xz(av), 'time', id_var_time_xz(av) )
	1715	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 173 )
	1716
	1717	nc_stat = NF90_INQUIRE_VARIABLE( id_set_xz(av), id_var_time_xz(av), &
	1718	dimids = id_dim_time_old )
	1719	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 174 )
	1720	id_dim_time_xz(av) = id_dim_time_old(1)
	1721
	1722	nc_stat = NF90_INQUIRE_DIMENSION( id_set_xz(av), id_dim_time_xz(av), &
	1723	len = do2d_xz_time_count(av) )
	1724	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 175 )
	1725
	1726	nc_stat = NF90_GET_VAR( id_set_xz(av), id_var_time_xz(av), &
	1727	last_time_coordinate, &
	1728	start = (/ do2d_xz_time_count(av) /), &
	1729	count = (/ 1 /) )
	1730	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 176 )
	1731
	1732	IF ( last_time_coordinate(1) >= simulated_time ) THEN
	1733	PRINT*, '+++ WARNING: NetCDF file for cross sections ' // &
	1734	TRIM( var ) // ' from previuos run found,'
	1735	PRINT*, ' but this file cannot be extended becaus' // &
	1736	'e the current output time'
	1737	PRINT*, ' is less or equal than the last output t' // &
	1738	'ime on this file.'
	1739	PRINT*, ' New file is created instead.'
	1740	do2d_xz_time_count(av) = 0
	1741	extend = .FALSE.
	1742	RETURN
	1743	ENDIF
	1744
	1745	!
	1746	!-- Dataset seems to be extendable.
	1747	!-- Now get the variable ids.
	1748	i = 1
	1749	DO WHILE ( do2d(av,i)(1:1) /= ' ' )
	1750	IF ( INDEX( do2d(av,i), 'xz' ) /= 0 ) THEN
	1751	netcdf_var_name = do2d(av,i)
	1752	CALL clean_netcdf_varname( netcdf_var_name )
	1753	nc_stat = NF90_INQ_VARID( id_set_xz(av), netcdf_var_name, &
	1754	id_var_do2d(av,i) )
	1755	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 177 )
	1756	ENDIF
	1757	i = i + 1
	1758	ENDDO
	1759
	1760	!
	1761	!-- Change the titel attribute on file
	1762	nc_stat = NF90_PUT_ATT( id_set_xz(av), NF90_GLOBAL, 'title', &
	1763	TRIM( run_description_header ) )
	1764	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 178 )
	1765	PRINT, '** NetCDF file for cross-sections ' // TRIM( var ) // &
	1766	' from previous run found.'
	1767	PRINT*, ' This file will be extended.'
	1768
	1769
	1770	CASE ( 'yz_new' )
	1771
	1772	!
	1773	!-- Define some global attributes of the dataset
	1774	nc_stat = NF90_PUT_ATT( id_set_yz(av), NF90_GLOBAL, 'Conventions', &
	1775	'COARDS' )
	1776	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 179 )
	1777
	1778	IF ( av == 0 ) THEN
	1779	time_average_text = ' '
	1780	ELSE
	1781	WRITE (time_average_text, '('', '',F7.1,'' s average'')') &
	1782	averaging_interval
	1783	ENDIF
	1784	nc_stat = NF90_PUT_ATT( id_set_yz(av), NF90_GLOBAL, 'title', &
	1785	TRIM( run_description_header ) // &
	1786	TRIM( time_average_text ) )
	1787	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 180 )
	1788	IF ( av == 1 ) THEN
	1789	WRITE ( time_average_text,'(F7.1,'' s avg'')' ) averaging_interval
	1790	nc_stat = NF90_PUT_ATT( id_set_yz(av), NF90_GLOBAL, 'time_avg', &
	1791	TRIM( time_average_text ) )
	1792	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 180 )
	1793	ENDIF
	1794
	1795	!
	1796	!-- Define time coordinate for yz sections (unlimited dimension)
	1797	nc_stat = NF90_DEF_DIM( id_set_yz(av), 'time', NF90_UNLIMITED, &
	1798	id_dim_time_yz(av) )
	1799	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 181 )
	1800
	1801	nc_stat = NF90_DEF_VAR( id_set_yz(av), 'time', NF90_DOUBLE, &
	1802	id_dim_time_yz(av), id_var_time_yz(av) )
	1803	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 182 )
	1804
	1805	nc_stat = NF90_PUT_ATT( id_set_yz(av), id_var_time_yz(av), 'units', &
	1806	'seconds')
	1807	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 183 )
	1808
	1809	!
	1810	!-- Define the spatial dimensions and coordinates for yz-sections.
	1811	!-- First, determine the number of vertical sections to be written.
	1812	IF ( section(1,3) == -9999 ) THEN
	1813	RETURN
	1814	ELSE
	1815	ns = 1
	1816	DO WHILE ( section(ns,3) /= -9999 .AND. ns <= 100 )
	1817	ns = ns + 1
	1818	ENDDO
	1819	ns = ns - 1
	1820	ENDIF
	1821
	1822	!
	1823	!-- Define x axis (for scalar position)
	1824	nc_stat = NF90_DEF_DIM( id_set_yz(av), 'x_yz', ns, id_dim_x_yz(av) )
	1825	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 184 )
	1826
	1827	nc_stat = NF90_DEF_VAR( id_set_yz(av), 'x_yz', NF90_DOUBLE, &
	1828	id_dim_x_yz(av), id_var_x_yz(av) )
	1829	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 185 )
	1830
	1831	nc_stat = NF90_PUT_ATT( id_set_yz(av), id_var_x_yz(av), 'units', &
	1832	'meters' )
	1833	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 186 )
	1834
	1835	!
	1836	!-- Define x axis (for u position)
	1837	nc_stat = NF90_DEF_DIM( id_set_yz(av), 'xu_yz', ns, id_dim_xu_yz(av) )
	1838	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 377 )
	1839
	1840	nc_stat = NF90_DEF_VAR( id_set_yz(av), 'xu_yz', NF90_DOUBLE, &
	1841	id_dim_xu_yz(av), id_var_xu_yz(av) )
	1842	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 378 )
	1843
	1844	nc_stat = NF90_PUT_ATT( id_set_yz(av), id_var_xu_yz(av), 'units', &
	1845	'meters' )
	1846	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 379 )
	1847
	1848	!
	1849	!-- Define a variable to store the layer indices of the vertical cross
	1850	!-- sections
	1851	nc_stat = NF90_DEF_VAR( id_set_yz(av), 'ind_x_yz', NF90_DOUBLE, &
	1852	id_dim_x_yz(av), id_var_ind_x_yz(av) )
	1853	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 187 )
	1854
	1855	nc_stat = NF90_PUT_ATT( id_set_yz(av), id_var_ind_x_yz(av), 'units', &
	1856	'gridpoints')
	1857	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 188 )
	1858
	1859	!
	1860	!-- Define y-axis (for scalar position)
	1861	nc_stat = NF90_DEF_DIM( id_set_yz(av), 'y', ny+2, id_dim_y_yz(av) )
	1862	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 189 )
	1863
	1864	nc_stat = NF90_DEF_VAR( id_set_yz(av), 'y', NF90_DOUBLE, &
	1865	id_dim_y_yz(av), id_var_y_yz(av) )
	1866	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 190 )
	1867
	1868	nc_stat = NF90_PUT_ATT( id_set_yz(av), id_var_y_yz(av), 'units', &
	1869	'meters' )
	1870	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 191 )
	1871
	1872	!
	1873	!-- Define y-axis (for v position)
	1874	nc_stat = NF90_DEF_DIM( id_set_yz(av), 'yv', ny+2, id_dim_yv_yz(av) )
	1875	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 380 )
	1876
	1877	nc_stat = NF90_DEF_VAR( id_set_yz(av), 'yv', NF90_DOUBLE, &
	1878	id_dim_yv_yz(av), id_var_yv_yz(av) )
	1879	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 381 )
	1880
	1881	nc_stat = NF90_PUT_ATT( id_set_yz(av), id_var_yv_yz(av), 'units', &
	1882	'meters' )
	1883	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 382 )
	1884
	1885	!
	1886	!-- Define the two z-axes (zu and zw)
	1887	nc_stat = NF90_DEF_DIM( id_set_yz(av), 'zu', nz+2, id_dim_zu_yz(av) )
	1888	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 192 )
	1889
	1890	nc_stat = NF90_DEF_VAR( id_set_yz(av), 'zu', NF90_DOUBLE, &
	1891	id_dim_zu_yz(av), id_var_zu_yz(av) )
	1892	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 193 )
	1893
	1894	nc_stat = NF90_PUT_ATT( id_set_yz(av), id_var_zu_yz(av), 'units', &
	1895	'meters' )
	1896	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 194 )
	1897
	1898	nc_stat = NF90_DEF_DIM( id_set_yz(av), 'zw', nz+2, id_dim_zw_yz(av) )
	1899	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 195 )
	1900
	1901	nc_stat = NF90_DEF_VAR( id_set_yz(av), 'zw', NF90_DOUBLE, &
	1902	id_dim_zw_yz(av), id_var_zw_yz(av) )
	1903	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 196 )
	1904
	1905	nc_stat = NF90_PUT_ATT( id_set_yz(av), id_var_zw_yz(av), 'units', &
	1906	'meters' )
	1907	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 197 )
	1908
	1909
	1910	!
	1911	!-- Define the variables
	1912	var_list = ';'
	1913	i = 1
	1914
	1915	DO WHILE ( do2d(av,i)(1:1) /= ' ' )
	1916
	1917	IF ( INDEX( do2d(av,i), 'yz' ) /= 0 ) THEN
	1918
	1919	!
	1920	!-- Check for the grid
	1921	found = .TRUE.
	1922	SELECT CASE ( do2d(av,i) )
	1923	!
	1924	!-- Most variables are defined on the zu grid
	1925	CASE ( 'e_yz', 'p_yz', 'pc_yz', 'pr_yz', 'pt_yz', 'q_yz', &
	1926	'ql_yz', 'ql_c_yz', 'ql_v_yz', 'ql_vp_yz', 'qv_yz', &
	1927	's_yz', 'vpt_yz' )
	1928
	1929	grid_x = 'x'
	1930	grid_y = 'y'
	1931	grid_z = 'zu'
	1932	!
	1933	!-- u grid
	1934	CASE ( 'u_yz' )
	1935
	1936	grid_x = 'xu'
	1937	grid_y = 'y'
	1938	grid_z = 'zu'
	1939	!
	1940	!-- v grid
	1941	CASE ( 'v_yz' )
	1942
	1943	grid_x = 'x'
	1944	grid_y = 'yv'
	1945	grid_z = 'zu'
	1946	!
	1947	!-- w grid
	1948	CASE ( 'w_yz' )
	1949
	1950	grid_x = 'x'
	1951	grid_y = 'y'
	1952	grid_z = 'zw'
	1953
	1954	CASE DEFAULT
	1955	!
	1956	!-- Check for user-defined quantities
	1957	CALL user_define_netcdf_grid( do2d(av,i), found, &
	1958	grid_x, grid_y, grid_z )
	1959
	1960	END SELECT
	1961
	1962	!
	1963	!-- Select the respective dimension ids
	1964	IF ( grid_x == 'x' ) THEN
	1965	id_x = id_dim_x_yz(av)
	1966	ELSEIF ( grid_x == 'xu' ) THEN
	1967	id_x = id_dim_xu_yz(av)
	1968	ENDIF
	1969
	1970	IF ( grid_y == 'y' ) THEN
	1971	id_y = id_dim_y_yz(av)
	1972	ELSEIF ( grid_y == 'yv' ) THEN
	1973	id_y = id_dim_yv_yz(av)
	1974	ENDIF
	1975
	1976	IF ( grid_z == 'zu' ) THEN
	1977	id_z = id_dim_zu_yz(av)
	1978	ELSEIF ( grid_z == 'zw' ) THEN
	1979	id_z = id_dim_zw_yz(av)
	1980	ENDIF
	1981
	1982	!
	1983	!-- Define the grid
	1984	nc_stat = NF90_DEF_VAR( id_set_yz(av), do2d(av,i), &
	1985	nc_precision(3), &
	1986	(/ id_x, id_y, id_z, id_dim_time_yz(av) /), &
	1987	id_var_do2d(av,i) )
	1988
	1989	IF ( .NOT. found ) THEN
	1990	PRINT*, '+++ define_netcdf_header: no grid defined for', &
	1991	' variable ', do2d(av,i)
	1992	ENDIF
	1993
	1994	var_list = TRIM( var_list ) // TRIM( do2d(av,i) ) // ';'
	1995
	1996	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 198 )
	1997	!
	1998	!-- Store the 'real' name of the variable (with *, for example)
	1999	!-- in the long_name attribute. This is evaluated by Ferret,
	2000	!-- for example.
	2001	nc_stat = NF90_PUT_ATT( id_set_yz(av), id_var_do2d(av,i), &
	2002	'long_name', do2d(av,i) )
	2003	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 199 )
	2004	!
	2005	!-- Define the variable's unit
	2006	nc_stat = NF90_PUT_ATT( id_set_yz(av), id_var_do2d(av,i), &
	2007	'units', TRIM( do2d_unit(av,i) ) )
	2008	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 356 )
	2009	ENDIF
	2010
	2011	i = i + 1
	2012
	2013	ENDDO
	2014
	2015	!
	2016	!-- No arrays to output. Close the netcdf file and return.
	2017	IF ( i == 1 ) RETURN
	2018
	2019	!
	2020	!-- Write the list of variables as global attribute (this is used by
	2021	!-- restart runs and by combine_plot_fields)
	2022	nc_stat = NF90_PUT_ATT( id_set_yz(av), NF90_GLOBAL, 'VAR_LIST', &
	2023	var_list )
	2024	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 200 )
	2025
	2026	!
	2027	!-- Leave NetCDF define mode
	2028	nc_stat = NF90_ENDDEF( id_set_yz(av) )
	2029	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 201 )
	2030
	2031	!
	2032	!-- Write axis data: x_yz, y, zu, zw
	2033	ALLOCATE( netcdf_data(1:ns) )
	2034
	2035	!
	2036	!-- Write x_yz data
	2037	DO i = 1, ns
	2038	IF( section(i,3) == -1 ) THEN
	2039	netcdf_data(i) = -1.0 ! section averaged along x
	2040	ELSE
	2041	netcdf_data(i) = section(i,3) * dx
	2042	ENDIF
	2043	ENDDO
	2044	nc_stat = NF90_PUT_VAR( id_set_yz(av), id_var_x_yz(av), netcdf_data, &
	2045	start = (/ 1 /), count = (/ ns /) )
	2046	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 202 )
	2047
	2048	!
	2049	!-- Write x_yz data (xu grid)
	2050	DO i = 1, ns
	2051	IF( section(i,3) == -1 ) THEN
	2052	netcdf_data(i) = -1.0 ! section averaged along x
	2053	ELSE
	2054	netcdf_data(i) = (section(i,3)-0.5) * dx
	2055	ENDIF
	2056	ENDDO
	2057	nc_stat = NF90_PUT_VAR( id_set_yz(av), id_var_xu_yz(av), netcdf_data, &
	2058	start = (/ 1 /), count = (/ ns /) )
	2059	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 383 )
	2060
	2061	!
	2062	!-- Write gridpoint number data
	2063	netcdf_data(1:ns) = section(1:ns,3)
	2064	nc_stat = NF90_PUT_VAR( id_set_yz(av), id_var_ind_x_yz(av), &
	2065	netcdf_data, start = (/ 1 /), &
	2066	count = (/ ns /) )
	2067	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 203 )
	2068
	2069	DEALLOCATE( netcdf_data )
	2070
	2071	!
	2072	!-- Write data for y and yv axis (shifted by -dy/2)
	2073	ALLOCATE( netcdf_data(0:ny+1) )
	2074
	2075	DO j = 0, ny+1
	2076	netcdf_data(j) = j * dy
	2077	ENDDO
	2078
	2079	nc_stat = NF90_PUT_VAR( id_set_yz(av), id_var_y_yz(av), netcdf_data, &
	2080	start = (/ 1 /), count = (/ ny+2 /) )
	2081	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 204 )
	2082
	2083	DO j = 0, ny+1
	2084	netcdf_data(j) = ( j - 0.5 ) * dy
	2085	ENDDO
	2086
	2087	nc_stat = NF90_PUT_VAR( id_set_yz(av), id_var_yv_yz(av), &
	2088	netcdf_data, start = (/ 1 /), &
	2089	count = (/ ny+2 /) )
	2090	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 384 )
	2091
	2092	DEALLOCATE( netcdf_data )
	2093
	2094	!
	2095	!-- Write zu and zw data (vertical axes)
	2096	ALLOCATE( netcdf_data(0:nz+1) )
	2097
	2098	netcdf_data(0:nz+1) = zu(nzb:nzt+1)
	2099	nc_stat = NF90_PUT_VAR( id_set_yz(av), id_var_zu_yz(av), &
	2100	netcdf_data, start = (/ 1 /), &
	2101	count = (/ nz+2 /) )
	2102	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 205 )
	2103
	2104	netcdf_data(0:nz+1) = zw(nzb:nzt+1)
	2105	nc_stat = NF90_PUT_VAR( id_set_yz(av), id_var_zw_yz(av), &
	2106	netcdf_data, start = (/ 1 /), &
	2107	count = (/ nz+2 /) )
	2108	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 206 )
	2109
	2110	DEALLOCATE( netcdf_data )
	2111
	2112
	2113	CASE ( 'yz_ext' )
	2114
	2115	!
	2116	!-- Get the list of variables and compare with the actual run.
	2117	!-- First var_list_old has to be reset, since GET_ATT does not assign
	2118	!-- trailing blanks.
	2119	var_list_old = ' '
	2120	nc_stat = NF90_GET_ATT( id_set_yz(av), NF90_GLOBAL, 'VAR_LIST', &
	2121	var_list_old )
	2122	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 207 )
	2123
	2124	var_list = ';'
	2125	i = 1
	2126	DO WHILE ( do2d(av,i)(1:1) /= ' ' )
	2127	IF ( INDEX( do2d(av,i), 'yz' ) /= 0 ) THEN
	2128	netcdf_var_name = do2d(av,i)
	2129	CALL clean_netcdf_varname( netcdf_var_name )
	2130	var_list = TRIM(var_list) // TRIM(netcdf_var_name) // ';'
	2131	ENDIF
	2132	i = i + 1
	2133	ENDDO
	2134
	2135	IF ( av == 0 ) THEN
	2136	var = '(yz)'
	2137	ELSE
	2138	var = '(yz_av)'
	2139	ENDIF
	2140
	2141	IF ( TRIM( var_list ) /= TRIM( var_list_old ) ) THEN
	2142	PRINT*, '+++ WARNING: NetCDF file for cross-sections ' // &
	2143	TRIM( var ) // ' from previuos run found,'
	2144	PRINT*, ' but this file cannot be extended due to' // &
	2145	' variable mismatch.'
	2146	PRINT*, ' New file is created instead.'
	2147	extend = .FALSE.
	2148	RETURN
	2149	ENDIF
	2150
	2151	!
	2152	!-- Calculate the number of current sections
	2153	ns = 1
	2154	DO WHILE ( section(ns,3) /= -9999 .AND. ns <= 100 )
	2155	ns = ns + 1
	2156	ENDDO
	2157	ns = ns - 1
	2158
	2159	!
	2160	!-- Get and compare the number of vertical cross sections
	2161	nc_stat = NF90_INQ_VARID( id_set_yz(av), 'x_yz', id_var_x_yz(av) )
	2162	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 208 )
	2163
	2164	nc_stat = NF90_INQUIRE_VARIABLE( id_set_yz(av), id_var_x_yz(av), &
	2165	dimids = id_dim_x_yz_old )
	2166	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 209 )
	2167	id_dim_x_yz(av) = id_dim_x_yz_old(1)
	2168
	2169	nc_stat = NF90_INQUIRE_DIMENSION( id_set_yz(av), id_dim_x_yz(av), &
	2170	len = ns_old )
	2171	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 210 )
	2172
	2173	IF ( ns /= ns_old ) THEN
	2174	PRINT*, '+++ WARNING: NetCDF file for cross-sections ' // &
	2175	TRIM( var ) // ' from previuos run found,'
	2176	PRINT*, ' but this file cannot be extended due to' // &
	2177	' mismatch in number of'
	2178	PRINT*, ' cross sections.'
	2179	PRINT*, ' New file is created instead.'
	2180	extend = .FALSE.
	2181	RETURN
	2182	ENDIF
	2183
	2184	!
	2185	!-- Get and compare the heights of the cross sections
	2186	ALLOCATE( netcdf_data(1:ns_old) )
	2187
	2188	nc_stat = NF90_GET_VAR( id_set_yz(av), id_var_x_yz(av), netcdf_data )
	2189	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 211 )
	2190
	2191	DO i = 1, ns
	2192	IF ( section(i,3) /= -1 ) THEN
	2193	IF ( ( section(i,3) * dx ) /= netcdf_data(i) ) THEN
	2194	PRINT*, '+++ WARNING: NetCDF file for cross-sections ' // &
	2195	TRIM( var ) // ' from previuos run found,'
	2196	PRINT*, ' but this file cannot be extended' // &
	2197	' due to mismatch in cross'
	2198	PRINT*, ' section indices.'
	2199	PRINT*, ' New file is created instead.'
	2200	extend = .FALSE.
	2201	RETURN
	2202	ENDIF
	2203	ELSE
	2204	IF ( -1.0 /= netcdf_data(i) ) THEN
	2205	PRINT*, '+++ WARNING: NetCDF file for cross-sections ' // &
	2206	TRIM( var ) // ' from previuos run found,'
	2207	PRINT*, ' but this file cannot be extended' // &
	2208	' due to mismatch in cross'
	2209	PRINT*, ' section indices.'
	2210	PRINT*, ' New file is created instead.'
	2211	extend = .FALSE.
	2212	RETURN
	2213	ENDIF
	2214	ENDIF
	2215	ENDDO
	2216
	2217	DEALLOCATE( netcdf_data )
	2218
	2219	!
	2220	!-- Get the id of the time coordinate (unlimited coordinate) and its
	2221	!-- last index on the file. The next time level is pl2d..count+1.
	2222	!-- The current time must be larger than the last output time
	2223	!-- on the file.
	2224	nc_stat = NF90_INQ_VARID( id_set_yz(av), 'time', id_var_time_yz(av) )
	2225	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 212 )
	2226
	2227	nc_stat = NF90_INQUIRE_VARIABLE( id_set_yz(av), id_var_time_yz(av), &
	2228	dimids = id_dim_time_old )
	2229	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 213 )
	2230	id_dim_time_yz(av) = id_dim_time_old(1)
	2231
	2232	nc_stat = NF90_INQUIRE_DIMENSION( id_set_yz(av), id_dim_time_yz(av), &
	2233	len = do2d_yz_time_count(av) )
	2234	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 214 )
	2235
	2236	nc_stat = NF90_GET_VAR( id_set_yz(av), id_var_time_yz(av), &
	2237	last_time_coordinate, &
	2238	start = (/ do2d_yz_time_count(av) /), &
	2239	count = (/ 1 /) )
	2240	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 215 )
	2241
	2242	IF ( last_time_coordinate(1) >= simulated_time ) THEN
	2243	PRINT*, '+++ WARNING: NetCDF file for cross sections ' // &
	2244	TRIM( var ) // ' from previuos run found,'
	2245	PRINT*, ' but this file cannot be extended becaus' // &
	2246	'e the current output time'
	2247	PRINT*, ' is less or equal than the last output t' // &
	2248	'ime on this file.'
	2249	PRINT*, ' New file is created instead.'
	2250	do2d_yz_time_count(av) = 0
	2251	extend = .FALSE.
	2252	RETURN
	2253	ENDIF
	2254
	2255	!
	2256	!-- Dataset seems to be extendable.
	2257	!-- Now get the variable ids.
	2258	i = 1
	2259	DO WHILE ( do2d(av,i)(1:1) /= ' ' )
	2260	IF ( INDEX( do2d(av,i), 'yz' ) /= 0 ) THEN
	2261	netcdf_var_name = do2d(av,i)
	2262	CALL clean_netcdf_varname( netcdf_var_name )
	2263	nc_stat = NF90_INQ_VARID( id_set_yz(av), netcdf_var_name, &
	2264	id_var_do2d(av,i) )
	2265	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 216 )
	2266	ENDIF
	2267	i = i + 1
	2268	ENDDO
	2269
	2270	!
	2271	!-- Change the titel attribute on file
	2272	nc_stat = NF90_PUT_ATT( id_set_yz(av), NF90_GLOBAL, 'title', &
	2273	TRIM( run_description_header ) )
	2274	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 217 )
	2275	PRINT, '** NetCDF file for cross-sections ' // TRIM( var ) // &
	2276	' from previous run found.'
	2277	PRINT*, ' This file will be extended.'
	2278
	2279
	2280	CASE ( 'pr_new' )
	2281
	2282	!
	2283	!-- Define some global attributes of the dataset
	2284	IF ( averaging_interval_pr /= 0.0 ) THEN
	2285	WRITE (time_average_text,'('', '',F7.1,'' s average'')') &
	2286	averaging_interval_pr
	2287	nc_stat = NF90_PUT_ATT( id_set_pr, NF90_GLOBAL, 'title', &
	2288	TRIM( run_description_header ) // &
	2289	TRIM( time_average_text ) )
	2290	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 218 )
	2291
	2292	WRITE ( time_average_text,'(F7.1,'' s avg'')' ) averaging_interval_pr
	2293	nc_stat = NF90_PUT_ATT( id_set_pr, NF90_GLOBAL, 'time_avg', &
	2294	TRIM( time_average_text ) )
	2295	ELSE
	2296	nc_stat = NF90_PUT_ATT( id_set_pr, NF90_GLOBAL, 'title', &
	2297	TRIM( run_description_header ) )
	2298	ENDIF
	2299	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 219 )
	2300
	2301	!
	2302	!-- Define time coordinate for profiles (unlimited dimension)
	2303	nc_stat = NF90_DEF_DIM( id_set_pr, 'time', NF90_UNLIMITED, &
	2304	id_dim_time_pr )
	2305	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 220 )
	2306
	2307	nc_stat = NF90_DEF_VAR( id_set_pr, 'time', NF90_DOUBLE, &
	2308	id_dim_time_pr, id_var_time_pr )
	2309	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 221 )
	2310
	2311	nc_stat = NF90_PUT_ATT( id_set_pr, id_var_time_pr, 'units', 'seconds')
	2312	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 222 )
	2313
	2314	!
	2315	!-- Define the variables
	2316	var_list = ';'
	2317	DO i = 1, dopr_n
	2318	!
	2319	!-- First, remove those characters not allowed by NetCDF
	2320	netcdf_var_name = data_output_pr(i)
	2321	CALL clean_netcdf_varname( netcdf_var_name )
	2322
	2323	IF ( statistic_regions == 0 ) THEN
	2324
	2325	!
	2326	!-- Define the z-axes (each variable gets its own z-axis)
	2327	nc_stat = NF90_DEF_DIM( id_set_pr, 'z'//TRIM(netcdf_var_name), &
	2328	nzt+2-nzb, id_dim_z_pr(i,0) )
	2329	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 223 )
	2330
	2331	nc_stat = NF90_DEF_VAR( id_set_pr, 'z'//TRIM(netcdf_var_name), &
	2332	NF90_DOUBLE, id_dim_z_pr(i,0), &
	2333	id_var_z_pr(i,0) )
	2334	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 224 )
	2335
	2336	nc_stat = NF90_PUT_ATT( id_set_pr, id_var_z_pr(i,0), 'units', &
	2337	'meters' )
	2338	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 225 )
	2339
	2340	!
	2341	!-- Define the variable
	2342	nc_stat = NF90_DEF_VAR( id_set_pr, netcdf_var_name, &
	2343	nc_precision(5), (/ id_dim_z_pr(i,0), &
	2344	id_dim_time_pr /), id_var_dopr(i,0) )
	2345	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 226 )
	2346
	2347	nc_stat = NF90_PUT_ATT( id_set_pr, id_var_dopr(i,0), &
	2348	'long_name', TRIM( data_output_pr(i) ) )
	2349	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 227 )
	2350
	2351	nc_stat = NF90_PUT_ATT( id_set_pr, id_var_dopr(i,0), &
	2352	'units', TRIM( dopr_unit(i) ) )
	2353	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 228 )
	2354
	2355	var_list = TRIM(var_list) // TRIM(netcdf_var_name) // ';'
	2356
	2357	ELSE
	2358	!
	2359	!-- If statistic regions are defined, add suffix _SR+#SR to the
	2360	!-- names
	2361	DO j = 0, statistic_regions
	2362	WRITE ( suffix, '(''_'',I1)' ) j
	2363
	2364	!
	2365	!-- Define the z-axes (each variable gets it own z-axis)
	2366	nc_stat = NF90_DEF_DIM( id_set_pr, &
	2367	'z'//TRIM(netcdf_var_name)//suffix, &
	2368	nzt+2-nzb, id_dim_z_pr(i,j) )
	2369	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 229 )
	2370
	2371	nc_stat = NF90_DEF_VAR( id_set_pr, &
	2372	'z'//TRIM(netcdf_var_name)//suffix, &
	2373	nc_precision(5), id_dim_z_pr(i,j), &
	2374	id_var_z_pr(i,j) )
	2375	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 230 )
	2376
	2377	nc_stat = NF90_PUT_ATT( id_set_pr, id_var_z_pr(i,j), &
	2378	'units', 'meters' )
	2379	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 231 )
	2380
	2381	!
	2382	!-- Define the variable
	2383	nc_stat = NF90_DEF_VAR( id_set_pr, &
	2384	TRIM( netcdf_var_name ) // suffix, &
	2385	nc_precision(5), &
	2386	(/ id_dim_z_pr(i,j), &
	2387	id_dim_time_pr /), id_var_dopr(i,j) )
	2388	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 232 )
	2389
	2390	nc_stat = NF90_PUT_ATT( id_set_pr, id_var_dopr(i,j), &
	2391	'long_name', &
	2392	TRIM( data_output_pr(i) ) // ' SR ' &
	2393	// suffix(2:2) )
	2394	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 233 )
	2395
	2396	nc_stat = NF90_PUT_ATT( id_set_pr, id_var_dopr(i,j), &
	2397	'units', TRIM( dopr_unit(i) ) )
	2398	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 234 )
	2399
	2400	var_list = TRIM(var_list) // TRIM(netcdf_var_name) // &
	2401	suffix // ';'
	2402
	2403	ENDDO
	2404
	2405	ENDIF
	2406
	2407	ENDDO
	2408
	2409	!
	2410	!-- Write the list of variables as global attribute (this is used by
	2411	!-- restart runs)
	2412	nc_stat = NF90_PUT_ATT( id_set_pr, NF90_GLOBAL, 'VAR_LIST', var_list )
	2413	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 235 )
	2414
	2415	!
	2416	!-- Define normalization variables (as time series)
	2417	DO i = 1, dopr_norm_num
	2418
	2419	nc_stat = NF90_DEF_VAR( id_set_pr, 'NORM_' // &
	2420	TRIM( dopr_norm_names(i) ), &
	2421	nc_precision(5), (/ id_dim_time_pr /), &
	2422	id_var_norm_dopr(i) )
	2423	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 236 )
	2424
	2425	nc_stat = NF90_PUT_ATT( id_set_pr, id_var_norm_dopr(i), &
	2426	'long_name', &
	2427	TRIM( dopr_norm_longnames(i) ) )
	2428	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 237 )
	2429
	2430	ENDDO
	2431
	2432	!
	2433	!-- Leave NetCDF define mode
	2434	nc_stat = NF90_ENDDEF( id_set_pr )
	2435	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 238 )
	2436
	2437	!
	2438	!-- Write z-axes data
	2439	DO i = 1, dopr_n
	2440	DO j = 0, statistic_regions
	2441
	2442	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_z_pr(i,j), &
	2443	hom(nzb:nzt+1,2,dopr_index(i),0), &
	2444	start = (/ 1 /), &
	2445	count = (/ nzt-nzb+2 /) )
	2446	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 239 )
	2447
	2448	ENDDO
	2449	ENDDO
	2450
	2451
	2452	CASE ( 'pr_ext' )
	2453
	2454	!
	2455	!-- Get the list of variables and compare with the actual run.
	2456	!-- First var_list_old has to be reset, since GET_ATT does not assign
	2457	!-- trailing blanks.
	2458	var_list_old = ' '
	2459	nc_stat = NF90_GET_ATT( id_set_pr, NF90_GLOBAL, 'VAR_LIST', &
	2460	var_list_old )
	2461	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 240 )
	2462
	2463	var_list = ';'
	2464	DO i = 1, dopr_n
	2465
	2466	netcdf_var_name = data_output_pr(i)
	2467	CALL clean_netcdf_varname( netcdf_var_name )
	2468
	2469	IF ( statistic_regions == 0 ) THEN
	2470	var_list = TRIM(var_list) // TRIM(netcdf_var_name) // ';'
	2471	ELSE
	2472	DO j = 0, statistic_regions
	2473	WRITE ( suffix, '(''_'',I1)' ) j
	2474	var_list = TRIM(var_list) // TRIM(netcdf_var_name) // &
	2475	suffix // ';'
	2476	ENDDO
	2477	ENDIF
	2478
	2479	ENDDO
	2480
	2481	IF ( TRIM( var_list ) /= TRIM( var_list_old ) ) THEN
	2482	PRINT*, '+++ WARNING: NetCDF file for vertical profiles from' // &
	2483	' previuos run found,'
	2484	PRINT*, ' but this file cannot be extended due to' // &
	2485	' variable mismatch.'
	2486	PRINT*, ' New file is created instead.'
	2487	extend = .FALSE.
	2488	RETURN
	2489	ENDIF
	2490
	2491	!
	2492	!-- Get the id of the time coordinate (unlimited coordinate) and its
	2493	!-- last index on the file. The next time level is dopr..count+1.
	2494	!-- The current time must be larger than the last output time
	2495	!-- on the file.
	2496	nc_stat = NF90_INQ_VARID( id_set_pr, 'time', id_var_time_pr )
	2497	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 241 )
	2498
	2499	nc_stat = NF90_INQUIRE_VARIABLE( id_set_pr, id_var_time_pr, &
	2500	dimids = id_dim_time_old )
	2501	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 242 )
	2502	id_dim_time_pr = id_dim_time_old(1)
	2503
	2504	nc_stat = NF90_INQUIRE_DIMENSION( id_set_pr, id_dim_time_pr, &
	2505	len = dopr_time_count )
	2506	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 243 )
	2507
	2508	nc_stat = NF90_GET_VAR( id_set_pr, id_var_time_pr, &
	2509	last_time_coordinate, &
	2510	start = (/ dopr_time_count /), &
	2511	count = (/ 1 /) )
	2512	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 244 )
	2513
	2514	IF ( last_time_coordinate(1) >= simulated_time ) THEN
	2515	PRINT*, '+++ WARNING: NetCDF file for vertical profiles from' // &
	2516	' previuos run found,'
	2517	PRINT*, ' but this file cannot be extended becaus' // &
	2518	'e the current output time'
	2519	PRINT*, ' is less or equal than the last output t' // &
	2520	'ime on this file.'
	2521	PRINT*, ' New file is created instead.'
	2522	dopr_time_count = 0
	2523	extend = .FALSE.
	2524	RETURN
	2525	ENDIF
	2526
	2527	!
	2528	!-- Dataset seems to be extendable.
	2529	!-- Now get the variable ids.
	2530	i = 1
	2531	DO i = 1, dopr_n
	2532
	2533	netcdf_var_name_base = data_output_pr(i)
	2534	CALL clean_netcdf_varname( netcdf_var_name_base )
	2535
	2536	IF ( statistic_regions == 0 ) THEN
	2537	nc_stat = NF90_INQ_VARID( id_set_pr, netcdf_var_name_base, &
	2538	id_var_dopr(i,0) )
	2539	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 245 )
	2540	ELSE
	2541	DO j = 0, statistic_regions
	2542	WRITE ( suffix, '(''_'',I1)' ) j
	2543	netcdf_var_name = TRIM( netcdf_var_name_base ) // suffix
	2544	nc_stat = NF90_INQ_VARID( id_set_pr, netcdf_var_name, &
	2545	id_var_dopr(i,j) )
	2546	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 246 )
	2547	ENDDO
	2548	ENDIF
	2549
	2550	ENDDO
	2551
	2552	!
	2553	!-- Get ids of the normalization variables
	2554	DO i = 1, dopr_norm_num
	2555	nc_stat = NF90_INQ_VARID( id_set_pr, &
	2556	'NORM_' // TRIM( dopr_norm_names(i) ), &
	2557	id_var_norm_dopr(i) )
	2558	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 247 )
	2559	ENDDO
	2560
	2561	!
	2562	!-- Change the title attribute on file
	2563	nc_stat = NF90_PUT_ATT( id_set_pr, NF90_GLOBAL, 'title', &
	2564	TRIM( run_description_header ) )
	2565	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 248 )
	2566	PRINT, '** NetCDF file for vertical profiles from previous run' // &
	2567	' found.'
	2568	PRINT*, ' This file will be extended.'
	2569
	2570
	2571	CASE ( 'ts_new' )
	2572
	2573	!
	2574	!-- Define some global attributes of the dataset
	2575	nc_stat = NF90_PUT_ATT( id_set_ts, NF90_GLOBAL, 'title', &
	2576	TRIM( run_description_header ) )
	2577	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 249 )
	2578
	2579	!
	2580	!-- Define time coordinate for time series (unlimited dimension)
	2581	nc_stat = NF90_DEF_DIM( id_set_ts, 'time', NF90_UNLIMITED, &
	2582	id_dim_time_ts )
	2583	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 250 )
	2584
	2585	nc_stat = NF90_DEF_VAR( id_set_ts, 'time', NF90_DOUBLE, &
	2586	id_dim_time_ts, id_var_time_ts )
	2587	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 251 )
	2588
	2589	nc_stat = NF90_PUT_ATT( id_set_ts, id_var_time_ts, 'units', 'seconds')
	2590	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 252 )
	2591
	2592	!
	2593	!-- Define the variables
	2594	var_list = ';'
	2595	DO i = 1, dots_num
	2596	!
	2597	!-- First, remove those characters not allowed by NetCDF
	2598	netcdf_var_name = dots_label(i)
	2599	CALL clean_netcdf_varname( netcdf_var_name )
	2600
	2601	IF ( statistic_regions == 0 ) THEN
	2602
	2603	nc_stat = NF90_DEF_VAR( id_set_ts, netcdf_var_name, &
	2604	nc_precision(6), (/ id_dim_time_ts /), &
	2605	id_var_dots(i,0) )
	2606	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 253 )
	2607
	2608	nc_stat = NF90_PUT_ATT( id_set_ts, id_var_dots(i,0), &
	2609	'long_name', TRIM( dots_label(i) ) )
	2610	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 254 )
	2611
	2612	nc_stat = NF90_PUT_ATT( id_set_ts, id_var_dots(i,0), &
	2613	'units', TRIM( dots_unit(i) ) )
	2614	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 255 )
	2615
	2616	var_list = TRIM(var_list) // TRIM(netcdf_var_name) // ';'
	2617
	2618	ELSE
	2619	!
	2620	!-- If statistic regions are defined, add suffix _SR+#SR to the
	2621	!-- names
	2622	DO j = 0, statistic_regions
	2623	WRITE ( suffix, '(''_'',I1)' ) j
	2624
	2625	nc_stat = NF90_DEF_VAR( id_set_ts, &
	2626	TRIM( netcdf_var_name ) // suffix, &
	2627	nc_precision(6), &
	2628	(/ id_dim_time_ts /), &
	2629	id_var_dots(i,j) )
	2630	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 256 )
	2631
	2632	nc_stat = NF90_PUT_ATT( id_set_ts, id_var_dots(i,j), &
	2633	'long_name', &
	2634	TRIM( dots_label(i) ) // ' SR ' // &
	2635	suffix(2:2) )
	2636	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 257 )
	2637
	2638	var_list = TRIM(var_list) // TRIM(netcdf_var_name) // &
	2639	suffix // ';'
	2640
	2641	ENDDO
	2642
	2643	ENDIF
	2644
	2645	ENDDO
	2646
	2647	!
	2648	!-- Write the list of variables as global attribute (this is used by
	2649	!-- restart runs)
	2650	nc_stat = NF90_PUT_ATT( id_set_ts, NF90_GLOBAL, 'VAR_LIST', var_list )
	2651	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 258 )
	2652
	2653	!
	2654	!-- Leave NetCDF define mode
	2655	nc_stat = NF90_ENDDEF( id_set_ts )
	2656	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 259 )
	2657
	2658
	2659	CASE ( 'ts_ext' )
	2660
	2661	!
	2662	!-- Get the list of variables and compare with the actual run.
	2663	!-- First var_list_old has to be reset, since GET_ATT does not assign
	2664	!-- trailing blanks.
	2665	var_list_old = ' '
	2666	nc_stat = NF90_GET_ATT( id_set_ts, NF90_GLOBAL, 'VAR_LIST', &
	2667	var_list_old )
	2668	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 260 )
	2669
	2670	var_list = ';'
	2671	i = 1
	2672	DO i = 1, dots_num
	2673
	2674	netcdf_var_name = dots_label(i)
	2675	CALL clean_netcdf_varname( netcdf_var_name )
	2676
	2677	IF ( statistic_regions == 0 ) THEN
	2678	var_list = TRIM(var_list) // TRIM(netcdf_var_name) // ';'
	2679	ELSE
	2680	DO j = 0, statistic_regions
	2681	WRITE ( suffix, '(''_'',I1)' ) j
	2682	var_list = TRIM(var_list) // TRIM(netcdf_var_name) // &
	2683	suffix // ';'
	2684	ENDDO
	2685	ENDIF
	2686
	2687	ENDDO
	2688
	2689	IF ( TRIM( var_list ) /= TRIM( var_list_old ) ) THEN
	2690	PRINT*, '+++ WARNING: NetCDF file for time series from' //&
	2691	' previuos run found,'
	2692	PRINT*, ' but this file cannot be extended due to' // &
	2693	' variable mismatch.'
	2694	PRINT*, ' New file is created instead.'
	2695	extend = .FALSE.
	2696	RETURN
	2697	ENDIF
	2698
	2699	!
	2700	!-- Get the id of the time coordinate (unlimited coordinate) and its
	2701	!-- last index on the file. The next time level is dots..count+1.
	2702	!-- The current time must be larger than the last output time
	2703	!-- on the file.
	2704	nc_stat = NF90_INQ_VARID( id_set_ts, 'time', id_var_time_ts )
	2705	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 261 )
	2706
	2707	nc_stat = NF90_INQUIRE_VARIABLE( id_set_ts, id_var_time_ts, &
	2708	dimids = id_dim_time_old )
	2709	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 262 )
	2710	id_dim_time_ts = id_dim_time_old(1)
	2711
	2712	nc_stat = NF90_INQUIRE_DIMENSION( id_set_ts, id_dim_time_ts, &
	2713	len = dots_time_count )
	2714	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 263 )
	2715
	2716	nc_stat = NF90_GET_VAR( id_set_ts, id_var_time_ts, &
	2717	last_time_coordinate, &
	2718	start = (/ dots_time_count /), &
	2719	count = (/ 1 /) )
	2720	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 264 )
	2721
	2722	IF ( last_time_coordinate(1) >= simulated_time ) THEN
	2723	PRINT*, '+++ WARNING: NetCDF file for time series from' // &
	2724	' previuos run found,'
	2725	PRINT*, ' but this file cannot be extended becaus' // &
	2726	'e the current output time'
	2727	PRINT*, ' is less or equal than the last output t' // &
	2728	'ime on this file.'
	2729	PRINT*, ' New file is created instead.'
	2730	dots_time_count = 0
	2731	extend = .FALSE.
	2732	RETURN
	2733	ENDIF
	2734
	2735	!
	2736	!-- Dataset seems to be extendable.
	2737	!-- Now get the variable ids
	2738	i = 1
	2739	DO i = 1, dots_num
	2740
	2741	netcdf_var_name_base = dots_label(i)
	2742	CALL clean_netcdf_varname( netcdf_var_name_base )
	2743
	2744	IF ( statistic_regions == 0 ) THEN
	2745	nc_stat = NF90_INQ_VARID( id_set_ts, netcdf_var_name_base, &
	2746	id_var_dots(i,0) )
	2747	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 265 )
	2748	ELSE
	2749	DO j = 0, statistic_regions
	2750	WRITE ( suffix, '(''_'',I1)' ) j
	2751	netcdf_var_name = TRIM( netcdf_var_name_base ) // suffix
	2752	nc_stat = NF90_INQ_VARID( id_set_ts, netcdf_var_name, &
	2753	id_var_dots(i,j) )
	2754	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 266 )
	2755	ENDDO
	2756	ENDIF
	2757
	2758	ENDDO
	2759
	2760	!
	2761	!-- Change the title attribute on file
	2762	nc_stat = NF90_PUT_ATT( id_set_ts, NF90_GLOBAL, 'title', &
	2763	TRIM( run_description_header ) )
	2764	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 267 )
	2765	PRINT, '** NetCDF file for time series from previous run found.'
	2766	PRINT*, ' This file will be extended.'
	2767
	2768
	2769	CASE ( 'sp_new' )
	2770
	2771	!
	2772	!-- Define some global attributes of the dataset
	2773	IF ( averaging_interval_sp /= 0.0 ) THEN
	2774	WRITE (time_average_text,'('', '',F7.1,'' s average'')') &
	2775	averaging_interval_sp
	2776	nc_stat = NF90_PUT_ATT( id_set_sp, NF90_GLOBAL, 'title', &
	2777	TRIM( run_description_header ) // &
	2778	TRIM( time_average_text ) )
	2779	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 268 )
	2780
	2781	WRITE ( time_average_text,'(F7.1,'' s avg'')' ) averaging_interval_sp
	2782	nc_stat = NF90_PUT_ATT( id_set_sp, NF90_GLOBAL, 'time_avg', &
	2783	TRIM( time_average_text ) )
	2784	ELSE
	2785	nc_stat = NF90_PUT_ATT( id_set_sp, NF90_GLOBAL, 'title', &
	2786	TRIM( run_description_header ) )
	2787	ENDIF
	2788	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 269 )
	2789
	2790	!
	2791	!-- Define time coordinate for spectra (unlimited dimension)
	2792	nc_stat = NF90_DEF_DIM( id_set_sp, 'time', NF90_UNLIMITED, &
	2793	id_dim_time_sp )
	2794	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 270 )
	2795
	2796	nc_stat = NF90_DEF_VAR( id_set_sp, 'time', NF90_DOUBLE, &
	2797	id_dim_time_sp, id_var_time_sp )
	2798	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 271 )
	2799
	2800	nc_stat = NF90_PUT_ATT( id_set_sp, id_var_time_sp, 'units', 'seconds')
	2801	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 272 )
	2802
	2803	!
	2804	!-- Define the spatial dimensions and coordinates for spectra.
	2805	!-- First, determine the number of vertical levels for which spectra
	2806	!-- are to be output.
	2807	ns = 1
	2808	DO WHILE ( comp_spectra_level(ns) /= 999999 .AND. ns <= 100 )
	2809	ns = ns + 1
	2810	ENDDO
	2811	ns = ns - 1
	2812
	2813	!
	2814	!-- Define vertical coordinate grid (zu grid)
	2815	nc_stat = NF90_DEF_DIM( id_set_sp, 'zu_sp', ns, id_dim_zu_sp )
	2816	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 273 )
	2817
	2818	nc_stat = NF90_DEF_VAR( id_set_sp, 'zu_sp', NF90_DOUBLE, &
	2819	id_dim_zu_sp, id_var_zu_sp )
	2820	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 274 )
	2821
	2822	nc_stat = NF90_PUT_ATT( id_set_sp, id_var_zu_sp, 'units', 'meters' )
	2823	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 275 )
	2824
	2825	!
	2826	!-- Define vertical coordinate grid (zw grid)
	2827	nc_stat = NF90_DEF_DIM( id_set_sp, 'zw_sp', ns, id_dim_zw_sp )
	2828	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 276 )
	2829
	2830	nc_stat = NF90_DEF_VAR( id_set_sp, 'zw_sp', NF90_DOUBLE, &
	2831	id_dim_zw_sp, id_var_zw_sp )
	2832	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 277 )
	2833
	2834	nc_stat = NF90_PUT_ATT( id_set_sp, id_var_zw_sp, 'units', 'meters' )
	2835	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 278 )
	2836
	2837	!
	2838	!-- Define x-axis
	2839	nc_stat = NF90_DEF_DIM( id_set_sp, 'k_x', nx/2, id_dim_x_sp )
	2840	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 279 )
	2841
	2842	nc_stat = NF90_DEF_VAR( id_set_sp, 'k_x', NF90_DOUBLE, id_dim_x_sp, &
	2843	id_var_x_sp )
	2844	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 280 )
	2845
	2846	nc_stat = NF90_PUT_ATT( id_set_sp, id_var_x_sp, 'units', 'm-1' )
	2847	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 281 )
	2848
	2849	!
	2850	!-- Define y-axis
	2851	nc_stat = NF90_DEF_DIM( id_set_sp, 'k_y', ny/2, id_dim_y_sp )
	2852	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 282 )
	2853
	2854	nc_stat = NF90_DEF_VAR( id_set_sp, 'k_y', NF90_DOUBLE, id_dim_y_sp, &
	2855	id_var_y_sp )
	2856	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 283 )
	2857
	2858	nc_stat = NF90_PUT_ATT( id_set_sp, id_var_y_sp, 'units', 'm-1' )
	2859	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 284 )
	2860
	2861	!
	2862	!-- Define the variables
	2863	var_list = ';'
	2864	i = 1
	2865	DO WHILE ( data_output_sp(i) /= ' ' .AND. i <= 10 )
	2866
	2867	IF ( INDEX( spectra_direction(i), 'x' ) /= 0 ) THEN
	2868
	2869	!
	2870	!-- Define the variable
	2871	netcdf_var_name = TRIM( data_output_sp(i) ) // '_x'
	2872	IF ( data_output_sp(i) == 'w' ) THEN
	2873	nc_stat = NF90_DEF_VAR( id_set_sp, netcdf_var_name, &
	2874	nc_precision(7), (/ id_dim_x_sp, &
	2875	id_dim_zw_sp, id_dim_time_sp /), &
	2876	id_var_dospx(i) )
	2877	ELSE
	2878	nc_stat = NF90_DEF_VAR( id_set_sp, netcdf_var_name, &
	2879	nc_precision(7), (/ id_dim_x_sp, &
	2880	id_dim_zu_sp, id_dim_time_sp /), &
	2881	id_var_dospx(i) )
	2882	ENDIF
	2883	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 285 )
	2884
	2885	nc_stat = NF90_PUT_ATT( id_set_sp, id_var_dospx(i), &
	2886	'long_name', netcdf_var_name )
	2887	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 286 )
	2888
	2889	nc_stat = NF90_PUT_ATT( id_set_sp, id_var_dospx(i), &
	2890	'units', 'unknown' )
	2891	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 287 )
	2892
	2893	var_list = TRIM( var_list ) // TRIM( netcdf_var_name ) // ';'
	2894
	2895	ENDIF
	2896
	2897	IF ( INDEX( spectra_direction(i), 'y' ) /= 0 ) THEN
	2898
	2899	!
	2900	!-- Define the variable
	2901	netcdf_var_name = TRIM( data_output_sp(i) ) // '_y'
	2902	IF ( data_output_sp(i) == 'w' ) THEN
	2903	nc_stat = NF90_DEF_VAR( id_set_sp, netcdf_var_name, &
	2904	nc_precision(7), (/ id_dim_y_sp, &
	2905	id_dim_zw_sp, id_dim_time_sp /), &
	2906	id_var_dospy(i) )
	2907	ELSE
	2908	nc_stat = NF90_DEF_VAR( id_set_sp, netcdf_var_name, &
	2909	nc_precision(7), (/ id_dim_y_sp, &
	2910	id_dim_zu_sp, id_dim_time_sp /), &
	2911	id_var_dospy(i) )
	2912	ENDIF
	2913	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 288 )
	2914
	2915	nc_stat = NF90_PUT_ATT( id_set_sp, id_var_dospy(i), &
	2916	'long_name', netcdf_var_name )
	2917	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 289 )
	2918
	2919	nc_stat = NF90_PUT_ATT( id_set_sp, id_var_dospy(i), &
	2920	'units', 'unknown' )
	2921	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 290 )
	2922
	2923	var_list = TRIM( var_list ) // TRIM( netcdf_var_name ) // ';'
	2924
	2925	ENDIF
	2926
	2927	i = i + 1
	2928
	2929	ENDDO
	2930
	2931	!
	2932	!-- Write the list of variables as global attribute (this is used by
	2933	!-- restart runs)
	2934	nc_stat = NF90_PUT_ATT( id_set_sp, NF90_GLOBAL, 'VAR_LIST', var_list )
	2935	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 291 )
	2936
	2937	!
	2938	!-- Leave NetCDF define mode
	2939	nc_stat = NF90_ENDDEF( id_set_sp )
	2940	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 292 )
	2941
	2942	!
	2943	!-- Write axis data: zu_sp, zw_sp, k_x, k_y
	2944	ALLOCATE( netcdf_data(1:ns) )
	2945
	2946	!
	2947	!-- Write zu data
	2948	netcdf_data(1:ns) = zu( comp_spectra_level(1:ns) )
	2949	nc_stat = NF90_PUT_VAR( id_set_sp, id_var_zu_sp, netcdf_data, &
	2950	start = (/ 1 /), count = (/ ns /) )
	2951	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 293 )
	2952
	2953	!
	2954	!-- Write zw data
	2955	netcdf_data(1:ns) = zw( comp_spectra_level(1:ns) )
	2956	nc_stat = NF90_PUT_VAR( id_set_sp, id_var_zw_sp, netcdf_data, &
	2957	start = (/ 1 /), count = (/ ns /) )
	2958	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 294 )
	2959
	2960	DEALLOCATE( netcdf_data )
	2961
	2962	!
	2963	!-- Write data for x and y axis (wavenumbers)
	2964	ALLOCATE( netcdf_data(nx/2) )
	2965	DO i = 1, nx/2
	2966	netcdf_data(i) = 2.0 * pi * i / ( dx * ( nx + 1 ) )
	2967	ENDDO
	2968
	2969	nc_stat = NF90_PUT_VAR( id_set_sp, id_var_x_sp, netcdf_data, &
	2970	start = (/ 1 /), count = (/ nx/2 /) )
	2971	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 295 )
	2972
	2973	DEALLOCATE( netcdf_data )
	2974
	2975	ALLOCATE( netcdf_data(ny/2) )
	2976	DO i = 1, ny/2
	2977	netcdf_data(i) = 2.0 * pi * i / ( dy * ( ny + 1 ) )
	2978	ENDDO
	2979
	2980	nc_stat = NF90_PUT_VAR( id_set_sp, id_var_y_sp, netcdf_data, &
	2981	start = (/ 1 /), count = (/ ny/2 /) )
	2982	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 296 )
	2983
	2984	DEALLOCATE( netcdf_data )
	2985
	2986
	2987	CASE ( 'sp_ext' )
	2988
	2989	!
	2990	!-- Get the list of variables and compare with the actual run.
	2991	!-- First var_list_old has to be reset, since GET_ATT does not assign
	2992	!-- trailing blanks.
	2993	var_list_old = ' '
	2994	nc_stat = NF90_GET_ATT( id_set_sp, NF90_GLOBAL, 'VAR_LIST', &
	2995	var_list_old )
	2996	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 297 )
	2997
	2998	var_list = ';'
	2999	i = 1
	3000	DO WHILE ( data_output_sp(i) /= ' ' .AND. i <= 10 )
	3001
	3002	IF ( INDEX( spectra_direction(i), 'x' ) /= 0 ) THEN
	3003	netcdf_var_name = TRIM( data_output_sp(i) ) // '_x'
	3004	var_list = TRIM( var_list ) // TRIM( netcdf_var_name ) // ';'
	3005	ENDIF
	3006
	3007	IF ( INDEX( spectra_direction(i), 'y' ) /= 0 ) THEN
	3008	netcdf_var_name = TRIM( data_output_sp(i) ) // '_y'
	3009	var_list = TRIM( var_list ) // TRIM( netcdf_var_name ) // ';'
	3010	ENDIF
	3011
	3012	i = i + 1
	3013
	3014	ENDDO
	3015
	3016	IF ( TRIM( var_list ) /= TRIM( var_list_old ) ) THEN
	3017	PRINT*, '+++ WARNING: NetCDF file for spectra from previous ' //&
	3018	'run found,'
	3019	PRINT*, ' but this file cannot be extended due to' // &
	3020	' variable mismatch.'
	3021	PRINT*, ' New file is created instead.'
	3022	extend = .FALSE.
	3023	RETURN
	3024	ENDIF
	3025
	3026	!
	3027	!-- Determine the number of current vertical levels for which spectra
	3028	!-- shall be output
	3029	ns = 1
	3030	DO WHILE ( comp_spectra_level(ns) /= 999999 .AND. ns <= 100 )
	3031	ns = ns + 1
	3032	ENDDO
	3033	ns = ns - 1
	3034
	3035	!
	3036	!-- Get and compare the number of vertical levels
	3037	nc_stat = NF90_INQ_VARID( id_set_sp, 'zu_sp', id_var_zu_sp )
	3038	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 298 )
	3039
	3040	nc_stat = NF90_INQUIRE_VARIABLE( id_set_sp, id_var_zu_sp, &
	3041	dimids = id_dim_zu_sp_old )
	3042	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 299 )
	3043	id_dim_zu_sp = id_dim_zu_sp_old(1)
	3044
	3045	nc_stat = NF90_INQUIRE_DIMENSION( id_set_sp, id_dim_zu_sp, &
	3046	len = ns_old )
	3047	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 300 )
	3048
	3049	IF ( ns /= ns_old ) THEN
	3050	PRINT*, '+++ WARNING: NetCDF file for spectra from previous ' //&
	3051	'run found,'
	3052	PRINT*, ' but this file cannot be extended due to' // &
	3053	' mismatch in number of'
	3054	PRINT*, ' vertical levels.'
	3055	PRINT*, ' New file is created instead.'
	3056	extend = .FALSE.
	3057	RETURN
	3058	ENDIF
	3059
	3060	!
	3061	!-- Get and compare the heights of the cross sections
	3062	ALLOCATE( netcdf_data(1:ns_old) )
	3063
	3064	nc_stat = NF90_GET_VAR( id_set_sp, id_var_zu_sp, netcdf_data )
	3065	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 301 )
	3066
	3067	DO i = 1, ns
	3068	IF ( zu(comp_spectra_level(i)) /= netcdf_data(i) ) THEN
	3069	PRINT*, '+++ WARNING: NetCDF file for spectra from previou' // &
	3070	's run found,'
	3071	PRINT*, ' but this file cannot be extended due' // &
	3072	' to mismatch in heights'
	3073	PRINT*, ' of vertical levels.'
	3074	PRINT*, ' New file is created instead.'
	3075	extend = .FALSE.
	3076	RETURN
	3077	ENDIF
	3078	ENDDO
	3079
	3080	DEALLOCATE( netcdf_data )
	3081
	3082	!
	3083	!-- Get the id of the time coordinate (unlimited coordinate) and its
	3084	!-- last index on the file. The next time level is plsp..count+1.
	3085	!-- The current time must be larger than the last output time
	3086	!-- on the file.
	3087	nc_stat = NF90_INQ_VARID( id_set_sp, 'time', id_var_time_sp )
	3088	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 302 )
	3089
	3090	nc_stat = NF90_INQUIRE_VARIABLE( id_set_sp, id_var_time_sp, &
	3091	dimids = id_dim_time_old )
	3092	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 303 )
	3093	id_dim_time_sp = id_dim_time_old(1)
	3094
	3095	nc_stat = NF90_INQUIRE_DIMENSION( id_set_sp, id_dim_time_sp, &
	3096	len = dosp_time_count )
	3097	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 304 )
	3098
	3099	nc_stat = NF90_GET_VAR( id_set_sp, id_var_time_sp, &
	3100	last_time_coordinate, &
	3101	start = (/ dosp_time_count /), &
	3102	count = (/ 1 /) )
	3103	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 305 )
	3104
	3105	IF ( last_time_coordinate(1) >= simulated_time ) THEN
	3106	PRINT*, '+++ WARNING: NetCDF file for spectra from previous ' // &
	3107	'run found,'
	3108	PRINT*, ' but this file cannot be extended becaus' // &
	3109	'e the current output time'
	3110	PRINT*, ' is less or equal than the last output t' // &
	3111	'ime on this file.'
	3112	PRINT*, ' New file is created instead.'
	3113	dosp_time_count = 0
	3114	extend = .FALSE.
	3115	RETURN
	3116	ENDIF
	3117
	3118	!
	3119	!-- Dataset seems to be extendable.
	3120	!-- Now get the variable ids.
	3121	i = 1
	3122	DO WHILE ( data_output_sp(i) /= ' ' .AND. i <= 10 )
	3123
	3124	IF ( INDEX( spectra_direction(i), 'x' ) /= 0 ) THEN
	3125	netcdf_var_name = TRIM( data_output_sp(i) ) // '_x'
	3126	nc_stat = NF90_INQ_VARID( id_set_sp, netcdf_var_name, &
	3127	id_var_dospx(i) )
	3128	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 306 )
	3129	ENDIF
	3130
	3131	IF ( INDEX( spectra_direction(i), 'y' ) /= 0 ) THEN
	3132	netcdf_var_name = TRIM( data_output_sp(i) ) // '_y'
	3133	nc_stat = NF90_INQ_VARID( id_set_sp, netcdf_var_name, &
	3134	id_var_dospy(i) )
	3135	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 307 )
	3136	ENDIF
	3137
	3138	i = i + 1
	3139
	3140	ENDDO
	3141
	3142	!
	3143	!-- Change the titel attribute on file
	3144	IF ( averaging_interval_sp /= 0.0 ) THEN
	3145	WRITE (time_average_text,'('', '',F7.1,'' s average'')') &
	3146	averaging_interval_sp
	3147	nc_stat = NF90_PUT_ATT( id_set_sp, NF90_GLOBAL, 'title', &
	3148	TRIM( run_description_header ) // &
	3149	TRIM( time_average_text ) )
	3150	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 308 )
	3151
	3152	WRITE ( time_average_text,'(F7.1,'' s avg'')' ) averaging_interval_sp
	3153	nc_stat = NF90_PUT_ATT( id_set_sp, NF90_GLOBAL, 'time_avg', &
	3154	TRIM( time_average_text ) )
	3155	ELSE
	3156	nc_stat = NF90_PUT_ATT( id_set_sp, NF90_GLOBAL, 'title', &
	3157	TRIM( run_description_header ) )
	3158	ENDIF
	3159	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 309 )
	3160
	3161	PRINT, '** NetCDF file for spectra from previous run found.'
	3162	PRINT*, ' This file will be extended.'
	3163
	3164
	3165	CASE ( 'pt_new' )
	3166
	3167	!
	3168	!-- Define some global attributes of the dataset
	3169	nc_stat = NF90_PUT_ATT( id_set_prt, NF90_GLOBAL, 'title', &
	3170	TRIM( run_description_header ) )
	3171	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 310 )
	3172
	3173	!
	3174	!-- Define time coordinate for particles (unlimited dimension)
	3175	nc_stat = NF90_DEF_DIM( id_set_prt, 'time', NF90_UNLIMITED, &
	3176	id_dim_time_prt )
	3177	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 311 )
	3178
	3179	nc_stat = NF90_DEF_VAR( id_set_prt, 'time', NF90_DOUBLE, &
	3180	id_dim_time_prt, id_var_time_prt )
	3181	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 312 )
	3182
	3183	nc_stat = NF90_PUT_ATT( id_set_prt, id_var_time_prt, 'units', &
	3184	'seconds' )
	3185	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 313 )
	3186
	3187	!
	3188	!-- Define particle coordinate (maximum particle number)
	3189	nc_stat = NF90_DEF_DIM( id_set_prt, 'particle_number', &
	3190	maximum_number_of_particles, id_dim_prtnum )
	3191	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 314 )
	3192
	3193	nc_stat = NF90_DEF_VAR( id_set_prt, 'particle_number', NF90_DOUBLE, &
	3194	id_dim_prtnum, id_var_prtnum )
	3195	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 315 )
	3196
	3197	nc_stat = NF90_PUT_ATT( id_set_prt, id_var_prtnum, 'units', &
	3198	'particle number' )
	3199	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 316 )
	3200
	3201	!
	3202	!-- Define variable which contains the real number of particles in use
	3203	nc_stat = NF90_DEF_VAR( id_set_prt, 'real_num_of_prt', NF90_INT, &
	3204	id_dim_time_prt, id_var_rnop_prt )
	3205	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 317 )
	3206
	3207	nc_stat = NF90_PUT_ATT( id_set_prt, id_var_rnop_prt, 'units', &
	3208	'particle number' )
	3209	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 318 )
	3210
	3211	!
	3212	!-- Define the variables
	3213	DO i = 1, 17
	3214
	3215	nc_stat = NF90_DEF_VAR( id_set_prt, prt_var_names(i), &
	3216	nc_precision(8), &
	3217	(/ id_dim_prtnum, id_dim_time_prt /), &
	3218	id_var_prt(i) )
	3219	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 319 )
	3220
	3221	nc_stat = NF90_PUT_ATT( id_set_prt, id_var_prt(i), &
	3222	'long_name', TRIM( prt_var_names(i) ) )
	3223	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 320 )
	3224
	3225	nc_stat = NF90_PUT_ATT( id_set_prt, id_var_prt(i), &
	3226	'units', TRIM( prt_var_units(i) ) )
	3227	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 321 )
	3228
	3229	ENDDO
	3230
	3231	!
	3232	!-- Leave NetCDF define mode
	3233	nc_stat = NF90_ENDDEF( id_set_prt )
	3234	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 322 )
	3235
	3236
	3237	CASE ( 'pt_ext' )
	3238
	3239	!
	3240	!-- Get the id of the time coordinate (unlimited coordinate) and its
	3241	!-- last index on the file. The next time level is prt..count+1.
	3242	!-- The current time must be larger than the last output time
	3243	!-- on the file.
	3244	nc_stat = NF90_INQ_VARID( id_set_prt, 'time', id_var_time_prt )
	3245	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 323 )
	3246
	3247	nc_stat = NF90_INQUIRE_VARIABLE( id_set_prt, id_var_time_prt, &
	3248	dimids = id_dim_time_old )
	3249	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 324 )
	3250	id_dim_time_prt = id_dim_time_old(1)
	3251
	3252	nc_stat = NF90_INQUIRE_DIMENSION( id_set_prt, id_dim_time_prt, &
	3253	len = prt_time_count )
	3254	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 325 )
	3255
	3256	nc_stat = NF90_GET_VAR( id_set_prt, id_var_time_prt, &
	3257	last_time_coordinate, &
	3258	start = (/ prt_time_count /), &
	3259	count = (/ 1 /) )
	3260	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 326 )
	3261
	3262	IF ( last_time_coordinate(1) >= simulated_time ) THEN
	3263	PRINT*, '+++ WARNING: NetCDF file for particles from previous ' //&
	3264	'run found,'
	3265	PRINT*, ' but this file cannot be extended becaus' // &
	3266	'e the current output time'
	3267	PRINT*, ' is less or equal than the last output t' // &
	3268	'ime on this file.'
	3269	PRINT*, ' New file is created instead.'
	3270	prt_time_count = 0
	3271	extend = .FALSE.
	3272	RETURN
	3273	ENDIF
	3274
	3275	!
	3276	!-- Dataset seems to be extendable.
	3277	!-- Now get the variable ids.
	3278	nc_stat = NF90_INQ_VARID( id_set_prt, 'real_num_of_prt', &
	3279	id_var_rnop_prt )
	3280	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 327 )
	3281
	3282	DO i = 1, 17
	3283
	3284	nc_stat = NF90_INQ_VARID( id_set_prt, prt_var_names(i), &
	3285	id_var_prt(i) )
	3286	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 328 )
	3287
	3288	ENDDO
	3289
	3290	IF ( myid == 0 ) THEN
	3291	PRINT, '** NetCDF file for particles from previous run found.'
	3292	PRINT*, ' This file will be extended.'
	3293	ENDIF
	3294
	3295
	3296	CASE ( 'ps_new' )
	3297
	3298	!
	3299	!-- Define some global attributes of the dataset
	3300	nc_stat = NF90_PUT_ATT( id_set_pts, NF90_GLOBAL, 'title', &
	3301	TRIM( run_description_header ) )
	3302	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 396 )
	3303
	3304	!
	3305	!-- Define time coordinate for particle time series (unlimited dimension)
	3306	nc_stat = NF90_DEF_DIM( id_set_pts, 'time', NF90_UNLIMITED, &
	3307	id_dim_time_pts )
	3308	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 397 )
	3309
	3310	nc_stat = NF90_DEF_VAR( id_set_pts, 'time', NF90_DOUBLE, &
	3311	id_dim_time_pts, id_var_time_pts )
	3312	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 398 )
	3313
	3314	nc_stat = NF90_PUT_ATT( id_set_pts, id_var_time_pts, 'units', &
	3315	'seconds')
	3316	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 399 )
	3317
	3318	!
	3319	!-- Define the variables. If more than one particle group is defined,
	3320	!-- define seperate variables for each group
	3321	var_list = ';'
	3322	DO i = 1, dopts_num
	3323
	3324	!
	3325	!-- First, remove those characters not allowed by NetCDF
	3326	netcdf_var_name = dopts_label(i)
	3327	CALL clean_netcdf_varname( netcdf_var_name )
	3328
	3329	DO j = 0, number_of_particle_groups
	3330
	3331	IF ( j == 0 ) THEN
	3332	suffix1 = ''
	3333	ELSE
	3334	WRITE ( suffix1, '(''_'',I2.2)' ) j
	3335	ENDIF
	3336
	3337	nc_stat = NF90_DEF_VAR( id_set_pts, &
	3338	TRIM( netcdf_var_name ) // suffix1, &
	3339	nc_precision(6), &
	3340	(/ id_dim_time_pts /), &
	3341	id_var_dopts(i,j) )
	3342	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 400 )
	3343
	3344	IF ( j == 0 ) THEN
	3345	nc_stat = NF90_PUT_ATT( id_set_pts, id_var_dopts(i,j), &
	3346	'long_name', &
	3347	TRIM( dopts_label(i) ) )
	3348	ELSE
	3349	nc_stat = NF90_PUT_ATT( id_set_pts, id_var_dopts(i,j), &
	3350	'long_name', &
	3351	TRIM( dopts_label(i) ) // ' PG ' // &
	3352	suffix1(2:3) )
	3353	ENDIF
	3354	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 401 )
	3355
	3356	nc_stat = NF90_PUT_ATT( id_set_pts, id_var_dopts(i,j), &
	3357	'units', TRIM( dopts_unit(i) ) )
	3358	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 402 )
	3359
	3360	var_list = TRIM(var_list) // TRIM(netcdf_var_name) // &
	3361	suffix1 // ';'
	3362
	3363	IF ( number_of_particle_groups == 1 ) EXIT
	3364
	3365	ENDDO
	3366
	3367	ENDDO
	3368
	3369	!
	3370	!-- Write the list of variables as global attribute (this is used by
	3371	!-- restart runs)
	3372	nc_stat = NF90_PUT_ATT( id_set_pts, NF90_GLOBAL, 'VAR_LIST', &
	3373	var_list )
	3374	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 403 )
	3375
	3376	!
	3377	!-- Leave NetCDF define mode
	3378	nc_stat = NF90_ENDDEF( id_set_pts )
	3379	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 404 )
	3380
	3381
	3382	CASE ( 'ps_ext' )
	3383
	3384	!
	3385	!-- Get the list of variables and compare with the actual run.
	3386	!-- First var_list_old has to be reset, since GET_ATT does not assign
	3387	!-- trailing blanks.
	3388	var_list_old = ' '
	3389	nc_stat = NF90_GET_ATT( id_set_pts, NF90_GLOBAL, 'VAR_LIST', &
	3390	var_list_old )
	3391	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 405 )
	3392
	3393	var_list = ';'
	3394	i = 1
	3395	DO i = 1, dopts_num
	3396
	3397	netcdf_var_name = dopts_label(i)
	3398	CALL clean_netcdf_varname( netcdf_var_name )
	3399
	3400	DO j = 0, number_of_particle_groups
	3401
	3402	IF ( j == 0 ) THEN
	3403	suffix1 = ''
	3404	ELSE
	3405	WRITE ( suffix1, '(''_'',I2.2)' ) j
	3406	ENDIF
	3407
	3408	var_list = TRIM(var_list) // TRIM(netcdf_var_name) // &
	3409	suffix1 // ';'
	3410
	3411	IF ( number_of_particle_groups == 1 ) EXIT
	3412
	3413	ENDDO
	3414
	3415	ENDDO
	3416
	3417	IF ( TRIM( var_list ) /= TRIM( var_list_old ) ) THEN
	3418	PRINT*, '+++ WARNING: NetCDF file for particle time series ' //&
	3419	'from previuos run found,'
	3420	PRINT*, ' but this file cannot be extended due to' // &
	3421	' variable mismatch.'
	3422	PRINT*, ' New file is created instead.'
	3423	extend = .FALSE.
	3424	RETURN
	3425	ENDIF
	3426
	3427	!
	3428	!-- Get the id of the time coordinate (unlimited coordinate) and its
	3429	!-- last index on the file. The next time level is dots..count+1.
	3430	!-- The current time must be larger than the last output time
	3431	!-- on the file.
	3432	nc_stat = NF90_INQ_VARID( id_set_pts, 'time', id_var_time_pts )
	3433	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 406 )
	3434
	3435	nc_stat = NF90_INQUIRE_VARIABLE( id_set_pts, id_var_time_pts, &
	3436	dimids = id_dim_time_old )
	3437	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 407 )
	3438	id_dim_time_pts = id_dim_time_old(1)
	3439
	3440	nc_stat = NF90_INQUIRE_DIMENSION( id_set_pts, id_dim_time_pts, &
	3441	len = dopts_time_count )
	3442	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 408 )
	3443
	3444	nc_stat = NF90_GET_VAR( id_set_pts, id_var_time_pts, &
	3445	last_time_coordinate, &
	3446	start = (/ dopts_time_count /), &
	3447	count = (/ 1 /) )
	3448	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 409 )
	3449
	3450	IF ( last_time_coordinate(1) >= simulated_time ) THEN
	3451	PRINT*, '+++ WARNING: NetCDF file for time series from' // &
	3452	' previuos run found,'
	3453	PRINT*, ' but this file cannot be extended becaus' // &
	3454	'e the current output time'
	3455	PRINT*, ' is less or equal than the last output t' // &
	3456	'ime on this file.'
	3457	PRINT*, ' New file is created instead.'
	3458	dopts_time_count = 0
	3459	extend = .FALSE.
	3460	RETURN
	3461	ENDIF
	3462
	3463	!
	3464	!-- Dataset seems to be extendable.
	3465	!-- Now get the variable ids
	3466	i = 1
	3467	DO i = 1, dopts_num
	3468
	3469	netcdf_var_name_base = dopts_label(i)
	3470	CALL clean_netcdf_varname( netcdf_var_name_base )
	3471
	3472	DO j = 0, number_of_particle_groups
	3473
	3474	IF ( j == 0 ) THEN
	3475	suffix1 = ''
	3476	ELSE
	3477	WRITE ( suffix1, '(''_'',I2.2)' ) j
	3478	ENDIF
	3479
	3480	netcdf_var_name = TRIM( netcdf_var_name_base ) // suffix1
	3481
	3482	nc_stat = NF90_INQ_VARID( id_set_pts, netcdf_var_name, &
	3483	id_var_dopts(i,j) )
	3484	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 410 )
	3485
	3486	IF ( number_of_particle_groups == 1 ) EXIT
	3487
	3488	ENDDO
	3489
	3490	ENDDO
	3491
	3492	!
	3493	!-- Change the title attribute on file
	3494	nc_stat = NF90_PUT_ATT( id_set_pts, NF90_GLOBAL, 'title', &
	3495	TRIM( run_description_header ) )
	3496	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 411 )
	3497	PRINT, '** NetCDF file for particle time series from previous ', &
	3498	'run found.'
	3499	PRINT*, ' This file will be extended.'
	3500
	3501
	3502	CASE DEFAULT
	3503
	3504	PRINT*, '+++ define_netcdf_header: mode "', mode, '" not supported'
	3505
	3506	END SELECT
	3507
	3508	#endif
	3509	END SUBROUTINE define_netcdf_header
	3510
	3511
	3512
	3513	SUBROUTINE handle_netcdf_error( errno )
	3514	#if defined( __netcdf )
	3515
	3516	!------------------------------------------------------------------------------!
	3517	!
	3518	! Description:
	3519	! ------------
	3520	! Prints out a text message corresponding to the current status.
	3521	!------------------------------------------------------------------------------!
	3522
	3523	USE netcdf
	3524	USE netcdf_control
	3525	USE pegrid
	3526
	3527	IMPLICIT NONE
	3528
	3529	INTEGER :: errno
	3530
	3531	IF ( nc_stat /= NF90_NOERR ) THEN
	3532	PRINT*, '+++ netcdf error ', errno,': ', TRIM( NF90_STRERROR( nc_stat ) )
	3533	#if defined( __parallel )
	3534	CALL MPI_ABORT( comm2d, 9999, ierr )
	3535	#else
	3536	CALL local_stop
	3537	#endif
	3538	ENDIF
	3539
	3540	#endif
	3541	END SUBROUTINE handle_netcdf_error
	3542
	3543
	3544
	3545	SUBROUTINE clean_netcdf_varname( string )
	3546	#if defined( __netcdf )
	3547
	3548	!------------------------------------------------------------------------------!
	3549	!
	3550	! Description:
	3551	! ------------
	3552	! Replace those characters in string which are not allowed by NetCDF.
	3553	!------------------------------------------------------------------------------!
	3554
	3555	USE netcdf_control
	3556
	3557	IMPLICIT NONE
	3558
	3559	CHARACTER (LEN=10), INTENT(INOUT) :: string
	3560
	3561	INTEGER :: i, ic
	3562
	3563	DO i = 1, replace_num
	3564	DO
	3565	ic = INDEX( string, replace_char(i) )
	3566	IF ( ic /= 0 ) THEN
	3567	string(ic:ic) = replace_by(i)
	3568	ELSE
	3569	EXIT
	3570	ENDIF
	3571	ENDDO
	3572	ENDDO
	3573
	3574	#endif
	3575	END SUBROUTINE clean_netcdf_varname

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

Context Navigation

source: palm/trunk/SOURCE/netcdf.f90 @ 77

Download in other formats: