I am currently doing a transient simulation of a glacier with 60 timesteps of 1 year.

I defined the year of reference (1958) and extract the results at different timesteps (1959, 1969, 1980, 1990, 1999, 2017).

Code: Select all

```
Simulation
Coordinate System = Cartesian 3D
! Transient simulation
Simulation Type = Transient
Timestepping Method = "bdf"
BDF Order = 2
Timestep Intervals(6) = 1 10 11 10 9 18
Output Intervals(6) = 1 1 1 1 1 1 1
Timestep Sizes(1) = 1
Extruded Mesh Levels = Integer 10
Extruded Mesh Ratio = Real 2.0
Steady State Max Iterations = 1
Steady State Min Iterations = 1
Post File = "#year#.vtu"
Restart File = "../INIT/OUTPUT/RUN_0_.result"
Restart Position = 0
Restart Time = Real 1958
Restart Before Initial Conditions = Logical True
max output level = 4
End
```

Code: Select all

```
Solver 9
Exec Solver = After Saving
Equation = SaveGrid
Procedure = "SaveGridData" "SaveGridData"
Filename Prefix = File "Grid_#year#"
! Results will be saved only on the top boundary where MaskSurf=Logical True
Mask Name = String "MaskSurf"
Check for Duplicates = Logical True
! We don't need a z-coordinate
Suppress Extra Dimension = Logical True
Single Precision = Logical True
!Grid definition
! same grid as basal.tif
Grid Origin At Corner = Logical True
Min Coordinate 1 = Real 647830.67578125
Max Coordinate 1 = Real 652120.67578125
Min Coordinate 2 = Real 7535416.9921875
Max Coordinate 2 = Real 7537846.9921875
!Number of gridspace (Dim-1)
Grid nx = Integer 429
Grid ny = Integer 243
!Variables to be saved
Scalar Field 1 = String "surface"
Scalar Field 2 = String "bed"
Scalar Field 3 = String "thickness"
Scalar Field 4 = String "bottom Copt"
Scalar Field 5 = String "mass balance"
Vector Field 1 = String "Velocity"
Vector Field 2 = String "int Velocity"
Vector Field 3 = String "bottom Velocity"
NetCDF Format = Logical True
End
```

Thanks in advance for your help,

Olivier