Changeset 4645 for palm/trunk/SOURCE/virtual_measurement_mod.f90

Timestamp:

Aug 24, 2020 1:55:58 PM (4 years ago)

Author:

suehring

Message:

Bugfix in output of E_UTM_soil coordinate

File:

• palm/trunk/SOURCE/virtual_measurement_mod.f90 (modified) (3 diffs)

palm/trunk/SOURCE/virtual_measurement_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix in output of E_UTM_soil coordinate
+!
+! 4642 2020-08-13 15:47:33Z suehring
 ! Do not set attribute bounds for time variable, as it refers to time_bounds which is not defined
 ! for non-aggregated quantities (according to data standard)
@@ -1143 +1146 @@
                          DO  k = kl, ku
                             meas_flag(k,j,i) = MERGE( IBSET( meas_flag(k,j,i), 0 ), 0,             &
-                                               BTEST( wall_flags_total_0(k,j,i), 0 ) )
+                                                      BTEST( wall_flags_total_0(k,j,i), 0 ) )
                          ENDDO
                       ENDDO
@@ -2083 +2086 @@
 !--          Output of Easting coordinate. Before output, recalculate EUTM.
              output_values_1d_target = init_model%origin_x                                         &
-               + REAL( vmea(l)%i(1:vmea(l)%ns_soil) + 0.5_wp, KIND = wp ) * dx                     &
+               + REAL( vmea(l)%i_soil(1:vmea(l)%ns_soil) + 0.5_wp, KIND = wp ) * dx                &
                * COS( init_model%rotation_angle * pi / 180.0_wp )                                  &
-               + REAL( vmea(l)%j(1:vmea(l)%ns_soil) + 0.5_wp, KIND = wp ) * dy                     &
+               + REAL( vmea(l)%j_soil(1:vmea(l)%ns_soil) + 0.5_wp, KIND = wp ) * dy                &
                * SIN( init_model%rotation_angle * pi / 180.0_wp )
              output_values_1d_pointer => output_values_1d_target