Where communities thrive


  • Join over 1.5M+ people
  • Join over 100K+ communities
  • Free without limits
  • Create your own community
People
Repo info
Activity
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?
JACK LONG
@longbiao107995_twitter
the 1/16 VIC parameter dataset is really important.
Saadi4469
@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?
Saadi4469
@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

Saadi4469
@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.)