[Chris G.]

Return-Path: Received: from compute1.internal (compute1.nyi.internal [10.202.2.41])
by sloti25t09 (Cyrus fastmail-fmjessie44472-15312-git-fastmail-15312) with LMTPA;
Wed, 12 Jul 2017 11:38:15 -0400
X-Cyrus-Session-Id: sloti25t09-1376490-1499873895-2-4042420544672845074
X-Sieve: CMU Sieve 3.0
X-Spam-known-sender: no
X-Spam-score: 0.0
X-Spam-hits: BAYES_00 -1.9, FREEMAIL_FROM 0.001, HTML_MESSAGE 0.001,
RCVD_IN_DNSWL_NONE -0.0001, RCVD_IN_MSPIKE_H3 -0.01,
RCVD_IN_MSPIKE_WL -0.01, SPF_PASS -0.001, URI_HEX 1.122, LANGUAGES en,
BAYES_USED global, SA_VERSION 3.4.0
X-Spam-source: IP=’209.85.217.170′, Host=’mail-ua0-f170.google.com’, Country=’US’,
FromHeader=’com’, MailFrom=’com’
X-Spam-charsets: to=’UTF-8′, plain=’UTF-8′, html=’UTF-8′
X-Resolved-to: [email protected]
X-Delivered-to: [email protected]
X-Mail-from: [email protected]
Received: from mx4 ([10.202.2.203])
by compute1.internal (LMTPProxy); Wed, 12 Jul 2017 11:38:15 -0400
Received: from mx4.messagingengine.com (localhost [127.0.0.1])
by mailmx.nyi.internal (Postfix) with ESMTP id 1400AC8547
for ; Wed, 12 Jul 2017 11:38:15 -0400 (EDT)
Received: from mx4.messagingengine.com (localhost [127.0.0.1])
by mx4.messagingengine.com (Authentication Milter) with ESMTP
id E67116435F4;
Wed, 12 Jul 2017 11:38:15 -0400
Authentication-Results: mx4.messagingengine.com;
dkim=pass (2048-bit rsa key) header.d=gmail.com [email protected] header.b=ppujFgKq;
dmarc=pass header.from=gmail.com;
spf=pass [email protected] smtp.helo=mail-ua0-f170.google.com;
x-google-dkim=pass (2048-bit rsa key) header.d=1e100.net [email protected] header.b=G6KFJayR
Received-SPF: pass
(gmail.com … _spf.google.com: Sender is authorized to use ‘[email protected]’ in ‘mfrom’ identity (mechanism ‘include:_netblocks.google.com’ matched))
receiver=mx4.messagingengine.com;
identity=mailfrom;
envelope-from=”[email protected]”;
helo=mail-ua0-f170.google.com;
client-ip=209.85.217.170
Received: from mail-ua0-f170.google.com (mail-ua0-f170.google.com [209.85.217.170])
(using TLSv1.2 with cipher ECDHE-RSA-AES128-GCM-SHA256 (128/128 bits))
(No client certificate requested)
by mx4.messagingengine.com (Postfix) with ESMTPS
for ; Wed, 12 Jul 2017 11:38:14 -0400 (EDT)
Received: by mail-ua0-f170.google.com with SMTP id g40so16774256uaa.3
for ; Wed, 12 Jul 2017 08:38:14 -0700 (PDT)
DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed;
d=gmail.com; s=20161025;
h=mime-version:in-reply-to:references:from:date:message-id:subject:to;
bh=pIiv7fXzq4peO+8BNKlGQ2pUpqtGBns7TgyuiWAHF5M=;
b=ppujFgKqw4vdQnaTDkc3WwretQ1SrRmgTCE1VH7mVkuBZWZoVGJe0N2UFJ2aXLJ14o
P10ApBBq5Byc+nmo0ZzVOz03L6miCmfJqfbDPU9yRnHYIy9xfPv7NeeyQX4b2vb+opzv
cItQdFqrW/LsrGOhAodXXhvULjZnfBSckYHdyKa2GT7otWeN6Mt6gPYgr293KeMkdFTD
fztlY9uFYxCg3tXZsa3gYh5gImFTW2qzAzaoZkeLFStKEHqUErsniqrRST7XYlrZVfma
rtynie0wjHcU5Xvcb1HfJGYkAwNzay8d23lwJqzmfvNa0Z02rh8KvbA4MYLqak1Gi+rN
P/Cw==
X-Google-DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed;
d=1e100.net; s=20161025;
h=x-gm-message-state:mime-version:in-reply-to:references:from:date
:message-id:subject:to;
bh=pIiv7fXzq4peO+8BNKlGQ2pUpqtGBns7TgyuiWAHF5M=;
b=G6KFJayRW8ehNfeW3x5ZyopF3Wn06a0WgD6aDl1vYdZVTuIVVCiymfoUUl47KHzoVp
8bnU/QvoinFBDzGreLLzNd9ozA+xcYVB2s/PLZeKyswc1uGumqMs+P5at6z6r0vvhdj2
kLItPc80CgWpTQv0X0rs1FJbfi7+lGEtLkLwjfaBKGHKHSPQZguZ7nT64rW3CLpD+lFJ
ECoc3BU/3cJGaBVNnlFj3MQRyb8NgWsEBess2sB47rHoJucugyp9b2OPT/gOeye2TLGc
jq8M+I3s5UCN2M5teFQI3KMhgnKoAc80eeAhx7c6vnFsjF7A4tlhTtPKFVNhHx/M1aTq
qZfA==
X-Gm-Message-State: AIVw112KyMs/o/eu3aYM+Enyk/hEyIUvCOmSgNVnB0pRqrW5IgGgE+cx
ycdD0lFbtIx/6wbkDJa/5lUeRGvJQd3u
X-Received: by 10.159.36.215 with SMTP id 81mr3685508uar.48.1499873893189;
Wed, 12 Jul 2017 08:38:13 -0700 (PDT)
MIME-Version: 1.0
Received: by 10.176.80.162 with HTTP; Wed, 12 Jul 2017 08:38:12 -0700 (PDT)
In-Reply-To: <[email protected]>
References: <[email protected]> <[email protected]>
<[email protected]>
From: Chris Goodell Date: Wed, 12 Jul 2017 08:38:12 -0700
Message-ID:
Subject: Re: Output Variables for Water Quality in VBA
To: =?UTF-8?B?Sm9zw6kgRMOtYXogW3ZpYSBIRUMtUkFTIEhlbHBd?=
Content-Type: multipart/alternative; boundary=”001a113e15c89c8ec90554209a46″

–001a113e15c89c8ec90554209a46
Content-Type: text/plain; charset=”UTF-8″
Content-Transfer-Encoding: quoted-printable

Jose-

Sorry, but HEC has not included water quality output extraction in the
HECRASController for Version 5.0.3.

Chris

On Tue, Jul 11, 2017 at 8:02 AM, Jos=C3=A9 D=C3=ADaz [via HEC-RAS Help] < [email protected]> wrote:

> Hello Chris,
>
> I am in the same question as Horace, is the 5.0.3 HecRASController now
> able to handle water quality data?
>
> ——————————
> If you reply to this email, your message will be added to the discussion
> below:
> http://hec-ras-help.1091112.n5.nabble.com/Output-
> Variables-for-Water-Quality-in-VBA-tp2913p4189.html
> To start a new topic under HECRAS Controller, email
> [email protected]
> To unsubscribe from HECRAS Controller, click here
> <http://hec-ras-help.1091112.n5.nabble.com/template/NamlServlet.jtp?macro=
=3Dunsubscribe_by_code&node=3D1184&code=3DdGhlcmFzc29sdXRpb25AZ21haWwuY29tf=
DExODR8MTE3NDg5NDgxNQ=3D=3D>
> .
> NAML
> <http://hec-ras-help.1091112.n5.nabble.com/template/NamlServlet.jtp?macro=
=3Dmacro_viewer&id=3Dinstant_html%21nabble%3Aemail.naml&base=3Dnabble.naml.=
namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.n=
aml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabb=
le.view.web.template.NodeNamespace&breadcrumbs=3Dnotify_subscribers%21nabbl=
e%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21na=
bble%3Aemail.naml>
>

–001a113e15c89c8ec90554209a46
Content-Type: text/html; charset=”UTF-8″
Content-Transfer-Encoding: quoted-printable

Jose-

Sorry, but HEC has not included w=
ater quality output extraction in the HECRASController for Version 5.0.3. =
=C2=A0

Chris

On Tue, Jul 11, 2017 at 8:02 AM, Jos=C3=A9 =
D=C3=ADaz [via HEC-RAS Help] <[email protected]=
e.com> wrote:

argin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

Hello Chris,

I am in =C2=A0the same question as Horace, is the 5.0.3 HecRASContr=
oller now able to handle water quality data?

=09
=09
=09

---

If you reply to this email, your message =
will be added to the discussion below:

http://hec=
-ras-help.1091112.n5.nabble.com/Output-Variables-for-Water-Qualit=
y-in-VBA-tp2913p4189.html

To start a new topic under HECRAS Controller, email ml+s1091112n1184h88@n5.=
nabble.com

To unsubscribe from HECRAS Controller, click here.
NAML

–001a113e15c89c8ec90554209a46–

[Chris G.]

Thanks for sharing that Jennifer. It’s common (and likely) to see differences from 1D to 2D, if for nothing else, because the 2D equations account for more of the losses that we tend to wrap into the n value for 1D. Assuming both 1D and 2D models are stable, I would tend to believe the 2D results more. But there are other things to consider as well. Cell size can play an important part in this. If your cells are too large across a conveyance path, this can cause a smoothing of the velocity profile, resulting in lower discrete velocities and higher water surface elevations. I generally target 5 to 7 cells (minimum) across an important conveyance path (e.g. across a main channel, or even just a concentration of flow in an overbank area). Of course, stability is the first thing to work out. Make sure you are targeting a Courant number of around 1 as much as possible in your 2D area. This will usually take care of most of the stability issues, including the pulsing of flow you mentioned.

[Chris G.]

Return-Path: X-Original-To: [email protected]
Delivered-To: [email protected]
Received: from mail-pg0-f47.google.com (mail-pg0-f47.google.com [74.125.83.47])
by mwork.nabble.com (Postfix) with ESMTP id DE2E832255673
for ; Fri, 10 Mar 2017 20:24:28 -0700 (MST)
Received: by mail-pg0-f47.google.com with SMTP id b129so45461355pgc.2
for ; Fri, 10 Mar 2017 19:24:28 -0800 (PST)
DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed;
d=gmail.com; s=20161025;
h=from:content-transfer-encoding:mime-version:date:subject:message-id
:references:in-reply-to:to;
bh=e34ZJabSrSQpcOrJfiTutJ+P417K6j+XtEN27qnIt+k=;
b=fQCBIfKCm4fzw9extOtOVnqeNp1RystPsu/IX1ljuWJQSRrw4EzySLx6xuBH06VL11
YmVfVocDxAnqKQM2EV+8cwPnwX6ji4biaelFX/88gX1JDf8HNpeoPlC0LfTZFHzjC/CO
1lk9C3lpLju4gI+F+43GB8cEWXpa+e+OKIeVtyfkrVBnArjmbs/6f5KpyTLAa2V7BhK6
tro/YoA7wSrsWNpTO8fwdgP51v7DnvRbOQ/N41OnXK3JY+/F8a+1buOHR3MZ/Orak/rq
vD0LmsZ1BmMUX+2+sFIo4pMbLnw/KLNN2rlkxRm7s1E5TKwRlV8oNrrsaa8m+cuMwvJs
c+4w==
X-Google-DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed;
d=1e100.net; s=20161025;
h=x-gm-message-state:from:content-transfer-encoding:mime-version:date
:subject:message-id:references:in-reply-to:to;
bh=e34ZJabSrSQpcOrJfiTutJ+P417K6j+XtEN27qnIt+k=;
b=nRtIaNOjU8ENyQKfAYwLeb0VzXi3e9869PqYA0FNXWPAEdBiSPk/P4Q6wpG/jmlsoy
olcG/NV7vJWIhYg1ooaY97GkmOfn6SDA5VXZKPAJFYI35lZnYflluNMvBfWbnMRQrkSb
JVtj3seI6AopWaZ8DNvd/NoY1sZFpS/EkytTQfglsAFir4mSFy+vgA9uEAaxGfLSm4lR
qPd/+Ds67dPvBoaB12cA21x1JterNlDFUXq/hakYq5QtP7qCN8/Jy2xrmTlB24E49AfY
NNdCWiuvGqunIAIECovP0SRzMRc0SiREKYUZSdhsqk1pHLnZpG5l972FL3OMFr195tlr
EJ4g==
X-Gm-Message-State: AMke39lS6W/xFULtxBRkW7BI3XrvRkTOcWTfZ2iMzaepLYitEAY7uzCPiV2Z2HD/nLJQhA==
X-Received: by 10.99.125.68 with SMTP id m4mr24541538pgn.13.1489202661879;
Fri, 10 Mar 2017 19:24:21 -0800 (PST)
Received: from [192.168.1.5] (c-73-37-84-115.hsd1.or.comcast.net. [73.37.84.115])
by smtp.gmail.com with ESMTPSA id v9sm20980350pfg.133.2017.03.10.19.24.21
for
(version=TLS1 cipher=ECDHE-RSA-AES128-SHA bits=128/128);
Fri, 10 Mar 2017 19:24:21 -0800 (PST)
From: The RAS Solution Content-Type: multipart/alternative;
boundary=Apple-Mail-F9A7553B-4BEC-4C2B-BD85-465486B090A0
Content-Transfer-Encoding: 7bit
Mime-Version: 1.0 (1.0)
Date: Fri, 10 Mar 2017 19:24:20 -0800
Subject: Re: save new plan through controller
Message-Id: <[email protected]>
References: <[email protected]> <[email protected]> <[email protected]>
In-Reply-To: <[email protected]>
To: “Bharath [via HEC-RAS Help]”
X-Mailer: iPhone Mail (14A456)

–Apple-Mail-F9A7553B-4BEC-4C2B-BD85-465486B090A0
Content-Type: text/plain;
charset=us-ascii
Content-Transfer-Encoding: quoted-printable

Once you bring in the newly created Geometry file, delete the associated HDF=
file. Then save the geometry in RAS. That will force RAS to rewrite the HDF=
file for the new geometry.=20

Sent from my iPhone

> On Mar 10, 2017, at 1:52 PM, Bharath [via HEC-RAS Help] wrote:
>=20
> Once again, Thanks a lot Chris. Does doing so affect the hdf file associat=
ed with the geometry file? Like am i supposed to create a new hdf file for a=
new geometry file?=20
>=20
>=20
> If you reply to this email, your message will be added to the discussion b=
elow:
> http://hec-ras-help.1091112.n5.nabble.com/save-new-plan-through-controller=
-tp3736p3739.html
> To start a new topic under HECRAS Controller, email ml-node+s1091112n1184h=
[email protected]=20
> To unsubscribe from HECRAS Controller, click here.
> NAML

–Apple-Mail-F9A7553B-4BEC-4C2B-BD85-465486B090A0
Content-Type: text/html;
charset=utf-8
Content-Transfer-Encoding: 7bit

Once you bring in the newly created Geometry file, delete the associated HDF file. Then save the geometry in RAS. That will force RAS to rewrite the HDF file for the new geometry. 

Sent from my iPhone

On Mar 10, 2017, at 1:52 PM, Bharath [via HEC-RAS Help] <[email protected]> wrote:

Once again, Thanks a lot Chris. Does doing so affect the hdf file associated with the geometry file? Like am i supposed to create a new hdf file for a new geometry file?

---

If you reply to this email, your message will be added to the discussion below:

http://hec-ras-help.1091112.n5.nabble.com/save-new-plan-through-controller-tp3736p3739.html

To start a new topic under HECRAS Controller, email [email protected]
To unsubscribe from HECRAS Controller, click here.
NAML

–Apple-Mail-F9A7553B-4BEC-4C2B-BD85-465486B090A0–

[Chris G.]

Yes, just give it a each new file a new name. If you want to stay consistent with RAS convention, keep the project name the same, just change the file extension (e.g. Project.g01, Project.g02, etc). You may have to rewrite the project and plan files accordingly so that RAS recognizes the new files.

[Chris G.]

I don’t know MatLab well enough to say what the problem is.

[Chris G.]

This from Arturo Leon: “The error says: No method ‘Compute_CurrentPlan’ with matching signature found for class ‘COM.RAS503_HECRASCONTROLLER’. It appears that the Matlab class was not defined. In our paper I defined the class as “h”, so to perform any operation the “h. XXXXXX” needs to be used, including defining each instance of HEC-RAS. For using HEC-RAS in parallel computing or to do multiple computations in series, h is replaced by a vector h(i) where “i” is the ID of the simulation (e.g., 1 to 100). For instance, to open a single instance of HEC-RAS 5.0, h=actxserver(‘RAS500.HECRASCONTROLLER’) is used in the paper. I would suggest the student to check our paper available at: http://www2.egr.uh.edu/~aleon3/papers_PDF/Journal/HEC_RAS_Controller.pdf&#8221;

[Chris G.]

As mentioned above, the arguments are:

‘Compute the current plan.
Dim lngMessages As Long
Dim strMessages() As String
Dim blnDidItCompute As Boolean

blnDidItCompute = RC.Compute_CurrentPlan(lngMessages, strMessages(), True)

With the 3rd argument set to “True”, the code will pause while HEC-RAS is computing.

Let me know if this works. I’ll check around to see if there is a MatLab compatibility issue with this command.

[Chris G.]

Bharath-

There are three arguments to include with the Compute_CurrentPlan Function.

1. nmsg (long). the number of returned messages
2. msg() (string). messages returned from the HECRASController during computations.
3. BlockingMode (boolean). True or False. If blocking mode is set to false, then code will continue to be readh while HEC-RAS is computing. So you want this to be set to True. You omitted the BlockingMode argument from your Compute_CurrentPlan call, so by default it is set to false. Your call should look more like this (in VBA syntax):

‘Compute the current plan.
Dim lngMessages As Long
Dim strMessages() As String
Dim blnDidItCompute As Boolean
blnDidItCompute = RC.Compute_CurrentPlan(lngMessages, strMessages(), True)

With the 3rd argument set to “True”, the code will pause while HEC-RAS is computing.

[Chris G.]

In version 5 and later, you need to include the QuitRAS subroutine at the end of your procedure. See example below…

Sub QuitRAS()

‘**********************************************************
‘Demonstrates the QuitRAS subroutine.

‘Written by Christopher Goodell
‘May 7, 2013

‘Closes HEC-RAS. _
Every procedure that instantiates an HECRASController _
class must conclude with a call to the QuitRAS _
subroutine.

‘Tested with the BEAVCREK.prj data set. Can be _
used with any HEC-RAS data set.
‘**********************************************************

‘Open a new HEC-RAS Project
Dim RC As New RAS503.HECRASController
mSamples.OpenRASProjectByRef RC

‘Show HEC-RAS
RC.ShowRAS

‘Message the user to click to close HEC-RAS
MsgBox “Click OK to close HEC-RAS.”

‘Close HEC-RAS
RC.QuitRAS ‘This is what you have to add to close HEC-RAS in Version 5 or later.

End Sub

Sub OpenRASProjectByRef(ByRef RC As RAS503.HECRASController)

‘**********************************************************
‘Demonstrates Project_Open subroutine and how to pass the _
HECRASController as a reference.

‘Written by Christopher Goodell
‘July 3, 2012

‘Opens a RAS project by being called from another _
subroutine
‘**********************************************************

‘Open the HEC-RAS Project
Dim strFilename As String
Sheets(“RASProjects”).Select
strFilename = Range(“C4”).Value
RC.Project_Open strFilename

End Sub

[Chris G.]

Lauren-

Not directly. However, you could “build” the cross section incrementally in the HD/Uniform Flow calculator in RAS (Run…Hydraulic Design Functions, Type…Uniform Flow). With each incremental addition of a station elevation point, run the uniform flow computations to solve for discharge. Save the station and the resulting discharge in an excel table. As you “build” the cross section you should have a cumulative discharge going across the cross section. Then don’t forget to divide all of those discharges by the square root of the energy slope to get conveyance. I haven’t tried this myself, but in principle it should work. If you already have the cross section in RAS, it would be easier to “deconstruct” it as you collect normal discharges, rather than build it (i.e. compute normal discharge, then remove the last sta-el pair from the table then repeat).

Chris

[Chris G.]

Yes, it’s for 5.0. Keep in mind that in 5.0, you now have to use the “QuitRAS” procedure to close HEC-RAS from your Task Manager at the end of your procedures. Version 4.1 did not require QuitRAS.

[Chris G.]

Horace-

Thanks for getting my book. Unfortunately, the HECRASController does not have any procedures that retrieve water quality output. And the water quality output is stored in a binary file (*.wq##), so it’s not readable/writable as a text file. I put in a request to HEC that they add this to the HECRASController for the 5.0.2 release. Cross your fingers.

Chris

[Chris G.]

You might explore using python to write a script to create a plan hdf file. Unfortunately, there’s no real documentation on how to do this. You can, however, write an output text file that can be read into GIS using GeoRAS. There you should be able to automate the process of creating shapefiles. I’m interested in hearing more about your probabilistic approach to floodplain modeling.

[Chris G.]

Thanks Paolo!

[Chris G.]

When you declare the HECRASController as public, RAS will remain open as an application outside the scope of a given procedure-QuitRAS only hides it. If you want RAS to close at the conclusion of the procedure, you have to declare the HECRASController within the scope of that procedure. Then at the end of the procedure, RAS will close as an application.

It’s possible that you could find some third-party code to manually close RAS as an application. I suggest Googling “vba close an application” and see if you can find a way. If you do, please share with the group.

Thanks-
Chris

[Author]

Posts