Home

Context Navigation

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

Last change on this file since 513 was 494, checked in by raasch, 15 years ago
last commit documented; configuration example file for netcdf4 added
Property svn:keywords set to `Id`
File size: 53.9 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |