Changeset 3913 for palm/trunk/SOURCE/virtual_measurement_mod.f90
- Timestamp:
- Apr 17, 2019 3:12:28 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/virtual_measurement_mod.f90
r3910 r3913 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Bugfix: rotate positions of measurements before writing them into file 28 ! 29 ! 3910 2019-04-17 11:46:56Z suehring 27 30 ! Bugfix in rotation of UTM coordinates 28 31 ! … … 962 965 IMPLICIT NONE 963 966 964 INTEGER(iwp) :: i 965 INTEGER(iwp) :: l 966 INTEGER(iwp) :: n 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 967 970 ! 968 971 !-- Header output on each PE … … 994 997 WRITE ( 27 ) vmea(l)%ns 995 998 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 ) 998 1004 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 ) 1001 1010 WRITE ( 27 ) 'Z_AG ' 1002 1011 WRITE ( 27 ) vmea(l)%z_ag(1:vmea(l)%ns) … … 1012 1021 WRITE ( 27 ) vmea(l)%ns_soil 1013 1022 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 ) 1017 1030 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 ) 1021 1038 WRITE ( 27 ) 'DEPTH ' 1022 1039 WRITE ( 27 ) vmea(l)%depth(1:vmea(l)%ns_soil)
Note: See TracChangeset
for help on using the changeset viewer.