search

MITgcm-contrib / llc_hires

Code customizations and scripts related to LLC4320 family of configurations.
17 stars 12 forks source link

LLC270 STDERR：** WARNING ** S.Edge of facet # 1 disconnected (facet_link= 0.00) #3

Closed liu-ran closed 3 years ago

liu-ran commented 3 years ago

STDERR.0000 (PID.TID 0000.0001) WARNING S.Edge of facet # 1 disconnected (facet_link= 0.00) (PID.TID 0000.0001) WARNING S.Edge of facet # 2 disconnected (facet_link= 0.00) (PID.TID 0000.0001) WARNING E.Edge of facet # 4 disconnected (facet_link= 0.00) (PID.TID 0000.0001) WARNING E.Edge of facet # 5 disconnected (facet_link= 0.00)

When running the LLC270 simulation, the STDERR files show the edge of a facet is disconnected from the other facet. Are these warnings the reason that leads to the discontinuous patterns at the junction of two facets (as shown in the picture) in the simulation? zonal velocity

menemenlis commented 3 years ago

Hi Liu-Ran, the warnings are normal. We use a 6-facet (cube-sphere) topology, with southern facet removed, to set-up llc topology. All the disconnects happen over land (Anatarctica) the southern edges of facets 1-2 and eastern edges of facets 4-5, which are rotated 90 degrees relative to 1-2. So that’s OK.

Your discontinuities are worrisome, however. I have never seen them (but may have missed them). Are you using one of the configurations described here: https://github.com/MITgcm-contrib/llc_hires/tree/master/llc_270 https://github.com/MITgcm-contrib/llc_hires/tree/master/llc_270 or here: https://github.com/MITgcm-contrib/ecco_darwin/tree/master/v04 https://github.com/MITgcm-contrib/ecco_darwin/tree/master/v04 ?

D.

On Aug 24, 2020, at 7:13 AM, liu-ran notifications@github.com wrote:

STDERR.0000 (PID.TID 0000.0001) WARNING S.Edge of facet # 1 disconnected (facet_link= 0.00) (PID.TID 0000.0001) WARNING S.Edge of facet # 2 disconnected (facet_link= 0.00) (PID.TID 0000.0001) WARNING E.Edge of facet # 4 disconnected (facet_link= 0.00) (PID.TID 0000.0001) WARNING E.Edge of facet # 5 disconnected (facet_link= 0.00)

When running the LLC270 simulation, the STDERR files show the edge of a facet is disconnected from the other facet. Are these warnings the reason that leads to the discontinuous patterns at the junction of two facets (as shown in the picture) in the simulation? https://user-images.githubusercontent.com/48266323/91055187-e9a98b80-e656-11ea-8538-a354cef67785.png — You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/MITgcm-contrib/llc_hires/issues/3, or unsubscribe https://github.com/notifications/unsubscribe-auth/ABH7CSOID44BON7VORCYN7LSCJYRZANCNFSM4QJRQKMQ.

liu-ran commented 3 years ago

Hi Liu-Ran, the warnings are normal. We use a 6-facet (cube-sphere) topology, with southern facet removed, to set-up llc topology. All the disconnects happen over land (Anatarctica) the southern edges of facets 1-2 and eastern edges of facets 4-5, which are rotated 90 degrees relative to 1-2. So that’s OK.

Your discontinuities are worrisome, however. I have never seen them (but may have missed them). Are you using one of the configurations described here: https://github.com/MITgcm-contrib/llc_hires/tree/master/llc_270 https://github.com/MITgcm-contrib/llc_hires/tree/master/llc_270 or here: https://github.com/MITgcm-contrib/ecco_darwin/tree/master/v04 https://github.com/MITgcm-contrib/ecco_darwin/tree/master/v04 ?

D. … On Aug 24, 2020, at 7:13 AM, liu-ran @.***> wrote:

STDERR.0000 (PID.TID 0000.0001) WARNING S.Edge of facet # 1 disconnected (facet_link= 0.00) (PID.TID 0000.0001) WARNING S.Edge of facet # 2 disconnected (facet_link= 0.00) (PID.TID 0000.0001) WARNING E.Edge of facet # 4 disconnected (facet_link= 0.00) (PID.TID 0000.0001) WARNING E.Edge of facet # 5 disconnected (facet_link= 0.00)

When running the LLC270 simulation, the STDERR files show the edge of a facet is disconnected from the other facet. Are these warnings the reason that leads to the discontinuous patterns at the junction of two facets (as shown in the picture) in the simulation? https://user-images.githubusercontent.com/48266323/91055187-e9a98b80-e656-11ea-8538-a354cef67785.png — You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub <#3>, or unsubscribe https://github.com/notifications/unsubscribe-auth/ABH7CSOID44BON7VORCYN7LSCJYRZANCNFSM4QJRQKMQ.

Hi Menemenlis, Thanks for your kind answer. The LLC270 configures (e.g. data, data.exch2 …) are from the https://github.com/MITgcm-contrib/llc_hires/tree/master/llc_270. Now I check the zonal velocity in the pickup which is downloaded from the website: https://data.nas.nasa.gov/ecco/data.php?dir=/eccodata/llc_270/iter42/input. The zonal velocity between 1 and 5 facet also exists the discontinuities. test_U

menemenlis commented 3 years ago

How are you plotting? Are you aware that facets 1-2 are rotated 90 degrees from facets 4-5? The model’s U/V direction is relative to the grid cells in the facets, not relative to geographical East/North directions. In particular U in the East edge of facet 2 matches with -V in the South edge of facet 4. In attached, I have indicated the five facets of llc grid: facet 1 = tiles 1-3 facet 2 = tiles 4-6 facet 3 = tiles 7 facet 4 = tiles 8-10 facet 5 = tiles 11-13

On Aug 24, 2020, at 8:44 PM, liu-ran notifications@github.com wrote:

The zonal velocity between 1 and 5 facet also exists the discontinuities.

menemenlis commented 3 years ago

with attachment

On Aug 24, 2020, at 9:21 PM, Dimitris Menemenlis dmenemenlis@gmail.com wrote:

How are you plotting? Are you aware that facets 1-2 are rotated 90 degrees from facets 4-5? The model’s U/V direction is relative to the grid cells in the facets, not relative to geographical East/North directions. In particular U in the East edge of facet 2 matches with -V in the South edge of facet 4. In attached, I have indicated the five facets of llc grid: facet 1 = tiles 1-3 facet 2 = tiles 4-6 facet 3 = tiles 7 facet 4 = tiles 8-10 facet 5 = tiles 11-13

On Aug 24, 2020, at 8:44 PM, liu-ran <notifications@github.com mailto:notifications@github.com> wrote:

The zonal velocity between 1 and 5 facet also exists the discontinuities.

liu-ran commented 3 years ago

Thank you very much! I try it again as you said and see the correct figure. I knew the 4,5 facets are rotated through the downloaded llc270 bathymetry. But I didn't realize that the model’s U/V direction is not relative to geographical East/North directions. Thanks again. test_U

How are you plotting? Are you aware that facets 1-2 are rotated 90 degrees from facets 4-5? The model’s U/V direction is relative to the grid cells in the facets, not relative to geographical East/North directions. In particular U in the East edge of facet 2 matches with -V in the South edge of facet 4. In attached, I have indicated the five facets of llc grid: facet 1 = tiles 1-3 facet 2 = tiles 4-6 facet 3 = tiles 7 facet 4 = tiles 8-10 facet 5 = tiles 11-13 …