Home

Context Navigation

source: palm/trunk/SOURCE/read_3d_binary.f90 @ 29

Last change on this file since 29 was 19, checked in by raasch, 18 years ago
preliminary version of modified boundary conditions at top
Property svn:keywords set to `Id`
File size: 16.7 KB

Rev	Line
[1]	1	SUBROUTINE read_3d_binary
	2
	3	!------------------------------------------------------------------------------!
	4	! Actual revisions:
	5	! -----------------
[19]	6	! +qswst, qswst_m, tswst, tswst_m
[1]	7	!
	8	! Former revisions:
	9	! -----------------
[3]	10	! $Id: read_3d_binary.f90 19 2007-02-23 04:53:48Z raasch $
	11	! RCS Log replace by Id keyword, revision history cleaned up
	12	!
[1]	13	! Revision 1.4 2006/08/04 15:02:32 raasch
	14	! +iran, iran_part
	15	!
	16	! Revision 1.1 2004/04/30 12:47:27 raasch
	17	! Initial revision
	18	!
	19	!
	20	! Description:
	21	! ------------
	22	! Binary input of variables and arrays from restart file
	23	!------------------------------------------------------------------------------!
	24
	25	USE arrays_3d
	26	USE averaging
	27	USE control_parameters
	28	USE cpulog
	29	USE indices
	30	USE interfaces
	31	USE particle_attributes
	32	USE pegrid
	33	USE profil_parameter
	34	USE random_function_mod
	35	USE statistics
	36
	37	IMPLICIT NONE
	38
	39	CHARACTER (LEN=10) :: binary_version, version_on_file
	40	CHARACTER (LEN=20) :: field_chr
	41	CHARACTER (LEN=10), DIMENSION(:), ALLOCATABLE :: chdum10
	42	CHARACTER (LEN=40), DIMENSION(:), ALLOCATABLE :: chdum40
	43	CHARACTER (LEN=100), DIMENSION(:), ALLOCATABLE :: chdum100
	44
	45	INTEGER :: idum1, myid_on_file, numprocs_on_file, nxl_on_file, &
	46	nxr_on_file, nyn_on_file, nys_on_file, nzb_on_file, nzt_on_file
	47
	48	INTEGER, DIMENSION(:), ALLOCATABLE :: idum
	49
	50	REAL, DIMENSION(:), ALLOCATABLE :: rdum
	51
	52	!
	53	!-- Read data from previous model run. unit 13 already opened in parin
	54	CALL cpu_log( log_point_s(14), 'read_3d_binary', 'start' )
	55
	56	!
	57	!-- First compare the version numbers
	58	READ ( 13 ) version_on_file
	59	binary_version = '3.0'
	60	IF ( TRIM( version_on_file ) /= TRIM( binary_version ) ) THEN
	61	IF ( myid == 0 ) THEN
	62	PRINT*, '+++ init_3d_model: version mismatch concerning data ', &
	63	'from prior run'
	64	PRINT*, ' version on file = "', TRIM( version_on_file ),&
	65	'"'
	66	PRINT*, ' version in program = "', TRIM( binary_version ), &
	67	'"'
	68	ENDIF
	69	CALL local_stop
	70	ENDIF
	71
	72	!
	73	!-- Read and compare number of processors, processor-id and array ranges
	74	READ ( 13 ) numprocs_on_file, myid_on_file, nxl_on_file, nxr_on_file, &
	75	nys_on_file, nyn_on_file, nzb_on_file, nzt_on_file
	76
	77	IF ( numprocs_on_file /= numprocs ) THEN
	78	PRINT*, '+++ init_3d_model: mismatch between actual data and data '
	79	PRINT*, ' from prior run on PE ', myid
	80	PRINT*, ' numprocs on file = ', numprocs_on_file
	81	PRINT*, ' numprocs = ', numprocs
	82	CALL local_stop
	83	ENDIF
	84
	85	IF ( myid_on_file /= myid ) THEN
	86	PRINT*, '+++ init_3d_model: mismatch between actual data and data '
	87	PRINT*, ' from prior run'
	88	PRINT*, ' myid_on_file = ', myid_on_file
	89	PRINT*, ' myid = ', myid
	90	#if defined( __parallel )
	91	CALL MPI_ABORT( comm2d, 9999, ierr )
	92	#else
	93	CALL local_stop
	94	#endif
	95	ENDIF
	96
	97	IF ( nxl_on_file /= nxl ) THEN
	98	PRINT*, '+++ init_3d_model: mismatch between actual data and data '
	99	PRINT*, ' from prior run on PE ', myid
	100	PRINT*, ' nxl on file = ', nxl_on_file
	101	PRINT*, ' nxl = ', nxl
	102	#if defined( __parallel )
	103	CALL MPI_ABORT( comm2d, 9999, ierr )
	104	#else
	105	CALL local_stop
	106	#endif
	107	ENDIF
	108
	109	IF ( nxr_on_file /= nxr ) THEN
	110	PRINT*, '+++ init_3d_model: mismatch between actual data and data '
	111	PRINT*, ' from prior run on PE ', myid
	112	PRINT*, ' nxr on file = ', nxr_on_file
	113	PRINT*, ' nxr = ', nxr
	114	#if defined( __parallel )
	115	CALL MPI_ABORT( comm2d, 9999, ierr )
	116	#else
	117	CALL local_stop
	118	#endif
	119	ENDIF
	120
	121	IF ( nys_on_file /= nys ) THEN
	122	PRINT*, '+++ init_3d_model: mismatch between actual data and data '
	123	PRINT*, ' from prior run on PE ', myid
	124	PRINT*, ' nys on file = ', nys_on_file
	125	PRINT*, ' nys = ', nys
	126	#if defined( __parallel )
	127	CALL MPI_ABORT( comm2d, 9999, ierr )
	128	#else
	129	CALL local_stop
	130	#endif
	131	ENDIF
	132
	133	IF ( nyn_on_file /= nyn ) THEN
	134	PRINT*, '+++ init_3d_model: mismatch between actual data and data '
	135	PRINT*, ' from prior run on PE ', myid
	136	PRINT*, ' nyn on file = ', nyn_on_file
	137	PRINT*, ' nyn = ', nyn
	138	#if defined( __parallel )
	139	CALL MPI_ABORT( comm2d, 9999, ierr )
	140	#else
	141	CALL local_stop
	142	#endif
	143	ENDIF
	144
	145	IF ( nzb_on_file /= nzb ) THEN
	146	PRINT*, '+++ init_3d_model: mismatch between actual data and data '
	147	PRINT*, ' from prior run on PE ', myid
	148	PRINT*, ' nzb on file = ', nzb_on_file
	149	PRINT*, ' nzb = ', nzb
	150	CALL local_stop
	151	ENDIF
	152
	153	IF ( nzt_on_file /= nzt ) THEN
	154	PRINT*, '+++ init_3d_model: mismatch between actual data and data '
	155	PRINT*, ' from prior run on PE ', myid
	156	PRINT*, ' nzt on file = ', nzt_on_file
	157	PRINT*, ' nzt = ', nzt
	158	CALL local_stop
	159	ENDIF
	160
	161	!
	162	!-- Local arrays that may be required for possible temporary information
	163	!-- storage in the following
	164	ALLOCATE( chdum10(crmax), chdum40(crmax), chdum100(crmax), &
	165	idum(100crmax), rdum(100crmax) )
	166
	167	!
	168	!-- Initialize spectra (for the case of just starting spectra computation
	169	!-- in continuation runs)
	170	IF ( dt_dosp /= 9999999.9 ) THEN
	171	spectrum_x = 0.0
	172	spectrum_y = 0.0
	173	ENDIF
	174
	175	!
	176	!-- Read arrays
	177	!-- ATTENTION: If the following read commands have been altered, the
	178	!-- ---------- version number of the variable binary_version must be altered,
	179	!-- too. Furthermore, the output list of arrays in write_3d_binary
	180	!-- must also be altered accordingly.
	181	READ ( 13 ) field_chr
	182	DO WHILE ( TRIM( field_chr ) /= '* end *' )
	183
	184	SELECT CASE ( TRIM( field_chr ) )
	185
	186	CASE ( 'e' )
	187	READ ( 13 ) e
	188	CASE ( 'e_av' )
	189	ALLOCATE( e_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	190	READ ( 13 ) e_av
	191	CASE ( 'e_m' )
	192	READ ( 13 ) e_m
	193	CASE ( 'iran' )
	194	READ ( 13 ) iran, iran_part
	195	CASE ( 'kh' )
	196	READ ( 13 ) kh
	197	CASE ( 'kh_m' )
	198	READ ( 13 ) kh_m
	199	CASE ( 'km' )
	200	READ ( 13 ) km
	201	CASE ( 'km_m' )
	202	READ ( 13 ) km_m
	203	CASE ( 'lwp_av' )
	204	ALLOCATE( lwp_av(nys-1:nyn+1,nxl-1:nxr+1) )
	205	READ ( 13 ) lwp_av
	206	CASE ( 'p' )
	207	READ ( 13 ) p
	208	CASE ( 'p_av' )
	209	ALLOCATE( p_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	210	READ ( 13 ) p_av
	211	CASE ( 'pc_av' )
	212	ALLOCATE( pc_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	213	READ ( 13 ) pc_av
	214	CASE ( 'pr_av' )
	215	ALLOCATE( pr_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	216	READ ( 13 ) pr_av
	217	CASE ( 'pt' )
	218	READ ( 13 ) pt
	219	CASE ( 'pt_av' )
	220	ALLOCATE( pt_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	221	READ ( 13 ) pt_av
	222	CASE ( 'pt_m' )
	223	READ ( 13 ) pt_m
	224	CASE ( 'q' )
	225	READ ( 13 ) q
	226	CASE ( 'q_av' )
	227	ALLOCATE( q_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	228	READ ( 13 ) q_av
	229	CASE ( 'q_m' )
	230	READ ( 13 ) q_m
	231	CASE ( 'ql' )
	232	READ ( 13 ) ql
	233	CASE ( 'ql_av' )
	234	ALLOCATE( ql_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	235	READ ( 13 ) ql_av
	236	CASE ( 'ql_c_av' )
	237	ALLOCATE( ql_c_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	238	READ ( 13 ) ql_c_av
	239	CASE ( 'ql_v_av' )
	240	ALLOCATE( ql_v_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	241	READ ( 13 ) ql_v_av
	242	CASE ( 'ql_vp_av' )
	243	ALLOCATE( ql_vp_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	244	READ ( 13 ) ql_vp_av
	245	CASE ( 'qs' )
	246	READ ( 13 ) qs
	247	CASE ( 'qsws' )
	248	READ ( 13 ) qsws
	249	CASE ( 'qsws_m' )
	250	READ ( 13 ) qsws_m
[19]	251	CASE ( 'qswst' )
	252	READ ( 13 ) qswst
	253	CASE ( 'qswst_m' )
	254	READ ( 13 ) qswst_m
[1]	255	CASE ( 'qv_av' )
	256	ALLOCATE( qv_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	257	READ ( 13 ) qv_av
	258	CASE ( 'random_iv' )
	259	READ ( 13 ) random_iv
	260	READ ( 13 ) random_iy
	261	CASE ( 'rif' )
	262	READ ( 13 ) rif
	263	CASE ( 'rif_m' )
	264	READ ( 13 ) rif_m
	265	CASE ( 's_av' )
	266	ALLOCATE( s_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	267	READ ( 13 ) s_av
	268	CASE ( 'shf' )
	269	READ ( 13 ) shf
	270	CASE ( 'shf_m' )
	271	READ ( 13 ) shf_m
[19]	272	CASE ( 'tswst' )
	273	READ ( 13 ) tswst
	274	CASE ( 'tswst_m' )
	275	READ ( 13 ) tswst_m
[1]	276	CASE ( 'spectrum_x' )
	277	READ ( 13 ) spectrum_x
	278	CASE ( 'spectrum_y' )
	279	READ ( 13 ) spectrum_y
	280	CASE ( 'ts' )
	281	READ ( 13 ) ts
	282	CASE ( 'ts_av' )
	283	ALLOCATE( ts_av(nys-1:nyn+1,nxl-1:nxr+1) )
	284	READ ( 13 ) ts_av
	285	CASE ( 'u' )
	286	READ ( 13 ) u
	287	CASE ( 'u_av' )
	288	ALLOCATE( u_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	289	READ ( 13 ) u_av
	290	CASE ( 'u_m' )
	291	READ ( 13 ) u_m
	292	CASE ( 'us' )
	293	READ ( 13 ) us
	294	CASE ( 'usws' )
	295	READ ( 13 ) usws
	296	CASE ( 'usws_m' )
	297	READ ( 13 ) usws_m
	298	CASE ( 'us_av' )
	299	ALLOCATE( us_av(nys-1:nyn+1,nxl-1:nxr+1) )
	300	READ ( 13 ) us_av
	301	CASE ( 'v' )
	302	READ ( 13 ) v
	303	CASE ( 'volume_flow_area' )
	304	READ ( 13 ) volume_flow_area
	305	CASE ( 'volume_flow_initial' )
	306	READ ( 13 ) volume_flow_initial
	307	CASE ( 'v_av' )
	308	ALLOCATE( v_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	309	READ ( 13 ) v_av
	310	CASE ( 'v_m' )
	311	READ (13 ) v_m
	312	CASE ( 'vpt' )
	313	READ ( 13 ) vpt
	314	CASE ( 'vpt_av' )
	315	ALLOCATE( vpt_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	316	READ ( 13 ) vpt_av
	317	CASE ( 'vpt_m' )
	318	READ ( 13 ) vpt_m
	319	CASE ( 'vsws' )
	320	READ ( 13 ) vsws
	321	CASE ( 'vsws_m' )
	322	READ ( 13 ) vsws_m
	323	CASE ( 'w' )
	324	READ ( 13 ) w
	325	CASE ( 'w_av' )
	326	ALLOCATE( w_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
	327	READ ( 13 ) w_av
	328	CASE ( 'w_m' )
	329	READ ( 13 ) w_m
	330	CASE ( 'z0' )
	331	READ ( 13 ) z0
	332
	333	CASE ( 'cross_linecolors' )
	334	IF ( use_prior_plot1d_parameters ) THEN
	335	READ ( 13 ) cross_linecolors
	336	ELSE
	337	READ ( 13 ) idum
	338	ENDIF
	339	CASE ( 'cross_linestyles' )
	340	IF ( use_prior_plot1d_parameters ) THEN
	341	READ ( 13 ) cross_linestyles
	342	ELSE
	343	READ ( 13 ) idum
	344	ENDIF
	345	CASE ( 'cross_normalized_x' )
	346	IF ( use_prior_plot1d_parameters ) THEN
	347	READ ( 13 ) cross_normalized_x
	348	ELSE
	349	READ ( 13 ) chdum10
	350	ENDIF
	351	CASE ( 'cross_normalized_y' )
	352	IF ( use_prior_plot1d_parameters ) THEN
	353	READ ( 13 ) cross_normalized_y
	354	ELSE
	355	READ ( 13 ) chdum10
	356	ENDIF
	357	CASE ( 'cross_normx_factor' )
	358	IF ( use_prior_plot1d_parameters ) THEN
	359	READ ( 13 ) cross_normx_factor
	360	ELSE
	361	READ ( 13 ) rdum
	362	ENDIF
	363	CASE ( 'cross_normy_factor' )
	364	IF ( use_prior_plot1d_parameters ) THEN
	365	READ ( 13 ) cross_normy_factor
	366	ELSE
	367	READ ( 13 ) rdum
	368	ENDIF
	369	CASE ( 'cross_profiles' )
	370	IF ( use_prior_plot1d_parameters ) THEN
	371	READ ( 13 ) cross_profiles
	372	ELSE
	373	READ ( 13 ) chdum100
	374	ENDIF
	375	CASE ( 'cross_profile_n_coun' )
	376	IF ( use_prior_plot1d_parameters ) THEN
	377	READ ( 13 ) cross_profile_number_count
	378	ELSE
	379	READ ( 13 ) idum(1:crmax)
	380	ENDIF
	381	CASE ( 'cross_profile_number' )
	382	IF ( use_prior_plot1d_parameters ) THEN
	383	READ ( 13 ) cross_profile_numbers
	384	ELSE
	385	READ ( 13 ) idum
	386	ENDIF
	387	CASE ( 'cross_uxmax' )
	388	IF ( use_prior_plot1d_parameters ) THEN
	389	READ ( 13 ) cross_uxmax
	390	ELSE
	391	READ ( 13 ) rdum(1:crmax)
	392	ENDIF
	393	CASE ( 'cross_uxmax_computed' )
	394	IF ( use_prior_plot1d_parameters ) THEN
	395	READ ( 13 ) cross_uxmax_computed
	396	ELSE
	397	READ ( 13 ) rdum(1:crmax)
	398	ENDIF
	399	CASE ( 'cross_uxmax_normaliz' )
	400	IF ( use_prior_plot1d_parameters ) THEN
	401	READ ( 13 ) cross_uxmax_normalized
	402	ELSE
	403	READ ( 13 ) rdum(1:crmax)
	404	ENDIF
	405	CASE ( 'cross_uxmax_norm_com' )
	406	IF ( use_prior_plot1d_parameters ) THEN
	407	READ ( 13 ) cross_uxmax_normalized_computed
	408	ELSE
	409	READ ( 13 ) rdum(1:crmax)
	410	ENDIF
	411	CASE ( 'cross_uxmin' )
	412	IF ( use_prior_plot1d_parameters ) THEN
	413	READ ( 13 ) cross_uxmin
	414	ELSE
	415	READ ( 13 ) rdum(1:crmax)
	416	ENDIF
	417	CASE ( 'cross_uxmin_computed' )
	418	IF ( use_prior_plot1d_parameters ) THEN
	419	READ ( 13 ) cross_uxmin_computed
	420	ELSE
	421	READ ( 13 ) rdum(1:crmax)
	422	ENDIF
	423	CASE ( 'cross_uxmin_normaliz' )
	424	IF ( use_prior_plot1d_parameters ) THEN
	425	READ ( 13 ) cross_uxmin_normalized
	426	ELSE
	427	READ ( 13 ) rdum(1:crmax)
	428	ENDIF
	429	CASE ( 'cross_uxmin_norm_com' )
	430	IF ( use_prior_plot1d_parameters ) THEN
	431	READ ( 13 ) cross_uxmin_normalized_computed
	432	ELSE
	433	READ ( 13 ) rdum(1:crmax)
	434	ENDIF
	435	CASE ( 'cross_uymax' )
	436	IF ( use_prior_plot1d_parameters ) THEN
	437	READ ( 13 ) cross_uymax
	438	ELSE
	439	READ ( 13 ) rdum(1:crmax)
	440	ENDIF
	441	CASE ( 'cross_uymin' )
	442	IF ( use_prior_plot1d_parameters ) THEN
	443	READ ( 13 ) cross_uymin
	444	ELSE
	445	READ ( 13 ) rdum(1:crmax)
	446	ENDIF
	447	CASE ( 'cross_xtext' )
	448	IF ( use_prior_plot1d_parameters ) THEN
	449	READ ( 13 ) cross_xtext
	450	ELSE
	451	READ ( 13 ) chdum40
	452	ENDIF
	453	CASE ( 'dopr_crossindex' )
	454	IF ( use_prior_plot1d_parameters ) THEN
	455	READ ( 13 ) dopr_crossindex
	456	ELSE
	457	READ ( 13 ) idum(1:100)
	458	ENDIF
	459	CASE ( 'dopr_time_count' )
	460	IF ( use_prior_plot1d_parameters ) THEN
	461	READ ( 13 ) dopr_time_count
	462	ELSE
	463	READ ( 13 ) idum1
	464	ENDIF
	465	CASE ( 'hom_sum' )
	466	READ ( 13 ) hom_sum
	467	CASE ( 'profile_columns' )
	468	IF ( use_prior_plot1d_parameters ) THEN
	469	READ ( 13 ) profile_columns
	470	ELSE
	471	READ ( 13 ) idum1
	472	ENDIF
	473	CASE ( 'profile_number' )
	474	IF ( use_prior_plot1d_parameters ) THEN
	475	READ ( 13 ) profile_number
	476	ELSE
	477	READ ( 13 ) idum1
	478	ENDIF
	479	CASE ( 'profile_rows' )
	480	IF ( use_prior_plot1d_parameters ) THEN
	481	READ ( 13 ) profile_rows
	482	ELSE
	483	READ ( 13 ) idum1
	484	ENDIF
	485
	486	CASE DEFAULT
	487	PRINT*, '+++ init_3d_model: unknown field named "', &
	488	TRIM( field_chr ), '" found in'
	489	PRINT*, ' data from prior run on PE ', myid
	490	CALL local_stop
	491
	492	END SELECT
	493	!
	494	!-- Read next character string
	495	READ ( 13 ) field_chr
	496
	497	ENDDO
	498
	499	DEALLOCATE( chdum10, chdum40, chdum100, idum, rdum )
	500
	501	!
	502	!-- End of time measuring for reading binary data
	503	CALL cpu_log( log_point_s(14), 'read_3d_binary', 'stop' )
	504
	505	END SUBROUTINE read_3d_binary

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |