fgmax_grids.data for Hurricane Irma

[siamak daneshvaran]

Hi

I ran a test of hurricane Irma and added two locations using fgmax_grids.data and got the below results after the code completed the analysis: 

-0.81600000000E+02   0.24130000000E+02   0  -0.99999000E+99  -0.99999000E+99  -0.99999000E+99  -0.99999000E+99  -0.99999000E+99  -0.99999000E+99

  -0.81788500000E+02   0.24700000000E+02   7  -0.89079590E+01   0.98420946E+01   0.58911453E+00   0.43624195E+04   0.33754113E+04   0.33754113E+04

Any idea why i am getting  -0.99999000E+99  for the first point? and for the second point a negative number and then positive, which is the maximum height proably.

Many thanks

Siamak Daneshvaran

[siamak daneshvaran]

Hi,

I want to add some 20 points to estimate storm surge at these points. Used the below code and input file that has the format needed to input the points.the code ran and did not give any error.  But, I get 0.9*E99 every where. Not sure what is wrong, but all points are supposed to have at least 0.5 meter storm surge above mean see level.  I appreciate your help, as I do not have much experience with this.  Thanks, Siamak.   

-0.81719732793E+02   0.23999060220E+02   0  -0.99999000E+99  -0.99999000E+99  -0.99999000E+99  -0.99999000E+99  -0.99999000E+99  -0.99999000E+99

   

 fg = fgmax_tools.FGmaxGrid()
    fgmax_grids = rundata.fgmax_data.fgmax_grids
    tstart_max = 8000.
    fg.point_style = 0
    fg.min_level_check = amrdata.amr_levels_max # which levels to monitor max on
    fg.tstart_max = tstart_max  # just before wave arrives
    fg.tend_max = 1.e10    # when to stop monitoring max values
    fg.dt_check = 20.      # how often to update max values
    # can specify that list of points is in a different file:
    fg.npts = 0
    fg.xy_fname = 'fgmax_points_list.data'
    fgmax_grids.append(fg)  # written to fgmax_grids.data

    return rundata

[siamak daneshvaran]

Ultimately I need to add coastal locations and extract water depth for these locations (storm elevation - topo). The number of locations are around 1000 points. 

I really appreciate help on what is best way to do this. 

Thanks 

Siamak D 

--
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/Xml33vjbElo/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/5b4247d6-22f2-445b-a09f-05a5b1c46afdn%40googlegroups.com.

[Kyle Mandli]

A couple of thoughts:

• Have you looked at the 2D output to make sure everything is working as you expect?  Those values look like they are either not being filled in correctly or the simulation is blowing up at some point.
• Gauges may also work for you depending on your use case.  That would give a time trace as to what is going on that may be helpful, at least until you trust the fgmax data.

Kyle

--
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/69b4395b-7b89-40e3-bf84-4d4b37e60095n%40googlegroups.com.

[Randall J LeVeque]

Hi Siamak,

You should take a look at

  https://www.clawpack.org/fgmax.html#format-of-the-output-files

for information about the columns in the fgmax file.  In particular, the first 5 columns are

  longitude, latitude, AMR level, B, h

where B is the topography and h the max inundation depth.

The value -0.99999000E+99 is stored before an fgmax point gets updated for the first time.  I notice you have 

    fg.tstart_max = tstart_max = tstart_max = 8000.

so nothing gets updated before this time.  If your run had a smaller tfinal then nothing would be updated in fgmax.

Also you have

   fg.min_level_check = amrdata.amr_levels_max 

which means it only monitors fgmax on the finest level, so any fgmax point that lies outside the union of the finest level grids at all times will never get updated.  If you want to monitor some points close to shore where you have a very fine grid and other points far offshore where there are only coarser grids, you may want to define a couple different fgmax grids with different values of fg.min_level_check.

See 

   https://www.clawpack.org/fgmax.html

for some more details.

Hope that helps,

  Randy

On Fri, Nov 24, 2023 at 12:30 PM siamak daneshvaran <siamak.da...@gmail.com> wrote:

Hi, I want to add some 20 points to estimate storm surge at these points. Used the below code and input file that has the format needed to input the points. the code ran and did not give any error. But, I get 0. 9*E99 every where. Not sure what

ZjQcmQRYFpfptBannerStart

This Message Is From an Untrusted Sender

You have not previously corresponded with this sender.

See https://itconnect.uw.edu/email-tags for additional information. Please contact the UW-IT Service Center, he...@uw.edu 206.221.5000, for assistance.

 

ZjQcmQRYFpfptBannerEnd

--

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/5b4247d6-22f2-445b-a09f-05a5b1c46afdn%40googlegroups.com.

[siamak daneshvaran]

Thanks Kyle. My application is estimating losses in catastrophe modeling. Ultimately, I need max water depth (elevation- topo) for each point. I thought putting all the locations (homes) in a fgmax_grids.data file will be good. 

Referring to your email, how do I look at the 2D output?  I have looked at the animations for the whole storms and it looks ok. 

I ran another case with larger fgmax_grids.data file and 2 of the points worked and rest have 99…. 

Would love to find a correct way to do this for loss modeling. 

Many thanks!

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

[siamak daneshvaran]

Dear Randall,

Many thanks for your reply. It was really useful to see the header for the gmax file. So the fifth column is actually the depth (elevation -topt) as I need. I made a  10000  points input file with option zero (reading the .data file). YEs, unfortunately a large number of points (location) give no values (after I  changed the tstart_max to 4 instead of 8000). Not sure what I can change next or what is my best option to get depth data for my locations. 

Many thanks

Siamak 

On Sat, Nov 25, 2023 at 3:03 PM Randall J LeVeque <r...@uw.edu> wrote:

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