Page 1 of 1

SyntSMB( Model, Solver, dt, TransientSimulation )

Posted: 09 Dec 2021, 17:44
by jeffjoy
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