Home

Context Navigation

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

Last change on this file since 357 was 345, checked in by heinze, 15 years ago
In case of restart runs without extension, initial profiles are not written to NetCDF-file anymore.
Property svn:keywords set to `Id`
File size: 25.7 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |