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

Timestamp:

Nov 6, 2018 2:51:27 PM (5 years ago)

Author:

suehring

Message:

Surface output revised and some bugs are fixed + new post-processing tool to convert binary surface output to Paraview readable VTK files

Location:

palm/trunk/SOURCE

Files:

• . (modified) (1 prop)
• virtual_measurement_mod.f90 (modified) (9 diffs)

palm/trunk/SOURCE/virtual_measurement_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfixing
+! 
+! 3473 2018-10-30 20:50:15Z suehring
 ! Initial revision
 ! 
@@ -31 +34 @@
 ! Authors:
 ! --------
-! @author Matthias Suehring
+! @author Matthias Suehring and Klaus Ketelsen
 ! 
 !
@@ -554 +557 @@
 !--          lower-left origin of the model domain, which is the difference 
 !--          betwen UTM coordinates.
-!              write(9,*) l, "before Eutm", e_utm(t,1:vmea(l)%dim_t(t)), "orig", init_model%origin_x
-!              write(9,*) l, "before Nutm", n_utm(t,1:vmea(l)%dim_t(t)), "orig", init_model%origin_y
-             
+            
              e_utm(t,1:vmea(l)%dim_t(t)) = e_utm(t,1:vmea(l)%dim_t(t))         &
                                          - init_model%origin_x
              n_utm(t,1:vmea(l)%dim_t(t)) = n_utm(t,1:vmea(l)%dim_t(t))         &
                                          - init_model%origin_y
-!              write(9,*) l, "Eutm", e_utm(t,1:vmea(l)%dim_t(t)) 
-!              write(9,*) l, "Nutm", n_utm(t,1:vmea(l)%dim_t(t)) 
-!              write(9,*)
 !
 !--          Compute grid indices relative to origin and check if these are
@@ -585 +583 @@
                 ks = k_prev
                 IF ( on_pe )  THEN
-!                    write(9,*) " I am onpe" 
                    ksurf = get_topography_top_index_ji( js, is, 's' )
                    ks = MINLOC( ABS( zu - zw(ksurf) - z_ag(t,n) ), DIM = 1 ) - 1
@@ -655 +652 @@
                    vmea(l)%j(ns)  = js
                    vmea(l)%k(ns)  = ks
-!                    write(9,*) TRIM( vmea(l)%feature_type), l, "ns", ns, "ijk", vmea(l)%i(ns), vmea(l)%j(ns), vmea(l)%k(ns)
                 ENDIF
 !
@@ -713 +709 @@
      INTEGER(iwp) ::  mm, j, i
 
-!      write(9,*) "sampling"
-!      flush(9)
 !
 !--  Loop over all stations. For each possible variable loop over all 
@@ -723 +717 @@
 !--     the same indices for scalar and velocity components are used.
 !--     ToDo: Revise this later.
-!         DO  m = 1, vmea(l)%ns
-!            j = vmea(l)%j(m)
-!            i = vmea(l)%i(m)
+        DO  m = 1, vmea(l)%ns
+           j = vmea(l)%j(m)
+           i = vmea(l)%i(m)
 !            
 !            IF ( i >= nxl  .AND.  i <= nxr  .AND.                               &
@@ -731 +725 @@
 !               IF ( surf_def_h(0)%start_index(j,i) <= &
 !                    surf_def_h(0)%end_index(j,i) )  THEN
-!                 
-!                  write(9,*) "sampled"
-!                  flush(9)
 !                  mm = surf_def_h(0)%end_index(j,i)
 !               
@@ -763 +754 @@
 !               ENDDO
 !            ENDIF
-!         ENDDO
+        ENDDO
 
      ENDDO