NoDAta when calculating time series in the subcatchments

[crisprie...@gmail.com]

Dear Wilko,

 I am getting 1e31 in all the time series

i am getting this for water balance, white is NoData plase, see the imagen attached

might it be due to something in the algorithm that sphy uses to compute ldd?

it would be great if someone can help on this as it is very important

basically i am not able to obtain the results in the sub-catchments, all i get is NoData , ie 1e31

I am also attaching the configuration file, i use the cmd, just in case someone can have a look

i am happy to send the .map and forcings too, i couldn't via here as it is too big

if someone could provide a hint about what is happening it would be wonderful, i am stucked on this 

looking forward to hearing from you 

thanks a lot

Cris

[WilcoTerink]

Hi Cris,

What I suspect is that the kc of at least one of your land use classes (as defined in the land_use_nominal.map) cannot be found in the kc.tbl lookup table. You need to make sure that all classes in your land use map have an associated kc value specified in that lookup table. Otherwise a missing value for ET will be generated for those cells, and as a result, the resulting fluxes will have NaNs as well.

If that is not the problem, then your ldd might not be sound.

Hope this fixes the problem.

Cheers,

Wilco

[crisprie...@gmail.com]

Hi Wilko,

i think is this ldd might not be sound.

how we can repair the ldd?

[crisprie...@gmail.com]

I also generated this kc table and land_use_nominal.map

could you please have a quick look at least to be sure this is not

it looks fine to me

[crisprie...@gmail.com]

I use the following in the configuration file

ETREF]

ETREF_FLAG      = 0

ETref           = 

Lat             = latitude.map

Gsc             = 0.0820

Tmin            = forcing\tmin\tmin

Tmax            = forcing\tmax\tmax

[crisprie...@gmail.com]

also, how do you generate the kc.tbl?

i use just a .tbl file - which is a .txt

and then the funciton pcrcalc luseCrisnominal.map=nominal(r1.map)

where r1 is the land use for 1 land use

i should use something like the below?

pcrcalc root_field.map = lookupscalar(root_field.tbl, soil.map)

[crisprie...@gmail.com]

ETA is 0 in some parts of the map, but the only Nans are out of the catchment - sea

[crisprie...@gmail.com]

i see that where ETA in the .map is NoData there is not ToRR, but i cannot see the problem with the kc table and the landuse nominal

maybe that i not used the looupscalar to create the table?

could you please have a quick look a the 3 files ? it would be very helpful

[WilcoTerink]

you can do this using the pcraster command: lddrepair

[WilcoTerink]

You only have one land use class? If so, then it looks fine to me

[WilcoTerink]

You make the *.tbl file manually. It's basically a txt file with a *.tbl extension. The IDs in your nominal land use map should all be present in this *.tbl file (values in the left column). The corresponding kc value for each class is then specified in the right column. See attached an example for a basin I am currently working at..

[WilcoTerink]

looks ok to me

[WilcoTerink]

see my comment in one of the posts below. The kc.tbl is not generated, but you have to make it yourself with a text editor.