Commit 4ad96c54 authored by Sabine's avatar Sabine

BUGFIX: wrong Releasepointcoordinates in netcdf file

parent c7d1052e
......@@ -112,6 +112,8 @@ module netcdf_output_mod
logical, parameter :: write_vol = .false.
logical, parameter :: write_area = .false.
! coordinate transformation from internal to world coord
real :: xp1,yp1,xp2,yp2
contains
!****************************************************************
......@@ -667,10 +669,14 @@ subroutine writeheader_netcdf(lnest)
call nf90_err(nf90_put_var(ncid, relstartID, ireleasestart(i), (/i/)))
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/)))
call nf90_err(nf90_put_var(ncid, relpartID, npart(i), (/i/)))
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment