Feel free to join the next Helmholtz Hacky Hour #30 on Wednesday, June 16, 2021 from 2PM to 3PM!

Commit 0cae46cc authored by Friedrich Boeing's avatar Friedrich Boeing
Browse files

correction of wrong date allocation

parent fa8d3201
......@@ -45,7 +45,7 @@ program ED_Kriging
character(256) :: fname
character(256) :: author_name
character(256) :: vname_data
integer(i4) :: i, j, k
integer(i4) :: i, j, k, t
integer(i4) :: icell ! loop varaible for cells
integer(i4) :: jDay ! loop variable - current julian day
integer(i4) :: itimer
......@@ -154,8 +154,10 @@ program ED_Kriging
tmp_array(i, gridMeteo%nrows - j + 1, :) = cell(k)%z
end do
end do
do i = 1, jEnd - jStart + 1
tmp_time(i) = i
t = 0
do i = 1, jEnd - jStart + 1
tmp_time(i) = t
t = t + 1
end do
call nc_time%setData(tmp_time)
......
......@@ -32,14 +32,14 @@ CONTAINS
! create dimensions
dim_x = nc%setDimension("x", gridMeteo%ncols)
dim_y = nc%setDimension("y", gridMeteo%nrows)
dim_time = nc%setDimension("time", 0)
dim_time = nc%setDimension("time", -1)
! create variables
var_time = nc%setVariable('time', "i32", (/dim_time/))
! add some variable attributes
call var_time%setAttribute("units", "days since " // trim(num2str(yStart, form='(I4)')) // "-"// &
trim(num2str(mStart, form='(I0.2)')) // "-" // &
trim(num2str(dStart, form='(I0.2)')) // "-" // "12:00:00")
trim(num2str(dStart, form='(I0.2)')) // "-" // "00:00:00")
allocate(dummy(gridMeteo%ncols, gridMeteo%nrows))
do i = 1, gridMeteo%nrows
......
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