Changeset 3904 for palm/trunk/SOURCE

Timestamp:

Apr 16, 2019 6:22:51 PM (6 years ago)

Author:

gronemeier

Message:

rotate coordinates of stations by given rotation_angle

File:

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

palm/trunk/SOURCE/virtual_measurement_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Rotate coordinates of stations by given rotation_angle
+! 
+! 3876 2019-04-08 18:41:49Z knoop
 ! Remove print statement
 ! 
@@ -82 +85 @@
     USE arrays_3d,                                                             &
         ONLY:  q, pt, u, v, w, zu, zw
+
+    USE basic_constants_and_equations_mod,                                     &
+        ONLY:  pi
 
     USE chem_gasphase_mod,                                                     &
@@ -696 +702 @@
 !--          First, compute relative x- and y-coordinates with respect to the 
 !--          lower-left origin of the model domain, which is the difference 
-!--          betwen UTM coordinates. Note, if the origin is not correct, the
+!--          between UTM coordinates. Note, if the origin is not correct, the
 !--          virtual sites will be misplaced. 
              e_utm(t,1:dim_ntime) = e_utm(t,1:dim_ntime) - init_model%origin_x
              n_utm(t,1:dim_ntime) = n_utm(t,1:dim_ntime) - init_model%origin_y
+             e_utm(t,1:dim_ntime) = COS( init_model%rotation_angle * pi / 180.0_wp ) &
+                                    * e_utm(t,1:dim_ntime)                           &
+                                  - SIN( init_model%rotation_angle * pi / 180.0_wp ) &
+                                    * n_utm(t,1:dim_ntime)
+             n_utm(t,1:dim_ntime) = SIN( init_model%rotation_angle * pi / 180.0_wp ) &
+                                    * e_utm(t,1:dim_ntime)                           &
+                                  + COS( init_model%rotation_angle * pi / 180.0_wp ) &
+                                    * n_utm(t,1:dim_ntime)
 !
 !--          Determine the individual time coordinate length for each station and
@@ -1078 +1092 @@
         ONLY:  exner, pt, q, u, v, w
 
-    USE basic_constants_and_equations_mod,                                     &
-        ONLY:  pi
-    
     USE radiation_model_mod,                                                   &
         ONLY:  radiation