Home

Context Navigation

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

Last change on this file since 1 was 1, checked in by raasch, 17 years ago
Initial repository layout and content
File size: 25.7 KB

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

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |