Crashing with WATERTRANSPORTMODEL_SNOW = NIED
This issue was reported by Diogo Costa.
Hi,
I am setting up a model for Svalbard with WATERTRANSPORTMODEL_SNOW = NEID, but it's crashing, and the crashing doesn't always happen in the same place. Everything works well when WATERTRANSPORTMODEL_SNOW = BUCKET. I am attaching the *ini file. Are there any key points to pay attention to when using NEID that I should be aware of?
Below is the console output with error info.
Thanks, Digoo
[i] [2020-10-08T11:51:36] ---> Start SNOWPACK in RESEARCH mode [i] [] /media/dcosta/data/megasync/ec_main/models/pulse/code/code_snowpack/snowpack_pulse/bin/snowpack compiled on Oct 6 2020 at 11:51:54 [i] [] Experiment : Precip_Neum [i] [] Output dir : ./bin/output
[i] [] Run on meteo station svalbard
[i] [] Reading snow cover data for station svalbard
[i] [] Finished initializing station svalbard
[i] [] No *.haz file(s) to erase
[i] [] No *.sno file(s) to erase
[i] [] No file ./bin/output/svalbard_Precip_Neum.ini to erase
[i] [2020-10-08T11:51:37] ---> Start simulation for svalbard on 1995-07-31T23:45:00+01:00
[i] [] End date specified by user: 2017-12-31T23:00:00+01:00
[i] [] Integration step length: 15.000000 min
[i] [1995-08-16T00:00:00] ---> Station svalbard (1 slope(s)): advanced to 16.08.1995 00:00:00 (2449945.500000) station time
[i] [1995-08-31T00:00:00] ---> Station svalbard (1 slope(s)): advanced to 31.08.1995 00:00:00 (2449960.500000) station time
[i] [1995-09-15T00:00:00] ---> Station svalbard (1 slope(s)): advanced to 15.09.1995 00:00:00 (2449975.500000) station time
[i] [1995-09-30T00:00:00] ---> Station svalbard (1 slope(s)): advanced to 30.09.1995 00:00:00 (2449990.500000) station time
[i] [1995-10-15T00:00:00] ---> Station svalbard (1 slope(s)): advanced to 15.10.1995 00:00:00 (2450005.500000) station time
[i] [1995-10-30T00:00:00] ---> Station svalbard (1 slope(s)): advanced to 30.10.1995 00:00:00 (2450020.500000) station time
[i] [1995-11-14T00:00:00] ---> Station svalbard (1 slope(s)): advanced to 14.11.1995 00:00:00 (2450035.500000) station time
[i] [1995-11-29T00:00:00] ---> Station svalbard (1 slope(s)): advanced to 29.11.1995 00:00:00 (2450050.500000) station time
[i] [1995-12-14T00:00:00] ---> Station svalbard (1 slope(s)): advanced to 14.12.1995 00:00:00 (2450065.500000) station time
[i] [1995-12-29T00:00:00] ---> Station svalbard (1 slope(s)): advanced to 29.12.1995 00:00:00 (2450080.500000) station time
[i] [1996-01-13T00:00:00] ---> Station svalbard (1 slope(s)): advanced to 13.01.1996 00:00:00 (2450095.500000) station time
[i] [1996-01-28T00:00:00] ---> Station svalbard (1 slope(s)): advanced to 28.01.1996 00:00:00 (2450110.500000) station time
[i] [1996-02-12T00:00:00] ---> Station svalbard (1 slope(s)): advanced to 12.02.1996 00:00:00 (2450125.500000) station time
[i] [1996-02-27T00:00:00] ---> Station svalbard (1 slope(s)): advanced to 27.02.1996 00:00:00 (2450140.500000) station time
[W] [2020-10-08T11:51:51] [/media/dcosta/data/megasync/ec_main/models/pulse/code/code_snowpack/snowpack_pulse/snowpack/DataClasses.cc:995] Negative AIR volumetric content: -0.0351
[I] [1996-03-07T03:45:00] [/media/dcosta/data/megasync/ec_main/models/pulse/code/code_snowpack/snowpack_pulse/snowpack/snowpackCore/PhaseChange.cc:441] Phase Change Begin: Element=80, nE=81 ICE 0.998368, Water 0.036690, Air -0.035057 Soil 0.000000
[W] [2020-10-08T11:51:51] [/media/dcosta/data/megasync/ec_main/models/pulse/code/code_snowpack/snowpack_pulse/snowpack/DataClasses.cc:995] Negative AIR volumetric content: -0.0351
[W] [2020-10-08T11:51:51] [/media/dcosta/data/megasync/ec_main/models/pulse/code/code_snowpack/snowpack_pulse/snowpack/snowpackCore/PhaseChange.cc:128] wrong volumetric content
[W] [2020-10-08T11:51:51] [/media/dcosta/data/megasync/ec_main/models/pulse/code/code_snowpack/snowpack_pulse/snowpack/DataClasses.cc:979] SUM of volumetric contents = 1.0351
[E] [1996-03-07T03:45:00] [/media/dcosta/data/megasync/ec_main/models/pulse/code/code_snowpack/snowpack_pulse/snowpack/snowpackCore/PhaseChange.cc:202] Sum theta[I,W,A,S] > 1
[i] [] Ice: 0.998368, Water: 0.036690, Air: 0.000000 Soil: 0.000000
[i] [] in compSubSurfaceMelt at element 80 of 81
[E] [1996-03-07T03:45:00] [/media/dcosta/data/megasync/ec_main/models/pulse/code/code_snowpack/snowpack_pulse/snowpack/snowpackCore/Snowpack.cc:1865] Snowpack computation not completed
terminate called after throwing an instance of 'mio::IOException'
what():
**** backtrace ****
(1) in PhaseChange::compPhaseChange(SnowStation&, mio::Date const&, bool const&) from /media/dcosta/data/megasync/ec_main/models/pulse/code/code_snowpack/snowpack_pulse/lib/libsnowpack.so.3
(2) in Snowpack::runSnowpackModel(CurrentMeteo, SnowStation&, double&, BoundCond&, SurfaceFluxes&) from /media/dcosta/data/megasync/ec_main/models/pulse/code/code_snowpack/snowpack_pulse/lib/libsnowpack.so.3
(3) in from /media/dcosta/data/megasync/ec_main/models/pulse/code/code_snowpack/snowpack_pulse/bin/snowpack
[PhaseChange.cc:460] Run-time error in compPhaseChange()