Changeset 4808 for palm/trunk/UTIL/inifor/src/inifor_transform.f90

Timestamp:

Dec 3, 2020 1:08:23 PM (3 years ago)

Author:

eckhard

Message:

Small error handling improvements

• Added offending height values to model-top error message
• Temporarily disabling dataset validation, enabling compatibility for HHL fields without time dimension

File:

• palm/trunk/UTIL/inifor/src/inifor_transform.f90 (modified) (3 diffs)

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

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Add offending height values to model-top error message
+!
+!
+! 4756 2020-10-26 10:05:58Z eckhard
 ! Fixed an error in surface pressure extrapolation where the cosmo grid was
 !    misinterpreted as palm grid
@@ -1105 +1109 @@
              palm_grid%w_verti(i,j,k,1:2) = - 2.0_wp
 
-             message = "PALM-4U grid extends above COSMO-DE model top."
+             message =                                                         &
+                "PALM-4U grid extends above COSMO model top. " //              &
+                "Height of COSMO model top: " //                               &
+                TRIM( real_to_str( column_top, format = '(F11.3)') ) //        &
+                " m, height of current PALM point: " //                        &
+                TRIM( real_to_str( current_height , format = '(F11.3)') ) //   &
+                " m, height of PALM model top: " //                            &
+                TRIM( real_to_str( MAXVAL(palm_grid%z) + palm_grid%z0, '(F11.3)') ) // " m."
              CALL inifor_abort('find_vertical_neighbours_and_weights', message)
 
@@ -1246 +1257 @@
              avg_grid%w(l,k_palm,1:2) = - 2.0_wp
 
-             message = "PALM-4U grid extends above COSMO-DE model top."
+             message =                                                         &
+                "PALM-4U grid extends above COSMO model top. " //              &
+                "Height of COSMO model top: " //                               &
+                TRIM( real_to_str( column_top, format = '(F11.3)') ) //        &
+                " m, height of current PALM point: " //                        &
+                TRIM( real_to_str( current_height , format = '(F11.3)') ) //   &
+                " m, height of PALM model top: " //                            &
+                TRIM( real_to_str( MAXVAL(avg_grid%z) + avg_grid%z0, '(F11.3)') ) // " m."
              CALL inifor_abort('find_vertical_neighbours_and_weights_average', message)