Home

Context Navigation

source: palm/trunk/SOURCE/data_output_profiles.f90 @ 285

Last change on this file since 285 was 274, checked in by heinze, 15 years ago
Indentation of the message calls corrected
Property svn:keywords set to `Id`
File size: 24.9 KB

Rev	Line
[1]	1	SUBROUTINE data_output_profiles
	2
	3	!------------------------------------------------------------------------------!
[254]	4	! Current revisions:
[1]	5	! -----------------
[263]	6	! Output of NetCDF messages with aid of message handling routine.
[254]	7	! Output of messages replaced by message handling routine.
[1]	8	!
[254]	9	!
[1]	10	! Former revisions:
	11	! -----------------
[3]	12	! $Id: data_output_profiles.f90 274 2009-03-26 15:11:21Z raasch $
[198]	13	!
	14	! 197 2008-09-16 15:29:03Z raasch
	15	! Time coordinate t=0 stored on netcdf-file only if an output is required for
	16	! this time for at least one of the profiles
	17	!
	18	! February 2007
[3]	19	! RCS Log replace by Id keyword, revision history cleaned up
	20	!
[90]	21	! 87 2007-05-22 15:46:47Z raasch
	22	! var_hom renamed pr_palm
	23	!
[1]	24	! Revision 1.18 2006/08/16 14:27:04 raasch
	25	! PRINT* statements for testing removed
	26	!
	27	! Revision 1.1 1997/09/12 06:28:48 raasch
	28	! Initial revision
	29	!
	30	!
	31	! Description:
	32	! ------------
	33	! Plot output of 1D-profiles for PROFIL
	34	!------------------------------------------------------------------------------!
	35
	36	USE control_parameters
	37	USE cpulog
	38	USE indices
	39	USE interfaces
	40	USE netcdf_control
	41	USE pegrid
	42	USE profil_parameter
	43	USE statistics
	44
	45	IMPLICIT NONE
	46
	47
	48	INTEGER :: i, id, ilc, ils, j, k, sr
[197]	49	LOGICAL :: output_for_t0
[1]	50	REAL :: uxma, uxmi
	51
	52
	53	!
	54	!-- If required, compute statistics
	55	IF ( .NOT. flow_statistics_called ) CALL flow_statistics
	56
	57	!
	58	!-- Flow_statistics has its own CPU time measurement
	59	CALL cpu_log( log_point(15), 'data_output_profiles', 'start' )
	60
	61	!
	62	!-- If required, compute temporal average
	63	IF ( averaging_interval_pr == 0.0 ) THEN
	64	hom_sum(:,:,:) = hom(:,1,:,:)
	65	ELSE
	66	IF ( average_count_pr > 0 ) THEN
	67	hom_sum = hom_sum / REAL( average_count_pr )
	68	ELSE
	69	!
	70	!-- This case may happen if dt_dopr is changed in the d3par-list of
	71	!-- a restart run
	72	RETURN
	73	ENDIF
	74	ENDIF
	75
	76
	77	IF ( myid == 0 ) THEN
	78
	79	!
	80	!-- Plot-output for each (sub-)region
	81
	82	!
	83	!-- Open file for profile output in NetCDF format
	84	IF ( netcdf_output ) THEN
	85	CALL check_open( 104 )
	86	ENDIF
	87
	88	!
	89	!-- Open PROFIL-output files for each (sub-)region
	90	IF ( profil_output ) THEN
	91	DO sr = 0, statistic_regions
	92	CALL check_open( 40 + sr )
	93	ENDDO
	94	ENDIF
	95
	96	!
	97	!-- Increment the counter for number of output times
	98	dopr_time_count = dopr_time_count + 1
	99
	100	!
	101	!-- Re-set to zero the counter for the number of profiles already written
	102	!-- at the current output time into the respective crosses
	103	cross_pnc_local = 0
	104
	105	!
	106	!-- Output of initial profiles
	107	IF ( dopr_time_count == 1 ) THEN
	108
	109	IF ( netcdf_output ) THEN
	110	#if defined( __netcdf )
	111	!
[197]	112	!-- Store initial time (t=0) to time axis, but only if an output
	113	!-- is required for at least one of the profiles
	114	output_for_t0 = .FALSE.
	115	DO i = 1, dopr_n
	116	IF ( dopr_initial_index(i) /= 0 ) THEN
	117	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_time_pr, &
	118	(/ 0.0 /), start = (/ 1 /), &
	119	count = (/ 1 /) )
[263]	120	CALL handle_netcdf_error( 'data_output_profiles', 329 )
[197]	121	output_for_t0 = .TRUE.
	122	EXIT
	123	ENDIF
	124	ENDDO
[1]	125
	126	!
	127	!-- Store normalization factors
	128	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(1), & ! wpt0
	129	(/ hom_sum(nzb,18,normalizing_region) /), &
	130	start = (/ 1 /), count = (/ 1 /) )
[263]	131	CALL handle_netcdf_error( 'data_output_profiles', 330 )
[1]	132
	133	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(2), & ! ws2
[87]	134	(/ hom_sum(nzb+8,pr_palm,normalizing_region)**2 /), &
[1]	135	start = (/ 1 /), count = (/ 1 /) )
[263]	136	CALL handle_netcdf_error( 'data_output_profiles', 331 )
[1]	137	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(3), & ! tsw2
[87]	138	(/ hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), &
[1]	139	start = (/ 1 /), count = (/ 1 /) )
[263]	140	CALL handle_netcdf_error( 'data_output_profiles', 332 )
[1]	141	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(4), & ! ws3
[87]	142	(/ hom_sum(nzb+8,pr_palm,normalizing_region)**3 /), &
[1]	143	start = (/ 1 /), count = (/ 1 /) )
[263]	144	CALL handle_netcdf_error( 'data_output_profiles', 333 )
[1]	145
	146	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(5), & ! ws2tsw
[87]	147	(/ hom_sum(nzb+8,pr_palm,normalizing_region)*3 &
	148	hom_sum(nzb+3,pr_palm,normalizing_region) /), &
[1]	149	start = (/ 1 /), count = (/ 1 /) )
[263]	150	CALL handle_netcdf_error( 'data_output_profiles', 334 )
[1]	151
	152	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(6), & ! wstsw2
[87]	153	(/ hom_sum(nzb+8,pr_palm,normalizing_region) * &
	154	hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), &
[1]	155	start = (/ 1 /), count = (/ 1 /) )
[263]	156	CALL handle_netcdf_error( 'data_output_profiles', 335 )
[1]	157
	158	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(7), & ! z_i
[87]	159	(/ hom_sum(nzb+6,pr_palm,normalizing_region) /), &
[1]	160	start = (/ 1 /), count = (/ 1 /) )
[263]	161	CALL handle_netcdf_error( 'data_output_profiles', 336 )
[1]	162	#endif
	163	ENDIF
	164	!
	165	!-- Loop over all 1D variables
	166	DO i = 1, dopr_n
	167
	168	IF ( dopr_initial_index(i) /= 0 ) THEN
	169
	170	!
	171	!-- Output for the individual (sub-)regions
	172	DO sr = 0, statistic_regions
	173
	174	IF ( profil_output ) THEN
	175	id = 40 + sr
	176	!
	177	!-- Write Label-Header
	178	WRITE ( id, 100 ) TRIM( data_output_pr(i) ), '(t=0)'
	179	!
	180	!-- Write total profile
	181	DO k = nzb, nzt+1
	182	WRITE ( id, 101 ) hom(k,2,dopr_initial_index(i),sr), &
	183	hom(k,1,dopr_initial_index(i),sr)
	184	ENDDO
	185	!
	186	!-- Write separation label
	187	WRITE ( id, 102 )
	188	ENDIF
	189
	190	IF ( netcdf_output ) THEN
	191	#if defined( __netcdf )
	192	!
	193	!-- Write data to netcdf file
	194	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_dopr(i,sr), &
	195	hom(nzb:nzt+1,1,dopr_initial_index(i),sr), &
	196	start = (/ 1, 1 /), &
	197	count = (/ nzt-nzb+2, 1 /) )
[263]	198	CALL handle_netcdf_error( 'data_output_profiles', 337 )
[1]	199	#endif
	200	ENDIF
	201
	202	ENDDO
	203
	204	IF ( profil_output ) THEN
	205	!
	206	!-- Determine indices for later NAMELIST-output (s. below)
	207	profile_number = profile_number + 1
	208	j = dopr_crossindex(i)
	209	IF ( j /= 0 ) THEN
	210	cross_profile_number_count(j) = &
	211	cross_profile_number_count(j) + 1
	212	k = cross_profile_number_count(j)
	213	cross_profile_numbers(k,j) = profile_number
	214	!
	215	!-- Initial profiles are always drawn as solid lines in
	216	!-- anti-background colour.
	217	cross_linecolors(k,j) = 1
	218	cross_linestyles(k,j) = 0
	219	!
	220	!-- If required, extend x-value range of the respective
	221	!-- cross, provided it has not been specified in &
	222	!-- check_parameters. Determination over all (sub-)regions.
	223	IF ( cross_uxmin(j) == 0.0 .AND. &
	224	cross_uxmax(j) == 0.0 ) THEN
	225
	226	DO sr = 0, statistic_regions
	227
	228	uxmi = &
	229	MINVAL( hom(:nz_do1d,1,dopr_initial_index(i),sr) )
	230
	231	uxma = &
	232	MAXVAL( hom(:nz_do1d,1,dopr_initial_index(i),sr) )
	233	!
	234	!-- When the value range of the first line in the
	235	!-- corresponding cross is determined, its value range
	236	!-- is simply adopted.
	237	IF ( cross_uxmin_computed(j) > &
	238	cross_uxmax_computed(j) ) THEN
	239	cross_uxmin_computed(j) = uxmi
	240	cross_uxmax_computed(j) = uxma
	241	ELSE
	242	cross_uxmin_computed(j) = &
	243	MIN( cross_uxmin_computed(j), uxmi )
	244	cross_uxmax_computed(j) = &
	245	MAX( cross_uxmax_computed(j), uxma )
	246	ENDIF
	247
	248	ENDDO
	249
	250	ENDIF
	251	!
	252	!-- If required, determine and note normalizing factors
	253	SELECT CASE ( cross_normalized_x(j) )
	254
	255	CASE ( 'ts2' )
	256	cross_normx_factor(k,j) = &
[87]	257	( hom_sum(nzb+3,pr_palm,normalizing_region) )**2
[1]	258	CASE ( 'wpt0' )
	259	cross_normx_factor(k,j) = &
	260	hom_sum(nzb,18,normalizing_region)
	261	CASE ( 'wsts2' )
	262	cross_normx_factor(k,j) = &
[87]	263	hom_sum(nzb+8,pr_palm,normalizing_region) &
	264	* ( hom_sum(nzb+3,pr_palm,normalizing_region) )**2
[1]	265	CASE ( 'ws2' )
	266	cross_normx_factor(k,j) = &
[87]	267	( hom_sum(nzb+8,pr_palm,normalizing_region) )**2
[1]	268	CASE ( 'ws2ts' )
	269	cross_normx_factor(k,j) = &
[87]	270	( hom_sum(nzb+8,pr_palm,normalizing_region) )**2 &
	271	* hom_sum(nzb+3,pr_palm,normalizing_region)
[1]	272	CASE ( 'ws3' )
	273	cross_normx_factor(k,j) = &
[87]	274	( hom_sum(nzb+8,pr_palm,normalizing_region) )**3
[1]	275
	276	END SELECT
	277
	278	SELECT CASE ( cross_normalized_y(j) )
	279
	280	CASE ( 'z_i' )
	281	cross_normy_factor(k,j) = &
[87]	282	hom_sum(nzb+6,pr_palm,normalizing_region)
[1]	283
	284	END SELECT
	285
	286	!
	287	!-- Check the normalizing factors for zeros and deactivate
	288	!-- the normalization, if required.
	289	IF ( cross_normx_factor(k,j) == 0.0 .OR. &
	290	cross_normy_factor(k,j) == 0.0 ) THEN
[274]	291	WRITE( message_string, * ) 'data_output_profiles: ', &
	292	'normalizing cross ',j, ' is not possible ', &
	293	'since one of the & normalizing factors ', &
	294	'is zero! & cross_normx_factor(',k,',',j,') = ', &
	295	cross_normx_factor(k,j), &
	296	' & cross_normy_factor(',k,',',j,') = ', &
	297	cross_normy_factor(k,j)
	298	CALL message( 'data_output_profiles', 'PA0185',&
	299	0, 1, 0, 6, 0 )
[1]	300	cross_normx_factor(k,j) = 1.0
	301	cross_normy_factor(k,j) = 1.0
	302	cross_normalized_x(j) = ' '
	303	cross_normalized_y(j) = ' '
	304	ENDIF
	305
	306	!
	307	!-- If required, extend normalized x-value range of the
	308	!-- respective cross, provided it has not been specified in
	309	!-- check_parameters. Determination over all (sub-)regions.
	310	IF ( cross_uxmin_normalized(j) == 0.0 .AND. &
	311	cross_uxmax_normalized(j) == 0.0 ) THEN
	312
	313	DO sr = 0, statistic_regions
	314
	315	uxmi = MINVAL( hom(:nz_do1d,1, &
	316	dopr_initial_index(i),sr) ) / &
	317	cross_normx_factor(k,j)
	318	uxma = MAXVAL( hom(:nz_do1d,1, &
	319	dopr_initial_index(i),sr) ) / &
	320	cross_normx_factor(k,j)
	321	!
	322	!-- When the value range of the first line in the
	323	!-- corresponding cross is determined, its value range
	324	!-- is simply adopted.
	325	IF ( cross_uxmin_normalized_computed(j) > &
	326	cross_uxmax_normalized_computed(j) ) THEN
	327	cross_uxmin_normalized_computed(j) = uxmi
	328	cross_uxmax_normalized_computed(j) = uxma
	329	ELSE
	330	cross_uxmin_normalized_computed(j) = &
	331	MIN( cross_uxmin_normalized_computed(j), uxmi )
	332	cross_uxmax_normalized_computed(j) = &
	333	MAX( cross_uxmax_normalized_computed(j), uxma )
	334	ENDIF
	335
	336	ENDDO
	337
	338	ENDIF
	339
	340	ENDIF ! Index determination
	341
	342	ENDIF ! profil output
	343
	344	ENDIF ! Initial profile available
	345
	346	ENDDO ! Loop over dopr_n for initial profiles
	347
[197]	348	IF ( netcdf_output .AND. output_for_t0 ) THEN
	349	dopr_time_count = dopr_time_count + 1
	350	ENDIF
[1]	351
	352	ENDIF ! Initial profiles
	353
	354	IF ( netcdf_output ) THEN
	355	#if defined( __netcdf )
	356	!
	357	!-- Store time to time axis
	358	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_time_pr, &
	359	(/ simulated_time /), &
	360	start = (/ dopr_time_count /), &
	361	count = (/ 1 /) )
[263]	362	CALL handle_netcdf_error( 'data_output_profiles', 338 )
[1]	363
	364	!
	365	!-- Store normalization factors
	366	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(1), & ! wpt0
	367	(/ hom_sum(nzb,18,normalizing_region) /), &
	368	start = (/ dopr_time_count /), &
	369	count = (/ 1 /) )
[263]	370	CALL handle_netcdf_error( 'data_output_profiles', 339 )
[1]	371
	372	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(2), & ! ws2
[87]	373	(/ hom_sum(nzb+8,pr_palm,normalizing_region)**2 /), &
[1]	374	start = (/ dopr_time_count /), &
	375	count = (/ 1 /) )
[263]	376	CALL handle_netcdf_error( 'data_output_profiles', 340 )
[1]	377
	378	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(3), & ! tsw2
[87]	379	(/ hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), &
[1]	380	start = (/ dopr_time_count /), &
	381	count = (/ 1 /) )
[263]	382	CALL handle_netcdf_error( 'data_output_profiles', 341 )
[1]	383
	384	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(4), & ! ws3
[87]	385	(/ hom_sum(nzb+8,pr_palm,normalizing_region)**3 /), &
[1]	386	start = (/ dopr_time_count /), &
	387	count = (/ 1 /) )
[263]	388	CALL handle_netcdf_error( 'data_output_profiles', 342 )
[1]	389
	390	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(5), & ! ws2tsw
[87]	391	(/ hom_sum(nzb+8,pr_palm,normalizing_region)*3 &
	392	hom_sum(nzb+3,pr_palm,normalizing_region) /), &
[1]	393	start = (/ dopr_time_count /), &
	394	count = (/ 1 /) )
[263]	395	CALL handle_netcdf_error( 'data_output_profiles', 343 )
	396
[1]	397	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(6), & ! wstsw2
[87]	398	(/ hom_sum(nzb+8,pr_palm,normalizing_region) * &
	399	hom_sum(nzb+3,pr_palm,normalizing_region)**2 /), &
[1]	400	start = (/ dopr_time_count /), &
	401	count = (/ 1 /) )
[263]	402	CALL handle_netcdf_error( 'data_output_profiles', 344 )
[1]	403
	404	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_norm_dopr(7), & ! z_i
[87]	405	(/ hom_sum(nzb+6,pr_palm,normalizing_region) /), &
[1]	406	start = (/ dopr_time_count /), &
	407	count = (/ 1 /) )
[263]	408	CALL handle_netcdf_error( 'data_output_profiles', 345 )
[1]	409	#endif
	410	ENDIF
	411
	412	!
	413	!-- Output of the individual (non-initial) profiles
	414	DO i = 1, dopr_n
	415
	416	!
	417	!-- Output for the individual (sub-)domains
	418	DO sr = 0, statistic_regions
	419
	420	IF ( profil_output ) THEN
	421	id = 40 + sr
	422	!
	423	!-- Write Label-Header
	424	WRITE ( id, 100 ) TRIM( dopr_label(i) ), simulated_time_chr
	425	!
	426	!-- Output of total profile
	427	DO k = nzb, nzt+1
	428	WRITE ( id, 101 ) hom(k,2,dopr_index(i),sr), &
	429	hom_sum(k,dopr_index(i),sr)
	430	ENDDO
	431	!
	432	!-- Write separation label
	433	WRITE ( id, 102 )
	434	ENDIF
	435
	436	IF ( netcdf_output ) THEN
	437	#if defined( __netcdf )
	438	!
	439	!-- Write data to netcdf file
	440	nc_stat = NF90_PUT_VAR( id_set_pr, id_var_dopr(i,sr), &
	441	hom_sum(nzb:nzt+1,dopr_index(i),sr),&
	442	start = (/ 1, dopr_time_count /), &
	443	count = (/ nzt-nzb+2, 1 /) )
[263]	444	CALL handle_netcdf_error( 'data_output_profiles', 346 )
[1]	445	#endif
	446	ENDIF
	447
	448	ENDDO
	449
	450	IF ( profil_output ) THEN
	451	!
	452	!-- Determine profile number on file and note the data for later
	453	!-- NAMELIST output, if the respective profile is to be drawn by
	454	!-- PROFIL (if it shall not be drawn, the variable dopr_crossindex has
	455	!-- the value 0, otherwise the number of the coordinate cross)
	456	profile_number = profile_number + 1
	457	j = dopr_crossindex(i)
	458
	459	IF ( j /= 0 ) THEN
	460	cross_profile_number_count(j) = cross_profile_number_count(j) +1
	461	k = cross_profile_number_count(j)
	462	cross_pnc_local(j) = cross_pnc_local(j) +1
	463	cross_profile_numbers(k,j) = profile_number
	464	ilc = MOD( dopr_time_count, 10 )
	465	IF ( ilc == 0 ) ilc = 10
	466	cross_linecolors(k,j) = linecolors(ilc)
	467	ils = MOD( cross_pnc_local(j), 11 )
	468	IF ( ils == 0 ) ils = 11
	469	cross_linestyles(k,j) = linestyles(ils)
	470	!
	471	!-- If required, extend x-value range of the respective coordinate
	472	!-- cross, provided it has not been specified in check_parameters.
	473	!-- Determination over all (sub-)regions.
	474	IF ( cross_uxmin(j) == 0.0 .AND. cross_uxmax(j) == 0.0 ) THEN
	475
	476	DO sr = 0, statistic_regions
	477
	478	uxmi = MINVAL( hom_sum(:nz_do1d,dopr_index(i),sr) )
	479	uxma = MAXVAL( hom_sum(:nz_do1d,dopr_index(i),sr) )
	480	!
	481	!-- When the value range of the first line in the
	482	!-- corresponding cross is determined, its value range is
	483	!-- simply adopted.
	484	IF ( cross_uxmin_computed(j) > cross_uxmax_computed(j) ) &
	485	THEN
	486	cross_uxmin_computed(j) = uxmi
	487	cross_uxmax_computed(j) = uxma
	488	ELSE
	489	cross_uxmin_computed(j) = &
	490	MIN( cross_uxmin_computed(j), uxmi )
	491	cross_uxmax_computed(j) = &
	492	MAX( cross_uxmax_computed(j), uxma )
	493	ENDIF
	494
	495	ENDDO
	496
	497	ENDIF
	498	!
	499	!-- If required, store the normalizing factors
	500	SELECT CASE ( cross_normalized_x(j) )
	501
	502	CASE ( 'tsw2' )
	503	cross_normx_factor(k,j) = &
[87]	504	( hom_sum(nzb+11,pr_palm,normalizing_region) )**2
[1]	505	CASE ( 'wpt0' )
	506	cross_normx_factor(k,j) = &
	507	hom_sum(nzb,18,normalizing_region)
	508	CASE ( 'wstsw2' )
	509	cross_normx_factor(k,j) = &
[87]	510	hom_sum(nzb+8,pr_palm,normalizing_region) &
	511	* ( hom_sum(nzb+11,pr_palm,normalizing_region) )**2
[1]	512	CASE ( 'ws2' )
	513	cross_normx_factor(k,j) = &
[87]	514	( hom_sum(nzb+8,pr_palm,normalizing_region) )**2
[1]	515	CASE ( 'ws2tsw' )
	516	cross_normx_factor(k,j) = &
[87]	517	( hom_sum(nzb+8,pr_palm,normalizing_region) )**2&
	518	* hom_sum(nzb+11,pr_palm,normalizing_region)
[1]	519	CASE ( 'ws3' )
	520	cross_normx_factor(k,j) = &
[87]	521	( hom_sum(nzb+8,pr_palm,normalizing_region) )**3
[1]	522
	523	END SELECT
	524	SELECT CASE ( cross_normalized_y(j) )
	525
	526	CASE ( 'z_i' )
	527	cross_normy_factor(k,j) = &
[87]	528	hom_sum(nzb+6,pr_palm,normalizing_region)
[1]	529
	530	END SELECT
	531
	532	!
[254]	533	!-- Check the normalizing factors for zeros and deactivate
	534	!-- the normalization, if required.
[1]	535	IF ( cross_normx_factor(k,j) == 0.0 .OR. &
	536	cross_normy_factor(k,j) == 0.0 ) THEN
[274]	537	WRITE( message_string, * ) 'data_output_profiles: ', &
	538	'normalizing cross ',j, ' is not possible ', &
	539	'since one of the & normalizing factors ', &
	540	'is zero! & cross_normx_factor(',k,',',j,') = ', &
	541	cross_normx_factor(k,j), &
	542	' & cross_normy_factor(',k,',',j,') = ', &
	543	cross_normy_factor(k,j)
	544	CALL message( 'data_output_profiles', 'PA0185',&
	545	0, 1, 0, 6, 0 )
[254]	546	cross_normx_factor(k,j) = 1.0
	547	cross_normy_factor(k,j) = 1.0
	548	cross_normalized_x(j) = ' '
	549	cross_normalized_y(j) = ' '
[1]	550	ENDIF
	551
	552	!
	553	!-- If required, extend normalized x-value range of the respective
	554	!-- cross, provided it has not been specified in check_parameters.
	555	!-- Determination over all (sub-)regions.
	556	IF ( cross_uxmin_normalized(j) == 0.0 .AND. &
	557	cross_uxmax_normalized(j) == 0.0 ) THEN
	558
	559	DO sr = 0, statistic_regions
	560
	561	uxmi = MINVAL( hom(:nz_do1d,1,dopr_index(i),sr) ) / &
	562	cross_normx_factor(k,j)
	563	uxma = MAXVAL( hom(:nz_do1d,1,dopr_index(i),sr) ) / &
	564	cross_normx_factor(k,j)
	565	!
	566	!-- When the value range of the first line in the
	567	!-- corresponding cross is determined, its value range is
	568	!-- simply adopted.
	569	IF ( cross_uxmin_normalized_computed(j) > &
	570	cross_uxmax_normalized_computed(j) ) THEN
	571	cross_uxmin_normalized_computed(j) = uxmi
	572	cross_uxmax_normalized_computed(j) = uxma
	573	ELSE
	574	cross_uxmin_normalized_computed(j) = &
	575	MIN( cross_uxmin_normalized_computed(j), uxmi )
	576	cross_uxmax_normalized_computed(j) = &
	577	MAX( cross_uxmax_normalized_computed(j), uxma )
	578	ENDIF
	579
	580	ENDDO
	581
	582	ENDIF
	583
	584	ENDIF ! Index determination
	585
	586	ENDIF ! profil output
	587
	588	ENDDO ! Loop over dopr_n
	589
	590	ENDIF ! Output on PE0
	591
	592	!
	593	!-- If averaging has been done above, the summation counter must be re-set.
	594	IF ( averaging_interval_pr /= 0.0 ) THEN
	595	average_count_pr = 0
	596	ENDIF
	597
	598	CALL cpu_log( log_point(15), 'data_output_profiles','stop', 'nobarrier' )
	599
	600	!
	601	!-- Formats
	602	100 FORMAT ('#1 ',A,1X,A)
	603	101 FORMAT (E15.7,1X,E15.7)
	604	102 FORMAT ('NEXT')
	605
	606	END SUBROUTINE data_output_profiles

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |