Changeset 4ad96c5 in flexpart.git for src/netcdf_output_mod.f90

Timestamp:

Dec 19, 2018, 2:09:04 PM (5 years ago)

Author:

Sabine <sabine.eckhardt@…>

Branches:

master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug

Children:

3d7eebf

Parents:

c7d1052

Message:

BUGFIX: wrong Releasepointcoordinates in netcdf file

File:

• src/netcdf_output_mod.f90 (modified) (2 diffs)

src/netcdf_output_mod.f90

@@ -113 +113 @@
   logical, parameter :: write_area = .false.
 
+  ! coordinate transformation from internal to world coord
+  real :: xp1,yp1,xp2,yp2
 contains
 
@@ -668 +670 @@
        call nf90_err(nf90_put_var(ncid, relendID, ireleaseend(i), (/i/)))
        call nf90_err(nf90_put_var(ncid, relkindzID, kindz(i), (/i/)))
-       call nf90_err(nf90_put_var(ncid, rellng1ID, xpoint1(i), (/i/)))
-       call nf90_err(nf90_put_var(ncid, rellng2ID, xpoint2(i), (/i/)))
-       call nf90_err(nf90_put_var(ncid, rellat1ID, ypoint1(i), (/i/)))
-       call nf90_err(nf90_put_var(ncid, rellat2ID, ypoint2(i), (/i/)))
+       xp1=xpoint1(i)*dx+xlon0
+       yp1=ypoint1(i)*dy+ylat0
+       xp2=xpoint2(i)*dx+xlon0
+       yp2=ypoint2(i)*dy+ylat0
+       call nf90_err(nf90_put_var(ncid, rellng1ID, xp1, (/i/)))
+       call nf90_err(nf90_put_var(ncid, rellng2ID, xp2, (/i/)))
+       call nf90_err(nf90_put_var(ncid, rellat1ID, yp1, (/i/)))
+       call nf90_err(nf90_put_var(ncid, rellat2ID, yp2, (/i/)))
        call nf90_err(nf90_put_var(ncid, relzz1ID, zpoint1(i), (/i/)))
        call nf90_err(nf90_put_var(ncid, relzz2ID, zpoint2(i), (/i/)))