search

NOAA-ORR-ERD / PyGnome

The General NOAA Operational Modeling Environment
https://gnome.orr.noaa.gov/doc/pygnome/index.html
Other
54 stars 44 forks source link

Output display format #147

Open larissafranklin23 opened 1 year ago

larissafranklin23 commented 1 year ago

Hello, I'm having some trouble to understand the format of my weathering output file. I'm using this as my script: spill1 = surface_point_line_spill(num_elements=2000, ## qtd de particulas start_position=(-26.21,-11.99), release_time=start_time, end_position=(-26.64, -11.70), end_release_time= start_time + gs.hours(5.1), amount = 20464, ## default volume_units=m^3 substance=gs.GnomeOil(filename='bachequero-heavy_AD00100.json'), units='m^3',
windage_range= (0.01,0.02), windage_persist = 1)

And when I open it to check the number is really different: image

The output format of the data is displaying too many decimal places. How do you handle these weathering readings?

ChrisBarker-NOAA commented 1 year ago

The output is using all the precision of the internal floating point numbers: about 15 decimal digits, which is the default behaviour.

It is a CSV file, so if you load it into, e.g. Excel, you can display as many or as few digits as you like.

That being said, it is a bit silly to store this many digits -- I'll make an issue to limit the CSV output to a more reasonable number.

larissafranklin23 commented 1 year ago

Hello Chris, thank you. I think I wasn't clear in my question.

I'm confused about the result because the amount that was spilled was 20464m^3, but when I look at my results, I find weathering results as evaporated, dispersed, sedimentation and other quantities that are larger than the originally spilled amount.

amount released = 20464, ## default volume_units=m^3 units='m^3',

For example, the result of one of the models run below: evaporated: 2321608.9919537366 dispersed: 2178.419994947777 sedimentation: 111.26444704240757 beached: 25655.915887327097 floating: 3112104.1564391125

Is there any conversion that I need to do ?

dylanrighi commented 1 year ago

The disconnect is that GNOME thinks in mass (kg) internally (correct me if I'm wrong Chris). So it has used your input m^3 and translated to kg using the oil density and then works in kg from there.

On Thu, Jun 15, 2023 at 12:42 PM larissafranklin23 @.***> wrote:

Hello Chris, thank you. I think I wasn't clear in my question.

I'm confused about the result because the amount that was spilled was 20464m^3, but when I look at my results, I find weathering results as evaporated, dispersed, sedimentation and other quantities that are larger than the originally spilled amount.

amount released = 20464, ## default volume_units=m^3 units='m^3',

For example, the result of one of the models run below: evaporated: 2321608.9919537366 dispersed: 2178.419994947777 sedimentation: 111.26444704240757 beached: 25655.915887327097 floating: 3112104.1564391125

Is there any conversion that I need to do ?

— Reply to this email directly, view it on GitHub https://github.com/NOAA-ORR-ERD/PyGnome/issues/147#issuecomment-1593622412, or unsubscribe https://github.com/notifications/unsubscribe-auth/ABQA2MAYRQYK25QBNUFR4QDXLNQQDANCNFSM6AAAAAAZIFBGEE . You are receiving this because you are subscribed to this thread.Message ID: @.***>

-- Dylan Righi, Oceanographer (he/him) @.***

Emergency Response Div (Genwest Systems) Office of Response and Restoration NOAA / National Ocean Service (NOS) 7600 Sand Point Way NE Seattle, WA 98115

ChrisBarker-NOAA commented 1 year ago

Thanks @dylanrighi -- that is correct -- we should put units in the head of that CSV file.