Home

Context Navigation

init_dvrp.f90 @ 149

Last change on this file since 149 was 139, checked in by raasch, 17 years ago

New:
---

Plant canopy model of Watanabe (2004,BLM 112,307-341) added.
It can be switched on by the inipar parameter plant_canopy.
The inipar parameter canopy_mode can be used to prescribe a
plant canopy type. The default case is a homogeneous plant
canopy. Heterogeneous distributions of the leaf area
density and the canopy drag coefficient can be defined in the
new routine user_init_plant_canopy (user_interface).
The inipar parameters lad_surface, lad_vertical_gradient and
lad_vertical_gradient_level can be used in order to
prescribe the vertical profile of leaf area density. The
inipar parameter drag_coefficient determines the canopy
drag coefficient.
Finally, the inipar parameter pch_index determines the
index of the upper boundary of the plant canopy.

Allow new case bc_uv_t = 'dirichlet_0' for channel flow.

For unknown variables (CASE DEFAULT) call new subroutine user_data_output_dvrp

Pressure boundary conditions for vertical walls added to the multigrid solver.
They are applied using new wall flag arrays (wall_flags_..) which are defined
for each grid level. New argument gls added to routine user_init_grid
(user_interface).

Frequence of sorting particles can be controlled with new particles_par
parameter dt_sort_particles. Sorting is moved from the SGS timestep loop in
advec_particles after the end of this loop.

advec_particles, check_parameters, data_output_dvrp, header, init_3d_model, init_grid, init_particles, init_pegrid, modules, package_parin, parin, plant_canopy_model, read_var_list, read_3d_binary, user_interface, write_var_list, write_3d_binary

Changed:

Redefine initial nzb_local as the actual total size of topography (later the
extent of topography in nzb_local is reduced by 1dx at the E topography walls
and by 1dy at the N topography walls to form the basis for nzb_s_inner);
for consistency redefine 'single_building' case.

Vertical profiles now based on nzb_s_inner; they are divided by
ngp_2dh_s_inner (scalars, procucts of scalars) and ngp_2dh (staggered velocity
components and their products, procucts of scalars and velocity components),
respectively.

Allow two instead of one digit to specify isosurface and slicer variables.

Status of 3D-volume NetCDF data file only depends on switch netcdf_64bit_3d (check_open)

prognostic_equations include the respective wall_*flux in the parameter list of
calls of diffusion_s. Same as before, only the values of wall_heatflux(0:4)
can be assigned. At present, wall_humidityflux, wall_qflux, wall_salinityflux,
and wall_scalarflux are kept zero. diffusion_s uses the respective wall_*flux
instead of wall_heatflux. This update serves two purposes:

it avoids errors in calculations with humidity/scalar/salinity and prescribed

non-zero wall_heatflux,

it prepares PALM for a possible assignment of wall fluxes of

humidity/scalar/salinity in a future release.

buoyancy, check_open, data_output_dvrp, diffusion_s, diffusivities, flow_statistics, header, init_3d_model, init_dvrp, init_grid, modules, prognostic_equations

Errors:

Bugfix: summation of sums_l_l in diffusivities.

Several bugfixes in the ocean part: Initial density rho is calculated
(init_ocean). Error in initializing u_init and v_init removed
(check_parameters). Calculation of density flux now starts from
nzb+1 (production_e).

Bugfix: pleft/pright changed to pnorth/psouth in sendrecv of particle tail
numbers along y, small bugfixes in the SGS part (advec_particles)

Bugfix: model_string needed a default value (combine_plot_fields)

Bugfix: wavenumber calculation for even nx in routines maketri (poisfft)

Bugfix: assignment of fluxes at walls

Bugfix: absolute value of f must be used when calculating the Blackadar mixing length (init_1d_model)

advec_particles, check_parameters, combine_plot_fields, diffusion_s, diffusivities, init_ocean, init_1d_model, poisfft, production_e

Property svn:keywords set to Id

File size: 20.6 KB

Rev	Line
[1]	1	SUBROUTINE init_dvrp
	2
	3	!------------------------------------------------------------------------------!
	4	! Actual revisions:
[82]	5	! -----------------
[139]	6	!
[1]	7	! TEST: print* statements
	8	! ToDo: checking of mode_dvrp for legal values is not correct
	9	!
	10	! Former revisions:
	11	! -----------------
[3]	12	! $Id: init_dvrp.f90 139 2007-11-29 09:37:41Z letzel $
[39]	13	!
[139]	14	! 130 2007-11-13 14:08:40Z letzel
	15	! allow two instead of one digit to specify isosurface and slicer variables
	16	! Test output of isosurface on camera file
	17	!
[83]	18	! 82 2007-04-16 15:40:52Z raasch
	19	! Preprocessor strings for different linux clusters changed to "lc",
	20	! routine local_flush is used for buffer flushing
	21	!
[39]	22	! 17 2007-02-19 01:57:39Z raasch
[17]	23	! dvrp_output_local activated for all streams
	24	!
	25	! 13 2007-02-14 12:15:07Z raasch
[3]	26	! RCS Log replace by Id keyword, revision history cleaned up
	27	!
[1]	28	! Revision 1.12 2006/02/23 12:30:22 raasch
	29	! ebene renamed section, pl.. replaced by do..,
	30	!
	31	! Revision 1.1 2000/04/27 06:24:39 raasch
	32	! Initial revision
	33	!
	34	!
	35	! Description:
	36	! ------------
	37	! Initializing actions needed when using dvrp-software
	38	!------------------------------------------------------------------------------!
	39	#if defined( __dvrp_graphics )
	40
	41	USE arrays_3d
	42	USE DVRP
	43	USE dvrp_variables
	44	USE grid_variables
	45	USE indices
	46	USE pegrid
	47	USE control_parameters
	48
	49	IMPLICIT NONE
	50
	51	CHARACTER (LEN=2) :: section_chr
	52	CHARACTER (LEN=80) :: dvrp_file_local
	53	INTEGER :: i, j, k, l, m, pn, tv, vn
	54	LOGICAL :: allocated
	55	REAL :: center(3), distance
	56
	57	REAL, DIMENSION(:,:,:), ALLOCATABLE :: local_pf
	58
	59	TYPE(CSTRING), SAVE :: dvrp_directory_c, dvrp_file_c, &
	60	dvrp_file_local_c,dvrp_host_c, &
	61	dvrp_password_c, dvrp_username_c, name_c
	62
	63	!
	64	!-- Set the maximum time the program can be suspended on user request (by
	65	!-- dvrp steering). This variable is defined in module DVRP.
	66	DVRP_MAX_SUSPEND_TIME = 7200
	67
	68	!
	69	!-- Allocate array holding the names and limits of the steering variables
	70	!-- (must have the same number of elements as array mode_dvrp!)
	71	ALLOCATE( steering_dvrp(10) )
	72
	73	!
	74	!-- Check, if output parameters are given and/or allowed
	75	!-- and set default-values, where necessary
	76	IF ( dvrp_username == ' ' ) THEN
	77	IF ( myid == 0 ) THEN
	78	PRINT*, '+++ init_dvrp: dvrp_username is undefined'
	79	CALL local_stop
	80	ENDIF
	81	ENDIF
	82
	83	IF ( dvrp_output /= 'ftp' .AND. dvrp_output /= 'rtsp' .AND. &
	84	dvrp_output /= 'local' ) THEN
	85	IF ( myid == 0 ) THEN
	86	PRINT*, '+++ init_dvrp: dvrp_output="', dvrp_output, '" not allowed'
	87	CALL local_stop
	88	ENDIF
	89	ENDIF
	90
	91	IF ( dvrp_directory == 'default' ) THEN
	92	dvrp_directory = TRIM( dvrp_username ) // '/' // TRIM( run_identifier )
	93	ENDIF
	94
[13]	95	IF ( dvrp_output /= 'local' ) THEN
	96	IF ( dvrp_file /= 'default' .AND. dvrp_file /= '/dev/null' ) THEN
	97	IF ( myid == 0 ) THEN
	98	PRINT*, '+++ init_dvrp: dvrp_file="', dvrp_file, '" not allowed'
	99	CALL local_stop
	100	ENDIF
[1]	101	ENDIF
	102	ENDIF
	103
	104	!
	105	!-- Strings are assigned to strings of special type which have a CHAR( 0 )
	106	!-- (C end-of-character symbol) at their end. This is needed when strings are
	107	!-- passed to C routines.
	108	dvrp_directory_c = dvrp_directory
	109	dvrp_file_c = dvrp_file
	110	dvrp_host_c = dvrp_host
	111	dvrp_password_c = dvrp_password
	112	dvrp_username_c = dvrp_username
	113
	114	!
	115	!-- Loop over all output modes choosed
	116	m = 1
	117	allocated = .FALSE.
	118	DO WHILE ( mode_dvrp(m) /= ' ' )
	119
	120	!
	121	!-- Check, if mode is allowed
	122	IF ( mode_dvrp(m)(1:10) /= 'isosurface' .AND. &
	123	mode_dvrp(m)(1:6) /= 'slicer' .AND. &
	124	mode_dvrp(m)(1:9) /= 'particles' ) THEN
	125
	126	IF ( myid == 0 ) THEN
	127	PRINT*, '+++ init_dvrp: mode_dvrp="', mode_dvrp, '" not allowed'
	128	ENDIF
	129	CALL local_stop
	130
	131	ENDIF
	132
	133	!
	134	!-- Camera position must be computed and written on file when no dvrp-output
	135	!-- has been generated so far (in former runs)
	136	! IF ( dvrp_filecount == 0 ) THEN
	137	!
	138	!-- Compute center of domain and distance of camera from center
	139	center(1) = ( nx + 1.0 ) * dx * 0.5 * superelevation_x
	140	center(2) = ( ny + 1.0 ) * dy * 0.5 * superelevation_y
	141	center(3) = ( zu(nz_do3d) - zu(nzb) ) * 0.5 * superelevation
	142	distance = 1.5 * MAX( ( nx + 1.0 ) * dx * superelevation_x, &
	143	( ny + 1.0 ) * dy * superelevation_y, &
	144	( zu(nz_do3d) - zu(nzb) ) * superelevation )
	145
	146	!
	147	!-- Write camera position on file
	148	CALL DVRP_INIT( m-1, 0 )
	149
	150	!
	151	!-- Create filename for camera
	152	IF ( dvrp_output == 'rtsp' ) THEN
	153
	154	WRITE ( 9, * ) '*** vor dvrp_output_rtsp'
[82]	155	CALL local_flush( 9 )
[13]	156
[1]	157	dvrp_file = TRIM( mode_dvrp(m) ) // '/camera.dvr'
	158	dvrp_file_c = dvrp_file
	159	CALL DVRP_OUTPUT_RTSP( m-1, dvrp_host_c, dvrp_username_c, &
	160	dvrp_password_c, dvrp_directory_c, &
	161	dvrp_file_c )
	162	WRITE ( 9, * ) '*** nach dvrp_output_rtsp'
[82]	163	CALL local_flush( 9 )
[1]	164
	165	ELSEIF ( dvrp_output == 'ftp' ) THEN
	166
	167	dvrp_file = TRIM( mode_dvrp(m) ) // '.camera.dvr'
	168	dvrp_file_c = dvrp_file
	169	! CALL DVRP_OUTPUT_FTP( m-1, 0, dvrp_host_c, dvrp_username_c, &
	170	! dvrp_password_c, dvrp_directory_c, &
	171	! dvrp_file_c )
	172
	173	ELSE
	174
	175	IF ( dvrp_file(1:9) /= '/dev/null' ) THEN
	176	dvrp_file_local = TRIM( mode_dvrp(m) ) // '.camera.dvr'
	177	dvrp_file_local_c = dvrp_file_local
	178	ELSE
	179	dvrp_file_local_c = dvrp_file_c
	180	ENDIF
[13]	181	CALL DVRP_OUTPUT_LOCAL( m-1, 0, dvrp_file_local_c )
[1]	182
	183	ENDIF
	184
	185	CALL DVRP_CAMERA( m-1, center, distance )
	186	WRITE ( 9, * ) '*** #1'
[82]	187	CALL local_flush( 9 )
[1]	188
	189	!
	190	!-- Define bounding box material and create a bounding box
	191	CALL DVRP_MATERIAL_RGB( m-1, 1, 0.5, 0.5, 0.5, 0.0 )
	192	CALL DVRP_BOUNDINGBOX( m-1, 1, 0.01, 0.0, 0.0, 0.0, &
	193	(nx+1) * dx * superelevation_x, &
	194	(ny+1) * dy * superelevation_y, &
	195	zu(nz_do3d) * superelevation )
	196
	197	CALL DVRP_VISUALIZE( m-1, 0, 0 )
	198	CALL DVRP_EXIT( m-1 )
	199	WRITE ( 9, * ) '*** #2'
[82]	200	CALL local_flush( 9 )
[1]	201
	202
	203	!
	204	!-- Write topography isosurface on file
	205	CALL DVRP_INIT( m-1, 0 )
	206
	207	!
	208	!-- Create filename for buildings
	209	IF ( dvrp_output == 'rtsp' ) THEN
	210
	211	dvrp_file = TRIM( mode_dvrp(m) ) // '/buildings.dvr'
	212	dvrp_file_c = dvrp_file
	213	CALL DVRP_OUTPUT_RTSP( m-1, dvrp_host_c, dvrp_username_c, &
	214	dvrp_password_c, dvrp_directory_c, &
	215	dvrp_file_c )
	216	WRITE ( 9, * ) '*** #3'
[82]	217	CALL local_flush( 9 )
[1]	218
[17]	219	ELSEIF ( dvrp_output == 'ftp' ) THEN
	220
	221	dvrp_file = TRIM( mode_dvrp(m) ) // '.buildings.dvr'
	222	dvrp_file_c = dvrp_file
	223	! CALL DVRP_OUTPUT_FTP( m-1, 0, dvrp_host_c, dvrp_username_c, &
	224	! dvrp_password_c, dvrp_directory_c, &
	225	! dvrp_file_c )
	226
	227	ELSE
	228
	229	IF ( dvrp_file(1:9) /= '/dev/null' ) THEN
	230	dvrp_file_local = TRIM( mode_dvrp(m) ) // '.buildings.dvr'
	231	dvrp_file_local_c = dvrp_file_local
	232	ELSE
	233	dvrp_file_local_c = dvrp_file_c
	234	ENDIF
	235	CALL DVRP_OUTPUT_LOCAL( m-1, 0, dvrp_file_local_c )
	236
[1]	237	ENDIF
	238
	239	!
	240	!-- Determine local gridpoint coordinates
	241	IF ( .NOT. allocated ) THEN
	242	ALLOCATE( xcoor_dvrp(nxl:nxr+1), ycoor_dvrp(nys:nyn+1), &
	243	zcoor_dvrp(nzb:nz_do3d) )
	244	allocated = .TRUE.
	245
	246	DO i = nxl, nxr+1
	247	xcoor_dvrp(i) = i * dx * superelevation_x
	248	ENDDO
	249	DO j = nys, nyn+1
	250	ycoor_dvrp(j) = j * dy * superelevation_y
	251	ENDDO
	252	zcoor_dvrp = zu(nzb:nz_do3d) * superelevation
	253	nx_dvrp = nxr+1 - nxl + 1
	254	ny_dvrp = nyn+1 - nys + 1
	255	nz_dvrp = nz_do3d - nzb + 1
	256	ENDIF
	257
	258	!
	259	!-- Define the grid used by dvrp
	260	CALL DVRP_GRID( m-1, nx_dvrp, ny_dvrp, nz_dvrp, xcoor_dvrp, &
	261	ycoor_dvrp, zcoor_dvrp )
	262	CALL DVRP_MATERIAL_RGB( m-1, 1, 0.8, 0.7, 0.6, 0.0 )
	263	WRITE ( 9, * ) '*** #4'
[82]	264	CALL local_flush( 9 )
[1]	265
	266	!
	267	!-- Compute and plot isosurface in dvr-format
	268	ALLOCATE( local_pf(nxl:nxr+1,nys:nyn+1,nzb:nz_do3d) )
	269	local_pf = 0.0
	270	DO i = nxl, nxr+1
	271	DO j = nys, nyn+1
	272	IF ( nzb_s_inner(j,i) > 0 ) THEN
	273	local_pf(i,j,nzb:nzb_s_inner(j,i)) = 1.0
	274	ENDIF
	275	ENDDO
	276	ENDDO
	277	WRITE ( 9, * ) '*** #4.1'
[82]	278	CALL local_flush( 9 )
[1]	279	CALL DVRP_DATA( m-1, local_pf, 1, nx_dvrp, ny_dvrp, nz_dvrp, &
	280	cyclic_dvrp, cyclic_dvrp, cyclic_dvrp )
	281	WRITE ( 9, * ) '*** #4.2'
[82]	282	CALL local_flush( 9 )
[1]	283	CALL DVRP_THRESHOLD( m-1, 1.0 )
	284	WRITE ( 9, * ) '*** #4.3'
[82]	285	CALL local_flush( 9 )
[1]	286	CALL DVRP_VISUALIZE( m-1, 1, 0 )
	287	WRITE ( 9, * ) '*** #4.4'
[82]	288	CALL local_flush( 9 )
[1]	289
	290	DEALLOCATE( local_pf )
	291
	292	CALL DVRP_EXIT( m-1 )
	293	WRITE ( 9, * ) '*** #5'
[82]	294	CALL local_flush( 9 )
[1]	295
	296	!
	297	!-- Write the surface isosurface on file
	298	CALL DVRP_INIT( m-1, 0 )
	299
	300	!
	301	!-- Create filename for surface
	302	IF ( dvrp_output == 'rtsp' ) THEN
	303
	304	dvrp_file = TRIM( mode_dvrp(m) ) // '/surface.dvr'
	305	dvrp_file_c = dvrp_file
	306	CALL DVRP_OUTPUT_RTSP( m-1, dvrp_host_c, dvrp_username_c, &
	307	dvrp_password_c, dvrp_directory_c, &
	308	dvrp_file_c )
	309	WRITE ( 9, * ) '*** #6'
[82]	310	CALL local_flush( 9 )
[1]	311
[17]	312	ELSEIF ( dvrp_output == 'ftp' ) THEN
	313
	314	dvrp_file = TRIM( mode_dvrp(m) ) // '.surface.dvr'
	315	dvrp_file_c = dvrp_file
	316	! CALL DVRP_OUTPUT_FTP( m-1, 0, dvrp_host_c, dvrp_username_c, &
	317	! dvrp_password_c, dvrp_directory_c, &
	318	! dvrp_file_c )
	319
	320	ELSE
	321
	322	IF ( dvrp_file(1:9) /= '/dev/null' ) THEN
	323	dvrp_file_local = TRIM( mode_dvrp(m) ) // '.surface.dvr'
	324	dvrp_file_local_c = dvrp_file_local
	325	ELSE
	326	dvrp_file_local_c = dvrp_file_c
	327	ENDIF
	328	CALL DVRP_OUTPUT_LOCAL( m-1, 0, dvrp_file_local_c )
	329
[1]	330	ENDIF
	331
	332	!
	333	!-- Determine local gridpoint coordinates
	334	IF ( .NOT. allocated ) THEN
	335	ALLOCATE( xcoor_dvrp(nxl:nxr+1), ycoor_dvrp(nys:nyn+1), &
	336	zcoor_dvrp(nzb:nz_do3d) )
	337	allocated = .TRUE.
	338
	339	DO i = nxl, nxr+1
	340	xcoor_dvrp(i) = i * dx * superelevation_x
	341	ENDDO
	342	DO j = nys, nyn+1
	343	ycoor_dvrp(j) = j * dy * superelevation_y
	344	ENDDO
	345	zcoor_dvrp = zu(nzb:nz_do3d) * superelevation
	346	nx_dvrp = nxr+1 - nxl + 1
	347	ny_dvrp = nyn+1 - nys + 1
	348	nz_dvrp = nz_do3d - nzb + 1
	349	ENDIF
	350
	351	!
	352	!-- Define the grid used by dvrp
	353	CALL DVRP_GRID( m-1, nx_dvrp, ny_dvrp, nz_dvrp, xcoor_dvrp, &
	354	ycoor_dvrp, zcoor_dvrp )
	355	CALL DVRP_MATERIAL_RGB( m-1, 1, 0.0, 0.6, 0.0, 0.0 )
	356	WRITE ( 9, * ) '*** #7'
[82]	357	CALL local_flush( 9 )
[1]	358
	359	!
	360	!-- Compute and plot isosurface in dvr-format
	361	ALLOCATE( local_pf(nxl:nxr+1,nys:nyn+1,nzb:nz_do3d) )
	362	local_pf = 0.0
	363	local_pf(:,:,0) = 1.0
	364
	365	CALL DVRP_DATA( m-1, local_pf, 1, nx_dvrp, ny_dvrp, nz_dvrp, &
	366	cyclic_dvrp, cyclic_dvrp, cyclic_dvrp )
	367	CALL DVRP_THRESHOLD( m-1, 1.0 )
	368	CALL DVRP_VISUALIZE( m-1, 1, 0 )
	369
	370	DEALLOCATE( local_pf )
	371
	372	CALL DVRP_EXIT( m-1 )
	373	WRITE ( 9, * ) '*** #8'
[82]	374	CALL local_flush( 9 )
[1]	375
	376
	377	! ENDIF
	378
	379
	380	!
	381	!-- Initialize dvrp for all dvrp-calls during the run
	382	CALL DVRP_INIT( m-1, 0 )
	383
	384	!
	385	!-- Preliminary definition of filename for dvrp-output
	386	IF ( dvrp_output == 'rtsp' ) THEN
	387
	388	!
	389	!-- First initialize parameters for possible interactive steering.
	390	!-- Every parameter has to be passed to the respective stream.
	391	pn = 1
	392	!
	393	!-- Initialize threshold counter needed for initialization of the
	394	!-- isosurface steering variables
	395	tv = 0
	396
	397	DO WHILE ( mode_dvrp(pn) /= ' ' )
	398
	399	IF ( mode_dvrp(pn)(1:10) == 'isosurface' ) THEN
	400
[130]	401	READ ( mode_dvrp(pn), '(10X,I2)' ) vn
[1]	402	steering_dvrp(pn)%name = do3d(0,vn)
	403	tv = tv + 1
	404
	405	IF ( do3d(0,vn)(1:1) == 'w' ) THEN
	406	steering_dvrp(pn)%min = -4.0
	407	steering_dvrp(pn)%max = 5.0
	408	ELSE
	409	steering_dvrp(pn)%min = 288.0
	410	steering_dvrp(pn)%max = 292.0
	411	ENDIF
	412
	413	name_c = TRIM( do3d(0,vn) )
	414	WRITE ( 9, * ) '*** #9'
[82]	415	CALL local_flush( 9 )
[1]	416	CALL DVRP_STEERING_INIT( m-1, name_c, steering_dvrp(pn)%min, &
	417	steering_dvrp(pn)%max, threshold(tv) )
	418	WRITE ( 9, * ) '*** #10'
[82]	419	CALL local_flush( 9 )
[1]	420
	421	ELSEIF ( mode_dvrp(pn)(1:6) == 'slicer' ) THEN
	422
[130]	423	READ ( mode_dvrp(pn), '(6X,I2)' ) vn
[1]	424	steering_dvrp(pn)%name = do2d(0,vn)
	425	name_c = TRIM( do2d(0,vn) )
	426
	427	l = MAX( 2, LEN_TRIM( do2d(0,vn) ) )
	428	section_chr = do2d(0,vn)(l-1:l)
	429	SELECT CASE ( section_chr )
	430	CASE ( 'xy' )
	431	steering_dvrp(pn)%imin = 0
	432	steering_dvrp(pn)%imax = nz_do3d
	433	slicer_position_dvrp(pn) = section(1,1)
	434	CALL DVRP_STEERING_INIT( m-1, name_c, &
	435	steering_dvrp(pn)%imin, &
	436	steering_dvrp(pn)%imax, &
	437	slicer_position_dvrp(pn) )
	438	CASE ( 'xz' )
	439	steering_dvrp(pn)%imin = 0
	440	steering_dvrp(pn)%imax = ny
	441	slicer_position_dvrp(pn) = section(1,2)
	442	CALL DVRP_STEERING_INIT( m-1, name_c, &
	443	steering_dvrp(pn)%imin, &
	444	steering_dvrp(pn)%imax, &
	445	slicer_position_dvrp(pn) )
	446	CASE ( 'yz' )
	447	steering_dvrp(pn)%imin = 0
	448	steering_dvrp(pn)%imax = nx
	449	slicer_position_dvrp(pn) = section(1,3)
	450	CALL DVRP_STEERING_INIT( m-1, name_c, &
	451	steering_dvrp(pn)%imin, &
	452	steering_dvrp(pn)%imax, &
	453	slicer_position_dvrp(pn) )
	454	END SELECT
	455
	456	ENDIF
	457
	458	pn = pn + 1
	459
	460	ENDDO
	461
	462	WRITE ( 9, * ) '*** #11'
[82]	463	CALL local_flush( 9 )
[1]	464
	465	dvrp_file = TRIM( mode_dvrp(m) ) // '/*****.dvr'
	466	dvrp_file_c = dvrp_file
	467	CALL DVRP_OUTPUT_RTSP( m-1, dvrp_host_c, dvrp_username_c, &
	468	dvrp_password_c, dvrp_directory_c, &
	469	dvrp_file_c )
	470	WRITE ( 9, * ) '*** #12'
[82]	471	CALL local_flush( 9 )
[1]	472
	473	ELSEIF ( dvrp_output == 'ftp' ) THEN
	474
	475	dvrp_file = TRIM( mode_dvrp(m) ) // '.%05d.dvr'
	476	dvrp_file_c = dvrp_file
	477	! CALL DVRP_OUTPUT_FTP( m-1, 0, dvrp_host_c, dvrp_username_c, &
	478	! dvrp_password_c, dvrp_directory_c, dvrp_file_c )
	479
	480	ELSE
	481
	482	IF ( dvrp_file(1:9) /= '/dev/null' ) THEN
[13]	483	dvrp_file_local = TRIM( mode_dvrp(m) ) // '_%05d.dvr'
[1]	484	dvrp_file_local_c = dvrp_file_local
	485	ELSE
	486	dvrp_file_local_c = dvrp_file_c
	487	ENDIF
[13]	488	CALL DVRP_OUTPUT_LOCAL( m-1, 0, dvrp_file_local_c )
[1]	489
	490	ENDIF
	491
	492	! dvrp_file = TRIM( mode_dvrp(m) ) // '.%05d.dvr' // CHAR( 0 )
	493	! dvrp_file = TRIM( mode_dvrp(m) ) // '/*****.dvr' // CHAR( 0 )
	494	! dvrp_file = '/dev/null' // CHAR( 0 )
	495	! CALL DVRP_OUTPUT_FTP( m-1, 0, dvrp_host, dvrp_username, dvrp_password, &
	496	! dvrp_directory, dvrp_file )
	497	! CALL DVRP_OUTPUT_RTSP( m-1, dvrp_host, dvrp_username, dvrp_password, &
	498	! dvrp_directory, dvrp_file )
	499	! CALL DVRP_OUTPUT_LOCAL( m-1, 0, dvrp_file )
	500
	501	!
	502	!-- Determine local gridpoint coordinates
	503	IF ( .NOT. allocated ) THEN
	504	ALLOCATE( xcoor_dvrp(nxl:nxr+1), ycoor_dvrp(nys:nyn+1), &
	505	zcoor_dvrp(nzb:nz_do3d) )
	506	allocated = .TRUE.
	507
	508	DO i = nxl, nxr+1
	509	xcoor_dvrp(i) = i * dx * superelevation_x
	510	ENDDO
	511	DO j = nys, nyn+1
	512	ycoor_dvrp(j) = j * dy * superelevation_y
	513	ENDDO
	514	zcoor_dvrp = zu(nzb:nz_do3d) * superelevation
	515	nx_dvrp = nxr+1 - nxl + 1
	516	ny_dvrp = nyn+1 - nys + 1
	517	nz_dvrp = nz_do3d - nzb + 1
	518	ENDIF
	519
	520	!
	521	!-- Define the grid used by dvrp
	522	WRITE ( 9, * ) '*** #13'
[82]	523	CALL local_flush( 9 )
[1]	524
	525	CALL DVRP_GRID( m-1, nx_dvrp, ny_dvrp, nz_dvrp, xcoor_dvrp, ycoor_dvrp, &
	526	zcoor_dvrp )
	527	WRITE ( 9, * ) '*** #14'
[82]	528	CALL local_flush( 9 )
[1]	529
	530
	531	m = m + 1
	532
	533	ENDDO
	534
	535	#endif
	536	END SUBROUTINE init_dvrp
	537
	538
	539	SUBROUTINE init_dvrp_logging
	540
	541	!------------------------------------------------------------------------------!
	542	! Description:
	543	! ------------
	544	! Initializes logging events for time measurement with dvrp software
	545	! and splits one PE from the global communicator in case that dvrp output
	546	! shall be done by one single PE.
	547	!------------------------------------------------------------------------------!
	548	#if defined( __dvrp_graphics )
	549
	550	USE dvrp_variables
	551	USE pegrid
	552
	553	IMPLICIT NONE
	554
	555	CHARACTER (LEN=4) :: chr
	556	INTEGER :: idummy
	557
	558	!
	559	!-- Initialize logging of calls by DVRP graphic software
	560	WRITE ( 9, * ) '*** myid=', myid, ' vor DVRP_LOG_INIT'
[82]	561	CALL local_flush( 9 )
[1]	562	CALL DVRP_LOG_INIT( 'DVRP_LOG' // CHAR( 0 ), 0 )
	563	WRITE ( 9, * ) '*** myid=', myid, ' nach DVRP_LOG_INIT'
[82]	564	CALL local_flush( 9 )
[1]	565
	566	!
	567	!-- User-defined logging events: #1 (total time needed by PALM)
	568	WRITE ( 9, * ) '*** myid=', myid, ' vor DVRP_LOG_SYMBOL'
[82]	569	CALL local_flush( 9 )
[1]	570	CALL DVRP_LOG_SYMBOL( 1, 'PALM_total' // CHAR( 0 ) )
	571	WRITE ( 9, * ) '*** myid=', myid, ' nach DVRP_LOG_SYMBOL'
[82]	572	CALL local_flush( 9 )
[1]	573	CALL DVRP_LOG_SYMBOL( 2, 'PALM_timestep' // CHAR( 0 ) )
	574	WRITE ( 9, * ) '*** myid=', myid, ' vor DVRP_LOG_EVENT'
[82]	575	CALL local_flush( 9 )
[1]	576	CALL DVRP_LOG_EVENT( 1, 1 )
	577	WRITE ( 9, * ) '*** myid=', myid, ' nach DVRP_LOG_EVENT'
[82]	578	CALL local_flush( 9 )
[1]	579
	580	#if defined( __parallel )
	581	!
	582	!-- Find out, if dvrp output shall be done by a dedicated PE
	583	CALL local_getenv( 'use_seperate_pe_for_dvrp_output', 31, chr, idummy )
	584	IF ( chr == 'true' ) THEN
	585	use_seperate_pe_for_dvrp_output = .TRUE.
	586	WRITE ( 9, * ) '*** myid=', myid, ' vor DVRP_SPLIT'
[82]	587	CALL local_flush( 9 )
[1]	588	CALL DVRP_SPLIT( MPI_COMM_WORLD, comm_palm )
	589	WRITE ( 9, * ) '*** myid=', myid, ' nach DVRP_SPLIT'
[82]	590	CALL local_flush( 9 )
[1]	591	CALL MPI_COMM_SIZE( comm_palm, numprocs, ierr )
	592	ENDIF
	593	#endif
	594
	595	#endif
	596	END SUBROUTINE init_dvrp_logging
	597
	598
	599	SUBROUTINE close_dvrp
	600
	601	!------------------------------------------------------------------------------!
	602	! Description:
	603	! ------------
	604	! Exit of dvrp software and finish dvrp logging
	605	!------------------------------------------------------------------------------!
	606	#if defined( __dvrp_graphics )
	607
	608	USE control_parameters
	609	USE dvrp
	610	USE dvrp_variables
	611
	612	INTEGER :: m
	613
	614	!
	615	!-- If required, close dvrp-software and logging of dvrp-calls
	616	IF ( dt_dvrp /= 9999999.9 ) THEN
	617	m = 1
	618	DO WHILE ( mode_dvrp(m) /= ' ' )
	619	CALL DVRP_EXIT( m-1 )
	620	m = m + 1
	621	ENDDO
	622	CALL DVRP_LOG_EVENT( -1, 1 ) ! Logging of total cpu-time used by PALM
	623	IF ( use_seperate_pe_for_dvrp_output ) THEN
	624	CALL DVRP_SPLIT_EXIT( 1 ) ! Argument 0: reduced output
	625	ELSE
	626	CALL DVRP_LOG_EXIT( 1 ) ! Argument 0: reduced output
	627	ENDIF
	628	ENDIF
	629
	630	#endif
	631	END SUBROUTINE close_dvrp

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

Context Navigation

source: palm/trunk/SOURCE/init_dvrp.f90 @ 149

Download in other formats: