Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

try fixing testing #39

Open
wants to merge 75 commits into
base: main
Choose a base branch
from
Open

try fixing testing #39

wants to merge 75 commits into from

Conversation

bolliger32
Copy link
Member

No description provided.

rjleveque and others added 30 commits April 14, 2024 08:25
Previous value -9999 did not have enough digits since new etopo 2022
data has this as an actual topography value in the Mariana Trench
(when written with `Z_format='%.0f'`).
change nodata_value to -99999 in topotools.py
Fix bug related to tracking pressure at gauges
These do not need to be passed into update, unpacked, and repacked,
since the objects created in the script will be used into update.
If blit==False then they are not needed at all.
The default is now '%.3f', millimeter resolution, making smaller dtopo
files than previously.
Improve chile2010_fgmax-fgout example script for fgout animation
Add dZ_format parameter to DTopography.write function
…iemann/src

Note that rpn2_geoclaw.f and geoclaw_riemann_utils.f were update in riemann/src
to handle 5 waves and rpt2_geoclaw_sym.f was discarded in favor of
standard rpt2.f, which works fine with 5 waves.
Do not want to touch the old topotools file for reference.
Includes a couple of minor bugs related to gauge plotting.  Major change
involves how we now plot the gauge data.
mandli and others added 13 commits June 6, 2024 14:33
improvements to fgout_module for t arrays and D-Claw
Now that `output_style==2` is supported for fgout grids (an array of times,
see clawpack#617), the function
fgout_tools.FGoutGrid.read_fgout_grids_data needs to be fixed to properly
read in the new format of `fgout_grids.data`, also note that `nout`
now comes before `tstart` and `tend` for `output_style==1`.
fix fgout_tools.FGoutGrid.read_fgout_grids_data for time array
@bolliger32
Copy link
Member Author

Waiting on clawpack#610 before these new tests are functional

rjleveque and others added 16 commits July 13, 2024 18:17
Modifications to fgout_module.f90 and fgout_tools.py, primarily.
Can now set `fgout.q_out_vars` to a list of which components of q
to output for each fgout frame (Fortran indexing). The default is
`fgout.q_out_vars = [1,2,3,4]` for all components of q[1:3] and
also eta (4), consistent with the previous behavior.

The user could also/instead ask to output the topo B
as component 5.  (If two out of three of h, eta, B are output then the
other can be computed from these.)

Note that this list is written out to `fgout_grids.data` rather than a
boolean list of True/False values for each possible component as is done
in other places in GeoClaw.
Add Plotting Capabilities to Storm Object
… D-Claw

When instantiating fgout_tools.FGoutGrid object, new qmap parameter can be set
to 'geoclaw', 'geoclaw-bouss', 'dclaw', or a custom dictionary.
(Eventually add 'dclaw-bouss' option also).
This dictionary has the form `{'h':1, 'hu':2, etc}` with mapping from fgout
variable names to Fortran indices of q in the code creating the fgout output.
The `character` variables that were storing variable names and IDs for NetCDF topography reading were only length 10.  This just makes them 64 characters long.  Probably too long but not really a big deal.
Now need to call fgout_grid.read_fgout_grids_data() explicitly,
but do not need to set format.
this is ignored since it is read from fgout_grids.data.
Also properly handle output_format == 2 for binary32.
Allow specifying which components of q to output on each fgout grid
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants