Home

Context Navigation

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

Last change on this file since 320 was 291, checked in by raasch, 16 years ago
changes for coupling with independent precursor runs; z_i calculation with Sullivan criterion
Property svn:keywords set to `Id`
File size: 25.0 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |