Where communities thrive


  • Join over 1.5M+ people
  • Join over 100K+ communities
  • Free without limits
  • Create your own community
People
Repo info
Activity
Saadi4469
@Saadi4469
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.)
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

Kiran KG
@Kiran-KG

Hi all, I complied the image version of VIC with 'ROUT=rout_rvic' option in the make file. Now, when I run an example file (Stehekin) I get the following error.
"Traceback (most recent call last):
/lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf3) [0x7f61f27330b3]
./vic_image.exe(+0x4dc48) [0x55940e9a1c48]
./vic_image.exe(+0xb6753) [0x55940ea0a753]
[ERROR] ../../extensions/rout_rvic/src/rout_alloc.c:48: errno: No such file or directory:
Error opening /home/vvs/Downloads/VIC_sample_data-master/image/Stehekin/parameters/parameters_01.rvic.prm.Stehekin.20150727.nc "

When I compile it without the rvic option it runs without a problem. Does anyone know what the issue could be and how to resolve it?. Thanks in advance