Changeset 3395 for palm/trunk/UTIL/inifor/src/types.f90

Timestamp:

Oct 22, 2018 5:32:49 PM (6 years ago)

Author:

eckhard

Message:

inifor: Added computation of geostrophic winds from COSMO input

File:

• palm/trunk/UTIL/inifor/src/types.f90 (modified) (7 diffs)

palm/trunk/UTIL/inifor/src/types.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Added *_is_set LOGICALs to inifor_config type to indicate option invocation
+!     from the command-line 
+! Added 1D index vertical weights lists to support addressing averaging regions
+!     by list of columns instead of index bounds
+! 
+! 
+! 3183 2018-07-27 14:25:55Z suehring
 ! Introduced new PALM grid stretching:
 ! - Converted vertical grid_definition coordinte variables to pointers
@@ -66 +73 @@
 
     CHARACTER(LEN=SNAME) ::  flow_prefix          !< Prefix of flow input files, e.g. 'laf' for COSMO-DE analyses
+    CHARACTER(LEN=SNAME) ::  input_prefix         !< Prefix of all input files, e.g. 'laf' for COSMO-DE analyses
+    CHARACTER(LEN=SNAME) ::  radiation_prefix     !< Prefix of radiation input files, e.g 'laf' for COSMO-DE analyses
     CHARACTER(LEN=SNAME) ::  soil_prefix          !< Prefix of soil input files, e.g. 'laf' for COSMO-DE analyses
-    CHARACTER(LEN=SNAME) ::  radiation_prefix     !< Prefix of radiation input files, e.g 'laf' for COSMO-DE analyses
     CHARACTER(LEN=SNAME) ::  soilmoisture_prefix  !< Prefix of input files for soil moisture spin-up, e.g 'laf' for COSMO-DE analyses
 
+    CHARACTER(LEN=SNAME) ::  averaging_mode
     CHARACTER(LEN=SNAME) ::  bc_mode
     CHARACTER(LEN=SNAME) ::  ic_mode
@@ -77 +86 @@
     REAL(dp)             ::  ug
     REAL(dp)             ::  vg
-    REAL(dp)             ::  z0 !< Elevation of the PALM-4U domain above sea level [m]
+    REAL(dp)             ::  z0                   !< Elevation of the PALM-4U domain above sea level [m]
+    REAL(dp)             ::  averaging_angle      !< latitudal and longitudal width of averaging regions [deg]
+    
+
+    LOGICAL              ::  debug
+    LOGICAL              ::  p0_is_set
+    LOGICAL              ::  ug_is_set
+    LOGICAL              ::  vg_is_set
+    LOGICAL              ::  flow_prefix_is_set          !< 
+    LOGICAL              ::  input_prefix_is_set         !< 
+    LOGICAL              ::  radiation_prefix_is_set     !< 
+    LOGICAL              ::  soil_prefix_is_set          !< 
+    LOGICAL              ::  soilmoisture_prefix_is_set  !< 
  END TYPE inifor_config
 
  TYPE grid_definition
     CHARACTER(LEN=SNAME)  ::  name(3)       !< names of the grid dimensions, e.g. (/'x', 'y', 'z'/) or (/'latitude', 'longitude', 'height'/)
+    CHARACTER(LEN=SNAME)  ::  kind          !< names of the grid dimensions, e.g. (/'x', 'y', 'z'/) or (/'latitude', 'longitude', 'height'/)
+    INTEGER               ::  k_min         !< Index of lowest PALM grid level that is not cut by local COSMO orography; vertically separates interpolation and extrapolation region.
     INTEGER               ::  nx            !< number of gridpoints in the first dimension
     INTEGER               ::  ny            !< number of gridpoints in the second dimension
-    INTEGER               ::  nz            !< number of gridpoints in the third dimension
+    INTEGER               ::  nz            !< number of gridpoints in the third dimension, used for PALM points
+    INTEGER               ::  nlev          !< number of COSMO grid levels
+    INTEGER               ::  n_columns     !< number of averaging columns of the source grid 
     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.
     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.
     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.
+    INTEGER, ALLOCATABLE  ::  iii(:)        !< profile averaging neighbour indices 
+    INTEGER, ALLOCATABLE  ::  jjj(:)        !< profile averaging neighbour indices 
+    INTEGER, ALLOCATABLE  ::  kkk(:,:,:)    !< indices of vertical interpolation neightbours, kkk(<source column>, <PALM k level>, <neighbour index>)
     REAL(dp)              ::  lx            !< domain length in the first dimension [m]
     REAL(dp)              ::  ly            !< domain length in the second dimension [m]
@@ -97 +125 @@
     REAL(dp), POINTER     ::  z(:)          !< coordinates of cell centers in z direction [m]
     REAL(dp), ALLOCATABLE ::  h(:,:,:)      !< heights grid point for intermediate grids [m]
+    REAL(dp), POINTER     ::  cosmo_h(:,:,:)!< pointer to appropriate COSMO level heights (scalar/w) [m]
     REAL(dp), POINTER     ::  hhl(:,:,:)    !< heights of half layers (cell faces) above sea level in COSMO-DE, read in from 
     REAL(dp), POINTER     ::  hfl(:,:,:)    !< heights of full layers (cell centres) above sea level in COSMO-DE, computed as arithmetic average of hhl
@@ -115 +144 @@
     REAL(dp), ALLOCATABLE ::  w_horiz(:,:,:)   !< weights for bilinear horizontal interpolation
     REAL(dp), ALLOCATABLE ::  w_verti(:,:,:,:) !< weights for linear vertical interpolation
+    REAL(dp), ALLOCATABLE ::  w(:,:,:)      !< vertical interpolation weights, w(<source_column>, <PALM k level>, <neighbour index>) [-]
  END TYPE grid_definition
 
@@ -150 +180 @@
     CHARACTER(LEN=SNAME)                  ::  task                      !< Processing task that generates this variable, e.g. 'interpolate_2d' or 'average profile'
     LOGICAL                               ::  to_be_processed = .FALSE. !< INIFOR flag indicating whether variable shall be processed
+    LOGICAL                               ::  is_internal = .FALSE.     !< INIFOR flag indicating whether variable shall be written to netCDF file (.FALSE.) or kept for later (.TRUE.)
     LOGICAL                               ::  is_read = .FALSE.         !< INIFOR flag indicating whether variable has been read
-    LOGICAL                               ::  is_upside_down  = .FALSE. !< INIFOR flag indicating whether variable shall be processed
+    LOGICAL                               ::  is_upside_down  = .FALSE. !< INIFOR flag indicating whether vertical dimension is reversed (typically the case with COSMO-DE atmospheric fields)
     TYPE(grid_definition), POINTER        ::  grid                      !< Pointer to the corresponding output grid
     TYPE(grid_definition), POINTER        ::  intermediate_grid         !< Pointer to the corresponding intermediate grid
+    TYPE(grid_definition), POINTER        ::  averaging_grid         !< Pointer to the corresponding intermediate grid
  END TYPE nc_var
 
@@ -159 +191 @@
  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.
     INTEGER                          ::  nt             !< maximum number of output time steps across all output variables
-    INTEGER                          ::  nv             !< number of output variables
+    INTEGER                          ::  nv             !< number of netCDF output variables
+    INTEGER                          ::  n_inputs       !< number of input variables
+    INTEGER                          ::  n_output_quantities !< number of physical quantities required for computing netCDF output variables
     CHARACTER(LEN=SNAME)             ::  kind           !< kind of I/O group
     CHARACTER(LEN=PATH), ALLOCATABLE ::  in_files(:)    !< list of nt input files