Code: Select all
SUBROUTINE SyntSMB( Model, Solver, dt, TransientSimulation )
!------------------------------------------------------------------------------
!******************************************************************************
!
! For the Upper surface, creates and updates a mask which may be equal to -1, 0 or 1
!
! Icy Mask = + 1 if Icy (glacier)
! = - 1 if Ice Free
! = 0 on the contouur of the glacier (first node with H=Hmin)
! < -1 for Icy Isolated nodes (useful to move them back to initial elevation)
!
! ARGUMENTS:
Code: Select all
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! BODY FORCE
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Body Force 1
Flow BodyForce 1 = 0.0
Flow BodyForce 2 = 0.0
Flow BodyForce 3 = #gravity
Temp = Real 273.15
Temp Condition =Variable Coordinate 3
Real
0.0 -1.0
399.9999 -1.0
400. 1.0
1000.0 1.0
End
!Variable Coordinate 3
! Real lua "settemperate(tx[0])"
End