Elevation/Speed Timeseries

[Behzad Ebrahimi]

Hi,

Could you suggest me a procedure to create time series of water elevation and speed for some specific points after finishing the simulation and using the output files? Supposing that I don't know the points before run and cannot use the gauge option.

Thanks,

Behzad

[Randall J LeVeque]

Dear Behzad,

If you want a time series at a point where there was no gauge during the simulation, you could extract values from the output frames and if you have output the solution frequently enough this might be sufficient.

To extract the value at a point from the finest level AMR patch available there, you could use the utilities now in visclaw gridtools module (as of v5.7.0).  This isn't well documented yet, but there is a notebook with examples of how to extract a 2d uniform grid or a 1d transect in the gallery,

  http://www.clawpack.org/gallery/_static/apps/notebooks/visclaw/gridtools.html

To extract at a single point or list of points you could call gridtools.grid_output_2d with xout,yout being numpy arrays containing the x and y location(s).

 - Randy

--
You received this message because you are subscribed to the Google Groups "claw-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to claw-users+...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/claw-users/27eceb3d-d23a-4412-b821-7817f0674e58n%40googlegroups.com.

[Behzad Ebrahimi]

Hi Randy,

I truly appreciate your useful information. I'll go through this and let you know if I have any questions.

Best,

Behzad

You received this message because you are subscribed to a topic in the Google Groups "claw-users" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/claw-users/czBGrn59QtI/unsubscribe.
To unsubscribe from this group and all its topics, send an email to claw-users+...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/claw-users/CADwYU3%2B6NvQE%3Dw%2BccpKAEo27QCQ%2Bf%2Bit_iL8wZf5EVj7VsTFBg%40mail.gmail.com.

[Behzad Ebrahimi]

Dear Randy,

I was trying the gridtool Jupyter notebook you suggested in your email, but I faced a problem and couldn't solve it. The screenshot is attached. Could you please take a look at it?

Thanks,

Behzad

On Fri, Sep 25, 2020 at 10:12 AM Randall J LeVeque <r...@uw.edu> wrote:

You received this message because you are subscribed to a topic in the Google Groups "claw-users" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/claw-users/czBGrn59QtI/unsubscribe.
To unsubscribe from this group and all its topics, send an email to claw-users+...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/claw-users/CADwYU3%2B6NvQE%3Dw%2BccpKAEo27QCQ%2Bf%2Bit_iL8wZf5EVj7VsTFBg%40mail.gmail.com.

[Kyle Mandli]

Hi Behzad,

Can you share the entire Jupyter notebook you are using by chance?

Kyle

To view this discussion on the web visit https://groups.google.com/d/msgid/claw-users/CA%2BiN8-f-%3D7eW0XgSOn3M_S2YUHhMyjwqZV9C%3DZLB8Ppdw7eqng%40mail.gmail.com.

[Behzad Ebrahimi]

Hi Kyle,

Sure, I have attached it to this email. Although based on the documentation above the notebook was supposed to be located under "amrclaw" I found it in this directory(under "visclaw"):

$CLAW/apps/notebooks/amrclaw/gridtools.ipynb

The error I receive might be due to the changes I made in the run directory and desired output grid to apply the notebook functions on Chile2010 example.

Thanks for your help in advance.

Behzad

To view this discussion on the web visit https://groups.google.com/d/msgid/claw-users/1c9c1f33-6885-4b13-a0ff-1fb773cffe2b%40Spark.

[Randall J LeVeque]

The problem is that you specified

ylat = 0.002
yout = ylat * ones(xout.shape)

and are asking for a transect at a latitude above the equator.  But the Chile2010 problem as set up in the examples has -60 <= y <=0, so this transect is outside the domain.

When I change to a ylat value in this range it seems to work fine.

But it shouldn't give this error -- it should return a 2d array of nan values since you specified 2 output variables, rather than a 1d array of nan's.  I'll take a look.

 - Randy

To view this discussion on the web visit https://groups.google.com/d/msgid/claw-users/CA%2BiN8-fZR7k%3DAP%2BkmF4UEF5mCGB0e%2Bk7gquXfyd5h_bN8qPipQ%40mail.gmail.com.

[Behzad Ebrahimi]

Thanks for your prompt response Randy.

Behzad

To view this discussion on the web visit https://groups.google.com/d/msgid/claw-users/CADwYU3JZLKD%2Boxa313Pw7oLcUOSJPxRV8Wi46w3Z49ZQtechVw%40mail.gmail.com.