Author Topic: [VIC] routing energy balance outputs with subdaily fluxes  (Read 3326 times)

graggio

  • Team tributary
  • *
  • Karma: +0/-0
    • View Profile
  • Institute : CIMA CONICET-UBA
  • Programming language : Python
[VIC] routing energy balance outputs with subdaily fluxes
« on: February 02, 2018, 01:33:01 AM »
Hello!

So far I have had no problems applying the water balance on a local basin using the 4.2.d or 4.0.1 versions of VIC and the Lohmann's routing rout.1.0. However, since most of the basin's water originates from glacier lakes in the Andes, I have been trying to use the energy balance mode.

According to the documentation, energy balance mode should be run with a time step of less than 24h, so I tried using a time step of 8h to run the energy balance mode. Fluxes are created with no problem (VIC 4.2.d), but for what I understand Lohmann's routing is not made to understand sub-daily fluxes (from the documentation and from the script make_convolution.f used in the compilation of the routing model). Trying to do so raises the following error:

Code: [Select]
Error reading time-series data insufficient data or missing input file
I made sure that the fluxes are located where the rout control file reads them.


So how can I rout fluxes with, for example, an 8h time step?

Thanks in advance!


p.s. Even though it's against what the documentation say, someone told me you can run energy balance mode with 24h time step (VIC 4.2.d). I tried this but it resulted in an error a particular cell:


Code: [Select]
ERROR: snow_melt failed to converge to a solution in root_brent.  Variable values will be dumped to the screen, check for invalid values.

Diwan

  • Moderator
  • Team sea
  • *****
  • Karma: +12/-0
  • Be a learner or Teacher or helping to learn
    • View Profile
  • Institute : National Institute of Technology Karnataka, Surathkal
  • Programming language : Python, R, Matlab
Re: [VIC] routing energy balance outputs with subdaily fluxes
« Reply #1 on: February 09, 2018, 02:56:02 PM »
Hello!

Quote
p.s. Even though it's against what the documentation say, someone told me you can run energy balance mode with 24h time step (VIC 4.2.d)

Yes. You can run energy balance computation with daily time step in VIC earlier version. But it will take more computational time depending on your grids.

Though my response to your query is late.. can you share some more information? eg. global parameter, soil and other parameters. 
Research Scholar
Dept. of Applied Mechanics & Hydraulics
National Institute of Technology Karnataka
Surathkal Mangalore - 575025

Muhammad Riaz

  • Team tributary
  • *
  • Karma: +0/-0
    • View Profile
  • Institute : NARC
  • Programming language : R
Re: [VIC] routing energy balance outputs with subdaily fluxes
« Reply #2 on: February 14, 2018, 03:32:17 PM »
Hi Diwan
I want to calibrate the VIC and route model
please suggest the algorithm for automatic calibration. and please also tell the procedure for automatic calibration

Diwan

  • Moderator
  • Team sea
  • *****
  • Karma: +12/-0
  • Be a learner or Teacher or helping to learn
    • View Profile
  • Institute : National Institute of Technology Karnataka, Surathkal
  • Programming language : Python, R, Matlab
Re: [VIC] routing energy balance outputs with subdaily fluxes
« Reply #3 on: February 14, 2018, 04:15:59 PM »
Hi Riaz

I am not sure about auto calibration. you May refer @Ila Chawla 's recent paper in AWR http://wateradroit.com/forum/index.php/topic,1344.0.html  and also see the MOCOM-UA method from VIC website.. http://www.hydro.washington.edu/Lettenmaier/Models/VIC/Documentation/MOCOM.shtml

Hope this will help.

Note: Kindly post your question with relevant topics. If you are not find any existing topic relevant to your query you can start a new topic. You can use this thread for general queries and other information related to VIC http://wateradroit.com/forum/index.php/topic,205.75.html 
« Last Edit: February 14, 2018, 04:35:13 PM by Diwan »
Research Scholar
Dept. of Applied Mechanics & Hydraulics
National Institute of Technology Karnataka
Surathkal Mangalore - 575025

Muhammad Riaz

  • Team tributary
  • *
  • Karma: +0/-0
    • View Profile
  • Institute : NARC
  • Programming language : R
Re: [VIC] routing energy balance outputs with subdaily fluxes
« Reply #4 on: February 17, 2018, 09:37:45 PM »
Hi Diwan
I run the model and try the I change the b_inf lower and upper limit
Please tell how I calibrate the VIC . Orange line is for simulated flow by model

Diwan

  • Moderator
  • Team sea
  • *****
  • Karma: +12/-0
  • Be a learner or Teacher or helping to learn
    • View Profile
  • Institute : National Institute of Technology Karnataka, Surathkal
  • Programming language : Python, R, Matlab
Re: [VIC] routing energy balance outputs with subdaily fluxes
« Reply #5 on: February 20, 2018, 10:10:18 AM »
Hi Riaz

Ds, Dsmax, Ws, binf and soil layer depths are main calibration parameter for the model. (other variables can be calibrated eg., snow partitioning, energy balance etc.)

some points for understanding the parameter and typical range
 
1) Higher value of Ds,[Fraction of Dsmax(0-~30, in mm/day)] the baseflow will be higher at lower water
content in lowest soil layer ( 0 - 1, units in percentage)

2) Higher value of Ws will raise the water content required for rapidly
increasing, non-linear baseflow, which will tend to delay runoff
peaks (0 - 1, units in percentage)

3) Higher value of binf gives lower infiltration and yields higher
surface runoff (0 - 0.4)

4) Soil depth effects many model variables (0.1 to 1.5, units in meter)

5) For runoff considerations, thicker soil depths slow down (baseflow
dominated) seasonal peak flows and increase the loss due to
evapotranspiration

You may fix the value by trail and error based on Ksat for the particular grids and knowledge on the local condition

HTH

Note: Kindly post your question with relevant topics. If you are not find any existing topic relevant to your query you can start a new topic. You can use this thread for general queries and other information related to VIC http://wateradroit.com/forum/index.php/topic,205.75.html 
« Last Edit: February 20, 2018, 10:23:42 AM by Diwan »
Research Scholar
Dept. of Applied Mechanics & Hydraulics
National Institute of Technology Karnataka
Surathkal Mangalore - 575025

Muhammad Riaz

  • Team tributary
  • *
  • Karma: +0/-0
    • View Profile
  • Institute : NARC
  • Programming language : R
Re: [VIC] routing energy balance outputs with subdaily fluxes
« Reply #6 on: February 20, 2018, 08:22:44 PM »
Thanks Diwan
I tried every parameter and use the lower and higher limits of every parameter. But the peak does not low . How we can reduce the peak