Changeset 3997 for palm/trunk/UTIL/inifor/src/inifor_grid.f90

Timestamp:

May 23, 2019 12:35:57 PM (5 years ago)

Author:

eckhard

Message:

inifor: Read origin_z from static driver if given; alert user to warnings

File:

• palm/trunk/UTIL/inifor/src/inifor_grid.f90 (modified) (5 diffs)

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

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Read origin_z from static driver if given on the command line
+! Warn if origin_z is specified twice (by static driver and --elevation)
+! 
+! 
+! 3866 2019-04-05 14:25:01Z eckhard
 ! Use PALM's working precision
 ! Catch errors while reading namelists
@@ -140 +145 @@
                SNAME, LNAME, PATH, FORCING_STEP, WATER_ID, FILL_ITERATIONS,    &
                BETA, P_SL, T_SL, BETA, RD, RV, G, P_REF, RD_PALM, CP_PALM,     &
-               RHO_L, OMEGA, HECTO, wp, iwp
+               RHO_L, OMEGA, HECTO, wp, iwp,                                   &
+               PIDS_ORIGIN_LON,                                                &
+               PIDS_ORIGIN_LAT,                                                &
+               PIDS_ORIGIN_Z
     USE inifor_io,                                                             &
         ONLY:  get_cosmo_grid, get_netcdf_attribute, get_netcdf_dim_vector,    &
@@ -440 +448 @@
 
     output_file%name = cfg%output_file
-    z0 = cfg%z0
-    p0 = cfg%p0
 
     init_variables_required = .TRUE.
@@ -536 +542 @@
     IF (TRIM(cfg%static_driver_file) .NE. '')  THEN
 
-       origin_lon = get_netcdf_attribute(cfg%static_driver_file, 'origin_lon')
-       origin_lat = get_netcdf_attribute(cfg%static_driver_file, 'origin_lat')
+       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 '"                   &
@@ -556 +564 @@
 
     CALL report('setup_parameters', message)
+
+    IF ( cfg%z0_is_set .AND. 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('setup_parameters', message)
+
+       z0 = cfg%z0
+
+    END IF
+
+    p0 = cfg%p0
 
     CALL log_runtime('time', 'read')