[MITgcm-support] Strange output frequency

Mon May 15 08:48:01 EDT 2023

If you use the diagnostics pkg and specify snapshot every 5 seconds, then these will be, by default, in the middle of the 5s interval. Your timestep prevents that (can’t reach 2.5), so it uses the closest. If you want out put at 5s,10s, etc. you need to specify the timePhase = 0., or 5., 

I have no idea why it sometime goes wrong. Maybe the rounding sometimes implies .4 and sometimes .6 as the neares to .5?

Martin

> On 15. May 2023, at 14:41, Gregorio Sandy <Sandy.Gregorio at uqar.ca> wrote:
> 
> Hi everyone,
> 
> I am running some experiments with the MITgcm. It was working fine until I increase my resolution and decrease my timestep (0.2s). I save my diagnostics as snap-shot every 5 seconds. When I look at my output, I have snap-shots filled with NaN. I checked my time output and I find some inconsistencies. While I still have output every 5 s, I have, at some times, output at 0.2s. It happens at around 1 min and around 17 min of my simulation, then after that, all my output are correct (I run my simulation 60 min). In total, I have more output than I should have (due to this strange output at 0.2s). Here are my first time output (in bold red are the output at 0.2s):
> 
>  T = 2.4, 7.4, 12.4, 17.4, 22.4, 27.4, 32.4, 37.4, 42.4, 47.4, 52.4, 57.4,
>     62.4, 67.4, 67.6, 72.4, 72.6, 77.4, 82.4, 87.4, 92.4, 97.4, 102.4, 107.4,
>     112.4, 117.4, 122.4, 127.4, 132.4, 137.4, 142.4, 147.4, 152.4, 157.4,
>     162.4, 167.4, 172.4, 177.4, 182.4, 187.4, 192.4, 197.4, 202.4, 207.4,
>     212.4, 217.4, 222.4, 227.4, 232.4, 237.4, 242.4, 247.4, 252.4, 257.4,
>     262.4, 267.4, 272.4, 277.4, 282.4, 287.4, 292.4, 297.4, 302.4, 307.4,
>     312.4, 317.4, 322.4, 327.4, 332.4, 337.4, 342.4, 347.4, 352.4, 357.4,
>     362.4, 367.4, 372.4, 377.4, 382.4, 387.4, 392.4, 397.4, 402.4, 407.4,
>     412.4, 417.4, 422.4, 427.4, 432.4, 437.4, 442.4, 447.4, 452.4, 457.4,
>     462.4, 467.4, 472.4, 477.4, 482.4, 487.4, 492.4, 497.4, 502.4, 507.4,
>     512.4, 517.4, 522.4, 527.4, 532.4, 537.4, 542.4, 547.4, 552.4, 557.4,
>     562.4, 567.4, 572.4, 577.4, 582.4, 587.4, 592.4, 597.4, 602.4, 607.4,
>     612.4, 617.4, 622.4, 627.4, 632.4, 637.4, 642.4, 647.4, 652.4, 657.4,
>     662.4, 667.4, 672.4, 677.4, 682.4, 687.4, 692.4, 697.4, 702.4, 707.4,
>     712.4, 717.4, 722.4, 727.4, 732.4, 737.4, 742.4, 747.4, 752.4, 757.4,
>     762.4, 767.4, 772.4, 777.4, 782.4, 787.4, 792.4, 797.4, 802.4, 807.4,
>     812.4, 817.4, 822.4, 827.4, 832.4, 837.4, 842.4, 847.4, 852.4, 857.4,
>     862.4, 867.4, 872.4, 877.4, 882.4, 887.4, 892.4, 897.4, 902.4, 907.4,
>     912.4, 917.4, 922.4, 927.4, 932.4, 937.4, 942.4, 947.4, 952.4, 957.4,
>     962.4, 967.4, 972.4, 977.4, 982.4, 987.4, 992.4, 997.4, 1002.4, 1007.4,
>     1012.4, 1017.4, 1022.4, 1027.4, 1027.6, 1032.4, 1032.6, 1037.4, 1037.6, 
>     1042.4, 1042.6, 1047.4, 1047.6, 1052.4, 1052.6, 1057.4, 1057.6, 1062.4,
>     1062.6, 1067.4, 1067.6, 1072.4, 1072.6, 1077.4, 1077.6, 1082.4, 1082.6, 
>     1087.4, 1087.6, 1092.4, 1092.6, 1097.4, 1097.6, 1102.4, 1102.6, 1107.4,
>     1107.6, 1112.4, 1112.6, 1117.4, 1117.6, 1122.4, 1122.6, 1127.4, 1127.6, 
>     1132.4, 1132.6, 1137.4, 1137.6, 1142.4, 1142.6, 1147.4, 1147.6, 1152.4,
>     1152.6, 1157.4, 1157.6, 1162.4, 1162.6, 1167.4, 1167.6, 1172.4, 1172.6, 
>     1177.4, 1177.6, 1182.4, 1182.6, 1187.4, 1187.6, 1192.4, 1192.6, 1197.4,
>     1197.6, 1202.4, 1202.6, 1207.4, 1207.6, 1212.4, 1212.6, 1217.4, 1217.6, 
>     1222.4, 1222.6, 1227.4, 1227.6, 1232.4, 1237.4, 1242.4, 1247.4, 1252.4,
>     1257.4, 1262.4, 1267.4, 1272.4, 1277.4, 1282.4, 1287.4, 1292.4, 1297.4,
>     1302.4, 1307.4, 1312.4, 1317.4, 1322.4, 1327.4, 1332.4, 1337.4, 1342.4,
>     1347.4, 1352.4, 1357.4, 1362.4, 1367.4, 1372.4, 1377.4, 1382.4, 1387.4,
>     1392.4, 1397.4, 1402.4, 1407.4, 1412.4, 1417.4, 1422.4, 1427.4, 1432.4,
>     1437.4, 1442.4, 1447.4, 1452.4, 1457.4, 1462.4, 1467.4, 1472.4, 1477.4,
>     1482.4, 1487.4, 1492.4, 1497.4, 1502.4, 1507.4, 1512.4, 1517.4, 1522.4,
>     1527.4, 1532.4, 1537.4, 1542.4, 1547.4, 1552.4, 1557.4, 1562.4, 1567.4,
>     1572.4, 1577.4, 1582.4, 1587.4, 1592.4, 1597.4, 1602.4, 1607.4, 1612.4,
>     1617.4, 1622.4, 1627.4, 1632.4, 1637.4, 1642.4, 1647.4, 1652.4, 1657.4,
>     1662.4, 1667.4, 1672.4, 1677.4, 1682.4, 1687.4, 1692.4, 1697.4, 1702.4,
> 
> Do you have any idea about what is going on? Also I was wondering, as a naïve question, why the first output is at 2.4s?
> 
> Many thanks for your help.
> 
> Cheers,
> 
> Sandy
>  
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org <mailto:MITgcm-support at mitgcm.org>
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20230515/80dbc41e/attachment-0001.html>