1)
Hi @Saadi4469 , I am not familiar with the flux2nc.py script, but converting from NetCDF to geotiff is fairly straightforward in R or MATLAB using the raster library in R or ncread and geotiffwrite functions in MATLAB.
If you are still using the flux2nc.py script, you might be able to use the snow water equivalent output and make a snow-cover mask from that.
cell: 287909, lat: 44.4062, long: -110.7188
"/mnt/d/Spring_2020/VIC/Forcing/data44.40625-110.71875" has been
opened for reading.
"/mnt/d/Spring_2020/VIC/VIC_Output/fluxes44.40625-110.71875" has been
truncated or created for writing.
Initializing Forcing Data
Read meteorological forcing file
Model State Initialization
Running Model
Maximum Energy Error: 0 135529.9936 135529.9936
Maximum Energy Error: 1 -26465462969.9630 -26465327439.9694
Maximum Energy Error: 2 -38711394088.5370 -32588360764.2532
Maximum Energy Error: 13 -927981257130.5930 -82156948242.1886
Maximum Energy Error: 14 -945935082001.6515 -143855386367.8646
Maximum Energy Error: 16 -948483648921.6132 -244228530258.2426
Maximum Energy Error: 67 -949117638426.6837 -509439174432.8633
Maximum Energy Error: 70 -949668692328.4795 -528280520730.2501
Maximum Energy Error: 77 -949806779513.0524 -544033664709.3660
Maximum Energy Error: 121 -981496935027.9856 -682611183273.2134
Maximum Energy Error: 122 -1011284747723.4287 -685305228883.4609
Maximum Energy Error: 124 -1011819008255.3538 -690320868419.7130
Maximum Energy Error: 163 -1011862276433.5739 -754886686051.3418
Maximum Energy Error: 273 -1011870514014.7301 -849742114138.4972
Maximum Energy Error: 482 -1019504728169.0823 -701651311958.5828
Maximum Energy Error: 497 -1019538206516.0092 -709517131515.5535
Maximum Energy Error: 498 -1019573828377.6522 -710139735324.5134
Maximum Energy Error: 499 -1019742778834.6741 -710760182305.4957
Maximum Energy Error: 668 -1020623171278.1158 -779199444826.9565
Maximum Energy Error: 675 -1020624667737.6036 -781150896042.3878
Total Cumulative Water Error for Grid Cell = -0.0000
Total Cumulative Energy Error for Grid Cell = -745953788780.1587
Total number of fallbacks in Tfoliage: 0
Total number of fallbacks in Tcanopy: 0
Total number of fallbacks in Tsnowsurf: 0
Total number of fallbacks in Tsurf: 512878
Hi @Arunavcoder123 In below link there is an pkg for auto calibraton:
https://github.com/thanhiwer/VICResOpt/tree/master/Routing/OptCalib
For Kel Markert scripts, I've tested them, they work.
Hi, I am getting an error when I run VIC with the Image Driver.
[ERROR] ../shared_image/src/get_nc_field.c:65: errno: NetCDF: Index exceeds dimension bound: Error getting values for tas in ${PATH}/data/netcdf_forcings/forc_tuo.2009.nc
I have checked my forcing files against the Stehekin sample inputs, and they appear to be essentially the same format. The only difference is the "_FillValue" parameter, which is set to 9.96921e+36f in the Stehekin forcings and unset in my forcing file (these cells are masked out by run_cell in the parameter file, anyway, so I don't think it matters). I added some printf lines to VIC to see what was causing the error, and it looks like the variables start, count, and var are set to very large values (i.e. -55532992, -55533024, and 573470640, respectively) in the get_nc_field_double
function, which reads the forcing data. I've been trying to figure out how the "start" and "count" variables are set, without much luck. Is there anyone who might be able to put me on the right track?