SyntSMB( Model, Solver, dt, TransientSimulation )
Posted: 09 Dec 2021, 17:44
Hello,in the case MLB,i don't know what's the meaning,did the ELA is glacier is 400 amsl?or 400 is the minimum of the bedDEM?Who knows?and what's the meaning of the 1000?
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