Home

Context Navigation

source: palm/trunk/UTIL/inifor/src/inifor.f90 @ 3201

Last change on this file since 3201 was 3183, checked in by suehring, 6 years ago
last commit documented
Property svn:keywords set to `Id`
File size: 12.4 KB

Rev	Line
[2696]	1	!> @file src/inifor.f90
	2	!------------------------------------------------------------------------------!
[2718]	3	! This file is part of the PALM model system.
[2696]	4	!
[2718]	5	! PALM is free software: you can redistribute it and/or modify it under the
	6	! terms of the GNU General Public License as published by the Free Software
	7	! Foundation, either version 3 of the License, or (at your option) any later
[2696]	8	! version.
	9	!
[2718]	10	! PALM is distributed in the hope that it will be useful, but WITHOUT ANY
	11	! WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
	12	! A PARTICULAR PURPOSE. See the GNU General Public License for more details.
[2696]	13	!
	14	! You should have received a copy of the GNU General Public License along with
	15	! PALM. If not, see <http://www.gnu.org/licenses/>.
	16	!
[2718]	17	! Copyright 2017-2018 Leibniz Universitaet Hannover
	18	! Copyright 2017-2018 Deutscher Wetterdienst Offenbach
[2696]	19	!------------------------------------------------------------------------------!
	20	!
	21	! Current revisions:
	22	! -----------------
[3183]	23	!
	24	!
	25	! Former revisions:
	26	! -----------------
	27	! $Id: inifor.f90 3183 2018-07-27 14:25:55Z sward $
[3182]	28	! Introduced new PALM grid stretching
	29	! Renamend initial-condition mode variable 'mode' to 'ic_mode'
	30	! Improved log messages
[2696]	31	!
	32	!
[3183]	33	! 3182 2018-07-27 13:36:03Z suehring
[2696]	34	! Initial revision
	35	!
	36	!
	37	!
	38	! Authors:
	39	! --------
	40	! @author Eckhard Kadasch
	41	!
	42	! Description:
	43	! ------------
	44	!> INIFOR is an interpolation tool for generating meteorological initialization
	45	!> and forcing data for the urban climate model PALM-4U. The required
	46	!> meteorological fields are interpolated from output data of the mesoscale
	47	!> model COSMO-DE. This is the main program file.
	48	!------------------------------------------------------------------------------!
	49	PROGRAM inifor
	50
	51	USE control
	52	USE defs
	53	USE grid, &
	54	ONLY: setup_parameters, setup_grids, setup_variable_tables, &
	55	setup_io_groups, fini_grids, fini_variables, fini_io_groups, &
[3182]	56	fini_file_lists, preprocess, origin_lon, origin_lat, &
[2696]	57	output_file, io_group_list, output_var_table, &
[3182]	58	cosmo_grid, palm_grid, nx, ny, nz, ug, vg, p0, cfg, &
	59	average_imin, average_imax, average_jmin, average_jmax
[2696]	60
	61	USE io
	62	USE transform, &
	63	ONLY: average_profile, average_2d, interpolate_2d, interpolate_3d
	64	USE types
	65
	66	IMPLICIT NONE
	67
	68	INTEGER :: igroup
	69	INTEGER :: ivar
	70	INTEGER :: iter
	71	REAL(dp), ALLOCATABLE, DIMENSION(:,:,:) :: output_arr
	72	TYPE(nc_var), POINTER :: output_var
	73	TYPE(io_group), POINTER :: group
	74	TYPE(container), ALLOCATABLE :: input_buffer(:)
	75
	76	!> \mainpage About INIFOR
	77	!> ...
	78	!
	79	!------------------------------------------------------------------------------
	80	!- Section 1: Initialization
	81	!------------------------------------------------------------------------------
	82	CALL run_control('init', 'void')
	83
	84	! Initialize INIFOR's parameters from command-line interface and namelists
	85	CALL setup_parameters()
	86
	87	! Initialize all grids, including interpolation neighbours and weights
	88	CALL setup_grids()
	89	CALL run_control('time', 'init')
	90
	91	! Initialize the netCDF output file and define dimensions
[3182]	92	CALL setup_netcdf_dimensions(output_file, palm_grid, cfg % start_date, &
	93	origin_lon, origin_lat)
[2696]	94	CALL run_control('time', 'write')
	95
	96	! Set up the tables containing the input and output variables and set
	97	! the corresponding netCDF dimensions for each output variable
[3182]	98	CALL setup_variable_tables(cfg % ic_mode)
[2696]	99	CALL run_control('time', 'write')
	100
	101	! Add the output variables to the netCDF output file
	102	CALL setup_netcdf_variables(output_file % name, output_var_table)
	103
[3182]	104	CALL setup_io_groups()
[2696]	105	CALL run_control('time', 'init')
	106
	107	!------------------------------------------------------------------------------
	108	!- Section 2: Main loop
	109	!------------------------------------------------------------------------------
	110
	111	DO igroup = 1, SIZE(io_group_list)
	112
	113	group => io_group_list(igroup)
	114	IF ( group % to_be_processed ) THEN
	115
	116	DO iter = 1, group % nt
	117
	118	!------------------------------------------------------------------------------
	119	!- Section 2.1: Read and preprocess input data
	120	!------------------------------------------------------------------------------
	121	CALL read_input_variables(group, iter, input_buffer)
	122	CALL run_control('time', 'read')
	123
	124	CALL preprocess(group, input_buffer, cosmo_grid, iter)
	125	CALL run_control('time', 'comp')
	126
[3182]	127	!TODO: move this assertion into 'preprocess'.
[2696]	128	IF ( .NOT. ALL(input_buffer(:) % is_preprocessed .AND. .TRUE.) ) THEN
	129	message = "Input buffers for group '" // TRIM(group % kind) // &
	130	"' could not be preprocessed sucessfully."
	131	CALL abort('main loop', message)
	132	END IF
	133
	134	!------------------------------------------------------------------------------
	135	!- Section 2.2: Interpolate each output variable of the group
	136	!------------------------------------------------------------------------------
	137	DO ivar = 1, group % nv
	138
	139	output_var => group % out_vars( ivar )
	140
	141	IF ( output_var % to_be_processed .AND. &
	142	iter .LE. output_var % nt ) THEN
	143
	144	message = "Processing '" // TRIM(output_var % name) // &
	145	"' (" // TRIM(output_var % kind) // &
	146	"), iteration " // TRIM(str(iter)) //" of " // &
	147	TRIM(str(output_var % nt))
	148	CALL report('main loop', message)
	149
	150	SELECT CASE( TRIM(output_var % task) )
	151
	152	CASE( 'interpolate_2d' )
	153
	154	SELECT CASE( TRIM(output_var % kind) )
	155
	156	CASE( 'init soil' )
	157
	158	ALLOCATE( output_arr( 0:output_var % grid % nx, &
	159	0:output_var % grid % ny, &
	160	SIZE(output_var % grid % depths) ) )
	161
	162	CASE ( 'surface forcing' )
	163
	164	ALLOCATE( output_arr( 0:output_var % grid % nx, &
	165	0:output_var % grid % ny, 1 ) )
	166
	167	CASE DEFAULT
	168
[3182]	169	message = "'" // TRIM(output_var % kind) // "' is not a soil variable"
	170	CALL abort("main loop", message)
[2696]	171
	172	END SELECT
	173	CALL run_control('time', 'alloc')
	174
	175	CALL interpolate_2d(input_buffer(output_var % input_id) % array(:,:,:), &
	176	output_arr(:,:,:), output_var % intermediate_grid, output_var)
	177	CALL run_control('time', 'comp')
	178
	179
	180	CASE ( 'interpolate_3d' )
	181
	182	ALLOCATE( output_arr( 0:output_var % grid % nx, &
	183	0:output_var % grid % ny, &
[3182]	184	1:output_var % grid % nz ) )
[2696]	185
	186	CALL run_control('time', 'alloc')
	187	CALL interpolate_3d( &
	188	input_buffer(output_var % input_id) % array(:,:,:), &
	189	output_arr(:,:,:), &
	190	output_var % intermediate_grid, &
	191	output_var % grid)
	192	CALL run_control('time', 'comp')
	193
	194	CASE ( 'average profile' )
	195
	196	ALLOCATE( output_arr( 0:output_var % grid % nx, &
	197	0:output_var % grid % ny, &
[3182]	198	1:output_var % grid % nz ) )
[2696]	199	CALL run_control('time', 'alloc')
	200
	201
	202	CALL average_profile( &
	203	input_buffer(output_var % input_id) % array(:,:,:), &
[3182]	204	output_arr(:,:,:), average_imin, average_imax, &
	205	average_jmin, average_jmax, &
[2696]	206	output_var % intermediate_grid, &
	207	output_var % grid)
	208	CALL run_control('time', 'comp')
	209
	210	CASE ( 'average scalar' )
	211
	212	ALLOCATE( output_arr(1,1,1) )
	213	CALL run_control('time', 'alloc')
	214	output_arr(1,1,1) = p0
	215	CALL run_control('time', 'comp')
	216
[3182]	217	CASE ( 'set profile' )
[2696]	218
[3182]	219	ALLOCATE( output_arr( 1, 1, 1:nz ) )
[2696]	220	CALL run_control('time', 'alloc')
	221
	222	SELECT CASE (TRIM(output_var % name))
	223
	224	CASE('ls_forcing_ug')
	225	output_arr(1, 1, :) = ug
	226
	227	CASE('ls_forcing_vg')
	228	output_arr(1, 1, :) = vg
	229
[3182]	230	CASE('nudging_tau')
	231	output_arr(1, 1, :) = NUDGING_TAU
	232
[2696]	233	CASE DEFAULT
	234	message = "'" // TRIM(output_var % name) // &
	235	"' is not a valid '" // TRIM(output_var % kind) //&
	236	"' variable kind."
	237	CALL abort('main loop', message)
	238	END SELECT
	239	CALL run_control('time', 'comp')
	240
[3182]	241	CASE('average large-scale profile')
	242	message = "Averaging of large-scale forcing profiles " //&
	243	"has not been implemented, yet."
	244	CALL abort('main loop', message)
	245	!ALLOCATE( output_arr( 1, 1, 1:nz ) )
	246
[2696]	247	CASE DEFAULT
	248	message = "Processing task '" // TRIM(output_var % task) //&
	249	"' not recognized."
	250	CALL abort('', message)
	251
	252	END SELECT
	253	CALL run_control('time', 'comp')
	254
	255	!------------------------------------------------------------------------------
	256	!- Section 2.3: Write current time step of current variable
	257	!------------------------------------------------------------------------------
	258	message = "Writing variable '" // TRIM(output_var%name) // "'."
	259	CALL report('main loop', message)
	260	CALL update_output(output_var, output_arr, iter, output_file)
	261	CALL run_control('time', 'write')
	262
	263	DEALLOCATE(output_arr)
	264	CALL run_control('time', 'alloc')
	265
	266	END IF
	267
	268	END DO ! ouput variables
	269
	270	IF ( group % kind == 'running average' .OR. &
	271	group % kind == 'accumulated' ) THEN
	272	! Keep input buffer around for averaged (radiation) and
	273	! accumulated COSMO-DE quantities (precipitation).
	274	ELSE
	275	CALL report('main loop', 'Deallocating input buffer')
	276	DEALLOCATE(input_buffer)
	277	END IF
	278	CALL run_control('time', 'alloc')
	279
	280	END DO ! time steps / input files
	281
	282	IF (ALLOCATED(input_buffer)) THEN
	283	CALL report('main loop', 'Deallocating input buffer')
	284	DEALLOCATE(input_buffer)
	285	END IF
	286	CALL run_control('time', 'alloc')
	287
	288	ELSE
	289
[3182]	290	message = "Skipping IO group " // TRIM(str(igroup)) // " '" // TRIM(group % kind) // "'"
[2696]	291	IF ( ALLOCATED(group % in_var_list) ) THEN
	292	message = TRIM(message) // " with input variable '" // &
	293	TRIM(group % in_var_list(1) % name) // "'."
	294	END IF
	295
	296	CALL report('main loop', message)
	297
	298	END IF ! IO group % to_be_processed
	299
	300	END DO ! IO groups
	301
	302	!------------------------------------------------------------------------------
	303	!- Section 3: Clean up.
	304	!------------------------------------------------------------------------------
	305	CALL fini_file_lists()
	306	CALL fini_io_groups()
	307	CALL fini_variables()
	308	!CALL fini_grids()
	309	CALL run_control('time', 'alloc')
	310	CALL run_control('report', 'void')
	311
[3182]	312	message = "Finished writing dynamic driver '" // TRIM(output_file % name) // &
[2696]	313	"' successfully."
	314	CALL report('main loop', message)
	315
	316
	317	END PROGRAM inifor

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

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |