Where communities thrive


  • Join over 1.5M+ people
  • Join over 100K+ communities
  • Free without limits
  • Create your own community
People
Repo info
Activity
singhharsh12
@singhharsh12
Hello Everyone,
I am using vic model 4.2d for my basin. I have used Indian Meteorological Data for temperature and rainfall. I have run the model for 0.05 degree grid size.
I have also routed the model for my region. Now all I want to do the calibrations, in which I need the Soil Raster and Precipitation Raster.
Can I know to get these raster data along with the Leaf area index and albedo data for the region. My research region is India.
hoseinsalehi
@hoseinsalehi
Hello Everyone
who can Help me to calibrate VIC.4.2.d with MOCOM-UA optimizer?
Thanks.
Saadi4469
@Saadi4469
Good day everybody,
Is it possible to modify the VIC model algorithm, as such that the user could tell the model that these particular grid cells are covered with snow as input?
arivoli
@arivoliii_twitter
Hi everyone, i am trying to run VIC Routing in my PC via ubuntu OS, though i have installed all libraries required it is showing a message "rout" not found. Please help over this.
Jacob
@jschap1
Hi, @arivoliii_twitter, if you are getting the message "rout not found," you may have a path issue. Check that you are calling ./rout from within the directory containing the rout command. Also, make sure you are following the instructions in the VIC documentation (https://vic.readthedocs.io/en/vic.4.2.d/Documentation/Routing/RunRouting/).
@Saadi4469 You can use the state file as an input to the VIC model to specify an initial condition where certain grid cells are snow-covered.
markwestcott34
@markwestcott34
HI all! Does anyone know where I can find a dataset of global input parameters? There is link to https://www.hydro.washington.edu/SurfaceWaterGroup/Data/vic_global_0.5deg.html on https://vic.readthedocs.io/en/master/Datasets/Datasets/#global-datasets but sadly returns a 404..
Jacob
@jschap1
@markwestcott34 Hi, for that particular dataset, you could try getting in touch with the folks at the NASA LIS project. If you're interested, I have a preliminary global 1/16 VIC parameter dataset that I could share with you, but it will need calibration.
SaddamWaheed
@SaddamWaheed
@jschap1 Hi, I also need that dataset. If you have an access or give us a direction how to access them, i would be very thankful for you.
Hi, all The global forcing data of 0.5 resolution cannot be found. Can anyone help me navigating to it please?
Jacob
@jschap1
@SaddamWaheed For the 1/2 degree dataset, you can try contacting someone working on the NASA LIS project (https://lis.gsfc.nasa.gov/contacts). I have a higher resolution dataset that is still in development, but I hope to finish it soon and post it here!
Saadi4469
@Saadi4469
@jschap1 Thanks Jacob, can you please also guide as to how can I convert VIC output files to Geotiff? I want to view the output files in ArcGIS.
Saadi4469
@Saadi4469
I tried using the flux2nc.py script on a Linux shell on windows 10, but I get an error that module NetCDF4 not found in line 29, although I have installed the python package for NetCDF.
Saadi4469
@Saadi4469
Actually I am able to run flux2nc.py script, but I am getting the following error:python2.7 flux2nc.py /mnt/d/Spring_2020/VIC/VIC_Output/ mnt/d/Spring_2020/VIC/VIC_Output/netCDFS
IMPORTANT: /mnt/d/Spring_2020/VIC/VICOutput/ SHOULD CONTAIN ONLY FLUXES FILES!!!
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
Also how can I add a 7th parameter to the script for snow cover?
Jacob
@jschap1

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.

Saadi4469
@Saadi4469
Thank you for your reply @jschap1 so how can I first convert the VIC output files to NetCDF? Is it possible to have an individual NetCDF file for 1 - Precipitation, 2- ET, 3- Runoff, 4- Base flow, 5- SWE and 6- Snow Cover .
Jacob
@jschap1
If you are running VIC in Image Mode, they should already be in NetCDF format. Otherwise, if you are running VIC-4 or VIC-5 in Classic Mode, you will have to convert the outputs. I have written a MATLAB toolbox that reads in VIC outputs in ASCII text format and converts them to geotiffs. https://github.com/jschap1/VICMATLAB
Saadi4469
@Saadi4469
Hi @jschap1 thank you for your quick response, I am using VIC 4.2.d, can you please guide me as how I can use the VIC model in image mode, also I don't have access to complete version of MATLAB.
hoseinsalehi
@hoseinsalehi
Hi
I run VIC.4.2.d version in an area where there are two dams.
Developers have said we can consider them as lakes. Can anyone help me? how to get their data and how to set up a global.param file?
Saadi4469
@Saadi4469
Hello everyone
I am getting some errors in VIC output file
First of all I am getting zero values for snow cover and SWE, why is that? when I change the model to full energy mode, I get MAXIMUM ENERGY ERRORs and the total cumulative error for each grid is quite big and is negative.
Saadi4469
@Saadi4469

YEAR MONTH DAY OUT_SNOW_COVER OUT_SURF_TEMP OUT_RUNOFF OUT_BASEFLOW OUT_SWE OUT_EVAP OUT_PREC

2000 01 01 0.0000 64623.2627 0.4104 1.6648 0.0000 17.1867 9.0000
And this is the error I am getting while running the model in energy balance mode

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

Saadi4469
@Saadi4469
Actually ignore my last two messages, I've got it working now.
Arunav Nanda
@Arunavcoder123
Hello everyone,
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.
Saadi4469
@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.
Saadi4469
@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.
Saadi4469
@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
Saadi4469
@Saadi4469
@hoseinsalehi https://vic.readthedocs.io/en/vic.4.2.d/Documentation/LakeParam/ use this link to get some ideas!
Aydinbakhtar
@Aydinbakhtar
HI
how to execute vic in cygwin?
is there any learning resources for VIC-ASSIST
Saadi4469
@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!
Aydinbakhtar
@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
Aydinbakhtar
@Aydinbakhtar
Hi
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?