Home

Context Navigation

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

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |