Changeset 3913


Ignore:
Timestamp:
Apr 17, 2019 3:12:28 PM (5 years ago)
Author:
gronemeier
Message:

Bugfix: rotate positions of measurements before writing them into file (virtual_measurement_mod)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/virtual_measurement_mod.f90

    r3910 r3913  
    2525! -----------------
    2626! $Id$
     27! Bugfix: rotate positions of measurements before writing them into file
     28!
     29! 3910 2019-04-17 11:46:56Z suehring
    2730! Bugfix in rotation of UTM coordinates
    2831!
     
    962965     IMPLICIT NONE
    963966         
    964      INTEGER(iwp) ::  i         !< running index over IO blocks   
    965      INTEGER(iwp) ::  l         !< running index over all stations
    966      INTEGER(iwp) ::  n         !< running index over all measured variables at a station
     967     INTEGER(iwp) ::  i  !< running index over IO blocks   
     968     INTEGER(iwp) ::  l  !< running index over all stations
     969     INTEGER(iwp) ::  n  !< running index over all measured variables at a station
    967970!
    968971!--  Header output on each PE
     
    994997                 WRITE ( 27 )  vmea(l)%ns
    995998                 WRITE ( 27 )  'E_UTM                             '
    996                  WRITE ( 27 )  init_model%origin_x +                           &
    997                         REAL( vmea(l)%i(1:vmea(l)%ns) + 0.5_wp, KIND = wp ) * dx
     999                 WRITE ( 27 )  init_model%origin_x                             &
     1000                    + REAL( vmea(l)%i(1:vmea(l)%ns) + 0.5_wp, KIND = wp ) * dx &
     1001                    * COS( init_model%rotation_angle * pi / 180.0_wp )         &
     1002                    + REAL( vmea(l)%j(1:vmea(l)%ns) + 0.5_wp, KIND = wp ) * dy &
     1003                    * SIN( init_model%rotation_angle * pi / 180.0_wp )
    9981004                 WRITE ( 27 )  'N_UTM                             '
    999                  WRITE ( 27 )  init_model%origin_y +                           &
    1000                         REAL( vmea(l)%j(1:vmea(l)%ns) + 0.5_wp, KIND = wp ) * dy
     1005                 WRITE ( 27 )  init_model%origin_y                             &
     1006                    - REAL( vmea(l)%i(1:vmea(l)%ns) + 0.5_wp, KIND = wp ) * dx &
     1007                    * SIN( init_model%rotation_angle * pi / 180.0_wp )         &
     1008                    + REAL( vmea(l)%j(1:vmea(l)%ns) + 0.5_wp, KIND = wp ) * dy &
     1009                    * COS( init_model%rotation_angle * pi / 180.0_wp )
    10011010                 WRITE ( 27 )  'Z_AG                              '
    10021011                 WRITE ( 27 )  vmea(l)%z_ag(1:vmea(l)%ns)
     
    10121021                    WRITE ( 27 )  vmea(l)%ns_soil
    10131022                    WRITE ( 27 )  'E_UTM soil                        '
    1014                     WRITE ( 27 )  init_model%origin_x +                        &
    1015                            REAL( vmea(l)%i_soil(1:vmea(l)%ns_soil) + 0.5_wp,   &
    1016                                  KIND = wp ) * dx
     1023                    WRITE ( 27 )  init_model%origin_x                          &
     1024                       + REAL( vmea(l)%i_soil(1:vmea(l)%ns_soil) + 0.5_wp,     &
     1025                               KIND = wp ) * dx                                &
     1026                       * COS( init_model%rotation_angle * pi / 180.0_wp )      &
     1027                       + REAL( vmea(l)%j_soil(1:vmea(l)%ns_soil) + 0.5_wp,     &
     1028                               KIND = wp ) * dy                                &
     1029                       * SIN( init_model%rotation_angle * pi / 180.0_wp )
    10171030                    WRITE ( 27 )  'N_UTM soil                        '
    1018                     WRITE ( 27 )  init_model%origin_y +                        &
    1019                            REAL( vmea(l)%j_soil(1:vmea(l)%ns_soil) + 0.5_wp,   &
    1020                                  KIND = wp ) * dy
     1031                    WRITE ( 27 )  init_model%origin_y                          &
     1032                       - REAL( vmea(l)%i_soil(1:vmea(l)%ns_soil) + 0.5_wp,     &
     1033                               KIND = wp ) * dx                                &
     1034                       * SIN( init_model%rotation_angle * pi / 180.0_wp )      &
     1035                       + REAL( vmea(l)%j_soil(1:vmea(l)%ns_soil) + 0.5_wp,     &
     1036                               KIND = wp ) * dy                                &
     1037                       * COS( init_model%rotation_angle * pi / 180.0_wp )
    10211038                    WRITE ( 27 )  'DEPTH                             '
    10221039                    WRITE ( 27 )  vmea(l)%depth(1:vmea(l)%ns_soil)
Note: See TracChangeset for help on using the changeset viewer.