segrelab / comets

Computation of Microbial Ecosystems in Time and Space
19 stars 10 forks source link

Timestep length-related instability #17

Closed mquintin closed 4 years ago

mquintin commented 4 years ago

Attached are two layouts and their associated files. They have radically different biomass output despite the only difference being that one has a timestep of 10h, the other is 10.0001h

comets_instability.zip

mquintin commented 4 years ago

Resolved as of version 2.7.4 (or at least no longer duplicating this particular error)

Results follow from running the attached layouts with that version. Models 0 and 1 use v2.7.4 and are timesteps 10 and 10.0001 respectively, models 2 and 3 use the old code.

image

t x y z model biomass modelName 0 0 0 0 0 0.000135000000000000 'TS-10' 1 0 0 0 0 0.000135000000000000 'TS-10' 2 0 0 0 0 0.000135000000000000 'TS-10' 3 0 0 0 0 0.000135000000000000 'TS-10' 4 0 0 0 0 0.000135000000000000 'TS-10' 5 0 0 0 0 0.000135000000000000 'TS-10' 6 0 0 0 0 0.000135000000000000 'TS-10' 7 0 0 0 0 0.000135000000000000 'TS-10' 8 0 0 0 0 0.000135000000000000 'TS-10' 9 0 0 0 0 0.000135000000000000 'TS-10' 10 0 0 0 0 0.000135000000000000 'TS-10' 11 0 0 0 0 0.000135000000000000 'TS-10' 12 0 0 0 0 0.000135000000000000 'TS-10' 13 0 0 0 0 0.000135000000000000 'TS-10' 14 0 0 0 0 0.000135000000000000 'TS-10' 15 0 0 0 0 0.000135000000000000 'TS-10' 16 0 0 0 0 0.000135000000000000 'TS-10' 17 0 0 0 0 0.000135000000000000 'TS-10' 18 0 0 0 0 0.000135000000000000 'TS-10' 19 0 0 0 0 0.000135086478870000 'TS-10' 20 0 0 0 0 0.000135193714640000 'TS-10' 21 0 0 0 0 0.000135315065350000 'TS-10' 22 0 0 0 0 0.000135444380790000 'TS-10' 23 0 0 0 0 0.000135575737270000 'TS-10' 24 0 0 0 0 0.000135703889520000 'TS-10' 25 0 0 0 0 0.000135824153910000 'TS-10' 26 0 0 0 0 0.000135932091860000 'TS-10' 27 0 0 0 0 0.000136024466800000 'TS-10' 28 0 0 0 0 0.000136097254260000 'TS-10' 29 0 0 0 0 0.000136147300300000 'TS-10' 30 0 0 0 0 0.000136172555490000 'TS-10' 31 0 0 0 0 0.000136172555490000 'TS-10' 32 0 0 0 0 0.000136196091140000 'TS-10' 33 0 0 0 0 0.000136196091140000 'TS-10' 34 0 0 0 0 0.000136206866150000 'TS-10' 35 0 0 0 0 0.000136206866150000 'TS-10' 36 0 0 0 0 0.000136206866150000 'TS-10' 37 0 0 0 0 0.000136212955150000 'TS-10' 38 0 0 0 0 0.000136212955150000 'TS-10' 39 0 0 0 0 0.000136212955150000 'TS-10' 40 0 0 0 0 0.000136212955150000 'TS-10' 41 0 0 0 0 0.000136212955150000 'TS-10' 42 0 0 0 0 0.000136221604360000 'TS-10' 43 0 0 0 0 0.000136221604360000 'TS-10' 44 0 0 0 0 0.000136221604360000 'TS-10' 45 0 0 0 0 0.000136221604360000 'TS-10' 46 0 0 0 0 0.000136221604360000 'TS-10' 47 0 0 0 0 0.000136221604360000 'TS-10' 48 0 0 0 0 0.000136221604360000 'TS-10' 49 0 0 0 0 0.000136221604360000 'TS-10' 50 0 0 0 0 0.000136223242790000 'TS-10' 51 0 0 0 0 0.000136223242790000 'TS-10' 52 0 0 0 0 0.000136223242790000 'TS-10' 53 0 0 0 0 0.000136223242790000 'TS-10' 54 0 0 0 0 0.000136223242790000 'TS-10' 55 0 0 0 0 0.000136223242790000 'TS-10' 56 0 0 0 0 0.000136223242790000 'TS-10' 57 0 0 0 0 0.000136223242790000 'TS-10' 58 0 0 0 0 0.000136223242790000 'TS-10' 59 0 0 0 0 0.000136223242790000 'TS-10' 60 0 0 0 0 0.000136223242790000 'TS-10' 61 0 0 0 0 0.000136223242790000 'TS-10' 62 0 0 0 0 0.000136223242790000 'TS-10' 63 0 0 0 0 0.000136223242790000 'TS-10' 64 0 0 0 0 0.000136223242790000 'TS-10' 65 0 0 0 0 0.000136223242790000 'TS-10' 66 0 0 0 0 0.000136223242790000 'TS-10' 67 0 0 0 0 0.000136223242790000 'TS-10' 68 0 0 0 0 0.000136223242790000 'TS-10' 69 0 0 0 0 0.000136223242790000 'TS-10' 70 0 0 0 0 0.000136223242790000 'TS-10' 71 0 0 0 0 0.000136223242790000 'TS-10' 72 0 0 0 0 0.000136223242790000 'TS-10' 73 0 0 0 0 0.000136223242790000 'TS-10' 74 0 0 0 0 0.000136223242790000 'TS-10' 75 0 0 0 0 0.000136223242790000 'TS-10' 76 0 0 0 0 0.000136223242790000 'TS-10' 77 0 0 0 0 0.000136223242790000 'TS-10' 78 0 0 0 0 0.000136223242790000 'TS-10' 79 0 0 0 0 0.000136223242790000 'TS-10' 80 0 0 0 0 0.000136223242790000 'TS-10' 81 0 0 0 0 0.000136223242790000 'TS-10' 82 0 0 0 0 0.000136223242790000 'TS-10' 83 0 0 0 0 0.000136223242790000 'TS-10' 84 0 0 0 0 0.000136223242790000 'TS-10' 85 0 0 0 0 0.000136223242790000 'TS-10' 86 0 0 0 0 0.000136223242790000 'TS-10' 87 0 0 0 0 0.000136223242790000 'TS-10' 88 0 0 0 0 0.000136223242790000 'TS-10' 89 0 0 0 0 0.000136223242790000 'TS-10' 90 0 0 0 0 0.000136223242790000 'TS-10' 91 0 0 0 0 0.000136223242790000 'TS-10' 92 0 0 0 0 0.000136223242790000 'TS-10' 93 0 0 0 0 0.000136223242790000 'TS-10' 94 0 0 0 0 0.000136223242790000 'TS-10' 95 0 0 0 0 0.000136223242790000 'TS-10' 96 0 0 0 0 0.000136223242790000 'TS-10' 97 0 0 0 0 0.000136223242790000 'TS-10' 98 0 0 0 0 0.000136223242790000 'TS-10' 99 0 0 0 0 0.000136223242790000 'TS-10' 100 0 0 0 0 0.000136223242790000 'TS-10' 101 0 0 0 0 0.000136223242790000 'TS-10' 0 0 0 0 1 0.000135000000000000 'TS-10.0001' 1 0 0 0 1 0.000135000000000000 'TS-10.0001' 2 0 0 0 1 0.000135000000000000 'TS-10.0001' 3 0 0 0 1 0.000135000000000000 'TS-10.0001' 4 0 0 0 1 0.000135000000000000 'TS-10.0001' 5 0 0 0 1 0.000135000000000000 'TS-10.0001' 6 0 0 0 1 0.000135000000000000 'TS-10.0001' 7 0 0 0 1 0.000135000000000000 'TS-10.0001' 8 0 0 0 1 0.000135000000000000 'TS-10.0001' 9 0 0 0 1 0.000135000000000000 'TS-10.0001' 10 0 0 0 1 0.000135000000000000 'TS-10.0001' 11 0 0 0 1 0.000135000000000000 'TS-10.0001' 12 0 0 0 1 0.000135000000000000 'TS-10.0001' 13 0 0 0 1 0.000135000000000000 'TS-10.0001' 14 0 0 0 1 0.000135000000000000 'TS-10.0001' 15 0 0 0 1 0.000135000000000000 'TS-10.0001' 16 0 0 0 1 0.000135000000000000 'TS-10.0001' 17 0 0 0 1 0.000135000000000000 'TS-10.0001' 18 0 0 0 1 0.000135000000000000 'TS-10.0001' 19 0 0 0 1 0.000135086500710000 'TS-10.0001' 20 0 0 0 1 0.000135193751000000 'TS-10.0001' 21 0 0 0 1 0.000135315115680000 'TS-10.0001' 22 0 0 0 1 0.000135444422710000 'TS-10.0001' 23 0 0 0 1 0.000135575790110000 'TS-10.0001' 24 0 0 0 1 0.000135704001850000 'TS-10.0001' 25 0 0 0 1 0.000135824056130000 'TS-10.0001' 26 0 0 0 1 0.000135932009860000 'TS-10.0001' 27 0 0 0 1 0.000136024041530000 'TS-10.0001' 28 0 0 0 1 0.000136096830920000 'TS-10.0001' 29 0 0 0 1 0.000136146876410000 'TS-10.0001' 30 0 0 0 1 0.000136172634760000 'TS-10.0001' 31 0 0 0 1 0.000136172634760000 'TS-10.0001' 32 0 0 0 1 0.000136195473250000 'TS-10.0001' 33 0 0 0 1 0.000136195473250000 'TS-10.0001' 34 0 0 0 1 0.000136204948650000 'TS-10.0001' 35 0 0 0 1 0.000136204948650000 'TS-10.0001' 36 0 0 0 1 0.000136204948650000 'TS-10.0001' 37 0 0 0 1 0.000136208944370000 'TS-10.0001' 38 0 0 0 1 0.000136208944370000 'TS-10.0001' 39 0 0 0 1 0.000136208944370000 'TS-10.0001' 40 0 0 0 1 0.000136208944370000 'TS-10.0001' 41 0 0 0 1 0.000136208944370000 'TS-10.0001' 42 0 0 0 1 0.000136214520700000 'TS-10.0001' 43 0 0 0 1 0.000136214520700000 'TS-10.0001' 44 0 0 0 1 0.000136214520700000 'TS-10.0001' 45 0 0 0 1 0.000136214520700000 'TS-10.0001' 46 0 0 0 1 0.000136214520700000 'TS-10.0001' 47 0 0 0 1 0.000136214520700000 'TS-10.0001' 48 0 0 0 1 0.000136214520700000 'TS-10.0001' 49 0 0 0 1 0.000136214520700000 'TS-10.0001' 50 0 0 0 1 0.000136214520700000 'TS-10.0001' 51 0 0 0 1 0.000136216257940000 'TS-10.0001' 52 0 0 0 1 0.000136216257940000 'TS-10.0001' 53 0 0 0 1 0.000136216257940000 'TS-10.0001' 54 0 0 0 1 0.000136216257940000 'TS-10.0001' 55 0 0 0 1 0.000136216257940000 'TS-10.0001' 56 0 0 0 1 0.000136216257940000 'TS-10.0001' 57 0 0 0 1 0.000136216257940000 'TS-10.0001' 58 0 0 0 1 0.000136216257940000 'TS-10.0001' 59 0 0 0 1 0.000136216257940000 'TS-10.0001' 60 0 0 0 1 0.000136216257940000 'TS-10.0001' 61 0 0 0 1 0.000136216257940000 'TS-10.0001' 62 0 0 0 1 0.000136216257940000 'TS-10.0001' 63 0 0 0 1 0.000136216257940000 'TS-10.0001' 64 0 0 0 1 0.000136216257940000 'TS-10.0001' 65 0 0 0 1 0.000136216257940000 'TS-10.0001' 66 0 0 0 1 0.000136216257940000 'TS-10.0001' 67 0 0 0 1 0.000136216257940000 'TS-10.0001' 68 0 0 0 1 0.000136216257940000 'TS-10.0001' 69 0 0 0 1 0.000136216257940000 'TS-10.0001' 70 0 0 0 1 0.000136216257940000 'TS-10.0001' 71 0 0 0 1 0.000136216257940000 'TS-10.0001' 72 0 0 0 1 0.000136216257940000 'TS-10.0001' 73 0 0 0 1 0.000136216257940000 'TS-10.0001' 74 0 0 0 1 0.000136216257940000 'TS-10.0001' 75 0 0 0 1 0.000136216257940000 'TS-10.0001' 76 0 0 0 1 0.000136216257940000 'TS-10.0001' 77 0 0 0 1 0.000136216257940000 'TS-10.0001' 78 0 0 0 1 0.000136216257940000 'TS-10.0001' 79 0 0 0 1 0.000136216257940000 'TS-10.0001' 80 0 0 0 1 0.000136216257940000 'TS-10.0001' 81 0 0 0 1 0.000136216257940000 'TS-10.0001' 82 0 0 0 1 0.000136216257940000 'TS-10.0001' 83 0 0 0 1 0.000136216257940000 'TS-10.0001' 84 0 0 0 1 0.000136216257940000 'TS-10.0001' 85 0 0 0 1 0.000136216257940000 'TS-10.0001' 86 0 0 0 1 0.000136216257940000 'TS-10.0001' 87 0 0 0 1 0.000136216257940000 'TS-10.0001' 88 0 0 0 1 0.000136216257940000 'TS-10.0001' 89 0 0 0 1 0.000136216257940000 'TS-10.0001' 90 0 0 0 1 0.000136216257940000 'TS-10.0001' 91 0 0 0 1 0.000136216257940000 'TS-10.0001' 92 0 0 0 1 0.000136216257940000 'TS-10.0001' 93 0 0 0 1 0.000136216257940000 'TS-10.0001' 94 0 0 0 1 0.000136216257940000 'TS-10.0001' 95 0 0 0 1 0.000136216257940000 'TS-10.0001' 96 0 0 0 1 0.000136216257940000 'TS-10.0001' 97 0 0 0 1 0.000136216257940000 'TS-10.0001' 98 0 0 0 1 0.000136216257940000 'TS-10.0001' 99 0 0 0 1 0.000136216257940000 'TS-10.0001' 100 0 0 0 1 0.000136216257940000 'TS-10.0001' 101 0 0 0 1 0.000136216257940000 'TS-10.0001' 0 1 1 0 3 0.000135000000000000 [] 1 1 1 0 3 0.000135000000000000 [] 2 1 1 0 3 0.000135000000000000 [] 3 1 1 0 3 0.000135000000000000 [] 4 1 1 0 3 0.000135000000000000 [] 5 1 1 0 3 0.000135000000000000 [] 6 1 1 0 3 0.000135000000000000 [] 7 1 1 0 3 0.000135000000000000 [] 8 1 1 0 3 0.000135000000000000 [] 9 1 1 0 3 0.000135000000000000 [] 10 1 1 0 3 0.000135000000000000 [] 11 1 1 0 3 0.000135000000000000 [] 12 1 1 0 3 0.000135000000000000 [] 13 1 1 0 3 0.000135000000000000 [] 14 1 1 0 3 0.000135000000000000 [] 15 1 1 0 3 0.000135000000000000 [] 16 1 1 0 3 0.000135000000000000 [] 17 1 1 0 3 0.000135000000000000 [] 18 1 1 0 3 0.000135000000000000 [] 19 1 1 0 3 0.000135030663140000 [] 20 1 1 0 3 0.000135082826440000 [] 21 1 1 0 3 0.000135149803990000 [] 22 1 1 0 3 0.000135225230360000 [] 23 1 1 0 3 0.000135302388660000 [] 24 1 1 0 3 0.000135376382250000 [] 25 1 1 0 3 0.000135442263710000 [] 26 1 1 0 3 0.000135495766350000 [] 27 1 1 0 3 0.000135534191330000 [] 28 1 1 0 3 0.000135552444090000 [] 29 1 1 0 3 0.000135552444090000 [] 30 1 1 0 3 0.000135577763110000 [] 31 1 1 0 3 0.000135577763110000 [] 32 1 1 0 3 0.000135598330970000 [] 33 1 1 0 3 0.000135598330970000 [] 34 1 1 0 3 0.000135605089520000 [] 35 1 1 0 3 0.000135605089520000 [] 36 1 1 0 3 0.000135605089520000 [] 37 1 1 0 3 0.000135605247470000 [] 38 1 1 0 3 0.000135605247470000 [] 39 1 1 0 3 0.000135605247470000 [] 40 1 1 0 3 0.000135605247470000 [] 41 1 1 0 3 0.000135605247470000 [] 42 1 1 0 3 0.000135605467760000 [] 43 1 1 0 3 0.000135605467760000 [] 44 1 1 0 3 0.000135605467760000 [] 45 1 1 0 3 0.000135605467760000 [] 46 1 1 0 3 0.000135605467760000 [] 47 1 1 0 3 0.000135605467760000 [] 48 1 1 0 3 0.000135605467760000 [] 49 1 1 0 3 0.000135605467760000 [] 50 1 1 0 3 0.000135605467760000 [] 51 1 1 0 3 0.000135605467760000 [] 52 1 1 0 3 0.000135605467760000 [] 53 1 1 0 3 0.000135607451560000 [] 54 1 1 0 3 0.000135607451560000 [] 55 1 1 0 3 0.000135607451560000 [] 56 1 1 0 3 0.000135607451560000 [] 57 1 1 0 3 0.000135607451560000 [] 58 1 1 0 3 0.000135607451560000 [] 59 1 1 0 3 0.000135607451560000 [] 60 1 1 0 3 0.000135607451560000 [] 61 1 1 0 3 0.000135607451560000 [] 62 1 1 0 3 0.000135607451560000 [] 63 1 1 0 3 0.000135607451560000 [] 64 1 1 0 3 0.000135607451560000 [] 65 1 1 0 3 0.000135607451560000 [] 66 1 1 0 3 0.000135607451560000 [] 67 1 1 0 3 0.000135607451560000 [] 68 1 1 0 3 0.000135607451560000 [] 69 1 1 0 3 0.000135607451560000 [] 70 1 1 0 3 0.000135607451560000 [] 71 1 1 0 3 0.000135607451560000 [] 72 1 1 0 3 0.000135607451560000 [] 73 1 1 0 3 0.000135607451560000 [] 74 1 1 0 3 0.000135607451560000 [] 75 1 1 0 3 0.000135607451560000 [] 76 1 1 0 3 0.000135607451560000 [] 77 1 1 0 3 0.000135607451560000 [] 78 1 1 0 3 0.000135607451560000 [] 79 1 1 0 3 0.000135607451560000 [] 80 1 1 0 3 0.000135607451560000 [] 81 1 1 0 3 0.000135607451560000 [] 82 1 1 0 3 0.000135607451560000 [] 83 1 1 0 3 0.000135607451560000 [] 84 1 1 0 3 0.000135607451560000 [] 85 1 1 0 3 0.000135607451560000 [] 86 1 1 0 3 0.000135607451560000 [] 87 1 1 0 3 0.000135607451560000 [] 88 1 1 0 3 0.000135607451560000 [] 89 1 1 0 3 0.000135607451560000 [] 90 1 1 0 3 0.000135607451560000 [] 91 1 1 0 3 0.000135607451560000 [] 92 1 1 0 3 0.000135607451560000 [] 93 1 1 0 3 0.000135607451560000 [] 94 1 1 0 3 0.000135607451560000 [] 95 1 1 0 3 0.000135607451560000 [] 96 1 1 0 3 0.000135607451560000 [] 97 1 1 0 3 0.000135607451560000 [] 98 1 1 0 3 0.000135607451560000 [] 99 1 1 0 3 0.000135607451560000 [] 100 1 1 0 3 0.000135607451560000 [] 101 1 1 0 3 0.000135607451560000 []