source: palm/trunk/UTIL/inifor/src/types.f90 @ 3182

Last change on this file since 3182 was 3182, checked in by suehring, 3 years ago

New Inifor features: grid stretching, improved command-interface, support start dates in different formats in both YYYYMMDD and YYYYMMDDHH, Ability to manually control input file prefixes (--radiation-prefix, --soil-preifx, --flow-prefix, --soilmoisture-prefix) for compatiblity with DWD forcast naming scheme; GNU-style short and long option; Prepared output of large-scale forcing profiles (no computation yet); Added preprocessor flag netcdf4 to switch output format between netCDF 3 and 4; Updated netCDF variable names and attributes to comply with PIDS v1.9; Inifor bugfixes: Improved compatibility with older Intel Intel compilers by avoiding implicit array allocation; Added origin_lon/_lat values and correct reference time in dynamic driver global attributes; corresponding PALM changes: adjustments to revised Inifor; variables names in dynamic driver adjusted; enable geostrophic forcing also in offline nested mode; variable names in LES-LES and COSMO offline nesting changed; lateral boundary flags for nesting, in- and outflow conditions renamed

  • Property svn:keywords set to Id
File size: 12.1 KB
Line 
1!> @file src/types.f90
2!------------------------------------------------------------------------------!
3! This file is part of the PALM model system.
4!
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
8! version.
9!
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.
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!
17! Copyright 2017-2018 Leibniz Universitaet Hannover
18! Copyright 2017-2018 Deutscher Wetterdienst Offenbach
19!------------------------------------------------------------------------------!
20!
21! Current revisions:
22! -----------------
23! Introduced new PALM grid stretching:
24! - Converted vertical grid_definition coordinte variables to pointers
25! Improved command line interface:
26! - Moved INIFOR configuration into a new derived data type
27! Removed unnecessary variables
28!
29!
30! Former revisions:
31! -----------------
32! $Id: types.f90 3182 2018-07-27 13:36:03Z suehring $
33! Initial revision
34!
35!
36!
37! Authors:
38! --------
39! @author Eckhard Kadasch
40!
41! Description:
42! ------------
43!> The types module provides derived data types used in INIFOR.
44!------------------------------------------------------------------------------!
45 MODULE types
46 
47 USE defs,                                                                     &
48    ONLY:  dp, DATE, PATH, SNAME, LNAME
49 USE netcdf,                                                                   &
50    ONLY:  NF90_MAX_VAR_DIMS, NF90_MAX_NAME
51
52 IMPLICIT NONE
53
54 TYPE inifor_config
55    CHARACTER(LEN=DATE)  ::  start_date           !< String of the FORMAT YYYYMMDDHH indicating the start of the intended PALM-4U simulation
56
57    CHARACTER(LEN=PATH)  ::  input_path           !< Path to the input data file directory
58    CHARACTER(LEN=PATH)  ::  hhl_file             !< Path to the file containing the COSMO-DE HHL variable (height of half layers, i.e. vertical cell faces)
59    CHARACTER(LEN=PATH)  ::  namelist_file        !< Path to the PALM-4U namelist file
60    CHARACTER(LEN=PATH)  ::  output_file          !< Path to the INIFOR output file (i.e. PALM-4U dynamic driver')
61    CHARACTER(LEN=PATH)  ::  soiltyp_file         !< Path to the file containing the COSMO-DE SOILTYP variable (map of COSMO-DE soil types)
62    CHARACTER(LEN=PATH)  ::  static_driver_file   !< Path to the file containing the COSMO-DE SOILTYP variable (map of COSMO-DE soil types)
63
64    CHARACTER(LEN=SNAME) ::  flow_prefix          !< Prefix of flow input files, e.g. 'laf' for COSMO-DE analyses
65    CHARACTER(LEN=SNAME) ::  soil_prefix          !< Prefix of soil input files, e.g. 'laf' for COSMO-DE analyses
66    CHARACTER(LEN=SNAME) ::  radiation_prefix     !< Prefix of radiation input files, e.g 'laf' for COSMO-DE analyses
67    CHARACTER(LEN=SNAME) ::  soilmoisture_prefix  !< Prefix of input files for soil moisture spin-up, e.g 'laf' for COSMO-DE analyses
68
69    CHARACTER(LEN=SNAME) ::  bc_mode
70    CHARACTER(LEN=SNAME) ::  ic_mode
71    CHARACTER(LEN=SNAME) ::  rotation_method
72
73    REAL(dp)             ::  p0
74    REAL(dp)             ::  ug
75    REAL(dp)             ::  vg
76    REAL(dp)             ::  z0 !< Elevation of the PALM-4U domain above sea level [m]
77 END TYPE inifor_config
78
79 TYPE grid_definition
80    CHARACTER(LEN=SNAME)  ::  name(3)       !< names of the grid dimensions, e.g. (/'x', 'y', 'z'/) or (/'latitude', 'longitude', 'height'/)
81    INTEGER               ::  nx            !< number of gridpoints in the first dimension
82    INTEGER               ::  ny            !< number of gridpoints in the second dimension
83    INTEGER               ::  nz            !< number of gridpoints in the third dimension
84    INTEGER, ALLOCATABLE  ::  ii(:,:,:)     !< Given a point (i,j,k) in the PALM-4U grid, ii(i,j,l) gives the x index of the l'th horizontl neighbour on the COSMO-DE grid.
85    INTEGER, ALLOCATABLE  ::  jj(:,:,:)     !< Given a point (i,j,k) in the PALM-4U grid, jj(i,j,l) gives the y index of the l'th horizontl neighbour on the COSMO-DE grid.
86    INTEGER, ALLOCATABLE  ::  kk(:,:,:,:)   !< Given a point (i,j,k) in the PALM-4U grid, kk(i,j,k,l) gives the z index of the l'th vertical neighbour in the intermediate grid.
87    REAL(dp)              ::  lx            !< domain length in the first dimension [m]
88    REAL(dp)              ::  ly            !< domain length in the second dimension [m]
89    REAL(dp)              ::  x0            !< x coordinate of PALM-4U domain projection centre, i.e. location of zero distortion
90    REAL(dp)              ::  y0            !< y coordinate of PALM-4U domain projection centre, i.e. location of zwro distortion
91    REAL(dp)              ::  z0            !< displacement of the coordinate origin above sea level [m]
92    REAL(dp), ALLOCATABLE ::  x(:)          !< coordinates of cell centers in x direction [m]
93    REAL(dp), ALLOCATABLE ::  y(:)          !< coordinates of cell centers in y direction [m]
94    REAL(dp), POINTER     ::  z(:)          !< coordinates of cell centers in z direction [m]
95    REAL(dp), ALLOCATABLE ::  h(:,:,:)      !< heights grid point for intermediate grids [m]
96    REAL(dp), POINTER     ::  hhl(:,:,:)    !< heights of half layers (cell faces) above sea level in COSMO-DE, read in from
97    REAL(dp), POINTER     ::  hfl(:,:,:)    !< heights of full layers (cell centres) above sea level in COSMO-DE, computed as arithmetic average of hhl
98    REAL(dp), POINTER     ::  depths(:)     !< depths of output soil layers, equal the depths of the source model (e.g. COSMO-DE)
99    REAL(dp), ALLOCATABLE ::  xu(:)         !< coordinates of cell faces in x direction [m]
100    REAL(dp), ALLOCATABLE ::  yv(:)         !< coordinates of cell faces in y direction [m]
101    REAL(dp), POINTER     ::  zw(:)         !< coordinates of cell faces in z direction [m]
102    REAL(dp), ALLOCATABLE ::  lat(:)        !< rotated-pole latitudes of scalars (cell centers) of the COSMO-DE grid [rad]
103    REAL(dp), ALLOCATABLE ::  lon(:)        !< rotated-pole longitudes of scalars (cell centres) of the COSMO-DE grid [rad]
104    REAL(dp), ALLOCATABLE ::  latv(:)       !< rotated-pole latitudes of v winds (face centres in latitudal/y direction) [rad]
105    REAL(dp), ALLOCATABLE ::  lonu(:)       !< rotated-pole latitudes of u winds (face centres in longitudal/x direction) [rad]
106    REAL(dp), ALLOCATABLE ::  clat(:,:)     !< latitudes of PALM-4U cell centres in COSMO-DE's rotated-pole grid [rad]
107    REAL(dp), ALLOCATABLE ::  clon(:,:)     !< longitudes of PALM-4U scalars (cell centres) in COSMO-DE's rotated-pole grid [rad]
108    REAL(dp), ALLOCATABLE ::  clatu(:,:)    !< latitudes of PALM-4U u winds (cell faces in u direction) in COSMO-DE's rotated-pole grid [rad]
109    REAL(dp), ALLOCATABLE ::  clonu(:,:)    !< longitudes of PALM-4U u winds (cell faces in u direction) in COSMO-DE's rotated-pole grid [rad]
110    REAL(dp), ALLOCATABLE ::  clatv(:,:)    !< latitudes of PALM-4U v winds (cell faces in v direction) in COSMO-DE's rotated-pole grid [rad]
111    REAL(dp), ALLOCATABLE ::  clonv(:,:)    !< longitudes of PALM-4U v winds (cell faces in v direction) in COSMO-DE's rotated-pole grid [rad]
112    REAL(dp), ALLOCATABLE ::  w_horiz(:,:,:)   !< weights for bilinear horizontal interpolation
113    REAL(dp), ALLOCATABLE ::  w_verti(:,:,:,:) !< weights for linear vertical interpolation
114 END TYPE grid_definition
115
116
117 TYPE nc_file
118    CHARACTER(LEN=PATH)   ::  name          !< file name
119    INTEGER               ::  dimid_time    !< NetCDF IDs of the time dimension
120    INTEGER               ::  dimids_scl(3) !< NetCDF IDs of the grid dimensions for scalar points x, y, z
121    INTEGER               ::  dimids_vel(3) !< NetCDF IDs of the grid dimensions for velocity points xu, yu, zu
122    INTEGER               ::  dimids_soil(3)!< NetCDF IDs of the grid dimensions for soil points x, y, depth
123    INTEGER               ::  dimvarid_time !< NetCDF IDs of the time variable
124    INTEGER               ::  dimvarids_scl(3) !< NetCDF IDs of the grid coordinates of scalars x, y, z
125    INTEGER               ::  dimvarids_vel(3) !< NetCDF IDs of the grid coordinates of velocities xu, yu, zu. Note that velocities are located at mix of both coordinates, e.g. u(xu, y, z).
126    INTEGER               ::  dimvarids_soil(3)!< NetCDF IDs of the grid coordinates for soil points x, y, depth
127    REAL(dp), POINTER     ::  time(:)       ! vector of output time steps
128 END TYPE nc_file
129
130
131 TYPE nc_var
132    INTEGER                               ::  varid     !< NetCDF ID of the variable
133    INTEGER                               ::  input_id  !< ID of the correpsonding input variables, only valid for output variables
134    INTEGER                               ::  ndim      !< number of NetCDF dimensions
135    INTEGER                               ::  nt        !< number of output time steps
136    INTEGER                               ::  lod       !< NetCDF attribute indicating the PALM-4U level of detail
137    INTEGER, DIMENSION(NF90_MAX_VAR_DIMS) ::  dimids    !< NetCDF IDs of the dimensions
138    INTEGER, DIMENSION(NF90_MAX_VAR_DIMS) ::  dimvarids !< IDs of NetCDF dimension variables
139    INTEGER, DIMENSION(NF90_MAX_VAR_DIMS) ::  dimlen    !< length of NetCDF dimensions
140    CHARACTER(LEN=NF90_MAX_NAME), DIMENSION(NF90_MAX_VAR_DIMS) ::  dimname !< names of NetCDF dimensions
141    CHARACTER(LEN=SNAME)                  ::  name                      !< NetCDF short name of the variable
142    CHARACTER(LEN=LNAME)                  ::  standard_name             !< NetCDF standard name of the variable
143    CHARACTER(LEN=LNAME)                  ::  long_name                 !< NetCDF long name of the variable
144    CHARACTER(LEN=LNAME)                  ::  source                    !< NetCDF attribute indicating the data source for the output
145    CHARACTER(LEN=SNAME)                  ::  units                     !< NetCDF units of the variable
146    CHARACTER(LEN=SNAME)                  ::  kind                      !< Kind of grid
147    CHARACTER(LEN=SNAME)                  ::  task                      !< Processing task that generates this variable, e.g. 'interpolate_2d' or 'average profile'
148    LOGICAL                               ::  to_be_processed = .FALSE. !< INIFOR flag indicating whether variable shall be processed
149    LOGICAL                               ::  is_read = .FALSE.         !< INIFOR flag indicating whether variable has been read
150    LOGICAL                               ::  is_upside_down  = .FALSE. !< INIFOR flag indicating whether variable shall be processed
151    TYPE(grid_definition), POINTER        ::  grid                      !< Pointer to the corresponding output grid
152    TYPE(grid_definition), POINTER        ::  intermediate_grid         !< Pointer to the corresponding intermediate grid
153 END TYPE nc_var
154
155
156 TYPE io_group                                          !< Input/Output group, groups together output variabels that share their input variables. For instance, all boundary surfaces and initialization fields of the potential temperature are base on T and p.
157    INTEGER                          ::  nt             !< maximum number of output time steps across all output variables
158    INTEGER                          ::  nv             !< number of output variables
159    CHARACTER(LEN=SNAME)             ::  kind           !< kind of I/O group
160    CHARACTER(LEN=PATH), ALLOCATABLE ::  in_files(:)    !< list of nt input files
161    TYPE(nc_var), ALLOCATABLE        ::  out_vars(:)    !< list of output variables
162    TYPE(nc_var), ALLOCATABLE        ::  in_var_list(:) !< list of input variables
163    LOGICAL                          ::  to_be_processed = .FALSE. !< Inifor flag indicating whether I/O group shall be processed
164    LOGICAL                          ::  is_accumulated = .FALSE.  !< Flag indicating whether this I/O group contains accumulated variables
165    LOGICAL                          ::  is_preprocessed = .FALSE. !< Inifor flag indicating whether the I/O group has been preprocessed
166 END TYPE io_group
167
168
169 TYPE container
170   REAL(dp), ALLOCATABLE ::  array(:,:,:)
171   LOGICAL               ::  is_preprocessed = .FALSE.
172 END TYPE container
173
174 END MODULE types
175
Note: See TracBrowser for help on using the repository browser.