Standard Output
+++ Running case: Example-Greenland
+++ working dir: /Users/jenkins/workspace/macOS-Intel-Examples/nightlylog
| < M A T L A B (R) >
| Copyright 1984-2023 The MathWorks, Inc.
| R2023b Update 6 (23.2.0.2485118) 64-bit (maci64)
| December 28, 2023
|
| To get started, type doc.
| For product information, visit www.mathworks.com.
|
| ISSM development path correctly loaded
| Step 1: Mesh creation
| Anisotropic mesh adaptation
| WARNING: mesh present but no geometry found. Reconstructing...
| new number of triangles = 6345
| Step 2: Parameterization
| Loading SeaRISE data from NetCDF
| Interpolating surface and bedrock
| Constructing thickness
| Interpolating velocities
| Interpolating temperatures
| Interpolating surface mass balance
| Construct basal friction parameters
| Construct ice rheological properties
| [Warning: paterson is outdated, please consider using cuffey instead]
| [> In paterson (line 10)
| In TemporaryParameterFile57804 (line 54)
| In parameterize (line 29)
| In runme (line 41)]
| Set other boundary conditions
| Set geothermal heat flux
| Set Pressure
| Single point constraints
| Step 3: Control method friction
| checking model consistency
| INFO: the outlog will look better if only md.verbose.control is turned on
| marshalling file /Users/jenkins/workspace/macOS-Intel-Examples//execution/SeaRISEgreenland-06-30-2026-05-17-52-57804/SeaRISEgreenland.bin
| launching solution sequence
| ───────────────────────────────────────────────────────────────────
| Ice-sheet and Sea-level System Model (ISSM) version 2026.2
| GitHub: https://github.com/ISSMteam/ISSM/
| Documentation: https://issmteam.github.io/ISSM-Documentation/
| ───────────────────────────────────────────────────────────────────
| call computational core:
| preparing initial solution
| x | Cost function f(x) | List of contributions
| ====================== step 1/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 59172.79 | 7092.887 52079.9 6.085212e-32
| x = 1 | f(x) = 50954.87 | 3600.244 47354.58 0.04587451
| ====================== step 2/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 50953.45 | 3598.809 47354.6 0.04587451
| x = 1 | f(x) = 44978.74 | 2137.679 42840.82 0.2370654
| ====================== step 3/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 44987.02 | 2146.197 42840.58 0.2370654
| x = 1 | f(x) = 44425.8 | 2044.873 42380.66 0.2722667
| ====================== step 4/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 44432.36 | 2051.851 42380.24 0.2722667
| x = 1 | f(x) = 42775.97 | 1747.394 41028.16 0.4148691
| ====================== step 5/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 42783.98 | 1756.726 41026.84 0.4148691
| x = 1 | f(x) = 41876.1 | 1659.336 40216.23 0.5302616
| ====================== step 6/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 41885.98 | 1670.347 40215.1 0.5302616
| x = 1 | f(x) = 40869.39 | 1561.173 39307.52 0.6932322
| ====================== step 7/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 40870.91 | 1563.273 39306.95 0.6932322
| x = 1 | f(x) = 39697.76 | 1338.256 38358.58 0.9236775
| ====================== step 8/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 39697.43 | 1338.655 38357.85 0.9236775
| x = 1 | f(x) = 39040.22 | 1208.675 37830.47 1.071791
| ====================== step 9/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 39039.49 | 1208.558 37829.86 1.071791
| x = 1 | f(x) = 38556.09 | 1126.764 37428.13 1.191478
| ====================== step 10/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 38552.98 | 1124.507 37427.28 1.191478
| x = 1 | f(x) = 38121.56 | 1092.495 37027.74 1.3288
| ====================== step 11/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 38122.77 | 1094.958 37026.49 1.3288
| x = 1 | f(x) = 37729.41 | 1149.628 36578.29 1.493932
| ====================== step 12/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 37740.46 | 1159.142 36579.82 1.493932
| x = 1 | f(x) = 37487.92 | 1087.186 36399.14 1.592633
| x = 0.381966 | f(x) = 37611.19 | 1112.303 36497.36 1.529775
| x = 0.618034 | f(x) = 37561.1 | 1101.089 36458.46 1.55338
| x = 0.763932 | f(x) = 37530.94 | 1093.798 36435.58 1.568476
| x = 0.854102 | f(x) = 37512.65 | 1089.499 36421.58 1.577683
| x = 0.90983 | f(x) = 37501.45 | 1086.979 36412.88 1.583344
| ====================== step 13/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 37484.86 | 1084.592 36398.67 1.592633
| x = 1 | f(x) = 37022.46 | 1018.285 36002.41 1.768007
| ====================== step 14/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 37030.9 | 1026.085 36003.05 1.768007
| x = 1 | f(x) = 36814.89 | 1058.135 35754.9 1.859438
| x = 0.381966 | f(x) = 36898.13 | 989.5907 35906.74 1.799676
| x = 0.618034 | f(x) = 36853.98 | 1003.884 35848.28 1.821225
| x = 0.763932 | f(x) = 36838.62 | 1025.038 35811.75 1.835325
| x = 0.937069 | f(x) = 36824.08 | 1052.478 35769.75 1.852853
| x = 0.870937 | f(x) = 36839.21 | 1051.966 35785.4 1.846056
| ====================== step 15/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 36824.8 | 1068.389 35754.55 1.859438
| x = 1 | f(x) = 36559.07 | 954.0377 35603.09 1.936544
| x = 0.381966 | f(x) = 36711.5 | 1017.403 35692.21 1.887487
| x = 0.618034 | f(x) = 36653.71 | 995.7569 35656.05 1.905392
| x = 0.763932 | f(x) = 36614.76 | 978.0481 35634.8 1.916939
| x = 0.854102 | f(x) = 36589.85 | 965.9689 35621.95 1.92432
| x = 0.90983 | f(x) = 36574.06 | 958.1749 35613.95 1.928974
| ====================== step 16/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 36552.27 | 948.6767 35601.65 1.936544
| x = 1 | f(x) = 36368.74 | 909.4905 35457.24 2.009954
| x = 0.381966 | f(x) = 36450.66 | 912.2196 35536.48 1.964688
| x = 0.618034 | f(x) = 36415.28 | 905.8741 35507.42 1.981913
| x = 0.763932 | f(x) = 36394.68 | 904.2792 35488.41 1.992565
| x = 0.854102 | f(x) = 36384.31 | 906.0028 35476.31 1.999281
| x = 0.90983 | f(x) = 36379.73 | 908.9319 35468.79 2.003357
| ====================== step 17/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 36373.34 | 914.7542 35456.58 2.009954
| x = 1 | f(x) = 36229.12 | 957.6674 35269.34 2.111313
| x = 0.381966 | f(x) = 36289.06 | 905.1848 35381.82 2.046875
| x = 0.618034 | f(x) = 36257.86 | 916.5034 35339.28 2.070593
| x = 0.763932 | f(x) = 36247.49 | 932.6872 35312.72 2.085785
| x = 0.913021 | f(x) = 36239.38 | 951.6387 35285.64 2.101794
| x = 0.856074 | f(x) = 36249.62 | 951.977 35295.55 2.095621
| ====================== step 18/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 36238.23 | 966.1893 35269.93 2.111313
| x = 1 | f(x) = 36005.81 | 893.306 35110.31 2.18899
| x = 0.381966 | f(x) = 36142.54 | 929.6447 35210.76 2.139784
| x = 0.618034 | f(x) = 36093.23 | 916.1558 35174.92 2.158128
| x = 0.763932 | f(x) = 36058.58 | 905.5225 35150.89 2.16975
| x = 0.854102 | f(x) = 36036.9 | 898.7214 35136 2.177043
| x = 0.90983 | f(x) = 36023.73 | 894.6151 35126.93 2.18159
| ====================== step 19/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 36005.22 | 890.2107 35112.82 2.18899
| x = 1 | f(x) = 35868.42 | 869.9581 34996.18 2.287732
| x = 0.381966 | f(x) = 35911.39 | 860.1821 35048.98 2.225541
| x = 0.618034 | f(x) = 35893.04 | 858.0691 35032.72 2.249067
| x = 0.763932 | f(x) = 35881.94 | 859.7353 35019.94 2.263658
| x = 0.854102 | f(x) = 35877.04 | 863.763 35011.01 2.272799
| x = 0.90983 | f(x) = 35875.47 | 868.0771 35005.11 2.278502
| ====================== step 20/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 35873.35 | 875.6681 34995.4 2.287732
| x = 1 | f(x) = 35723.35 | 910.5513 34810.43 2.374027
| x = 0.381966 | f(x) = 35789.66 | 860.5594 34926.79 2.319061
| x = 0.618034 | f(x) = 35757.04 | 868.789 34885.91 2.339232
| x = 0.763932 | f(x) = 35744.03 | 883.7317 34857.95 2.352231
| x = 0.854102 | f(x) = 35739.94 | 897.5494 34840.03 2.360435
| x = 0.930997 | f(x) = 35737.52 | 910.28 34824.87 2.367548
| ====================== step 21/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 35735.84 | 922.2347 34811.23 2.374027
| x = 1 | f(x) = 35602.51 | 870.7784 34729.3 2.423398
| x = 0.381966 | f(x) = 35668.82 | 880.0971 34786.33 2.391628
| x = 0.618034 | f(x) = 35644.9 | 872.7159 34769.79 2.403223
| x = 0.763932 | f(x) = 35629.08 | 870.0702 34756.6 2.410606
| x = 0.854102 | f(x) = 35618.25 | 868.4465 34747.39 2.415368
| x = 0.90983 | f(x) = 35610.89 | 867.4925 34740.98 2.418388
| ====================== step 22/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 35599.85 | 867.5045 34729.92 2.423398
| x = 1 | f(x) = 35435.33 | 841.7746 34591.04 2.515761
| x = 0.381966 | f(x) = 35503.09 | 830.83 34669.8 2.457304
| x = 0.618034 | f(x) = 35463.09 | 825.8312 34634.77 2.479098
| x = 0.763932 | f(x) = 35440.81 | 828.0732 34610.24 2.492898
| x = 0.854102 | f(x) = 35432.56 | 834.5397 34595.51 2.501553
| x = 0.90983 | f(x) = 35438.35 | 841.9242 34593.91 2.50695
| ====================== step 23/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 35439.26 | 840.4865 34596.27 2.501553
| x = 1 | f(x) = 35377.76 | 827.9225 34547.29 2.55107
| x = 0.381966 | f(x) = 35416.76 | 819.8319 34594.41 2.518729
| x = 0.618034 | f(x) = 35399.79 | 819.4099 34577.85 2.530419
| x = 0.763932 | f(x) = 35393.8 | 822.2266 34569.04 2.538053
| x = 0.944313 | f(x) = 35386.99 | 828.3334 34556.1 2.547926
| x = 0.875413 | f(x) = 35393.31 | 829.0398 34561.73 2.544098
| ====================== step 24/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 35389.05 | 833.3283 34553.17 2.55107
| x = 1 | f(x) = 35252.08 | 814.3328 34435.09 2.649675
| x = 0.381966 | f(x) = 35317.08 | 812.978 34501.52 2.588143
| x = 0.618034 | f(x) = 35288.37 | 809.7173 34476.04 2.611631
| x = 0.763932 | f(x) = 35271.33 | 809.3561 34459.34 2.626304
| x = 0.854102 | f(x) = 35263.63 | 811.1878 34449.81 2.635229
| x = 0.90983 | f(x) = 35260.58 | 813.8445 34444.09 2.640723
| ====================== step 25/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 35256.61 | 818.8748 34435.09 2.649675
| x = 1 | f(x) = 35162.53 | 858.2923 34301.5 2.735676
| x = 0.381966 | f(x) = 35180.62 | 812.9346 34365.01 2.680864
| x = 0.618034 | f(x) = 35153.13 | 821.9691 34328.46 2.701321
| x = 0.763932 | f(x) = 35158.4 | 837.054 34318.64 2.714146
| x = 0.645763 | f(x) = 35167.63 | 833.394 34331.53 2.703727
| x = 0.527864 | f(x) = 35170.22 | 824.2448 34343.28 2.693583
| ====================== step 26/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 35160.91 | 827.091 34331.12 2.701321
| x = 1 | f(x) = 35121.46 | 799.6831 34319.05 2.727669
| x = 0.381966 | f(x) = 35148.31 | 811.7235 34333.88 2.710518
| x = 0.618034 | f(x) = 35139.87 | 807.5893 34329.56 2.716727
| x = 0.763932 | f(x) = 35133.07 | 803.7521 34326.59 2.720776
| x = 0.854102 | f(x) = 35127.96 | 801.0169 34324.22 2.723359
| x = 0.90983 | f(x) = 35124.32 | 799.2107 34322.38 2.724986
| ====================== step 27/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 35118.96 | 797.3025 34318.93 2.727669
| x = 1 | f(x) = 34992.01 | 809.2473 34179.93 2.832751
| x = 0.381966 | f(x) = 35054.02 | 790.9422 34260.31 2.766228
| x = 0.618034 | f(x) = 35029.18 | 793.4468 34232.94 2.790931
| x = 0.763932 | f(x) = 35015.1 | 799.223 34213.07 2.806635
| x = 0.854102 | f(x) = 35009.09 | 806.1057 34200.17 2.816507
| x = 0.90983 | f(x) = 35007.38 | 812.598 34191.96 2.822672
| ====================== step 28/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 35004.32 | 822.0413 34179.44 2.832751
| x = 1 | f(x) = 34934.87 | 806.5674 34125.43 2.869513
| x = 0.381966 | f(x) = 34947.17 | 786.7793 34157.54 2.846471
| x = 0.618034 | f(x) = 34937.13 | 788.8927 34145.38 2.855048
| x = 0.763932 | f(x) = 34935.7 | 795.0416 34137.8 2.860456
| x = 0.748008 | f(x) = 34937.94 | 796.546 34138.54 2.859861
| x = 0.854102 | f(x) = 34938.57 | 802.5726 34133.14 2.863862
| ====================== step 29/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 34941.3 | 812.7678 34125.66 2.869513
| x = 1 | f(x) = 34824.78 | 779.576 34042.26 2.935698
| x = 0.381966 | f(x) = 34880.03 | 788.1485 34088.99 2.894439
| x = 0.618034 | f(x) = 34858.6 | 784.1171 34071.57 2.91002
| x = 0.763932 | f(x) = 34843.88 | 780.8715 34060.09 2.919804
| x = 0.854102 | f(x) = 34835.32 | 779.3894 34053.01 2.925896
| x = 0.90983 | f(x) = 34830.74 | 779.0445 34048.76 2.929623
| ====================== step 30/30 ===============================
| x = 0 | computing velocities
| computing adjoint
| saving results
| f(x) = 34824.9 | 779.6846 34042.28 2.935698
| x = 1 | f(x) = 34721.73 | 773.8226 33944.9 3.012811
| x = 0.381966 | f(x) = 34777.78 | 771.9937 34002.82 2.964712
| x = 0.618034 | f(x) = 34756.76 | 772.5218 33981.26 2.983082
| x = 0.763932 | f(x) = 34743.45 | 772.8096 33967.65 2.994363
| x = 0.854102 | f(x) = 34735.33 | 773.2066 33959.12 3.001402
| x = 0.90983 | f(x) = 34730.47 | 773.6083 33953.86 3.005763
| preparing final solution
| computing new velocity
| write lock file:
| FemModel initialization elapsed time: 0.134429
| Total Core solution elapsed time: 97.4304
| Linear solver elapsed time: 58.3618 (60%)
| Total elapsed time: 0 hrs 1 min 37 sec
| loading results from cluster
| Step 4: Transient run
| checking model consistency
| marshalling file /Users/jenkins/workspace/macOS-Intel-Examples//execution/SeaRISEgreenland-06-30-2026-05-19-39-57804/SeaRISEgreenland.bin
| launching solution sequence
| ───────────────────────────────────────────────────────────────────
| Ice-sheet and Sea-level System Model (ISSM) version 2026.2
| GitHub: https://github.com/ISSMteam/ISSM/
| Documentation: https://issmteam.github.io/ISSM-Documentation/
| ───────────────────────────────────────────────────────────────────
| Input updates from constant
| Input updates from constant
| Renumbering degrees of freedom
| Renumbering degrees of freedom
| Renumbering degrees of freedom
| Renumbering degrees of freedom
| Renumbering degrees of freedom
| Renumbering degrees of freedom
| call computational core:
| Updating Mmes
| iteration 1/100 time [yr]: 0.20 (time step: 0.20)
| computing smb
| computing new velocity
| Updating constraints and active domain of analysis StressbalanceAnalysis for time: 0.2
| Get solution from inputs
| Reduce vector from g to f set
| Input updates from constant
| Updating inputs from solution for StressbalanceAnalysis
| Allocating matrices (Kff stiffness matrix size: 7398 x 7398)
| Assembling matrices
| Create nodal constraints
| Dirichlet lifting applied to load vector
| Solving matrix system
| solver residue: norm(KU-F)/norm(F)=2.05155e-16
| Merging solution vector from fset to gset
| checking convergence
| mechanical equilibrium convergence criterion 41.3709 > 1 %
| Convergence criterion: norm(du)/norm(u) 63.1692 > 10 %
| Convergence criterion: max(du) 0.000157031
| Input updates from constant
| Updating inputs from solution for StressbalanceAnalysis
| number of unstable constraints: 0
| Allocating matrices (Kff stiffness matrix size: 7398 x 7398)
| Assembling matrices
| Create nodal constraints
| Dirichlet lifting applied to load vector
| Solving matrix system
| solver residue: norm(KU-F)/norm(F)=2.08287e-16
| Merging solution vector from fset to gset
| checking convergence
| mechanical equilibrium convergence criterion 3.67351 > 1 %
| Convergence criterion: norm(du)/norm(u) 7.15747 < 10 %
| Convergence criterion: max(du) 8.86618e-06
| Input updates from constant
| Updating inputs from solution for StressbalanceAnalysis
| number of unstable constraints: 0
| Allocating matrices (Kff stiffness matrix size: 7398 x 7398)
| Assembling matrices
| Create nodal constraints
| Dirichlet lifting applied to load vector
| Solving matrix system
| solver residue: norm(KU-F)/norm(F)=2.21945e-16
| Merging solution vector from fset to gset
| checking convergence
| mechanical equilibrium convergence criterion 1.10917 > 1 %
| Convergence criterion: norm(du)/norm(u) 2.86495 < 10 %
| Convergence criterion: max(du) 5.14918e-06
| Input updates from constant
| Updating inputs from solution for StressbalanceAnalysis
| number of unstable constraints: 0
| Allocating matrices (Kff stiffness matrix size: 7398 x 7398)
| Assembling matrices
| Create nodal constraints
| Dirichlet lifting applied to load vector
| Solving matrix system
| solver residue: norm(KU-F)/norm(F)=2.01372e-16
| Merging solution vector from fset to gset
| checking convergence
| mechanical equilibrium convergence criterion 0.445157 < 1 %
| Convergence criterion: norm(du)/norm(u) 1.23492 < 10 %
| Convergence criterion: max(du) 2.12087e-06
| Input updates from constant
| Updating inputs from solution for StressbalanceAnalysis
| number of unstable constraints: 0
| total number of iterations: 4
| computing basal mass balance
| computing mass transport
| Updating constraints and active domain of analysis MasstransportAnalysis for time: 0.2
| Allocating matrices
| [1] ??? Error using ==> ./analyses/MasstransportAnalysis.cpp:720
| [1] CreatePVectorCG error message: Assertion "connectedtoocean_input" failed, please report bug at https://github.com/ISSMteam/ISSM/
| (Kff stiffness matrix size: 3699 x 3699)
| Assembling matrices
| [0] ??? Error using ==> ./analyses/MasstransportAnalysis.cpp:720
| [0] CreatePVectorCG error message: Assertion "connectedtoocean_input" failed, please report bug at https://github.com/ISSMteam/ISSM/
| loading results from cluster
| [Warning: Could not copy
| /Users/jenkins/workspace/macOS-Intel-Examples//execution/SeaRISEgreenland-06-30-2026-05-19-39-57804//SeaRISEgreenland.outbin]
| [> In issmscpin (line 22)
| In generic/Download (line 328)
| In loadresultsfromcluster (line 45)
| In solve (line 210)
| In runme (line 111)]
| Error using loadresultsfromdisk
| =========================================================================
| Binary file SeaRISEgreenland.outbin not found
| This typically results from an error encountered during the simulation
| Please check for error messages above or in the outlog
| =========================================================================
| Error in loadresultsfromcluster (line 48)
| md=loadresultsfromdisk(md,d.miscellaneous.name '.outbin']);
| Error in solve (line 210)
| md=loadresultsfromcluster(md);
| Error in runme (line 111)
| md=solve(md,'Transient');FAILURE
+++ exit code: 0