Home

Context Navigation

write_3d_binary.f90 @ 91

Last change on this file since 91 was 77, checked in by raasch, 18 years ago

New:
---

particle reflection from vertical walls implemented, particle SGS model adjusted to walls

Wall functions for vertical walls now include diabatic conditions. New subroutines wall_fluxes, wall_fluxes_e. New 4D-array rif_wall.

new d3par-parameter netcdf_64bit_3d to switch on 64bit offset only for 3D files

new d3par-parameter dt_max to define the maximum value for the allowed timestep

new inipar-parameter loop_optimization to control the loop optimization method

new inipar-parameter pt_refrence. If given, this value is used as the reference that in buoyancy terms (otherwise, the instantaneous horizontally averaged temperature is used).

new user interface user_advec_particles

new initializing action "by_user" calls user_init_3d_model and allows the initial setting of all 3d arrays

topography height informations are stored on arrays zu_s_inner and zw_w_inner and output to the 2d/3d NetCDF files

samples added to the user interface which show how to add user-define time series quantities.

calculation/output of precipitation amount, precipitation rate and z0 (by setting "pra*", "prr*", "z0*" with data_output). The time interval on which the precipitation amount is defined is set by new d3par-parameter precipitation_amount_interval

unit 9 opened for debug output (file DEBUG_<pe#>)

Makefile, advec_particles, average_3d_data, buoyancy, calc_precipitation, check_open, check_parameters, data_output_2d, diffusion_e, diffusion_u, diffusion_v, diffusion_w, diffusivities, header, impact_of_latent_heat, init_particles, init_3d_model, modules, netcdf, parin, production_e, read_var_list, read_3d_binary, sum_up_3d_data, user_interface, write_var_list, write_3d_binary

New: wall_fluxes

Changed:

General revision of non-cyclic horizontal boundary conditions: radiation boundary conditions are now used instead of Neumann conditions at the outflow (calculation needs velocity values for t-dt, which are stored on new arrays u_m_l, u_m_r, etc.), calculation of mean outflow is not needed any more, volume flow control is added for the outflow boundary (currently only for the north boundary!!), additional gridpoints along x and y (uxrp, vynp) are not needed any more, routine "boundary_conds" now operates on timelevel t+dt and is not split in two parts (main, uvw_outflow) any more, Neumann boundary conditions at inflow/outflow in case of non-cyclic boundary conditions for all 2d-arrays that are handled by exchange_horiz_2d

The FFT-method for solving the Poisson-equation is now working with Neumann boundary conditions both at the bottom and the top. This requires adjustments of the tridiagonal coefficients and subtracting the horizontally averaged mean from the vertical velocity field.

+age_m in particle_type

Particles-package is now part of the default code ("-p particles" is not needed any more).

Move call of user_actions( 'after_integration' ) below increment of times
and counters. user_actions is now called for each statistic region and has as an argument the number of the respective region (sr)

d3par-parameter data_output_ts removed. Timeseries output for "profil" removed. Timeseries are now switched on by dt_dots. Timeseries data is collected in flow_statistics.

Initial velocities at nzb+1 are regarded for volume flow control in case they have been set zero before (to avoid small timesteps); see new internal parameters u/v_nzb_p1_for_vfc.

q is not allowed to become negative (prognostic_equations).

poisfft_init is only called if fft-solver is switched on (init_pegrid).

d3par-parameter moisture renamed to humidity.

Subversion global revision number is read from mrun and added to the run description header and to the run control (_rc) file.

vtk directives removed from main program.

The uitility routine interpret_config reads PALM environment variables from NAMELIST instead using the system call GETENV.

advec_u_pw, advec_u_up, advec_v_pw, advec_v_up, asselin_filter, check_parameters, coriolis, data_output_dvrp, data_output_ptseries, data_output_ts, data_output_2d, data_output_3d, diffusion_u, diffusion_v, exchange_horiz, exchange_horiz_2d, flow_statistics, header, init_grid, init_particles, init_pegrid, init_rankine, init_pt_anomaly, init_1d_model, init_3d_model, modules, palm, package_parin, parin, poisfft, poismg, prandtl_fluxes, pres, production_e, prognostic_equations, read_var_list, read_3d_binary, sor, swap_timelevel, time_integration, write_var_list, write_3d_binary

Errors:

Bugfix: preset of tendencies te_em, te_um, te_vm in init_1d_model

Bugfix in sample for reading user defined data from restart file (user_init)

Bugfix in setting diffusivities for cases with the outflow damping layer extending over more than one subdomain (init_3d_model)

Check for possible negative humidities in the initial humidity profile.

in Makefile, default suffixes removed from the suffix list to avoid calling of m2c in
# case of .mod files

Makefile
check_parameters, init_1d_model, init_3d_model, user_interface

Property svn:keywords set to Id

File size: 11.4 KB

Rev	Line
[1]	1	SUBROUTINE write_3d_binary
	2
	3	!------------------------------------------------------------------------------!
	4	! Actual revisions:
	5	! -----------------
[77]	6	!
[1]	7	!
	8	! Former revisions:
	9	! -----------------
[3]	10	! $Id: write_3d_binary.f90 77 2007-03-29 04:26:56Z raasch $
[39]	11	!
[77]	12	! 75 2007-03-22 09:54:05Z raasch
	13	! +precipitation_amount, precipitation_rate_av, rif_wall, u_m_l, u_m_r, etc.,
	14	! z0_av, moisture renamed humidity
	15	!
[39]	16	! 19 2007-02-23 04:53:48Z raasch
	17	! +qswst, qswst_m, tswst, tswst_m
	18	!
[3]	19	! RCS Log replace by Id keyword, revision history cleaned up
	20	!
[1]	21	! Revision 1.21 2006/08/04 15:05:11 raasch
	22	! +iran, iran_part
	23	!
	24	! Revision 1.1 1998/03/18 20:20:21 raasch
	25	! Initial revision
	26	!
	27	!
	28	! Description:
	29	! ------------
	30	! Binary output of variables and arrays for restarts.
	31	!------------------------------------------------------------------------------!
	32
	33	USE arrays_3d
	34	USE averaging
[72]	35	USE cloud_parameters
[1]	36	USE control_parameters
	37	USE cpulog
	38	USE indices
	39	USE interfaces
	40	USE particle_attributes
	41	USE pegrid
	42	USE profil_parameter
	43	USE random_function_mod
	44	USE statistics
	45
	46	IMPLICIT NONE
	47
	48	CHARACTER (LEN=10) :: binary_version
	49
	50
	51	CALL cpu_log( log_point(22), 'write_3d_binary', 'start' )
	52
	53	CALL check_open( 14 )
	54
	55	!
	56	!-- Write control parameters and other variables for restart.
	57	CALL write_var_list
	58
	59	!
	60	!-- Write arrays.
	61	binary_version = '3.0'
	62
	63	WRITE ( 14 ) binary_version
	64
	65	WRITE ( 14 ) numprocs, myid, nxl, nxr, nys, nyn, nzb, nzt
	66
	67	!
	68	!-- Attention: After changes to the following output commands the version number
	69	!-- --------- of the variable binary_version must be changed!
	70	!-- Also, the list of arrays to be read in read_3d_binary must be
	71	!-- adjusted accordingly.
	72
	73	WRITE ( 14 ) 'e '; WRITE ( 14 ) e
	74	IF ( ALLOCATED( e_av ) ) THEN
	75	WRITE ( 14 ) 'e_av '; WRITE ( 14 ) e_av
	76	ENDIF
	77	WRITE ( 14 ) 'e_m '; WRITE ( 14 ) e_m
	78	WRITE ( 14 ) 'iran '; WRITE ( 14 ) iran, iran_part
	79	WRITE ( 14 ) 'kh '; WRITE ( 14 ) kh
	80	IF ( timestep_scheme(1:5) /= 'runge' ) THEN
	81	WRITE ( 14 ) 'kh_m '; WRITE ( 14 ) kh_m
	82	ENDIF
	83	WRITE ( 14 ) 'km '; WRITE ( 14 ) km
	84	IF ( timestep_scheme(1:5) /= 'runge' ) THEN
	85	WRITE ( 14 ) 'km_m '; WRITE ( 14 ) km_m
	86	ENDIF
	87	IF ( ALLOCATED( lwp_av ) ) THEN
	88	WRITE ( 14 ) 'lwp_av '; WRITE ( 14 ) lwp_av
	89	ENDIF
	90	WRITE ( 14 ) 'p '; WRITE ( 14 ) p
	91	IF ( ALLOCATED( p_av ) ) THEN
	92	WRITE ( 14 ) 'p_av '; WRITE ( 14 ) p_av
	93	ENDIF
	94	IF ( ALLOCATED( pc_av ) ) THEN
	95	WRITE ( 14 ) 'pc_av '; WRITE ( 14 ) pc_av
	96	ENDIF
	97	IF ( ALLOCATED( pr_av ) ) THEN
	98	WRITE ( 14 ) 'pr_av '; WRITE ( 14 ) pr_av
	99	ENDIF
[72]	100	IF ( ALLOCATED( precipitation_amount ) ) THEN
	101	WRITE ( 14 ) 'precipitation_amount'; WRITE ( 14 ) precipitation_amount
	102	ENDIF
	103	IF ( ALLOCATED( precipitation_rate_av ) ) THEN
	104	WRITE ( 14 ) 'precipitation_rate_a'; WRITE ( 14 ) &
	105	precipitation_rate_av
	106	ENDIF
[1]	107	WRITE ( 14 ) 'pt '; WRITE ( 14 ) pt
	108	IF ( ALLOCATED( pt_av ) ) THEN
	109	WRITE ( 14 ) 'pt_av '; WRITE ( 14 ) pt_av
	110	ENDIF
	111	WRITE ( 14 ) 'pt_m '; WRITE ( 14 ) pt_m
[75]	112	IF ( humidity .OR. passive_scalar ) THEN
[1]	113	WRITE ( 14 ) 'q '; WRITE ( 14 ) q
	114	IF ( ALLOCATED( q_av ) ) THEN
	115	WRITE ( 14 ) 'q_av '; WRITE ( 14 ) q_av
	116	ENDIF
	117	WRITE ( 14 ) 'q_m '; WRITE ( 14 ) q_m
	118	IF ( cloud_physics ) THEN
	119	WRITE ( 14 ) 'ql '; WRITE ( 14 ) ql
	120	IF ( ALLOCATED( ql_av ) ) THEN
	121	WRITE ( 14 ) 'ql_av '; WRITE ( 14 ) ql_av
	122	ENDIF
	123	ENDIF
	124	WRITE ( 14 ) 'qs '; WRITE ( 14 ) qs
	125	WRITE ( 14 ) 'qsws '; WRITE ( 14 ) qsws
	126	IF ( timestep_scheme(1:5) /= 'runge' ) THEN
	127	WRITE ( 14 ) 'qsws_m '; WRITE ( 14 ) qsws_m
	128	ENDIF
[19]	129	WRITE ( 14 ) 'qswst '; WRITE ( 14 ) qswst
	130	IF ( timestep_scheme(1:5) /= 'runge' ) THEN
	131	WRITE ( 14 ) 'qswst_m '; WRITE ( 14 ) qswst_m
	132	ENDIF
[1]	133	ENDIF
	134	IF ( ALLOCATED( ql_c_av ) ) THEN
	135	WRITE ( 14 ) 'ql_c_av '; WRITE ( 14 ) ql_c_av
	136	ENDIF
	137	IF ( ALLOCATED( ql_v_av ) ) THEN
	138	WRITE ( 14 ) 'ql_v_av '; WRITE ( 14 ) ql_v_av
	139	ENDIF
	140	IF ( ALLOCATED( ql_vp_av ) ) THEN
	141	WRITE ( 14 ) 'ql_vp_av '; WRITE ( 14 ) ql_vp_av
	142	ENDIF
	143	IF ( ALLOCATED( qv_av ) ) THEN
	144	WRITE ( 14 ) 'qv_av '; WRITE ( 14 ) qv_av
	145	ENDIF
	146	WRITE ( 14 ) 'random_iv '; WRITE ( 14 ) random_iv
	147	WRITE ( 14 ) random_iy
	148	WRITE ( 14 ) 'rif '; WRITE ( 14 ) rif
	149	IF ( timestep_scheme(1:5) /= 'runge' ) THEN
	150	WRITE ( 14 ) 'rif_m '; WRITE ( 14 ) rif_m
	151	ENDIF
[51]	152	IF ( topography /= 'flat' ) THEN
	153	WRITE ( 14 ) 'rif_wall '; WRITE ( 14 ) rif_wall
	154	ENDIF
[1]	155	IF ( ALLOCATED( s_av ) ) THEN
	156	WRITE ( 14 ) 's_av '; WRITE ( 14 ) s_av
	157	ENDIF
	158	WRITE ( 14 ) 'shf '; WRITE ( 14 ) shf
	159	IF ( timestep_scheme(1:5) /= 'runge' ) THEN
	160	WRITE ( 14 ) 'shf_m '; WRITE ( 14 ) shf_m
	161	ENDIF
	162	WRITE ( 14 ) 'ts '; WRITE ( 14 ) ts
	163	IF ( ALLOCATED( ts_av ) ) THEN
	164	WRITE ( 14 ) 'ts_av '; WRITE ( 14 ) ts_av
	165	ENDIF
[19]	166	WRITE ( 14 ) 'tswst '; WRITE ( 14 ) tswst
	167	IF ( timestep_scheme(1:5) /= 'runge' ) THEN
	168	WRITE ( 14 ) 'tswst_m '; WRITE ( 14 ) tswst_m
	169	ENDIF
[1]	170	WRITE ( 14 ) 'u '; WRITE ( 14 ) u
	171	IF ( ALLOCATED( u_av ) ) THEN
	172	WRITE ( 14 ) 'u_av '; WRITE ( 14 ) u_av
	173	ENDIF
	174	WRITE ( 14 ) 'u_m '; WRITE ( 14 ) u_m
[73]	175	IF ( ALLOCATED( u_m_l ) ) THEN
	176	WRITE ( 14 ) 'u_m_l '; WRITE ( 14 ) u_m_l
	177	ENDIF
	178	IF ( ALLOCATED( u_m_n ) ) THEN
	179	WRITE ( 14 ) 'u_m_n '; WRITE ( 14 ) u_m_n
	180	ENDIF
	181	IF ( ALLOCATED( u_m_r ) ) THEN
	182	WRITE ( 14 ) 'u_m_r '; WRITE ( 14 ) u_m_r
	183	ENDIF
	184	IF ( ALLOCATED( u_m_s ) ) THEN
	185	WRITE ( 14 ) 'u_m_s '; WRITE ( 14 ) u_m_s
	186	ENDIF
[1]	187	WRITE ( 14 ) 'us '; WRITE ( 14 ) us
	188	WRITE ( 14 ) 'usws '; WRITE ( 14 ) usws
	189	IF ( timestep_scheme(1:5) /= 'runge' ) THEN
	190	WRITE ( 14 ) 'usws_m '; WRITE ( 14 ) usws_m
	191	ENDIF
	192	IF ( ALLOCATED( us_av ) ) THEN
	193	WRITE ( 14 ) 'us_av '; WRITE ( 14 ) us_av
	194	ENDIF
	195	WRITE ( 14 ) 'v '; WRITE ( 14 ) v
	196	WRITE ( 14 ) 'volume_flow_area '; WRITE ( 14 ) volume_flow_area
	197	WRITE ( 14 ) 'volume_flow_initial '; WRITE ( 14 ) volume_flow_initial
	198	IF ( ALLOCATED( v_av ) ) THEN
	199	WRITE ( 14 ) 'v_av '; WRITE ( 14 ) v_av
	200	ENDIF
	201	WRITE ( 14 ) 'v_m '; WRITE ( 14 ) v_m
[73]	202	IF ( ALLOCATED( v_m_l ) ) THEN
	203	WRITE ( 14 ) 'v_m_l '; WRITE ( 14 ) v_m_l
	204	ENDIF
	205	IF ( ALLOCATED( v_m_n ) ) THEN
	206	WRITE ( 14 ) 'v_m_n '; WRITE ( 14 ) v_m_n
	207	ENDIF
	208	IF ( ALLOCATED( v_m_r ) ) THEN
	209	WRITE ( 14 ) 'v_m_r '; WRITE ( 14 ) v_m_r
	210	ENDIF
	211	IF ( ALLOCATED( v_m_s ) ) THEN
	212	WRITE ( 14 ) 'v_m_s '; WRITE ( 14 ) v_m_s
	213	ENDIF
[75]	214	IF ( humidity ) THEN
[1]	215	WRITE ( 14 ) 'vpt '; WRITE ( 14 ) vpt
	216	IF ( ALLOCATED( vpt_av ) ) THEN
	217	WRITE ( 14 ) 'vpt_av '; WRITE ( 14 ) vpt_av
	218	ENDIF
	219	IF ( timestep_scheme(1:5) /= 'runge' ) THEN
	220	WRITE ( 14 ) 'vpt_m '; WRITE ( 14 ) vpt_m
	221	ENDIF
	222	ENDIF
	223	WRITE ( 14 ) 'vsws '; WRITE ( 14 ) vsws
	224	IF ( timestep_scheme(1:5) /= 'runge' ) THEN
	225	WRITE ( 14 ) 'vsws_m '; WRITE ( 14 ) vsws_m
	226	ENDIF
	227	WRITE ( 14 ) 'w '; WRITE ( 14 ) w
	228	IF ( ALLOCATED( w_av ) ) THEN
	229	WRITE ( 14 ) 'w_av '; WRITE ( 14 ) w_av
	230	ENDIF
	231	WRITE ( 14 ) 'w_m '; WRITE ( 14 ) w_m
[73]	232	IF ( ALLOCATED( w_m_l ) ) THEN
	233	WRITE ( 14 ) 'w_m_l '; WRITE ( 14 ) w_m_l
	234	ENDIF
	235	IF ( ALLOCATED( w_m_n ) ) THEN
	236	WRITE ( 14 ) 'w_m_n '; WRITE ( 14 ) w_m_n
	237	ENDIF
	238	IF ( ALLOCATED( w_m_r ) ) THEN
	239	WRITE ( 14 ) 'w_m_r '; WRITE ( 14 ) w_m_r
	240	ENDIF
	241	IF ( ALLOCATED( w_m_s ) ) THEN
	242	WRITE ( 14 ) 'w_m_s '; WRITE ( 14 ) w_m_s
	243	ENDIF
[1]	244	WRITE ( 14 ) 'z0 '; WRITE ( 14 ) z0
[72]	245	IF ( ALLOCATED( z0_av ) ) THEN
	246	WRITE ( 14 ) 'z0_av '; WRITE ( 14 ) z0_av
	247	ENDIF
[1]	248
	249	WRITE ( 14 ) 'cross_linecolors '; WRITE ( 14 ) cross_linecolors
	250	WRITE ( 14 ) 'cross_linestyles '; WRITE ( 14 ) cross_linestyles
	251	WRITE ( 14 ) 'cross_normalized_x '; WRITE ( 14 ) cross_normalized_x
	252	WRITE ( 14 ) 'cross_normalized_y '; WRITE ( 14 ) cross_normalized_y
	253	WRITE ( 14 ) 'cross_normx_factor '; WRITE ( 14 ) cross_normx_factor
	254	WRITE ( 14 ) 'cross_normy_factor '; WRITE ( 14 ) cross_normy_factor
	255	WRITE ( 14 ) 'cross_profiles '; WRITE ( 14 ) cross_profiles
	256	WRITE ( 14 ) 'cross_profile_n_coun'
	257	WRITE ( 14 ) cross_profile_number_count
	258	WRITE ( 14 ) 'cross_profile_number'; WRITE ( 14 ) cross_profile_numbers
	259	WRITE ( 14 ) 'cross_uxmax '; WRITE ( 14 ) cross_uxmax
	260	WRITE ( 14 ) 'cross_uxmax_computed'; WRITE ( 14 ) cross_uxmax_computed
	261	WRITE ( 14 ) 'cross_uxmax_normaliz'; WRITE ( 14 ) cross_uxmax_normalized
	262	WRITE ( 14 ) 'cross_uxmax_norm_com'
	263	WRITE ( 14 ) cross_uxmax_normalized_computed
	264	WRITE ( 14 ) 'cross_uxmin '; WRITE ( 14 ) cross_uxmin
	265	WRITE ( 14 ) 'cross_uxmin_computed'; WRITE ( 14 ) cross_uxmin_computed
	266	WRITE ( 14 ) 'cross_uxmin_normaliz'; WRITE ( 14 ) cross_uxmin_normalized
	267	WRITE ( 14 ) 'cross_uxmin_norm_com'
	268	WRITE ( 14 ) cross_uxmin_normalized_computed
	269	WRITE ( 14 ) 'cross_uymax '; WRITE ( 14 ) cross_uymax
	270	WRITE ( 14 ) 'cross_uymin '; WRITE ( 14 ) cross_uymin
	271	WRITE ( 14 ) 'cross_xtext '; WRITE ( 14 ) cross_xtext
	272	WRITE ( 14 ) 'dopr_crossindex '; WRITE ( 14 ) dopr_crossindex
	273	WRITE ( 14 ) 'dopr_time_count '; WRITE ( 14 ) dopr_time_count
	274	WRITE ( 14 ) 'hom_sum '; WRITE ( 14 ) hom_sum
	275	WRITE ( 14 ) 'profile_columns '; WRITE ( 14 ) profile_columns
	276	WRITE ( 14 ) 'profile_number '; WRITE ( 14 ) profile_number
	277	WRITE ( 14 ) 'profile_rows '; WRITE ( 14 ) profile_rows
	278	IF ( ALLOCATED( spectrum_x ) ) THEN
	279	WRITE ( 14 ) 'spectrum_x '; WRITE ( 14 ) spectrum_x
	280	WRITE ( 14 ) 'spectrum_y '; WRITE ( 14 ) spectrum_y
	281	ENDIF
	282
	283	!
	284	!-- Write end label. Unit 14 is closed in the main program.
	285	WRITE ( 14 ) '* end * '
	286
	287
	288	CALL cpu_log( log_point(22), 'write_3d_binary', 'stop' )
	289
	290
	291	END SUBROUTINE write_3d_binary

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

Context Navigation

source: palm/trunk/SOURCE/write_3d_binary.f90 @ 91

Download in other formats: