GeoscienceAustralia / tcrm

A statistical-parametric model for assessing wind hazard from tropical cyclones
http://geoscienceaustralia.github.io/tcrm
Other
83 stars 53 forks source link

Problem in example simulation #151

Open busy67 opened 3 months ago

busy67 commented 3 months ago

Dear expert,

I met a problem in running the Tropical Cyclone Risk Model. 1

The output hazard maps of the example simulation: python tcrm.py -c example/port_hedland.ini (nothing modified, data downloaded from ibtracs manually) showed very different from https://geoscienceaustralia.github.io/tcrm/docs/examples.html.

On the right figure was my output, you can see scattered contors and obviously lower wind speed. I tried different regions, and tired to decrease the "Initial pressure difference too small" error in the log file by modifying "def validInitPressure(track)" in TrackGenerator.py to: if track.EnvPressure[0] - track.CentralPressure[0] < 0.5: but things didnt get better.

It is very kind of you if you can help me. What is probably wrong?

wcarthur commented 2 months ago

@busy67 the results you are getting are consistent with the latest model release. I've not updated the figures in the documentation since a much earlier release.

The patchiness results from the limited number of synthetic tracks used in the example. The calculation of the ARI wind speeds requires a minimum number of values, which may not be met by the number set in the example (1000 tracks). Changing the pressure difference in TrackGenerator.py will not directly influence this result.