Ignore:
Timestamp:
Nov 6, 2018 2:51:27 PM (3 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:
2 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE

  • palm/trunk/SOURCE/virtual_measurement_mod.f90

    r3473 r3494  
    2525! -----------------
    2626! $Id$
     27! Bugfixing
     28!
     29! 3473 2018-10-30 20:50:15Z suehring
    2730! Initial revision
    2831!
     
    3134! Authors:
    3235! --------
    33 ! @author Matthias Suehring
     36! @author Matthias Suehring and Klaus Ketelsen
    3437!
    3538!
     
    554557!--          lower-left origin of the model domain, which is the difference
    555558!--          betwen UTM coordinates.
    556 !              write(9,*) l, "before Eutm", e_utm(t,1:vmea(l)%dim_t(t)), "orig", init_model%origin_x
    557 !              write(9,*) l, "before Nutm", n_utm(t,1:vmea(l)%dim_t(t)), "orig", init_model%origin_y
    558              
     559           
    559560             e_utm(t,1:vmea(l)%dim_t(t)) = e_utm(t,1:vmea(l)%dim_t(t))         &
    560561                                         - init_model%origin_x
    561562             n_utm(t,1:vmea(l)%dim_t(t)) = n_utm(t,1:vmea(l)%dim_t(t))         &
    562563                                         - init_model%origin_y
    563 !              write(9,*) l, "Eutm", e_utm(t,1:vmea(l)%dim_t(t))
    564 !              write(9,*) l, "Nutm", n_utm(t,1:vmea(l)%dim_t(t))
    565 !              write(9,*)
    566564!
    567565!--          Compute grid indices relative to origin and check if these are
     
    585583                ks = k_prev
    586584                IF ( on_pe )  THEN
    587 !                    write(9,*) " I am onpe"
    588585                   ksurf = get_topography_top_index_ji( js, is, 's' )
    589586                   ks = MINLOC( ABS( zu - zw(ksurf) - z_ag(t,n) ), DIM = 1 ) - 1
     
    655652                   vmea(l)%j(ns)  = js
    656653                   vmea(l)%k(ns)  = ks
    657 !                    write(9,*) TRIM( vmea(l)%feature_type), l, "ns", ns, "ijk", vmea(l)%i(ns), vmea(l)%j(ns), vmea(l)%k(ns)
    658654                ENDIF
    659655!
     
    713709     INTEGER(iwp) ::  mm, j, i
    714710
    715 !      write(9,*) "sampling"
    716 !      flush(9)
    717711!
    718712!--  Loop over all stations. For each possible variable loop over all
     
    723717!--     the same indices for scalar and velocity components are used.
    724718!--     ToDo: Revise this later.
    725 !         DO  m = 1, vmea(l)%ns
    726 !            j = vmea(l)%j(m)
    727 !            i = vmea(l)%i(m)
     719        DO  m = 1, vmea(l)%ns
     720           j = vmea(l)%j(m)
     721           i = vmea(l)%i(m)
    728722!           
    729723!            IF ( i >= nxl  .AND.  i <= nxr  .AND.                               &
     
    731725!               IF ( surf_def_h(0)%start_index(j,i) <= &
    732726!                    surf_def_h(0)%end_index(j,i) )  THEN
    733 !                 
    734 !                  write(9,*) "sampled"
    735 !                  flush(9)
    736727!                  mm = surf_def_h(0)%end_index(j,i)
    737728!               
     
    763754!               ENDDO
    764755!            ENDIF
    765 !         ENDDO
     756        ENDDO
    766757
    767758     ENDDO
Note: See TracChangeset for help on using the changeset viewer.