|
|
@ -258,7 +258,7 @@ VYLETData::Data VYLETData::ReadL() const |
|
|
|
{ |
|
|
|
{ |
|
|
|
time = R2Time((*vylet)[tcol - 1][iy * lons.size() + ix]); |
|
|
|
time = R2Time((*vylet)[tcol - 1][iy * lons.size() + ix]); |
|
|
|
lambda = (*vylet)[lcol - 1][iy * lons.size() + ix]; |
|
|
|
lambda = (*vylet)[lcol - 1][iy * lons.size() + ix]; |
|
|
|
days = static_cast<real>(time - start) / MDateTime::secondsperday; |
|
|
|
days = (time - start).D(); |
|
|
|
out(ix, iy) = lambda / days; |
|
|
|
out(ix, iy) = lambda / days; |
|
|
|
} |
|
|
|
} |
|
|
|
|
|
|
|
|
|
|
@ -284,7 +284,7 @@ VYLETData::Data VYLETData::ReadT() const |
|
|
|
const real acc = vylet->ParameterRValue("accuracy", 1.0); |
|
|
|
const real acc = vylet->ParameterRValue("accuracy", 1.0); |
|
|
|
const real tstep = 2.0 * M_PI * 1000.0 / acc; |
|
|
|
const real tstep = 2.0 * M_PI * 1000.0 / acc; |
|
|
|
|
|
|
|
|
|
|
|
const real maxdays = static_cast<real>(end - start) / MDateTime::secondsperday; |
|
|
|
const real maxdays = (end - start).D(); |
|
|
|
|
|
|
|
|
|
|
|
MDateTime time; |
|
|
|
MDateTime time; |
|
|
|
real days; |
|
|
|
real days; |
|
|
@ -297,7 +297,7 @@ VYLETData::Data VYLETData::ReadT() const |
|
|
|
x = (*vylet)[2][iy * lons.size() + ix]; |
|
|
|
x = (*vylet)[2][iy * lons.size() + ix]; |
|
|
|
y = (*vylet)[3][iy * lons.size() + ix]; |
|
|
|
y = (*vylet)[3][iy * lons.size() + ix]; |
|
|
|
time = R2Time((*vylet)[tcol - 1][iy * lons.size() + ix]); |
|
|
|
time = R2Time((*vylet)[tcol - 1][iy * lons.size() + ix]); |
|
|
|
days = static_cast<real>(time - start) / MDateTime::secondsperday; |
|
|
|
days = (time - start).D(); |
|
|
|
if(days <= tstep * 1.5) days = 0.0; |
|
|
|
if(days <= tstep * 1.5) days = 0.0; |
|
|
|
inside = x > xl && x < xr && y > yd && y < yu; |
|
|
|
inside = x > xl && x < xr && y > yd && y < yu; |
|
|
|
maxtime = days >= maxdays - tstep * 0.5; |
|
|
|
maxtime = days >= maxdays - tstep * 0.5; |
|
|
@ -505,7 +505,7 @@ VYLETData::Data VYLETData::ReadTmask() const |
|
|
|
const real acc = vylet->ParameterRValue("accuracy", 1.0); |
|
|
|
const real acc = vylet->ParameterRValue("accuracy", 1.0); |
|
|
|
const real tstep = 2.0 * M_PI * 1000.0 / acc; |
|
|
|
const real tstep = 2.0 * M_PI * 1000.0 / acc; |
|
|
|
|
|
|
|
|
|
|
|
const real maxdays = static_cast<real>(end - start) / MDateTime::secondsperday; |
|
|
|
const real maxdays = (end - start).D(); |
|
|
|
|
|
|
|
|
|
|
|
MDateTime time; |
|
|
|
MDateTime time; |
|
|
|
real days; |
|
|
|
real days; |
|
|
@ -515,7 +515,7 @@ VYLETData::Data VYLETData::ReadTmask() const |
|
|
|
for(size_t ix = 0; ix < lons.size(); ix++) |
|
|
|
for(size_t ix = 0; ix < lons.size(); ix++) |
|
|
|
{ |
|
|
|
{ |
|
|
|
time = R2Time((*vylet)[tcol - 1][iy * lons.size() + ix]); |
|
|
|
time = R2Time((*vylet)[tcol - 1][iy * lons.size() + ix]); |
|
|
|
days = static_cast<real>(time - start) / MDateTime::secondsperday; |
|
|
|
days = (time - start).D(); |
|
|
|
maxtime = days >= maxdays - tstep * 0.5; |
|
|
|
maxtime = days >= maxdays - tstep * 0.5; |
|
|
|
out(ix, iy) = maxtime ? 1.0 : NAN; |
|
|
|
out(ix, iy) = maxtime ? 1.0 : NAN; |
|
|
|
} |
|
|
|
} |
|
|
|