Home

Context Navigation

init_dvrp.f90 @ 114

Last change on this file since 114 was 83, checked in by raasch, 18 years ago

New:
---

Changed:

PALM can be generally installed on any kind of Linux-, IBM-AIX-, or NEC-SX-system by adding appropriate settings to the configuration file.

Scripts are also running under the public domain ksh.

All system relevant compile and link options as well as the host identifier (local_host) are specified in the configuration file.

Filetransfer by ftp removed (options -f removed from mrun and mbuild).

Call of (system-)FLUSH routine moved to new routine local_flush.

return_addres and return_username are read from ENVPAR-NAMELIST-file instead of using local_getenv.

Preprocessor strings for different linux clusters changed to "lc", some preprocessor directives renamed (new: intel_openmp_bug), preprocessor directives for old systems removed

advec_particles, check_open, cpu_log, cpu_statistics, data_output_dvrp, flow_statistics, header, init_dvrp, init_particles, init_1d_model, init_dvrp, init_pegrid, local_getenv, local_system, local_tremain, local_tremain_ini, modules, palm, parin, run_control

new:
local_flush

mbuild, mrun

Errors:

Property svn:keywords set to Id

File size: 20.5 KB

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

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

Context Navigation

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

Download in other formats: