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

Timestamp:

May 18, 2020 1:45:35 PM (4 years ago)

Author:

eckhard

Message:

inifor: Fix issue where --elevation/-z option was ignored, make it mandatory

File:

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

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

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Make setting the vertical PALM origin mandatory
+! 
+! 
+! 4523 2020-05-07 15:58:16Z eckhard
 ! respect integer working precision (iwp) specified in inifor_defs.f90
 ! 
@@ -139 +143 @@
                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
+               NC_POLE_LONGITUDE_NAME, RHO_L, iwp, wp,                         &
+               PIDS_ORIGIN_LON, PIDS_ORIGIN_LAT, PIDS_ORIGIN_Z
     USE inifor_types
     USE inifor_util,                                                           &
@@ -431 +436 @@
     cfg%soil_prefix_is_set = .FALSE.
     cfg%soilmoisture_prefix_is_set = .FALSE.
+    cfg%static_driver_is_set = .FALSE.
     cfg%ug_defined_by_user = .FALSE.
     cfg%vg_defined_by_user = .FALSE.
@@ -503 +509 @@
 
              CASE( '-static', '-t', '--static-driver' )
+                cfg%static_driver_is_set = .TRUE.
                 CALL get_option_argument( i, arg )
                 cfg%static_driver_file = TRIM(arg)
@@ -762 +769 @@
        message = "You specified only one component of the geostrophic " // &
                  "wind. Please specify either both or none."
+       CALL inifor_abort( 'validate_config', message )
+    ENDIF
+
+    IF ( .NOT. cfg%static_driver_is_set .AND. .NOT. cfg%z0_is_set )  THEN
+       message =                                                               &
+          "The vertical origin of the PALM grid has not been defined. " // NEW_LINE(" ") // &
+          "Please specify the right value for your setup by either " // NEW_LINE(" ") // &
+          "  - using the command-line option --elevation <height above sea level>, or by" // NEW_LINE(" ") // &
+          "  - specifying a static driver file using --static <filename> in order to use " // NEW_LINE(" ") // &
+          "    use the value of origin_z (and origin_lon and origin_lat) specifed therein."
        CALL inifor_abort( 'validate_config', message )
     ENDIF
@@ -1512 +1529 @@
  END SUBROUTINE check
 
+
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> Setup the origin of the PALM coordinate system based on what is given in the
+!> INIFOR namelist file and via an optional static driver.
+!------------------------------------------------------------------------------!
+ SUBROUTINE set_palm_origin( cfg, namelist_longitude, namelist_latitude,       &
+                             origin_lon, origin_lat, z0 )
+
+    TYPE(inifor_config), INTENT(IN) ::  cfg
+    REAL(wp), INTENT(IN)            ::  namelist_longitude, namelist_latitude
+    REAL(wp), INTENT(OUT)           ::  origin_lon, origin_lat, z0
+
+    message = 'Reading PALM-4U origin from'
+    IF ( TRIM( cfg%static_driver_file ) .NE. '' )  THEN
+
+       origin_lon = get_netcdf_attribute( cfg%static_driver_file, TRIM( PIDS_ORIGIN_LON ) )
+       origin_lat = get_netcdf_attribute( cfg%static_driver_file, TRIM( PIDS_ORIGIN_LAT ) )
+       z0         = get_netcdf_attribute( cfg%static_driver_file, TRIM( PIDS_ORIGIN_Z ) )
+
+       message = TRIM(message) // " static driver file '"                      &
+                               // TRIM( cfg%static_driver_file ) // "'"
+
+
+    ELSE
+
+       origin_lon = namelist_longitude
+       origin_lat = namelist_latitude
+
+       message = TRIM( message ) // " namlist file '"                          &
+                                 // TRIM( cfg%namelist_file ) // "'"
+
+    ENDIF
+    origin_lon = origin_lon * TO_RADIANS
+    origin_lat = origin_lat * TO_RADIANS
+
+    CALL report('set_palm_origin', message)
+
+    IF ( cfg%z0_is_set )  THEN
+       z0 = cfg%z0
+       IF ( TRIM( cfg%static_driver_file ) .NE. '' )  THEN
+          message = 'You specified both --static-driver/-t and --elevation/-z0. ' // &
+                    'Using the command line value (' // TRIM( real_to_str_f( cfg%z0 ) ) // &
+                    ') instead of static driver value (' // TRIM( real_to_str_f( z0 ) ) // ').'
+          CALL warn( 'set_palm_origin', message )
+       ENDIF
+    ENDIF
+
+ END SUBROUTINE set_palm_origin
+
  END MODULE inifor_io
 #endif