Vector-based implementation of VIC

[pyested]

Hi! I'm planning to carry out a vector-based implementation of VIC for my study catchments. My idea is to use non-gridded HRUs as a way to explore the effect of other spatial configurations different from the classical gridded domains.

I would like to know is such feature is currently implemented in the Image driver. While for the Classic driver I find it easy to implement (every HRU would correspond to a line in the soil parameter file), I think it is not so simple for the Image driver as most of the variables and parameters use latitude and longitude as coordinates and not the HRU id.

As a possible workaround for the Image driver I thought that using a mask in the domain file indicating 1 only for the HRUs I'm interested in could do the job. This is still going to generate 3d outputs, but later I can mask the resulting netcdf.

Thanks a lot for your support.