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

Timestamp:

Jun 19, 2020 11:56:30 AM (4 years ago)

Author:

eckhard

Message:

Handle COSMO soil data with and without additional surface temperature

File:

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

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

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Handle COSMO soil data with and without additional surface temperature
+! 
+! 
+! 4553 2020-06-03 16:34:15Z eckhard
 ! Option --help now points the user to INIFOR's wiki page
 ! Remove deprecated command-line options -clon and -clat
@@ -146 +150 @@
     USE inifor_defs,                                                           &
         ONLY:  DATE, SNAME, PATH, PI, TO_RADIANS, TO_DEGREES, VERSION,         &
-               NC_DEPTH_NAME, NC_HHL_NAME, NC_RLAT_NAME, NC_RLON_NAME,         &
-               NC_ROTATED_POLE_NAME, NC_POLE_LATITUDE_NAME,                    &
+               NC_DEPTH_DIM_IDX, NC_DEPTH_NAME, NC_HHL_NAME, NC_RLAT_NAME,     &
+               NC_RLON_NAME, NC_ROTATED_POLE_NAME, NC_POLE_LATITUDE_NAME,      &
                NC_POLE_LONGITUDE_NAME, RHO_L, iwp, wp,                         &
                PIDS_ORIGIN_LON, PIDS_ORIGIN_LAT, PIDS_ORIGIN_Z
     USE inifor_types
     USE inifor_util,                                                           &
-        ONLY:  add_hours_to, reverse, str, real_to_str
+        ONLY:  add_hours_to, nearly_equal, reverse, str, real_to_str
     USE netcdf
 
@@ -355 +359 @@
 
     start = (/ 1, 1, 1 /)
-    IF ( TRIM(in_var%name) .EQ. 'T_SO' )  THEN
+    IF ( TRIM(in_var%name) .EQ. 'T_SO' .AND.                                &
+         in_var%has_redundant_first_level )  THEN
+       
 !
 !--    Skip depth = 0.0 for T_SO and reduce number of depths from 9 to 8
@@ -1578 +1584 @@
  END SUBROUTINE set_palm_origin
 
+
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+! This function is meant to check weather a COSMO soil variable has an
+! additional and redunant surface value at depth = 0.0. For instance operational
+! DWD COSMO output contains the surface temperature in T_SO as a copy of the
+! values in the first soil layer.
+!------------------------------------------------------------------------------!
+ LOGICAL FUNCTION has_surface_value( soil_var, filename ) 
+
+    TYPE(nc_var), INTENT(IN)     ::  soil_var
+    CHARACTER(LEN=*), INTENT(IN) ::  filename
+    
+    CHARACTER(LEN=NF90_MAX_NAME) ::  dimname
+    REAL(wp), ALLOCATABLE        ::  depths(:)
+
+    CALL get_dimension_vector_of_variable(                                     &
+       soil_var%name,                                                          &
+       dim_idx = NC_DEPTH_DIM_IDX,                                             &
+       filename = filename,                                                    &
+       dim_vector = depths                                                     &
+    )
+
+    has_surface_value = nearly_equal( depths(1), 0.0_wp, 10 * EPSILON(1.0_wp) )
+
+ END FUNCTION has_surface_value 
+
+
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+! This routine reads the dim_idx-th dimension vector of the variable varname
+! from netCDF file filename. It is used for finding the depth coordinate vector
+! of COSMO soil variables without knowing its name.
+!------------------------------------------------------------------------------!
+ SUBROUTINE get_dimension_vector_of_variable( varname, dim_idx, filename, dim_vector )
+    CHARACTER(LEN=*), INTENT(IN)              ::  varname, filename
+    INTEGER, INTENT(IN)                       ::  dim_idx
+
+    REAL(wp), INTENT(OUT), ALLOCATABLE ::  dim_vector(:)
+
+    INTEGER                            ::  dimids(NF90_MAX_VAR_DIMS)
+    INTEGER                            ::  varid
+    CHARACTER(LEN=NF90_MAX_NAME)       ::  dimname
+
+    INTEGER ::  ncid
+
+    IF ( nf90_open( TRIM( filename ), NF90_NOWRITE, ncid ) .EQ. NF90_NOERR )  THEN
+
+!
+!--    get id of variable varname
+       CALL check( nf90_inq_varid( ncid, TRIM( varname ), varid ) )
+       
+!
+!--    get dimension ids of variable with varid
+       CALL check( nf90_inquire_variable( ncid, varid, dimids = dimids ) )
+
+!
+!--    get name of dim_idx-th dimension variable
+       CALL check( nf90_inquire_dimension( ncid, dimids(dim_idx), name = dimname ) )
+       CALL check( nf90_close( ncid ) )
+
+    ELSE
+
+       message = "Failed to open file '" // TRIM(filename) // "'."
+       CALL inifor_abort('get_netcdf_variable', message)
+
+    ENDIF
+
+    ! get dimension vector with dimname
+    CALL get_netcdf_dim_vector( filename, dimname, dim_vector )
+
+ END SUBROUTINE get_dimension_vector_of_variable
+
+
  END MODULE inifor_io
 #endif