[Dhsvm-users] Re: Dhsvm-users Digest, Vol 25, Issue 17
Nathalie Voisin
nathalie at hydro.washington.edu
Mon Oct 30 13:36:39 PST 2006
Mark,
I will try to give brief answers.
(1) I believe that the tutorial provides you with the plot_modelmap.scr
script (GMT scripts) . If you
read it, you will see that it uses the myconvert to go back to ascii
format, then adds the previously saved header and there you are, you have
your arcGIS format for plotting. ( Just ignore the step with the GMT
command). I will put the script on the website as well, under the TOOLS
section.
(2) in the configuration file, you specify which statefile to use. If you
change your period run of the input file, you should change this as well.
The tutorial might not have specified this step.
(3.b)
- Units per cubic meter per time step, i.e. cubic meter per 3 hours if 3
hour is your time step.
- negative numbers. correct me if I am wrong but it looks like you have no
water ( precipitation) and this is a start of a run. Those negative
numbers are very small and could be round off errors. DO you still have
negative number at the end of your run?
(3.c) Stream.Flow. The documentation on the website at
http://www.hydro.washington.edu/Lettenmaier/Models/DHSVM/output_network.shtml
explains that "if the line is for a single segment the remaining fields
are as follows:
6 real ( columnn #6)
Change in segment storage, m, occurring within the previous routing time
step."
meaning that this can be negative during a recession period.
"totals" if for the entire network ( all stream segments).
(3.d to 3.f) you are correct. See earlier remarks on the units.
Nathalie
/******************************************************************\
* Nathalie Voisin *
* Hydrology Group *
* Civil and Environmental Engineering *
* University of Washington *
* Graduate Research Assistant *
* p:206-685-1796 *
\******************************************************************/
On Sat, 28 Oct 2006 dhsvm-users-request at mailman1.u.washington.edu wrote:
> Send Dhsvm-users mailing list submissions to
> dhsvm-users at u.washington.edu
>
> To subscribe or unsubscribe via the World Wide Web, visit
> http://mailman1.u.washington.edu/mailman/listinfo/dhsvm-users
> or, via email, send a message with subject or body 'help' to
> dhsvm-users-request at mailman1.u.washington.edu
>
> You can reach the person managing the list at
> dhsvm-users-owner at mailman1.u.washington.edu
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of Dhsvm-users digest..."
>
>
> Today's Topics:
>
> 1. Please help orient me in these output files. (meisenbi at vt.edu)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Fri, 27 Oct 2006 16:07:47 -0400
> From: meisenbi at vt.edu
> Subject: [Dhsvm-users] Please help orient me in these output files.
> To: dhsvm-users at u.washington.edu
> Message-ID: <1161979667.45426713578ac at webmail.vt.edu>
> Content-Type: text/plain; charset=ISO-8859-1
>
> Perhaps I am just tired today, but I haven't been able to figure a couple of
> things out. I also don't know whether to post this as one email or several.
>
>
>
> (1) On the second run of the tutorial, the model start and end are 1994 and
> 1996, but the model state date is 1995. For what purpose?
>
> Hopefully I was correct in entering a value of "4" for the number of map
> variables in the MODEL MAPS section based on the tutorial. Also, I assume I
> can choose to have the same map variable output for several dates.
>
> Last, how do I convert these back to ASCII so I can import them into ARCGIS? Is
> there a reverse program to mycovert.c that I missed?
>
>
> (2) When we copy the *1994* state files to the modelstate/ directory, How does
> DHSVM tell which state file to use? I didn't overwrite the original state
> files by that set of 5 commands.
>
>
> (3) I have been reading some of my model output dat files.
>
> (3a) Aggregated.values seems to be a straight forward table. The variable and
> units are described on the output variable ID table at
>
> http://www.hydro.washington.edu/Lettenmaier/Models/DHSVM/output_varid.shtml
>
> (3b) Mass.Balance. These variables are described at
>
> http://www.hydro.washington.edu/Lettenmaier/Models/DHSVM/output_default.shtml
>
> These units are cubic meters per time step?
>
> What is the meaning of negative numbers?
>
> -----------------------------------------------
> 12/02/1992-09:00:00 0.0362209 0.0277253 0.192406 0.163873 -6.12781e-05
> 6.07176e-05 5.65189e-07 3.33302e-07 1.27417e-06 0 -6.51353e-06 -1.38137e-05
> 0.420165 0 4.68232e-09 2.41709e-07
> 12/02/1992-12:00:00 0.0362151 0.0271805 0.192351 0.163877 -6.06088e-05
> 6.00584e-05 5.5501e-07 3.26566e-07 0.000523764 0 -5.38485e-07 -1.66305e-05
> 0.419563 0 4.60556e-09 -9.42536e-08
> 12/02/1992-15:00:00 0.0362306 0.0268704 0.192275 0.163873 -5.99511e-05
> 5.94106e-05 5.4503e-07 3.19953e-07 1.12927e-22 0 -4.20345e-06 -0.000310081
> 0.419189 0 4.53005e-09 1.82519e-07
> 12/02/1992-18:00:00 0.0362272 0.0265855 0.192219 0.163867 -5.93048e-05
> 5.8774e-05 5.35244e-07 3.13459e-07 9.50061e-23 0 -6.26889e-06 -0.000284931
> 0.418839 0 4.45577e-09 2.05322e-08
> 12/02/1992-21:00:00 0.0362205 0.0263901 0.192167 0.163858 -5.86698e-05
> 5.81485e-05 5.25652e-07 3.07084e-07 6.56981e-23 0 -8.36691e-06 -0.000195309
> 0.418577 0 4.38269e-09 -1.34228e-07
> -----------------------------------------------
>
> (3c) Stream.Flow
>
> The format of this file doesn't seem to match the description in the
> documentation. The second column is the total stream flow, but I thought where
> we placed the stream save indicator was the same location as the main outlet.
> So what are these columns, and the 0,1 rows. What are the negative values?
> Are the units cubic meters per time step?
>
> --------------------------------------------------
> 10/01/1990-03:00:00 0 84961 3.5516e+07 0
> 0 3.5431e+07 "Totals"
> 10/01/1990-06:00:00 1 2.8809e+06 812.38 3.2518e+06
> -3.7008e+05 "RAINY CREEK"
> 10/01/1990-06:00:00 0 2.3307e+05 3.2518e+06 1.5836e+05
> -3.0188e+06 0.96875 "Totals"
> 10/01/1990-09:00:00 1 5.3944e+05 3027.8 5.6868e+05
> -26206 "RAINY CREEK"
> 10/01/1990-09:00:00 0 4.5391e+05 5.6868e+05 43591
> -1.1476e+05 0.14062 "Totals"
> 10/01/1990-12:00:00 1 4.4188e+05 2585.5 4.4556e+05
> -1097.9 "RAINY CREEK"
> 10/01/1990-12:00:00 0 4.4164e+05 4.4556e+05 39674
> -3917 0.083984 "Totals"
> 10/01/1990-15:00:00 1 3.4923e+05 4445.7 3.547e+05
> -1017 "RAINY CREEK"
> 10/01/1990-15:00:00 0 3.4607e+05 3.547e+05 31050
> -8623.6 0.054688 "Totals"
> 10/01/1990-18:00:00 1 2.7998e+05 3140.5 2.8391e+05
> -790.42 "RAINY CREEK"
> 10/01/1990-18:00:00 0 2.7736e+05 2.8391e+05 24494
> -6555.9 0.087891 "Totals"
> 10/01/1990-21:00:00 1 2.3293e+05 4463.6 2.379e+05
> -512.31 "RAINY CREEK"
> 10/01/1990-21:00:00 0 2.3369e+05 2.379e+05 20280
> -4214.1 -0.023926 "Totals"
> 10/02/1990-00:00:00 1 2.341e+05 3411.1 2.3751e+05
> 1.4055 "RAINY CREEK"
> 10/02/1990-00:00:00 0 2.3754e+05 2.3751e+05 20302
> 22.298 0.0011063 "Totals"
> -------------------------------------------------------
>
> (3d) Stream.Only
>
> According to the documentation, this file contains the discarge for places where
> my SAVE indicators occured. However, these values appear to correspond to the
> second column, row "0" of Stream.Flow, not the one that was labled "Rainy
> Creek". Also, are the units cubic meters per time step?
>
> 10/01/1990-03:00:00
> 10/01/1990-06:00:00 3.2518e+06
> 10/01/1990-09:00:00 5.6868e+05
> 10/01/1990-12:00:00 4.4556e+05
> 10/01/1990-15:00:00 3.547e+05
> 10/01/1990-18:00:00 2.8391e+05
> 10/01/1990-21:00:00 2.379e+05
> 10/02/1990-00:00:00 2.3751e+05
>
> (3e) Road.Flow.
>
> According to the documentation, these columns refer to each culvert with a save
> marker. Nothing in my road.network.dat had a save indicator, so I am unsure
> what these columns refer to.
>
> 10/01/1990-03:00:00 0 86090 84865 0
> 0 -1225.5 "Totals"
> 10/01/1990-06:00:00 0 96293 96192 1326.5
> 101.01 0.0051804 "Totals"
> 10/01/1990-09:00:00 0 88086 88289 1123.1
> -203.42 -0.026489 "Totals"
> 10/01/1990-12:00:00 0 75111 75278 956.79
> -166.3 -0.0010376 "Totals"
> 10/01/1990-15:00:00 0 59079 59295 740.58
> -216.21 -0.0023651 "Totals"
> 10/01/1990-18:00:00 0 50495 50578 658.09
> -82.491 0.0049438 "Totals"
> 10/01/1990-21:00:00 0 41307 41455 510
> -148.09 0.004364 "Totals"
> 10/02/1990-00:00:00 0 41937 41937 509.71
> -0.29639 0.00048476 "Totals"
> 10/02/1990-03:00:00 0 42215 42221 503.69
> -6.0197 -0.0040379 "Totals"
> 10/02/1990-06:00:00 0 35521 35597 427.72
> -75.967 0.0019989 "Totals"
> 10/02/1990-09:00:00 0 45265 45157 535.01
> 107.29 -0.0036469 "Totals"
> 10/02/1990-12:00:00 0 52757 52648 644.09
> 109.09 0.0021439 "Totals"
> 10/02/1990-15:00:00 0 58367 58288 723.13
> 79.033 -0.0059052 "Totals"
> 10/02/1990-18:00:00 0 60199 60177 745.58
> 22.455 -0.006115 "Totals"
> 10/02/1990-21:00:00 0 63655 63614 787.2
> 41.617 -0.0043449 "Totals"
>
> (3f) Roadflow.Only
>
> Only has a column with a date, time step. Everything else is empty (I assume
> this is because I have no SAVE indicators).
>
>
>
> Sorry for the long post.
>
> Thanks for everyone's help.
>
>
>
>
> --
> Mark H. Eisenbies
> Postdoctoral Associate
> Virginia Tech Dept. of Forestry
> 228 Cheatham Hall
> Blacksburg, VA 24060
>
>
> ------------------------------
>
> _______________________________________________
> Dhsvm-users mailing list
> Dhsvm-users at mailman1.u.washington.edu
> http://mailman1.u.washington.edu/mailman/listinfo/dhsvm-users
>
>
> End of Dhsvm-users Digest, Vol 25, Issue 17
> *******************************************
>
More information about the Dhsvm-users
mailing list