Home

Context Navigation

check_open.f90 @ 394

Last change on this file since 394 was 392, checked in by raasch, 15 years ago

New:
---

Adapted for machine lck
(mrun, mbuild, subjob)

bc_lr/bc_ns in most subroutines replaced by LOGICAL variables bc_lr_cyc,
bc_ns_cyc for speed optimization
(check_parameters, diffusion_u, diffusion_v, diffusion_w, modules)

Additional timestep criterion in case of simulations with plant canopy (timestep)

Check for illegal entries in section_xy|xz|yz that exceed nz+1|ny+1|nx+1
(check_parameters)

Clipping of dvrp output implemented. Default colourtable for particles
implemented, particle attributes (color, dvrp_size) can be set with new
parameters particle_color, particle_dvrpsize, color_interval,
dvrpsize_interval (init_dvrp, data_output_dvrp, modules, user_data_output_dvrp).
Slicer attributes (dvrp) are set with new routine set_slicer_attributes_dvrp
and are controlled with existing parameters slicer_range_limits.
(set_slicer_attributes_dvrp)

Ocean atmosphere coupling allows to use independent precursor runs in order
to account for different spin-up times. The time when coupling has to be
started is given by new inipar parameter coupling_start_time. The precursor
ocean run has to be started using new mrun option "-y" in order to add
appendix "_O" to all output files.
(check_for_restart, check_parameters, data_output_2d, data_output_3d,
data_output_profiles, data_output_ptseries, data_output_spectra,
data_output_tseries, header, init_coupling, modules, mrun,
parin, read_var_list, surface_coupler, time_integration, write_var_list)

Polygon reduction for topography and ground plate isosurface. Reduction level
for buildings can be chosen with parameter cluster_size. (init_dvrp)

External pressure gradient (check_parameters, header, init_3d_model, modules,
parin, prognostic_equations, read_var_list, write_var_list)

New topography case 'single_street_canyon' (header, init_grid, modules, parin,
read_var_list, user_check_parameters, user_header, user_init_grid, write_var_list)

Option to predefine a target bulk velocity for conserve_volume_flow
(check_parameters, header, init_3d_model, modules, parin, read_var_list,
write_var_list)

Option for user defined 2D data output in xy cross sections at z=nzb+1
(data_output_2d, user_data_output_2d)

xy cross section output of surface heatfluxes (latent, sensible)
(average_3d_data, check_parameters, data_output_2d, modules, read_3d_binary,
sum_up_3d_data, write_3d_binary)

average_3d_data, check_for_restart, check_parameters, data_output_2d, data_output_3d, data_output_dvrp, data_output_profiles, data_output_ptseries, data_output_spectra, data_output_tseries, init_coupling, init_dvrp, init_grid, init_3d_model, header, mbuild, modules, mrun, package_parin, parin, prognostic_equations, read_3d_binary, read_var_list, subjob, surface_coupler, timestep, time_integration, user_check_parameters, user_data_output_2d, user_data_output_dvrp, user_header, user_init_grid, write_3d_binary, write_var_list

New: set_particle_attributes, set_slicer_attributes_dvrp

Changed:

lcmuk changed to lc to avoid problems with Intel compiler on sgi-ice
(poisfft)

For extended NetCDF files, the updated title attribute includes an update of
time_average_text where appropriate. (netcdf)

In case of restart runs without extension, initial profiles are not written
to NetCDF-file anymore. (data_output_profiles, modules, read_var_list, write_var_list)

Small change in formatting of the message handling routine concering the output in the
job protocoll. (message)

initializing_actions='read_data_for_recycling' renamed to 'cyclic_fill', now
independent of turbulent_inflow (check_parameters, header, init_3d_model)

2 NetCDF error numbers changed. (data_output_3d)

A Link to the website appendix_a.html is printed for further information
about the possible errors. (message)

Temperature gradient criterion for estimating the boundary layer height
replaced by the gradient criterion of Sullivan et al. (1998). (flow_statistics)

NetCDF unit attribute in timeseries output in case of statistic regions added
(netcdf)

Output of NetCDF messages with aid of message handling routine.
(check_open, close_file, data_output_2d, data_output_3d,
data_output_profiles, data_output_ptseries, data_output_spectra,
data_output_tseries, netcdf, output_particles_netcdf)

Output of messages replaced by message handling routine.
(advec_particles, advec_s_bc, buoyancy, calc_spectra, check_for_restart,
check_open, coriolis, cpu_log, data_output_2d, data_output_3d, data_output_dvrp,
data_output_profiles, data_output_spectra, fft_xy, flow_statistics, header,
init_1d_model, init_3d_model, init_dvrp, init_grid, init_particles, init_pegrid,
netcdf, parin, plant_canopy_model, poisfft_hybrid, poismg, read_3d_binary,
read_var_list, surface_coupler, temperton_fft, timestep, user_actions,
user_data_output_dvrp, user_dvrp_coltab, user_init_grid, user_init_plant_canopy,
user_parin, user_read_restart_data, user_spectra )

Maximum number of tails is calculated from maximum number of particles and
skip_particles_for_tail (init_particles)

Value of vertical_particle_advection may differ for each particle group
(advec_particles, header, modules)

First constant in array den also defined as type double. (eqn_state_seawater)

Parameter dvrp_psize moved from particles_par to dvrp_graphics_par. (package_parin)

topography_grid_convention moved from userpar to inipar (check_parameters,
header, parin, read_var_list, user_check_parameters, user_header,
user_init_grid, user_parin, write_var_list)

Default value of grid_matching changed to strict.

Adjustments for runs on lcxt4 (necessary due to an software update on CRAY) and
for coupled runs on ibmy (mrun, subjob)

advec_particles, advec_s_bc, buoyancy, calc_spectra, check_for_restart, check_open, check_parameters, close_file, coriolis, cpu_log, data_output_2d, data_output_3d, data_output_dvrp, data_output_profiles, data_output_ptseries, data_output_spectra, data_output_tseries, eqn_state_seawater, fft_xy, flow_statistics, header, init_1d_model, init_3d_model, init_dvrp, init_grid, init_particles, init_pegrid, message, mrun, netcdf, output_particles_netcdf, package_parin, parin, plant_canopy_model, poisfft, poisfft_hybrid, poismg, read_3d_binary, read_var_list, sort_particles, subjob, user_check_parameters, user_header, user_init_grid, user_parin, surface_coupler, temperton_fft, timestep, user_actions, user_data_output_dvrp, user_dvrp_coltab, user_init_grid, user_init_plant_canopy, user_parin, user_read_restart_data, user_spectra, write_var_list

Errors:

Bugfix: Initial hydrostatic pressure profile in case of ocean runs is now
calculated in 5 iteration steps. (init_ocean)

Bugfix: wrong sign in buoyancy production of ocean part in case of not using
the reference density (only in 3D routine production_e) (production_e)

Bugfix: output of averaged 2d/3d quantities requires that an avaraging
interval has been set, respective error message is included (check_parameters)

Bugfix: Output on unit 14 only if requested by write_binary.
(user_last_actions)

Bugfix to avoid zero division by km_neutral (production_e)

Bugfix for extended NetCDF files: In order to avoid 'data mode' errors if
updated attributes are larger than their original size, NF90_PUT_ATT is called
in 'define mode' enclosed by NF90_REDEF and NF90_ENDDEF calls. This implies a
possible performance loss; an alternative strategy would be to ensure equal
attribute size in a job chain. (netcdf)

Bugfix: correction of initial volume flow for non-flat topography (init_3d_model)
Bugfix: zero initialization of arrays within buildings for 'cyclic_fill' (init_3d_model)

Bugfix: to_be_resorted => s_av for time-averaged scalars (data_output_2d, data_output_3d)

Bugfix: error in formatting the output (message)

Bugfix: avoid that ngp_2dh_s_inner becomes zero (init_3_model)

Typographical error: unit of wpt in dots_unit (modules)

Bugfix: error in check, if particles moved further than one subdomain length.
This check must not be applied for newly released particles. (advec_particles)

Bugfix: several tail counters are initialized, particle_tail_coordinates is
only written to file if its third index is > 0, arrays for tails are allocated
with a minimum size of 10 tails if there is no tail initially (init_particles,
advec_particles)

Bugfix: pressure included for profile output (check_parameters)

Bugfix: Type of count and count_rate changed to default INTEGER on NEC machines
(cpu_log)

Bugfix: output if particle time series only if particle advection is switched
on. (time_integration)

Bugfix: qsws was calculated in case of constant heatflux = .FALSE. (prandtl_fluxes)

Bugfix: averaging along z is not allowed for 2d quantities (e.g. u* and z0) (data_output_2d)

Typographical errors (netcdf)

If the inversion height calculated by the prerun is zero, inflow_damping_height
must be explicitly specified (init_3d_model)

Small bugfix concerning 3d 64bit netcdf output format (header)

Bugfix: dt_fixed removed from the restart file, because otherwise, no change
from a fixed to a variable timestep would be possible in restart runs.
(read_var_list, write_var_list)

Bugfix: initial setting of time_coupling in coupled restart runs (time_integration)

advec_particles, check_parameters, cpu_log, data_output_2d, data_output_3d, header, init_3d_model, init_particles, init_ocean, modules, netcdf, prandtl_fluxes, production_e, read_var_list, time_integration, user_last_actions, write_var_list

Property svn:keywords set to Id

File size: 46.2 KB

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

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

Context Navigation

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

Download in other formats: