Inconsistencies between versions 5.9 and 5.10

[brey]

Dear All, running the same model with versions 5.9 and 5.10.1 gives different results. There are some small differences in the written maxelev data and some larger differences in the maxdahv values. I don't know if they grow in time.

The most important difference is that the model produces a singular high (negative) value on a single mesh point. See figure.

I am attaching the model for clarity. The only difference in the param.nml is the configuration of the output:

❯ diff param.nml ../schism59/param.nml
150,151c150,151
<     iof_hydro(1:31) = 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
<                       0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
---
>     iof_hydro(1:30) = 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
>                       0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0

I would appreciate your help in sorting this out as we are moving to the new IO.

test_model.zip

[josephzhang8]

Hi George: The differences may or may not be important, and may be related to some small changes in wet/dry rule in hydrologic regime. The default meth_sink=1 in newer version prevents ‘groundwater’ level from sinking below an amount. In looking at the elev results, always compare those with the local bottom elev. If it’s dry, the value of surface elev is irrelevant. In the case you showed below, there could be some topologic anomaly around that node that led to the water being trapped there. Or maybe init condition for elev plays a role… see the following chapter in manual: https://schism-dev.github.io/schism/master/case-study.html#hydrologic-flow Regards, ------------------------------- Joseph Zhang (804)684 7595 (office) SCHISM web: http://ccrm.vims.edu/schism/ From: George Breyiannis ***@***.***> Sent: Sunday, February 26, 2023 1:45 PM To: schism-dev/schism ***@***.***> Cc: Subscribed ***@***.***> Subject: [schism-dev/schism] Inconsistencies between versions 5.9 and 5.10 (Issue #99) [EXTERNAL to VIMS received message] Dear All, running the same model with versions 5.9 and 5.10.1 gives different results. There are some small differences in the written maxelev data and some larger differences in the maxdahv values. I don't know if they grow in time. The most important difference is that the model produces a singular high (negative) value on a single mesh point. See figure. I am attaching the model for clarity. The only difference in the param.nml is the configuration of the output: ❯ diff param.nml ../schism59/param.nml 150,151c150,151 < iof_hydro(1:31) = 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, < 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0

---

iof_hydro(1:30) = 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0

I would appreciate your help in sorting this out as we are moving to the new IO. [Figure 3]<https://nam11.safelinks.protection.outlook.com/?url=https%3A%2F%2Fuser-images.githubusercontent.com%2F5442433%2F221429959-a91a60f2-497d-428a-b4b2-905a645ef1e6.png&data=05%7C01%7Cyjzhang%40vims.edu%7Cafd10a780a55411ce85e08db182987f6%7C8cbcddd9588d4e3b9c1e2367dbdf1740%7C0%7C1%7C638130338878759202%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=88OYa16iV%2FDX%2Bus16OEIw0FYItz%2BL8BSXsCQ3bOKqGI%3D&reserved=0> test_model.zip<https://nam11.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2Fschism-dev%2Fschism%2Ffiles%2F10834438%2Ftest_model.zip&data=05%7C01%7Cyjzhang%40vims.edu%7Cafd10a780a55411ce85e08db182987f6%7C8cbcddd9588d4e3b9c1e2367dbdf1740%7C0%7C1%7C638130338878759202%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=B2dCj6B0xd7nO8QllngsuSS%2BKipAJBwKuojkWnw%2BKTo%3D&reserved=0> — Reply to this email directly, view it on GitHub<https://nam11.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2Fschism-dev%2Fschism%2Fissues%2F99&data=05%7C01%7Cyjzhang%40vims.edu%7Cafd10a780a55411ce85e08db182987f6%7C8cbcddd9588d4e3b9c1e2367dbdf1740%7C0%7C1%7C638130338878759202%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=aFb%2Bz4LP35QwZQNS78AW4XK0XsG6TICuXSWRr8gle7Q%3D&reserved=0>, or unsubscribe<https://nam11.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2Fnotifications%2Funsubscribe-auth%2FAFBKNZ7IWJB4TPCWZXKATXLWZOQB3ANCNFSM6AAAAAAVITTOEM&data=05%7C01%7Cyjzhang%40vims.edu%7Cafd10a780a55411ce85e08db182987f6%7C8cbcddd9588d4e3b9c1e2367dbdf1740%7C0%7C1%7C638130338878759202%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=V7aGTNWhgN87EhWp7vcrQ5ihYuUCenT22q1FPiwyshk%3D&reserved=0>. You are receiving this because you are subscribed to this thread.Message ID: ***@***.******@***.***>>