UW-Hydro/VIC

The Variable Infiltration Capacity (VIC) Macroscale Hydrologic Model

Arunav Nanda

@Arunavcoder123

I am trying to use VIC for my study area. The issue is manual calibration is taking a hell lot of time to be done. Is there a way to auto-calibrate my VIC model? Are there any documentations or resources for it? I tried the MOCOM-UA one, but the issue was it was only for the Wisconsin river. How do i change it for my study area? ANy insights and help would be appreciated.

Saad

@Saadi4469

Try using Kel Markert's auto calibration scripts from Github.

hoseinsalehi

@hoseinsalehi

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.

hoseinsalehi

@hoseinsalehi

Hi @Arunavcoder123 , Yes of course. It's my gmail account(hoseinsalehi111.hs@gmail.com). I will be happy to be able to help you.

Saad

@Saadi4469

Good day all,

Can anyone please share a source from where I can download a 500m resolution VIC input data-set for the Pacific North West in the US?

Jacob

@jschap1

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?

1 reply

hoseinsalehi

@hoseinsalehi

Hi
I want to consider the Dam in the simulation.
Can anyone help me prepare the parameters of the dam (which I should consider a lake)?
I use VIC.4.2.d version for simulation.

Saad

@Saadi4469

Is this group active? I need some help with the Tonic module!

I am using Tonic to convert VIC 4.2.d parameter file to NetCDF, but I am getting a format error, which I don't know how to fix as the parameter files read perfectly fine in VIC 4.2.d. basically the error is ValueError: invalid literal for int() with base 10:''.

This error has to do with line 1503 of grid_params.py

Saad

@Saadi4469

@hoseinsalehi https://vic.readthedocs.io/en/vic.4.2.d/Documentation/LakeParam/ use this link to get some ideas!

Aydin Bakhtar

@Aydinbakhtar

how to execute vic in cygwin?

is there any learning resources for VIC-ASSIST

Saad

@Saadi4469

@Aydinbakhtar instead of cygwin, use the Linux shell available in Windows 10. You can get the Linux shell from Windows store. Once you set that up, it's easy to make the VIC in the Linux shell and run the model. Hope this helps!

Aydin Bakhtar

@Aydinbakhtar

Hi is there anybody can help me about VIC-ASSIST?

Jacob

@jschap1

@Aydinbakhtar I don't know any learning resources for VIC-ASSIST, but you could try the R package VICmodel (https://github.com/Sibada/VICmodel) or my MATLAB toolbox, VICMATLAB (https://github.com/jschap1/VICMATLAB).

SOUBHIK MONDAL

@SOUBHIK86871142_twitter

WARNING: Root zone fractions sum to more than 1 ( = 1238.245850), normalizing fractions. If the sum is large, check that your vegetation parameter file is in the form - <zone 1 depth> <zone 1 fract> <zone 2 depth> <zone 2 fract> ...
Model run-time error... error while using stehekin.sample in VIC 4.0.6

Aydin Bakhtar

@Aydinbakhtar

Why everything in VIC is complicated I lost. :)
Plz give me some resources to learn how to run? and How to Make input files.
and which Is better VIC image, classic, R, python?

JACK LONG

@longbiao107995_twitter

the 1/16 VIC parameter dataset is really important.

Saad

@Saadi4469

When you try to create the elevation Bands in ArcGIS with the same interval as in the elevation band file, the distribution of the bands within the grid cells in ArcGIS does not match with the elevation band file. I am using the format_snow_params.py script to create the elevation band file. The purpose is to slice the gridcell according to the fractional areas of the Grid cells in the elevation band file.

RAJTANTRA LILHARE

@rlilhare

Hi @Saadi4469
I can see that you were working with the tonic module, any success in using that? I am actually working on that and wanted to do some discussion. Let me know your thoughts on this?

Saad

@Saadi4469

Hi @rlilhare, so I was able to successfully setup tonic and I did run it successfully on my data i.e. without any errors. Having said that the results were not accurate and it seems like that the library was created and tested for a particular dataset. Thus, I think the library needs tweaking in a sense that it can work with any dataset.

RAJTANTRA LILHARE

@rlilhare

Thank you @Saadi4469 for your reply. I just sent you a message (DM) on this app for further discussion.

Ehsan Jalilvand

@ehsan1325_twitter

Hello everyone,

I have created and defined both "veg_forcing" and "veg Param" in the global parameter file. Do you know which one is used by VIC for simulation?

I am using VIC 4.2.d

Saad

@Saadi4469

@ehsan1325_twitter "Veg Param"

Pritam Das

@pritamd47

Hi! I wanted to know if VIC can be built on the new M1 Chip? (Sorry if this has been discussed earlier, I did a search but didn't get any results.)

Alex Young

@alex-l-young

Hi, I'm wondering if anyone knows how the sample RVIC routing file here (UW-Hydro/VIC_sample_data/image/Stehekin/parameters/stehekin_parameters_01.rvic.prm.Stehekin.20150727.nc) was generated. Is this from the RVIC API?

singhharsh12

@singhharsh12

Hello Everyone,
I am running vic model for my region and want to calibrate in using optimization technique. Can anyone from you, who have done calibration help me regarding this. I will be very thankful for this. My email address is singhharsh127@gmail.com. Kindly help me regarding this.

singhharsh12

@singhharsh12

I am running VIC4.2d on water balance mode and getting some error. My fluxes have a lots of nan values. How can I fix it. Kindly help me in this.

"/home/rahul/Documents/vic_calib/fluxes/fluxes_23.000_77.250" has been
truncated or created for writing.
Initializing Forcing Data

Read meteorological forcing file
Model State Initialization
Running Model
Total Cumulative Water Error for Grid Cell = nan
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: 0
Total number of fallbacks in soil T profile: 0

cell: 171, lat: 23.0000, long: 77.5000
WARNING: Cv > 0.99 and Cv < 1.0 at grid cell 171, model assuming that bare soil is not to be run - fractions being adjusted to equal 1

Arunav Nanda

@Arunavcoder123

I also faced the similar issue of lot of 'nan' values. In order to resolve it I suggest you can write a simple python script that can read a flux file, go through all the elements in the file and replace the 'nan' values with zero. Then use the new flux files to rout the hydrological model.

farahmand-gh

@farahmand-gh

Hi, Kindly I have the same request for you please:
In a Ph.D. thesis, I'm trying to apply the Vic model but confronted a problem to access the Linux manual and documents. First, I examined the CYGWIN emulator as the VIC manual suggestion but It hasn't supported recently. Our model is ready but I have a problem running the model in Linux. I kindly request send me a manual or everything that helps us to finish the last part.

singhharsh12

@singhharsh12

When I am doing routing of the model following error is coming, How can i solve it?

/home/rahul/Documents/vic_calib/rout_fotran/src/fluxes/fluxes_23.0000_77.2500 NOT FOUND, INSERTING ZEROS
/home/rahul/Documents/vic_calib/rout_fotran/src/fluxes/fluxes_23.2500_77.2500 NOT FOUND, INSERTING ZEROS
/home/rahul/Documents/vic_calib/rout_fotran/src/fluxes/fluxes_23.5000_77.2500 NOT FOUND, INSERTING ZEROS
/home/rahul/Documents/vic_calib/rout_fotran/src/fluxes/fluxes_23.7500_77.2500 NOT FOUND, INSERTING ZEROS
/home/rahul/Documents/vic_calib/rout_fotran/src/fluxes/fluxes_24.0000_77.2500 NOT FOUND, INSERTING ZEROS
/home/rahul/Documents/vic_calib/rout_fotran/src/fluxes/fluxes_23.0000_77.0000 NOT FOUND, INSERTING ZEROS
/home/rahul/Documents/vic_calib/rout_fotran/src/fluxes/fluxes_23.2500_77.0000 NOT FOUND, INSERTING ZEROS

88 files not found, zero runoff/baseflow used
writing data...
Note: The following floating-point exceptions are signalling: IEEE_UNDERFLOW_FLAG IEEE_DENORMAL

mandrek7

@mandrek7

Hello @Saadi4469 how did you solve this error I am also stuck there : Choose output parameter
1 - Precipitation
2 - Evapotranspiration
3 - Runoff
4 - Base flow
5 - Snow Water Equivalent
6 - Soil moisture
Choose output (1 a 6)>3
Enter start year:2000
End year:2001
Traceback (most recent call last):
File "flux2nc.py", line 240, in <module>
main()
File "flux2nc.py", line 234, in main
flux2nc(sys.argv[1],sys.argv[2])
File "flux2nc.py", line 142, in flux2nc
latitude = float(str.split(f, sep="")[1])
TypeError: split() takes no keyword arguments

or can anyone assist with this error Choose output parameter
1 - Precipitation
2 - Evapotranspiration
3 - Runoff
4 - Base flow
5 - Snow Water Equivalent
6 - Soil moisture
Choose output (1 a 6)>3
Enter start year:2000
End year:2001
Traceback (most recent call last):
File "flux2nc.py", line 240, in <module>
main()
File "flux2nc.py", line 234, in main
flux2nc(sys.argv[1],sys.argv[2])
File "flux2nc.py", line 142, in flux2nc
latitude = float(str.split(f, sep="")[1])
TypeError: split() takes no keyword arguments

hoseinsalehi

@hoseinsalehi

Hello there,
I need help for preparing forcing files for VIC.5.
I download precipitation , mean temp, long and short wave radiation and surface pressure from ERA5 but I can't find a dataset for Vapor pressure.
Is there any way to calculate this or some database provide it?

Thanks in advance.

gia_nguyen

@giaMSU

If I ran my VIC model for 10years, how many yearly Land Cover layers do I need? Thanks

Arunav Nanda

@Arunavcoder123

@hoseinsalehi Hello there bro! I found this link in the VIC repository https://github.com/UW-Hydro/VIC_sample_data/tree/1554cb017c71e59e114b853af5808649f80ad836/image/Stehekin . Check this link.It contains the sample forcing datasets needed for running VIC-5.

hoseinsalehi

@hoseinsalehi

Dear @Arunavcoder123 , Thanks or your reply. I need a dataset for Vapor pressure, not a sample for this!

Shubham Goswami

@Shubh43_twitter

Can someone elaborate the flow_distance calculation for RVIC. I found this information on the RVIC page but the sample data given on github page makes it more confusing.

Saad

@Saadi4469

Hi can anyone can please tell if it's possible to convert the meteorological forcing data for VIC 4.2.d to readable text files?

hoseinsalehi

@hoseinsalehi

Dear @Saadi4469 , I think this pkg can halp you.
https://github.com/LSRathore/vic-auto/blob/master/force_auto.py

Saad

@Saadi4469

Hi Hosein, thanks but I already have the forcing dat a but it is in "bin" format, and right now I am just looking for a way to convert them to ASCII/text format files.

hoseinsalehi

@hoseinsalehi

Dear @Saadi4469 , you can convert your data to NetCDF and use that script. I think it's so easy.
https://www.researchgate.net/post/Convert-rainfall-dataset-from-binary-bin-to-netcdf-nc

Saad

@Saadi4469

@hoseinsalehi it sure is thanks so much for help brother.

ehsanj25

@ehsanj25

Hey Guys,

I am running a VIC simulation and finding an odd initialization error. Specifically, I am running an open loop VIC simulation from 01/01/2015-12/31/2017, and writing a statefile from this simulation on 2016/12/31. I then run an identical simulation (using the same parameter and forcing inputs), and initialize the model with the aforementioned 2016/12/31 statefile, and begin the simulation on 2017/1/1. Oddly enough, these two simulations do not produce identical soil moisture outputs!

Have you guys run into this issue before?

hoseinsalehi

@hoseinsalehi

Hi @ehsanj25
Could you send you error for me?
hoseinsalehi111.hs@gmail.com

Where communities thrive

People

Repo info

Activity