Changeset 4790 for palm/trunk/UTIL/inifor/src/inifor_io.f90

Timestamp:

Nov 18, 2020 1:05:58 PM (3 years ago)

Author:

eckhard

Message:

inifor: Validate netCDF dimensions of COSMO input files

File:

• palm/trunk/UTIL/inifor/src/inifor_io.f90 (modified) (4 diffs)

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

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Validate mesoscale input data by the number of dimensions and their lengths
+!
+!
+! 4675 2020-09-11 10:00:26Z eckhard
 ! New command-line options for soil profile initialization and homogeneous
 !    boundary conditions
@@ -409 +413 @@
  SUBROUTINE netcdf_define_variable(var, ncid)
 
-     TYPE(nc_var), INTENT(INOUT) ::  var
-     INTEGER, INTENT(IN)         ::  ncid
-
-     CALL check(nf90_def_var(ncid, var%name, NF90_FLOAT,       var%dimids(1:var%ndim), var%varid))
-     CALL check(nf90_put_att(ncid, var%varid, "long_name",     var%long_name))
-     CALL check(nf90_put_att(ncid, var%varid, "units",         var%units))
-     IF ( var%lod .GE. 0 )  THEN
-        CALL check(nf90_put_att(ncid, var%varid, "lod",           var%lod))
-     ENDIF 
-     CALL check(nf90_put_att(ncid, var%varid, "source",        var%source))
-     CALL check(nf90_put_att(ncid, var%varid, "_FillValue",    NF90_FILL_REAL))
+    TYPE(nc_var), INTENT(INOUT) ::  var
+    INTEGER, INTENT(IN)         ::  ncid
+
+    CALL check(nf90_def_var(ncid, var%name, NF90_FLOAT,       var%dimids(1:var%ndim), var%varid))
+    CALL check(nf90_put_att(ncid, var%varid, "long_name",     var%long_name))
+    CALL check(nf90_put_att(ncid, var%varid, "units",         var%units))
+    IF ( var%lod .GE. 0 )  THEN
+       CALL check(nf90_put_att(ncid, var%varid, "lod",           var%lod))
+    ENDIF 
+    CALL check(nf90_put_att(ncid, var%varid, "source",        var%source))
+    CALL check(nf90_put_att(ncid, var%varid, "_FillValue",    NF90_FILL_REAL))
 
  END SUBROUTINE netcdf_define_variable
@@ -433 +437 @@
  SUBROUTINE netcdf_get_dimensions(var, ncid)
 
-     TYPE(nc_var), INTENT(INOUT) ::  var
-     INTEGER, INTENT(IN)         ::  ncid
-     INTEGER                     ::  i
-     CHARACTER(SNAME)            ::  null
-
-     DO  i = 1, var%ndim
-        CALL check(nf90_inquire_dimension(ncid, var%dimids(i), &
-                                          name = null, &
-                                          len  = var%dimlen(i)  ) )
-     ENDDO
+    TYPE(nc_var), INTENT(INOUT) ::  var
+    INTEGER, INTENT(IN)         ::  ncid
+    INTEGER                     ::  i
+    CHARACTER(SNAME)            ::  null
+
+    DO  i = 1, var%ndim
+       CALL check(nf90_inquire_dimension(ncid, var%dimids(i), &
+                                         name = null, &
+                                         len  = var%dimlen(i)  ) )
+    ENDDO
 
  END SUBROUTINE netcdf_get_dimensions
@@ -819 +823 @@
 
 
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> Checks wheather the COSMO grid matches the shape of the meteorological input
+!> data by comparing the number of netCDF dimensions and their lengths in the
+!> hhl.nc and the first of the *-flow files.
+!------------------------------------------------------------------------------!
+ SUBROUTINE validate_dataset(flow_files, hhl_file)
+    CHARACTER(LEN=PATH), INTENT(IN) ::  flow_files(:) !< paths to files containing atmospheric variables
+    CHARACTER(LEN=PATH), INTENT(IN) ::  hhl_file      !< path to file containing the HHL variable (height of half layers)
+
+    CHARACTER(SNAME), PARAMETER ::  NC_W_NAME = 'W'
+    TYPE(nc_var)                ::  hhl_var, flow_var
+    INTEGER                     ::  dim_idx, ncid, ndims_flow, ndims_hhl, varid
+    REAL(wp), ALLOCATABLE       ::  hhl_dim_vector(:), flow_dim_vector(:)
+    LOGICAL                     ::  dims_have_same_length
+
+    IF ( nf90_open( TRIM( flow_files(1) ), NF90_NOWRITE, ncid ) .EQ. NF90_NOERR )  THEN
+
+       CALL check( nf90_inq_varid( ncid, NC_W_NAME, varid=varid ) )
+       CALL check( nf90_inquire_variable( ncid, varid, ndims=ndims_flow ) )
+       CALL netcdf_get_dimensions( flow_var, ncid )
+       CALL check( nf90_close( ncid ) )
+
+    ELSE
+
+       message = "Failed to read netCDF dimensions'" //                        &
+                 "' from file '" // TRIM( flow_files(1) ) // "'."
+       CALL inifor_abort( 'validate_dataset', message )
+
+    ENDIF
+
+    IF ( nf90_open( TRIM( hhl_file ), NF90_NOWRITE, ncid ) .EQ. NF90_NOERR )  THEN
+
+       CALL check( nf90_inq_varid( ncid, NC_HHL_NAME, varid=varid ) )
+       CALL check( nf90_inquire_variable( ncid, varid, ndims=ndims_hhl ) )
+       CALL netcdf_get_dimensions( hhl_var, ncid )
+       CALL check( nf90_close( ncid ) )
+
+    ELSE
+
+       message = "Failed to read netCDF dimensions'" //                        &
+                 "' from file '" // TRIM(hhl_file) // "'."
+       CALL inifor_abort( 'validate_dataset', message )
+
+    ENDIF
+
+!
+!-- Compare number dimensions of 'HHL' in hhl file and 'W' in first flow file
+    IF  ( .NOT. ndims_flow .EQ. ndims_hhl )  THEN
+       message = "Mesoscale data inconsistent. Number of dimensions in the " //&
+                 "hhl file does not match with the meteorologial fields " //   &
+                 "in the *-flow files (" //                                    &
+                 "HHL: ndims = " // TRIM( str( ndims_hhl ) )  // ", " //       &
+                 "W: ndims = "   // TRIM( str( ndims_flow ) ) // ")."
+       CALL inifor_abort( 'validate_dataset', message )
+    ENDIF
+
+
+!
+!-- Compare lengths of each dimension, ignoring time (dim_idx = 1)
+    DO dim_idx = 2, ndims_hhl
+
+       CALL get_dimension_vector_of_variable(                                  &
+          NC_HHL_NAME,                                                         &
+          dim_idx = dim_idx,                                                   &
+          filename = hhl_file,                                                 &
+          dim_vector = hhl_dim_vector                                          &
+       )
+
+       CALL get_dimension_vector_of_variable(                                  &
+          NC_W_NAME,                                                           &
+          dim_idx = dim_idx,                                                   &
+          filename = flow_files(1),                                            &
+          dim_vector = flow_dim_vector                                         &
+       )
+
+       dims_have_same_length = SIZE( flow_dim_vector ) .EQ. SIZE( hhl_dim_vector )
+       IF  ( .NOT. dims_have_same_length )  THEN
+          message = &
+             "Mesoscale data inconsistent. Number of grid points " //          &
+             "in dimension #" // TRIM( str( dim_idx ) ) //                     &
+             " do not match in the hhl and *-flow files (" //                  &
+             "HHL: len = " // TRIM( str( SIZE( hhl_dim_vector ) ) ) // ", " // &
+             "W: len = "   // TRIM( str( SIZE( flow_dim_vector ) ) )// ")."
+          CALL inifor_abort( 'validate_dataset', message )
+       ENDIF
+
+    ENDDO
+
+ END SUBROUTINE validate_dataset
+
  SUBROUTINE get_cosmo_grid( hhl_file, soil_file, rlon, rlat, hhl, hfl, depths, &
-                            d_depth, d_depth_rho_inv, phi_n, lambda_n,       &
-                            phi_equat,                                       &
-                            lonmin_cosmo, lonmax_cosmo,                      &
-                            latmin_cosmo, latmax_cosmo,                      &
+                            d_depth, d_depth_rho_inv, phi_n, lambda_n,         &
+                            phi_equat,                                         &
+                            lonmin_cosmo, lonmax_cosmo,                        &
+                            latmin_cosmo, latmax_cosmo,                        &
                             nlon, nlat, nlev, ndepths )