Home

Context Navigation

check_open.f90 @ 80

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

New:
---

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

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

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

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

new inipar-parameter loop_optimization to control the loop optimization method

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

new user interface user_advec_particles

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

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

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

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

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

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

New: wall_fluxes

Changed:

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

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

+age_m in particle_type

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

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

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

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

q is not allowed to become negative (prognostic_equations).

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

d3par-parameter moisture renamed to humidity.

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

vtk directives removed from main program.

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

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

Errors:

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

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

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

Check for possible negative humidities in the initial humidity profile.

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

Makefile
check_parameters, init_1d_model, init_3d_model, user_interface

Property svn:keywords set to Id

File size: 44.4 KB

Rev	Line
[1]	1	SUBROUTINE check_open( file_id )
	2
	3	!------------------------------------------------------------------------------!
	4	! Actual revisions:
	5	! -----------------
[77]	6	!
[1]	7	!
	8	! Former revisions:
	9	! -----------------
[3]	10	! $Id: check_open.f90 77 2007-03-29 04:26:56Z raasch $
[77]	11	!
	12	! 46 2007-03-05 06:00:47Z raasch
	13	! +netcdf_64bit_3d to switch on 64bit offset only for 3D files
	14	!
[3]	15	! RCS Log replace by Id keyword, revision history cleaned up
	16	!
[1]	17	! Revision 1.44 2006/08/22 13:48:34 raasch
	18	! xz and yz cross sections now up to nzt+1
	19	!
	20	! Revision 1.1 1997/08/11 06:10:55 raasch
	21	! Initial revision
	22	!
	23	!
	24	! Description:
	25	! ------------
	26	! Check if file unit is open. If not, open file and, if necessary, write a
	27	! header or start other initializing actions, respectively.
	28	!------------------------------------------------------------------------------!
	29
	30	USE array_kind
	31	USE arrays_3d
	32	USE control_parameters
	33	USE grid_variables
	34	USE indices
	35	USE netcdf_control
	36	USE particle_attributes
	37	USE pegrid
	38	USE profil_parameter
	39	USE statistics
	40
	41	IMPLICIT NONE
	42
	43	CHARACTER (LEN=2) :: suffix
	44	CHARACTER (LEN=20) :: return_addres, return_usern, xtext = 'time in s'
	45	CHARACTER (LEN=30) :: filename
	46	CHARACTER (LEN=40) :: avs_coor_file, avs_coor_file_localname, &
	47	avs_data_file_localname
	48	CHARACTER (LEN=80) :: rtext
	49	CHARACTER (LEN=100) :: avs_coor_file_catalog, avs_data_file_catalog, &
	50	batch_scp, zeile
	51	CHARACTER (LEN=400) :: command
	52
	53	INTEGER :: av, anzzeile = 1, cranz, file_id, i, iaddres, ierr1, iusern, &
	54	j, k, legpos = 1, timodex = 1
	55	INTEGER, DIMENSION(10) :: cucol, klist, lstyle
	56
	57	LOGICAL :: avs_coor_file_found = .FALSE., avs_data_file_found = .FALSE., &
	58	datleg = .TRUE., get_filenames, grid = .TRUE., netcdf_extend, &
	59	rand = .TRUE., swap = .TRUE., twoxa = .TRUE., twoya = .TRUE.
	60
	61	REAL :: ansx = -999.999, ansy = -999.999, gwid = 0.1, rlegfak = 1.5, &
	62	sizex = 250.0, sizey = 40.0, texfac = 1.5
	63
	64	REAL, DIMENSION(:), ALLOCATABLE :: eta, ho, hu
	65	REAL(spk), DIMENSION(:), ALLOCATABLE :: xkoor, ykoor, zkoor
	66
	67
	68	NAMELIST /RAHMEN/ anzzeile, cranz, datleg, rtext, swap
	69	NAMELIST /CROSS/ ansx, ansy, cucol, grid, gwid, klist, legpos, lstyle, &
	70	rand, rlegfak, sizex, sizey, texfac, &
	71	timodex, twoxa, twoya, xtext
	72
	73
	74	!
	75	!-- Immediate return if file already open
	76	IF ( openfile(file_id)%opened ) RETURN
	77
	78	!
	79	!-- Only certain files are allowed to be re-opened
	80	!-- NOTE: some of the other files perhaps also could be re-opened, but it
	81	!-- has not been checked so far, if it works!
	82	IF ( openfile(file_id)%opened_before ) THEN
	83	SELECT CASE ( file_id )
	84	CASE ( 14, 21, 22, 23, 80:85 )
	85	IF ( file_id == 14 .AND. openfile(file_id)%opened_before ) THEN
	86	IF ( myid == 0 ) PRINT*, '+++ check_open: re-open of unit ', &
	87	' 14 is not verified. Please check results!'
	88	ENDIF
	89	CONTINUE
	90	CASE DEFAULT
	91	IF ( myid == 0 ) THEN
	92	PRINT*, '+++ check_open: re-opening of file-id ', file_id, &
	93	' is not allowed'
	94	ENDIF
	95	RETURN
	96	END SELECT
	97	ENDIF
	98
	99	!
	100	!-- Check if file may be opened on the relevant PE
	101	SELECT CASE ( file_id )
	102
	103	CASE ( 15, 16, 17, 18, 19, 40:49, 50:59, 81:84, 101:107, 109, 111:113, &
	104	116 )
	105
	106	IF ( myid /= 0 ) THEN
	107	PRINT*,'+++ check_open: opening file-id ',file_id, &
	108	' not allowed for PE ',myid
	109	#if defined( __parallel )
	110	CALL MPI_ABORT( comm2d, 9999, ierr )
	111	#else
	112	CALL local_stop
	113	#endif
	114	ENDIF
	115
	116	CASE ( 21, 22, 23 )
	117
	118	IF ( .NOT. data_output_2d_on_each_pe ) THEN
	119	IF ( myid /= 0 ) THEN
	120	PRINT*,'+++ check_open: opening file-id ',file_id, &
	121	' not allowed for PE ',myid
	122	#if defined( __parallel )
	123	CALL MPI_ABORT( comm2d, 9999, ierr )
	124	#else
	125	CALL local_stop
	126	#endif
	127	ENDIF
	128	ENDIF
	129
	130	CASE ( 27, 28, 29, 31, 32, 33, 71:73, 90:99 )
	131
	132	!
	133	!-- File-ids that are used temporarily in other routines
	134	PRINT*,'+++ check_open: opening file-id ',file_id, &
	135	' is not allowed since it is used otherwise'
	136
	137	END SELECT
	138
	139	!
	140	!-- Open relevant files
	141	SELECT CASE ( file_id )
	142
	143	CASE ( 11 )
	144
	145	OPEN ( 11, FILE='PARIN', FORM='FORMATTED', STATUS='OLD' )
	146
	147	CASE ( 13 )
	148
	149	IF ( myid_char == '' ) THEN
	150	OPEN ( 13, FILE='BININ'//myid_char, FORM='UNFORMATTED', &
	151	STATUS='OLD' )
	152	ELSE
	153	#if defined( __t3eh ) \|\| defined( __t3eb ) \|\| defined( __t3ej5 )
	154	!
	155	!-- Declare a cache layer for faster I/O
	156	CALL ASNUNIT ( 13, '-F cachea:512:1:1', ierr1 )
	157	#endif
	158	OPEN ( 13, FILE='BININ/'//myid_char, FORM='UNFORMATTED', &
	159	STATUS='OLD' )
	160	ENDIF
	161
	162	CASE ( 14 )
	163
	164	IF ( myid_char == '' ) THEN
	165	OPEN ( 14, FILE='BINOUT'//myid_char, FORM='UNFORMATTED', &
	166	POSITION='APPEND' )
	167	ELSE
	168	IF ( myid == 0 .AND. .NOT. openfile(file_id)%opened_before ) THEN
	169	CALL local_system( 'mkdir BINOUT' )
	170	ENDIF
	171	#if defined( __parallel )
	172	!
	173	!-- Set a barrier in order to allow that all other processors in the
	174	!-- directory created by PE0 can open their file
	175	CALL MPI_BARRIER( comm2d, ierr )
	176	#endif
	177	#if defined( __t3eh ) \|\| defined( __t3eb ) \|\| defined( __t3ej5 )
	178	!
	179	!-- Declare a cache layer for faster I/O
	180	CALL ASNUNIT ( 14, '-F cachea:512:1:1', ierr1 )
	181	#endif
	182	OPEN ( 14, FILE='BINOUT/'//myid_char_14, FORM='UNFORMATTED', &
	183	POSITION='APPEND' )
	184	ENDIF
	185
	186	CASE ( 15 )
	187
	188	OPEN ( 15, FILE='RUN_CONTROL', FORM='FORMATTED' )
	189
	190	CASE ( 16 )
	191
	192	OPEN ( 16, FILE='LIST_PROFIL', FORM='FORMATTED' )
	193
	194	CASE ( 17 )
	195
	196	OPEN ( 17, FILE='LIST_PROFIL_1D', FORM='FORMATTED' )
	197
	198	CASE ( 18 )
	199
	200	OPEN ( 18, FILE='CPU_MEASURES', FORM='FORMATTED' )
	201
	202	CASE ( 19 )
	203
	204	OPEN ( 19, FILE='HEADER', FORM='FORMATTED' )
	205
	206	CASE ( 20 )
	207
	208	IF ( myid == 0 .AND. .NOT. openfile(file_id)%opened_before ) THEN
	209	CALL local_system( 'mkdir DATA_LOG' )
	210	ENDIF
	211	IF ( myid_char == '' ) THEN
	212	OPEN ( 20, FILE='DATA_LOG/_0000', FORM='UNFORMATTED', &
	213	POSITION='APPEND' )
	214	ELSE
	215	#if defined( __parallel )
	216	!
	217	!-- Set a barrier in order to allow that all other processors in the
	218	!-- directory created by PE0 can open their file
	219	CALL MPI_BARRIER( comm2d, ierr )
	220	#endif
	221	OPEN ( 20, FILE='DATA_LOG/'//myid_char, FORM='UNFORMATTED', &
	222	POSITION='APPEND' )
	223	ENDIF
	224
	225	CASE ( 21 )
	226
	227	IF ( data_output_2d_on_each_pe ) THEN
	228	OPEN ( 21, FILE='PLOT2D_XY'//myid_char, FORM='UNFORMATTED', &
	229	POSITION='APPEND' )
	230	ELSE
	231	OPEN ( 21, FILE='PLOT2D_XY', FORM='UNFORMATTED', &
	232	POSITION='APPEND' )
	233	ENDIF
	234
	235	IF ( myid == 0 .AND. .NOT. openfile(file_id)%opened_before ) THEN
	236	!
	237	!-- Output for combine_plot_fields
	238	IF ( data_output_2d_on_each_pe .AND. myid_char /= '' ) THEN
	239	WRITE (21) -1, nx+1, -1, ny+1 ! total array size
	240	WRITE (21) 0, nx+1, 0, ny+1 ! output part
	241	ENDIF
	242	!
	243	!-- Determine and write ISO2D coordiante header
	244	ALLOCATE( eta(0:ny+1), ho(0:nx+1), hu(0:nx+1) )
	245	hu = 0.0
	246	ho = (ny+1) * dy
	247	DO i = 1, ny
	248	eta(i) = REAL( i ) / ( ny + 1.0 )
	249	ENDDO
	250	eta(0) = 0.0
	251	eta(ny+1) = 1.0
	252
	253	WRITE (21) dx,eta,hu,ho
	254	DEALLOCATE( eta, ho, hu )
	255
	256	!
	257	!-- Create output file for local parameters
	258	IF ( iso2d_output ) THEN
	259	OPEN ( 27, FILE='PLOT2D_XY_LOCAL', FORM='FORMATTED', &
	260	DELIM='APOSTROPHE' )
	261	openfile(27)%opened = .TRUE.
	262	ENDIF
	263
	264	ENDIF
	265
	266	CASE ( 22 )
	267
	268	IF ( data_output_2d_on_each_pe ) THEN
	269	OPEN ( 22, FILE='PLOT2D_XZ'//myid_char, FORM='UNFORMATTED', &
	270	POSITION='APPEND' )
	271	ELSE
	272	OPEN ( 22, FILE='PLOT2D_XZ', FORM='UNFORMATTED', &
	273	POSITION='APPEND' )
	274	ENDIF
	275
	276	IF ( myid == 0 .AND. .NOT. openfile(file_id)%opened_before ) THEN
	277	!
	278	!-- Output for combine_plot_fields
	279	IF ( data_output_2d_on_each_pe .AND. myid_char /= '' ) THEN
	280	WRITE (22) -1, nx+1, 0, nz+1 ! total array size
	281	WRITE (22) 0, nx+1, 0, nz+1 ! output part
	282	ENDIF
	283	!
	284	!-- Determine and write ISO2D coordiante header
	285	ALLOCATE( eta(0:nz+1), ho(0:nx+1), hu(0:nx+1) )
	286	hu = 0.0
	287	ho = zu(nz+1)
	288	DO i = 1, nz
	289	eta(i) = REAL( zu(i) ) / zu(nz+1)
	290	ENDDO
	291	eta(0) = 0.0
	292	eta(nz+1) = 1.0
	293
	294	WRITE (22) dx,eta,hu,ho
	295	DEALLOCATE( eta, ho, hu )
	296	!
	297	!-- Create output file for local parameters
	298	OPEN ( 28, FILE='PLOT2D_XZ_LOCAL', FORM='FORMATTED', &
	299	DELIM='APOSTROPHE' )
	300	openfile(28)%opened = .TRUE.
	301
	302	ENDIF
	303
	304	CASE ( 23 )
	305
	306	IF ( data_output_2d_on_each_pe ) THEN
	307	OPEN ( 23, FILE='PLOT2D_YZ'//myid_char, FORM='UNFORMATTED', &
	308	POSITION='APPEND' )
	309	ELSE
	310	OPEN ( 23, FILE='PLOT2D_YZ', FORM='UNFORMATTED', &
	311	POSITION='APPEND' )
	312	ENDIF
	313
	314	IF ( myid == 0 .AND. .NOT. openfile(file_id)%opened_before ) THEN
	315	!
	316	!-- Output for combine_plot_fields
	317	IF ( data_output_2d_on_each_pe .AND. myid_char /= '' ) THEN
	318	WRITE (23) -1, ny+1, 0, nz+1 ! total array size
	319	WRITE (23) 0, ny+1, 0, nz+1 ! output part
	320	ENDIF
	321	!
	322	!-- Determine and write ISO2D coordiante header
	323	ALLOCATE( eta(0:nz+1), ho(0:ny+1), hu(0:ny+1) )
	324	hu = 0.0
	325	ho = zu(nz+1)
	326	DO i = 1, nz
	327	eta(i) = REAL( zu(i) ) / zu(nz+1)
	328	ENDDO
	329	eta(0) = 0.0
	330	eta(nz+1) = 1.0
	331
	332	WRITE (23) dx,eta,hu,ho
	333	DEALLOCATE( eta, ho, hu )
	334	!
	335	!-- Create output file for local parameters
	336	OPEN ( 29, FILE='PLOT2D_YZ_LOCAL', FORM='FORMATTED', &
	337	DELIM='APOSTROPHE' )
	338	openfile(29)%opened = .TRUE.
	339
	340	ENDIF
	341
	342	CASE ( 30 )
	343
	344
	345	#if defined( __t3eb ) \|\| defined( __t3eh ) \|\| defined( __t3ej5 )
	346	!
	347	!-- In case of active data compression 32Bit integer output
	348	IF ( do3d_compress ) THEN
	349	CALL ASNUNIT( 30, '-F f77 -N ieee_32', ierr1 )
	350	ENDIF
	351	#endif
	352
	353	OPEN ( 30, FILE='PLOT3D_DATA'//myid_char, FORM='UNFORMATTED' )
	354	!
	355	!-- Write coordinate file for AVS
	356	IF ( myid == 0 ) THEN
	357	#if defined( __parallel )
	358	!
	359	!-- Specifications for combine_plot_fields
	360	IF ( .NOT. do3d_compress ) THEN
	361	WRITE ( 30 ) -1,nx+1,-1,ny+1,0,nz_do3d
	362	WRITE ( 30 ) 0,nx+1,0,ny+1,0,nz_do3d
	363	ENDIF
	364	#endif
	365	!
	366	!-- Write coordinate file for AVS:
	367	!-- First determine file names (including cyle numbers) of AVS files on
	368	!-- target machine (to which the files are to be transferred).
	369	!-- Therefore path information has to be obtained first.
	370	IF ( avs_output ) THEN
	371	CALL local_getenv( 'return_addres', 13, return_addres, iaddres )
	372	CALL local_getenv( 'return_username', 15, return_usern, iusern )
	373
	374	OPEN ( 3, FILE='OUTPUT_FILE_CONNECTIONS', FORM='FORMATTED' )
	375	DO WHILE ( .NOT. avs_coor_file_found .OR. &
	376	.NOT. avs_data_file_found )
	377
	378	READ ( 3, '(A)', END=1 ) zeile
	379
	380	SELECT CASE ( zeile(1:11) )
	381
	382	CASE ( 'PLOT3D_COOR' )
	383	READ ( 3, '(A/A)' ) avs_coor_file_catalog, &
	384	avs_coor_file_localname
	385	avs_coor_file_found = .TRUE.
	386
	387	CASE ( 'PLOT3D_DATA' )
	388	READ ( 3, '(A/A)' ) avs_data_file_catalog, &
	389	avs_data_file_localname
	390	avs_data_file_found = .TRUE.
	391
	392	CASE DEFAULT
	393	READ ( 3, '(A/A)' ) zeile, zeile
	394
	395	END SELECT
	396
	397	ENDDO
	398	!
	399	!-- Now the cycle numbers on the remote machine must be obtained
	400	!-- using batch_scp
	401	1 CLOSE ( 3 )
	402	IF ( .NOT. avs_coor_file_found .OR. &
	403	.NOT. avs_data_file_found ) THEN
	404	PRINT*, '+++ check_open: no filename for AVS-data-file ', &
	405	'found in MRUN-config-file'
	406	PRINT*, ' filename in FLD-file set to ', &
	407	'"unknown"'
	408
	409	avs_coor_file = 'unknown'
	410	avs_data_file = 'unknown'
	411	ELSE
	412	get_filenames = .TRUE.
	413	IF ( TRIM( host ) == 'hpmuk' .OR. &
	414	TRIM( host ) == 'lcmuk' ) THEN
	415	batch_scp = '/home/raasch/pub/batch_scp'
	416	ELSEIF ( TRIM( host ) == 'nech' ) THEN
	417	batch_scp = '/ipf/b/b323011/pub/batch_scp'
	418	ELSEIF ( TRIM( host ) == 'ibmh' .OR. &
	419	TRIM( host ) == 'ibmb' ) THEN
	420	batch_scp = '/home/h/niksiraa/pub/batch_scp'
	421	ELSEIF ( TRIM( host ) == 't3eb' ) THEN
	422	batch_scp = '/home/nhbksira/pub/batch_scp'
	423	ELSE
	424	PRINT*,'+++ check_open: no path for batch_scp on host "',&
	425	TRIM( host ), '"'
	426	get_filenames = .FALSE.
	427	ENDIF
	428
	429	IF ( get_filenames ) THEN
	430	!
	431	!-- Determine the coordinate file name.
	432	!-- /etc/passwd serves as Dummy-Datei, because it is not
	433	!-- really transferred.
	434	command = TRIM( batch_scp ) // ' -n -u ' // &
	435	return_usern(1:iusern) // ' ' // &
	436	return_addres(1:iaddres) // ' /etc/passwd "' // &
	437	TRIM( avs_coor_file_catalog ) // '" ' // &
	438	TRIM( avs_coor_file_localname ) // ' > REMOTE_FILENAME'
	439
	440	CALL local_system( command )
	441	OPEN ( 3, FILE='REMOTE_FILENAME', FORM='FORMATTED' )
	442	READ ( 3, '(A)' ) avs_coor_file
	443	CLOSE ( 3 )
	444	!
	445	!-- Determine the data file name
	446	command = TRIM( batch_scp ) // ' -n -u ' // &
	447	return_usern(1:iusern) // ' ' // &
	448	return_addres(1:iaddres) // ' /etc/passwd "' // &
	449	TRIM( avs_data_file_catalog ) // '" ' // &
	450	TRIM( avs_data_file_localname ) // ' > REMOTE_FILENAME'
	451
	452	CALL local_system( command )
	453	OPEN ( 3, FILE='REMOTE_FILENAME', FORM='FORMATTED' )
	454	READ ( 3, '(A)' ) avs_data_file
	455	CLOSE ( 3 )
	456
	457	ELSE
	458
	459	avs_coor_file = 'unknown'
	460	avs_data_file = 'unknown'
	461
	462	ENDIF
	463
	464	ENDIF
	465
	466	!
	467	!-- Output of the coordinate file description for FLD-file
	468	OPEN ( 33, FILE='PLOT3D_FLD_COOR', FORM='FORMATTED' )
	469	openfile(33)%opened = .TRUE.
	470	WRITE ( 33, 3300 ) TRIM( avs_coor_file ), &
	471	TRIM( avs_coor_file ), (nx+2)*4, &
	472	TRIM( avs_coor_file ), (nx+2)4+(ny+2)4
	473
	474
	475	ALLOCATE( xkoor(0:nx+1), ykoor(0:ny+1), zkoor(0:nz_do3d) )
	476	DO i = 0, nx+1
	477	xkoor(i) = i * dx
	478	ENDDO
	479	DO j = 0, ny+1
	480	ykoor(j) = j * dy
	481	ENDDO
	482	DO k = 0, nz_do3d
	483	zkoor(k) = zu(k)
	484	ENDDO
	485
	486	!
	487	!-- Create and write on AVS coordinate file
	488	OPEN ( 31, FILE='PLOT3D_COOR', FORM='UNFORMATTED' )
	489	openfile(31)%opened = .TRUE.
	490
	491	WRITE (31) xkoor, ykoor, zkoor
	492	DEALLOCATE( xkoor, ykoor, zkoor )
	493
	494	!
	495	!-- Create FLD file (being written on in close_file)
	496	OPEN ( 32, FILE='PLOT3D_FLD', FORM='FORMATTED' )
	497	openfile(32)%opened = .TRUE.
	498
	499	!
	500	!-- Create flag file for compressed 3D output,
	501	!-- influences output commands in mrun
	502	IF ( do3d_compress ) THEN
	503	OPEN ( 3, FILE='PLOT3D_COMPRESSED', FORM='FORMATTED' )
	504	WRITE ( 3, '(1X)' )
	505	CLOSE ( 3 )
	506	ENDIF
	507
	508	ENDIF
	509
	510	ENDIF
	511
	512	!
	513	!-- In case of data compression output of the coordinates of the
	514	!-- corresponding partial array of a PE only once at the top of the file
	515	IF ( avs_output .AND. do3d_compress ) THEN
	516	WRITE ( 30 ) nxl-1, nxr+1, nys-1, nyn+1, nzb, nz_do3d
	517	ENDIF
	518
	519	CASE ( 40:49 )
	520
	521	IF ( statistic_regions == 0 .AND. file_id == 40 ) THEN
	522	suffix = ''
	523	ELSE
	524	WRITE ( suffix, '(''_'',I1)' ) file_id - 40
	525	ENDIF
	526	OPEN ( file_id, FILE='PLOT1D_DATA'//TRIM( suffix ), FORM='FORMATTED' )
	527	!
	528	!-- Write contents comments at the top of the file
	529	WRITE ( file_id, 4000 ) TRIM( run_description_header ) // ' ' // &
	530	TRIM( region( file_id - 40 ) )
	531
	532	CASE ( 50:59 )
	533
	534	IF ( statistic_regions == 0 .AND. file_id == 50 ) THEN
	535	suffix = ''
	536	ELSE
	537	WRITE ( suffix, '(''_'',I1)' ) file_id - 50
	538	ENDIF
	539	OPEN ( file_id, FILE='PLOTTS_DATA'//TRIM( suffix ), FORM='FORMATTED',&
	540	RECL=496 )
	541	!
	542	!-- Write PROFIL parameter file for output of time series
	543	!-- NOTE: To be on the safe side, this output is done at the beginning of
	544	!-- the model run (in case of collapse) and it is repeated in
	545	!-- close_file, then, however, with value ranges for the coordinate
	546	!-- systems
	547	!
	548	!-- Firstly determine the number of the coordinate systems to be drawn
	549	cranz = 0
	550	DO j = 1, 10
	551	IF ( cross_ts_number_count(j) /= 0 ) cranz = cranz+1
	552	ENDDO
	553	rtext = '\1.0 ' // TRIM( run_description_header ) // ' ' // &
	554	TRIM( region( file_id - 50 ) )
	555	!
	556	!-- Write RAHMEN parameter
	557	OPEN ( 90, FILE='PLOTTS_PAR'//TRIM( suffix ), FORM='FORMATTED', &
	558	DELIM='APOSTROPHE' )
	559	WRITE ( 90, RAHMEN )
	560	!
	561	!-- Determine and write CROSS parameters for the individual coordinate
	562	!-- systems
	563	DO j = 1, 10
	564	k = cross_ts_number_count(j)
	565	IF ( k /= 0 ) THEN
	566	!
	567	!-- Store curve numbers, colours and line style
	568	klist(1:k) = cross_ts_numbers(1:k,j)
	569	klist(k+1:10) = 999999
	570	cucol(1:k) = linecolors(1:k)
	571	lstyle(1:k) = linestyles(1:k)
	572	!
	573	!-- Write CROSS parameter
	574	WRITE ( 90, CROSS )
	575
	576	ENDIF
	577	ENDDO
	578
	579	CLOSE ( 90 )
	580	!
	581	!-- Write all labels at the top of the data file, but only during the
	582	!-- first run of a sequence of jobs. The following jobs copy the time
	583	!-- series data to the bottom of that file.
	584	IF ( runnr == 0 ) THEN
	585	WRITE ( file_id, 5000 ) TRIM( run_description_header ) // &
	586	' ' // TRIM( region( file_id - 50 ) )
	587	ENDIF
	588
	589
	590	CASE ( 80 )
	591
	592	IF ( myid_char == '' ) THEN
	593	OPEN ( 80, FILE='PARTICLE_INFOS'//myid_char, FORM='FORMATTED', &
	594	POSITION='APPEND' )
	595	ELSE
	596	IF ( myid == 0 .AND. .NOT. openfile(80)%opened_before ) THEN
	597	CALL local_system( 'mkdir PARTICLE_INFOS' )
	598	ENDIF
	599	#if defined( __parallel )
	600	!
	601	!-- Set a barrier in order to allow that thereafter all other
	602	!-- processors in the directory created by PE0 can open their file.
	603	!-- WARNING: The following barrier will lead to hanging jobs, if
	604	!-- check_open is first called from routine
	605	!-- allocate_prt_memory!
	606	IF ( .NOT. openfile(80)%opened_before ) THEN
	607	CALL MPI_BARRIER( comm2d, ierr )
	608	ENDIF
	609	#endif
	610	OPEN ( 80, FILE='PARTICLE_INFOS/'//myid_char, FORM='FORMATTED', &
	611	POSITION='APPEND' )
	612	ENDIF
	613
	614	IF ( .NOT. openfile(80)%opened_before ) THEN
	615	WRITE ( 80, 8000 ) TRIM( run_description_header )
	616	ENDIF
	617
	618	CASE ( 81 )
	619
	620	OPEN ( 81, FILE='PLOTSP_X_PAR', FORM='FORMATTED', &
	621	DELIM='APOSTROPHE', RECL=1500, POSITION='APPEND' )
	622
	623	CASE ( 82 )
	624
	625	OPEN ( 82, FILE='PLOTSP_X_DATA', FORM='FORMATTED', &
	626	POSITION = 'APPEND' )
	627
	628	CASE ( 83 )
	629
	630	OPEN ( 83, FILE='PLOTSP_Y_PAR', FORM='FORMATTED', &
	631	DELIM='APOSTROPHE', RECL=1500, POSITION='APPEND' )
	632
	633	CASE ( 84 )
	634
	635	OPEN ( 84, FILE='PLOTSP_Y_DATA', FORM='FORMATTED', &
	636	POSITION='APPEND' )
	637
	638	CASE ( 85 )
	639
	640	IF ( myid_char == '' ) THEN
	641	OPEN ( 85, FILE='PARTICLE_DATA'//myid_char, FORM='UNFORMATTED', &
	642	POSITION='APPEND' )
	643	ELSE
	644	IF ( myid == 0 .AND. .NOT. openfile(85)%opened_before ) THEN
	645	CALL local_system( 'mkdir PARTICLE_DATA' )
	646	ENDIF
	647	#if defined( __parallel )
	648	!
	649	!-- Set a barrier in order to allow that thereafter all other
	650	!-- processors in the directory created by PE0 can open their file
	651	CALL MPI_BARRIER( comm2d, ierr )
	652	#endif
	653	OPEN ( 85, FILE='PARTICLE_DATA/'//myid_char, FORM='UNFORMATTED', &
	654	POSITION='APPEND' )
	655	ENDIF
	656
	657	IF ( .NOT. openfile(85)%opened_before ) THEN
	658	WRITE ( 85 ) run_description_header
	659	!
	660	!-- Attention: change version number whenever the output format on
	661	!-- unit 85 is changed (see also in routine advec_particles)
	662	rtext = 'data format version 3.0'
	663	WRITE ( 85 ) rtext
	664	WRITE ( 85 ) number_of_particle_groups, &
	665	max_number_of_particle_groups
	666	WRITE ( 85 ) particle_groups
	667	ENDIF
	668
	669	#if defined( __netcdf )
	670	CASE ( 101, 111 )
	671	!
	672	!-- Set filename depending on unit number
	673	IF ( file_id == 101 ) THEN
	674	filename = 'DATA_2D_XY_NETCDF'
	675	av = 0
	676	ELSE
	677	filename = 'DATA_2D_XY_AV_NETCDF'
	678	av = 1
	679	ENDIF
	680	!
	681	!-- Inquire, if there is a NetCDF file from a previuos run. This should
	682	!-- be opened for extension, if its dimensions and variables match the
	683	!-- actual run.
	684	INQUIRE( FILE=filename, EXIST=netcdf_extend )
	685
	686	IF ( netcdf_extend ) THEN
	687	!
	688	!-- Open an existing NetCDF file for output
	689	nc_stat = NF90_OPEN( filename, NF90_WRITE, id_set_xy(av) )
	690	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 20 )
	691	!
	692	!-- Read header information and set all ids. If there is a mismatch
	693	!-- between the previuos and the actual run, netcdf_extend is returned
	694	!-- as .FALSE.
	695	CALL define_netcdf_header( 'xy', netcdf_extend, av )
	696
	697	!
	698	!-- Remove the local file, if it can not be extended
	699	IF ( .NOT. netcdf_extend ) THEN
	700	nc_stat = NF90_CLOSE( id_set_xy(av) )
	701	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 21 )
	702	CALL local_system( 'rm ' // TRIM( filename ) )
	703	ENDIF
	704
	705	ENDIF
	706
	707	IF ( .NOT. netcdf_extend ) THEN
	708	!
	709	!-- Create a new NetCDF output file
	710	IF ( netcdf_64bit ) THEN
	711	#if defined( __netcdf_64bit )
	712	nc_stat = NF90_CREATE( filename, &
	713	OR( NF90_NOCLOBBER, NF90_64BIT_OFFSET ),&
	714	id_set_xy(av) )
	715	#else
	716	IF ( myid == 0 ) PRINT*, '+++ WARNING: NetCDF: no 64-bit ', &
	717	'offset allowed on this machine'
	718	nc_stat = NF90_CREATE( filename, NF90_NOCLOBBER, id_set_xy(av) )
	719	#endif
	720	ELSE
	721	nc_stat = NF90_CREATE( filename, NF90_NOCLOBBER, id_set_xy(av) )
	722	ENDIF
	723	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 22 )
	724	!
	725	!-- Define the header
	726	CALL define_netcdf_header( 'xy', netcdf_extend, av )
	727
	728	ENDIF
	729
	730	CASE ( 102, 112 )
	731	!
	732	!-- Set filename depending on unit number
	733	IF ( file_id == 102 ) THEN
	734	filename = 'DATA_2D_XZ_NETCDF'
	735	av = 0
	736	ELSE
	737	filename = 'DATA_2D_XZ_AV_NETCDF'
	738	av = 1
	739	ENDIF
	740	!
	741	!-- Inquire, if there is a NetCDF file from a previuos run. This should
	742	!-- be opened for extension, if its dimensions and variables match the
	743	!-- actual run.
	744	INQUIRE( FILE=filename, EXIST=netcdf_extend )
	745
	746	IF ( netcdf_extend ) THEN
	747	!
	748	!-- Open an existing NetCDF file for output
	749	nc_stat = NF90_OPEN( filename, NF90_WRITE, id_set_xz(av) )
	750	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 23 )
	751	!
	752	!-- Read header information and set all ids. If there is a mismatch
	753	!-- between the previuos and the actual run, netcdf_extend is returned
	754	!-- as .FALSE.
	755	CALL define_netcdf_header( 'xz', netcdf_extend, av )
	756
	757	!
	758	!-- Remove the local file, if it can not be extended
	759	IF ( .NOT. netcdf_extend ) THEN
	760	nc_stat = NF90_CLOSE( id_set_xz(av) )
	761	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 24 )
	762	CALL local_system( 'rm ' // TRIM( filename ) )
	763	ENDIF
	764
	765	ENDIF
	766
	767	IF ( .NOT. netcdf_extend ) THEN
	768	!
	769	!-- Create a new NetCDF output file
	770	IF ( netcdf_64bit ) THEN
	771	#if defined( __netcdf_64bit )
	772	nc_stat = NF90_CREATE( filename, &
	773	OR( NF90_NOCLOBBER, NF90_64BIT_OFFSET ),&
	774	id_set_xz(av) )
	775	#else
	776	IF ( myid == 0 ) PRINT*, '+++ WARNING: NetCDF: no 64-bit ', &
	777	'offset allowed on this machine'
	778	nc_stat = NF90_CREATE( filename, NF90_NOCLOBBER, id_set_xz(av) )
	779	#endif
	780	ELSE
	781	nc_stat = NF90_CREATE( filename, NF90_NOCLOBBER, id_set_xz(av) )
	782	ENDIF
	783	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 25 )
	784	!
	785	!-- Define the header
	786	CALL define_netcdf_header( 'xz', netcdf_extend, av )
	787
	788	ENDIF
	789
	790	CASE ( 103, 113 )
	791	!
	792	!-- Set filename depending on unit number
	793	IF ( file_id == 103 ) THEN
	794	filename = 'DATA_2D_YZ_NETCDF'
	795	av = 0
	796	ELSE
	797	filename = 'DATA_2D_YZ_AV_NETCDF'
	798	av = 1
	799	ENDIF
	800	!
	801	!-- Inquire, if there is a NetCDF file from a previuos run. This should
	802	!-- be opened for extension, if its dimensions and variables match the
	803	!-- actual run.
	804	INQUIRE( FILE=filename, EXIST=netcdf_extend )
	805
	806	IF ( netcdf_extend ) THEN
	807	!
	808	!-- Open an existing NetCDF file for output
	809	nc_stat = NF90_OPEN( filename, NF90_WRITE, id_set_yz(av) )
	810	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 26 )
	811	!
	812	!-- Read header information and set all ids. If there is a mismatch
	813	!-- between the previuos and the actual run, netcdf_extend is returned
	814	!-- as .FALSE.
	815	CALL define_netcdf_header( 'yz', netcdf_extend, av )
	816
	817	!
	818	!-- Remove the local file, if it can not be extended
	819	IF ( .NOT. netcdf_extend ) THEN
	820	nc_stat = NF90_CLOSE( id_set_yz(av) )
	821	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 27 )
	822	CALL local_system( 'rm ' // TRIM( filename ) )
	823	ENDIF
	824
	825	ENDIF
	826
	827	IF ( .NOT. netcdf_extend ) THEN
	828	!
	829	!-- Create a new NetCDF output file
	830	IF ( netcdf_64bit ) THEN
	831	#if defined( __netcdf_64bit )
	832	nc_stat = NF90_CREATE( filename, &
	833	OR( NF90_NOCLOBBER, NF90_64BIT_OFFSET), &
	834	id_set_yz(av) )
	835	#else
	836	IF ( myid == 0 ) PRINT*, '+++ WARNING: NetCDF: no 64-bit ', &
	837	'offset allowed on this machine'
	838	nc_stat = NF90_CREATE( filename, NF90_NOCLOBBER, id_set_yz(av) )
	839	#endif
	840	ELSE
	841	nc_stat = NF90_CREATE( filename, NF90_NOCLOBBER, id_set_yz(av) )
	842	ENDIF
	843	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 28 )
	844	!
	845	!-- Define the header
	846	CALL define_netcdf_header( 'yz', netcdf_extend, av )
	847
	848	ENDIF
	849
	850	CASE ( 104 )
	851	!
	852	!-- Inquire, if there is a NetCDF file from a previuos run. This should
	853	!-- be opened for extension, if its variables match the actual run.
	854	INQUIRE( FILE='DATA_1D_PR_NETCDF', EXIST=netcdf_extend )
	855
	856	IF ( netcdf_extend ) THEN
	857	!
	858	!-- Open an existing NetCDF file for output
	859	nc_stat = NF90_OPEN( 'DATA_1D_PR_NETCDF', NF90_WRITE, id_set_pr )
	860	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 29 )
	861	!
	862	!-- Read header information and set all ids. If there is a mismatch
	863	!-- between the previuos and the actual run, netcdf_extend is returned
	864	!-- as .FALSE.
	865	CALL define_netcdf_header( 'pr', netcdf_extend, 0 )
	866
	867	!
	868	!-- Remove the local file, if it can not be extended
	869	IF ( .NOT. netcdf_extend ) THEN
	870	nc_stat = NF90_CLOSE( id_set_pr )
	871	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 30 )
	872	CALL local_system( 'rm DATA_1D_PR_NETCDF' )
	873	ENDIF
	874
	875	ENDIF
	876
	877	IF ( .NOT. netcdf_extend ) THEN
	878	!
	879	!-- Create a new NetCDF output file
	880	IF ( netcdf_64bit ) THEN
	881	#if defined( __netcdf_64bit )
	882	nc_stat = NF90_CREATE( 'DATA_1D_PR_NETCDF', &
	883	OR( NF90_NOCLOBBER, NF90_64BIT_OFFSET ),&
	884	id_set_pr )
	885	#else
	886	IF ( myid == 0 ) PRINT*, '+++ WARNING: NetCDF: no 64-bit ', &
	887	'offset allowed on this machine'
	888	nc_stat = NF90_CREATE( 'DATA_1D_PR_NETCDF', NF90_NOCLOBBER, &
	889	id_set_pr )
	890	#endif
	891	ELSE
	892	nc_stat = NF90_CREATE( 'DATA_1D_PR_NETCDF', NF90_NOCLOBBER, &
	893	id_set_pr )
	894	ENDIF
	895	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 31 )
	896	!
	897	!-- Define the header
	898	CALL define_netcdf_header( 'pr', netcdf_extend, 0 )
	899
	900	ENDIF
	901
	902	CASE ( 105 )
	903	!
	904	!-- Inquire, if there is a NetCDF file from a previuos run. This should
	905	!-- be opened for extension, if its variables match the actual run.
	906	INQUIRE( FILE='DATA_1D_TS_NETCDF', EXIST=netcdf_extend )
	907
	908	IF ( netcdf_extend ) THEN
	909	!
	910	!-- Open an existing NetCDF file for output
	911	nc_stat = NF90_OPEN( 'DATA_1D_TS_NETCDF', NF90_WRITE, id_set_ts )
	912	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 32 )
	913	!
	914	!-- Read header information and set all ids. If there is a mismatch
	915	!-- between the previuos and the actual run, netcdf_extend is returned
	916	!-- as .FALSE.
	917	CALL define_netcdf_header( 'ts', netcdf_extend, 0 )
	918
	919	!
	920	!-- Remove the local file, if it can not be extended
	921	IF ( .NOT. netcdf_extend ) THEN
	922	nc_stat = NF90_CLOSE( id_set_ts )
	923	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 33 )
	924	CALL local_system( 'rm DATA_1D_TS_NETCDF' )
	925	ENDIF
	926
	927	ENDIF
	928
	929	IF ( .NOT. netcdf_extend ) THEN
	930	!
	931	!-- Create a new NetCDF output file
	932	IF ( netcdf_64bit ) THEN
	933	#if defined( __netcdf_64bit )
	934	nc_stat = NF90_CREATE( 'DATA_1D_TS_NETCDF', &
	935	OR( NF90_NOCLOBBER, NF90_64BIT_OFFSET ),&
	936	id_set_ts )
	937	#else
	938	IF ( myid == 0 ) PRINT*, '+++ WARNING: NetCDF: no 64-bit ', &
	939	'offset allowed on this machine'
	940	nc_stat = NF90_CREATE( 'DATA_1D_TS_NETCDF', NF90_NOCLOBBER, &
	941	id_set_ts )
	942	#endif
	943	ELSE
	944	nc_stat = NF90_CREATE( 'DATA_1D_TS_NETCDF', NF90_NOCLOBBER, &
	945	id_set_ts )
	946	ENDIF
	947	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 34 )
	948	!
	949	!-- Define the header
	950	CALL define_netcdf_header( 'ts', netcdf_extend, 0 )
	951
	952	ENDIF
	953
	954
	955	CASE ( 106, 116 )
	956	!
	957	!-- Set filename depending on unit number
	958	IF ( file_id == 106 ) THEN
	959	filename = 'DATA_3D_NETCDF'
	960	av = 0
	961	ELSE
	962	filename = 'DATA_3D_AV_NETCDF'
	963	av = 1
	964	ENDIF
	965	!
	966	!-- Inquire, if there is a NetCDF file from a previuos run. This should
	967	!-- be opened for extension, if its dimensions and variables match the
	968	!-- actual run.
	969	INQUIRE( FILE=filename, EXIST=netcdf_extend )
	970
	971	IF ( netcdf_extend ) THEN
	972	!
	973	!-- Open an existing NetCDF file for output
	974	nc_stat = NF90_OPEN( filename, NF90_WRITE, id_set_3d(av) )
	975	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 35 )
	976	!
	977	!-- Read header information and set all ids. If there is a mismatch
	978	!-- between the previuos and the actual run, netcdf_extend is returned
	979	!-- as .FALSE.
	980	CALL define_netcdf_header( '3d', netcdf_extend, av )
	981
	982	!
	983	!-- Remove the local file, if it can not be extended
	984	IF ( .NOT. netcdf_extend ) THEN
	985	nc_stat = NF90_CLOSE( id_set_3d(av) )
	986	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 36 )
	987	CALL local_system('rm ' // TRIM( filename ) )
	988	ENDIF
	989
	990	ENDIF
	991
	992	IF ( .NOT. netcdf_extend ) THEN
	993	!
	994	!-- Create a new NetCDF output file
[46]	995	IF ( netcdf_64bit .AND. netcdf_64bit_3d ) THEN
[1]	996	#if defined( __netcdf_64bit )
	997	nc_stat = NF90_CREATE( filename, &
	998	OR( NF90_NOCLOBBER, NF90_64BIT_OFFSET ),&
	999	id_set_3d(av) )
	1000	#else
	1001	IF ( myid == 0 ) PRINT*, '+++ WARNING: NetCDF: no 64-bit ', &
	1002	'offset allowed on this machine'
	1003	nc_stat = NF90_CREATE( filename, NF90_NOCLOBBER, id_set_3d(av) )
	1004	#endif
	1005	ELSE
	1006	nc_stat = NF90_CREATE( filename, NF90_NOCLOBBER, id_set_3d(av) )
	1007	ENDIF
	1008	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 37 )
	1009	!
	1010	!-- Define the header
	1011	CALL define_netcdf_header( '3d', netcdf_extend, av )
	1012
	1013	ENDIF
	1014
	1015
	1016	CASE ( 107 )
	1017	!
	1018	!-- Inquire, if there is a NetCDF file from a previuos run. This should
	1019	!-- be opened for extension, if its variables match the actual run.
	1020	INQUIRE( FILE='DATA_1D_SP_NETCDF', EXIST=netcdf_extend )
	1021
	1022	IF ( netcdf_extend ) THEN
	1023	!
	1024	!-- Open an existing NetCDF file for output
	1025	nc_stat = NF90_OPEN( 'DATA_1D_SP_NETCDF', NF90_WRITE, id_set_sp )
	1026	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 38 )
	1027	!
	1028	!-- Read header information and set all ids. If there is a mismatch
	1029	!-- between the previuos and the actual run, netcdf_extend is returned
	1030	!-- as .FALSE.
	1031	CALL define_netcdf_header( 'sp', netcdf_extend, 0 )
	1032
	1033	!
	1034	!-- Remove the local file, if it can not be extended
	1035	IF ( .NOT. netcdf_extend ) THEN
	1036	nc_stat = NF90_CLOSE( id_set_sp )
	1037	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 39 )
	1038	CALL local_system( 'rm DATA_1D_SP_NETCDF' )
	1039	ENDIF
	1040
	1041	ENDIF
	1042
	1043	IF ( .NOT. netcdf_extend ) THEN
	1044	!
	1045	!-- Create a new NetCDF output file
	1046	IF ( netcdf_64bit ) THEN
	1047	#if defined( __netcdf_64bit )
	1048	nc_stat = NF90_CREATE( 'DATA_1D_SP_NETCDF', &
	1049	OR( NF90_NOCLOBBER, NF90_64BIT_OFFSET ),&
	1050	id_set_sp )
	1051	#else
	1052	IF ( myid == 0 ) PRINT*, '+++ WARNING: NetCDF: no 64-bit ', &
	1053	'offset allowed on this machine'
	1054	nc_stat = NF90_CREATE( 'DATA_1D_SP_NETCDF', NF90_NOCLOBBER, &
	1055	id_set_sp )
	1056	#endif
	1057	ELSE
	1058	nc_stat = NF90_CREATE( 'DATA_1D_SP_NETCDF', NF90_NOCLOBBER, &
	1059	id_set_sp )
	1060	ENDIF
	1061	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 40 )
	1062	!
	1063	!-- Define the header
	1064	CALL define_netcdf_header( 'sp', netcdf_extend, 0 )
	1065
	1066	ENDIF
	1067
	1068
	1069	CASE ( 108 )
	1070
	1071	IF ( myid_char == '' ) THEN
	1072	filename = 'DATA_PRT_NETCDF'
	1073	ELSE
	1074	filename = 'DATA_PRT_NETCDF/' // myid_char
	1075	ENDIF
	1076	!
	1077	!-- Inquire, if there is a NetCDF file from a previuos run. This should
	1078	!-- be opened for extension, if its variables match the actual run.
	1079	INQUIRE( FILE=filename, EXIST=netcdf_extend )
	1080
	1081	IF ( netcdf_extend ) THEN
	1082	!
	1083	!-- Open an existing NetCDF file for output
	1084	nc_stat = NF90_OPEN( filename, NF90_WRITE, id_set_prt )
	1085	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 41 )
	1086	!
	1087	!-- Read header information and set all ids. If there is a mismatch
	1088	!-- between the previuos and the actual run, netcdf_extend is returned
	1089	!-- as .FALSE.
	1090	CALL define_netcdf_header( 'pt', netcdf_extend, 0 )
	1091
	1092	!
	1093	!-- Remove the local file, if it can not be extended
	1094	IF ( .NOT. netcdf_extend ) THEN
	1095	nc_stat = NF90_CLOSE( id_set_prt )
	1096	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 42 )
	1097	CALL local_system( 'rm ' // filename )
	1098	ENDIF
	1099
	1100	ENDIF
	1101
	1102	IF ( .NOT. netcdf_extend ) THEN
	1103
	1104	!
	1105	!-- For runs on multiple processors create the subdirectory
	1106	IF ( myid_char /= '' ) THEN
	1107	IF ( myid == 0 .AND. .NOT. openfile(file_id)%opened_before ) &
	1108	THEN ! needs modification in case of non-extendable sets
	1109	CALL local_system( 'mkdir DATA_PRT_NETCDF/' )
	1110	ENDIF
	1111	#if defined( __parallel )
	1112	!
	1113	!-- Set a barrier in order to allow that all other processors in the
	1114	!-- directory created by PE0 can open their file
	1115	CALL MPI_BARRIER( comm2d, ierr )
	1116	#endif
	1117	ENDIF
	1118
	1119	!
	1120	!-- Create a new NetCDF output file
	1121	IF ( netcdf_64bit ) THEN
	1122	#if defined( __netcdf_64bit )
	1123	nc_stat = NF90_CREATE( filename, &
	1124	OR( NF90_NOCLOBBER, NF90_64BIT_OFFSET ),&
	1125	id_set_prt )
	1126	#else
	1127	IF ( myid == 0 ) PRINT*, '+++ WARNING: NetCDF: no 64-bit ', &
	1128	'offset allowed on this machine'
	1129	nc_stat = NF90_CREATE( filename, NF90_NOCLOBBER, id_set_prt )
	1130	#endif
	1131	ELSE
	1132	nc_stat = NF90_CREATE( filename, NF90_NOCLOBBER, id_set_prt )
	1133	ENDIF
	1134	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 43 )
	1135
	1136	!
	1137	!-- Define the header
	1138	CALL define_netcdf_header( 'pt', netcdf_extend, 0 )
	1139
	1140	ENDIF
	1141
	1142	CASE ( 109 )
	1143	!
	1144	!-- Inquire, if there is a NetCDF file from a previuos run. This should
	1145	!-- be opened for extension, if its variables match the actual run.
	1146	INQUIRE( FILE='DATA_1D_PTS_NETCDF', EXIST=netcdf_extend )
	1147
	1148	IF ( netcdf_extend ) THEN
	1149	!
	1150	!-- Open an existing NetCDF file for output
	1151	nc_stat = NF90_OPEN( 'DATA_1D_PTS_NETCDF', NF90_WRITE, id_set_pts )
	1152	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 393 )
	1153	!
	1154	!-- Read header information and set all ids. If there is a mismatch
	1155	!-- between the previuos and the actual run, netcdf_extend is returned
	1156	!-- as .FALSE.
	1157	CALL define_netcdf_header( 'ps', netcdf_extend, 0 )
	1158
	1159	!
	1160	!-- Remove the local file, if it can not be extended
	1161	IF ( .NOT. netcdf_extend ) THEN
	1162	nc_stat = NF90_CLOSE( id_set_pts )
	1163	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 394 )
	1164	CALL local_system( 'rm DATA_1D_PTS_NETCDF' )
	1165	ENDIF
	1166
	1167	ENDIF
	1168
	1169	IF ( .NOT. netcdf_extend ) THEN
	1170	!
	1171	!-- Create a new NetCDF output file
	1172	IF ( netcdf_64bit ) THEN
	1173	#if defined( __netcdf_64bit )
	1174	nc_stat = NF90_CREATE( 'DATA_1D_PTS_NETCDF', &
	1175	OR( NF90_NOCLOBBER, NF90_64BIT_OFFSET ),&
	1176	id_set_pts )
	1177	#else
	1178	IF ( myid == 0 ) PRINT*, '+++ WARNING: NetCDF: no 64-bit ', &
	1179	'offset allowed on this machine'
	1180	nc_stat = NF90_CREATE( 'DATA_1D_PTS_NETCDF', NF90_NOCLOBBER, &
	1181	id_set_pts )
	1182	#endif
	1183	ELSE
	1184	nc_stat = NF90_CREATE( 'DATA_1D_PTS_NETCDF', NF90_NOCLOBBER, &
	1185	id_set_pts )
	1186	ENDIF
	1187	IF ( nc_stat /= NF90_NOERR ) CALL handle_netcdf_error( 395 )
	1188	!
	1189	!-- Define the header
	1190	CALL define_netcdf_header( 'ps', netcdf_extend, 0 )
	1191
	1192	ENDIF
	1193	#else
	1194
	1195	CASE ( 101:109, 111:113, 116 )
	1196
	1197	!
	1198	!-- Nothing is done in case of missing netcdf support
	1199	RETURN
	1200
	1201	#endif
	1202
	1203	CASE DEFAULT
	1204
	1205	PRINT*,'+++ check_open: no OPEN-statement for file-id ',file_id
	1206	#if defined( __parallel )
	1207	CALL MPI_ABORT( comm2d, 9999, ierr )
	1208	#else
	1209	CALL local_stop
	1210	#endif
	1211
	1212	END SELECT
	1213
	1214	!
	1215	!-- Set open flag
	1216	openfile(file_id)%opened = .TRUE.
	1217
	1218	!
	1219	!-- Formats
	1220	3300 FORMAT ('#'/ &
	1221	'coord 1 file=',A,' filetype=unformatted'/ &
	1222	'coord 2 file=',A,' filetype=unformatted skip=',I6/ &
	1223	'coord 3 file=',A,' filetype=unformatted skip=',I6/ &
	1224	'#')
	1225	4000 FORMAT ('# ',A)
	1226	5000 FORMAT ('# ',A/ &
	1227	'#1 E'/'#2 E'/'#3 dt'/'#4 u'/'#5 th*'/'#6 umax'/'#7 vmax'/ &
	1228	'#8 wmax'/'#9 div_new'/'#10 div_old'/'#11 z_i_wpt'/'#12 z_i_pt'/ &
	1229	'#13 w*'/'#14 w''pt''0'/'#15 w''pt'''/'#16 wpt'/'#17 pt(0)'/ &
	1230	'#18 pt(zp)'/'#19 splptx'/'#20 splpty'/'#21 splptz')
	1231	8000 FORMAT (A/ &
	1232	' step time # of parts lPE sent/recv rPE sent/recv ',&
	1233	'sPE sent/recv nPE sent/recv max # of parts'/ &
	1234	103('-'))
	1235
	1236	END SUBROUTINE check_open

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

Context Navigation

source: palm/trunk/SOURCE/check_open.f90 @ 80

Download in other formats: