SyntSMB( Model, Solver, dt, TransientSimulation )

Extension of Elmer in computational glaciology
Post Reply
jeffjoy
Posts: 20
Joined: 10 Nov 2020, 04:01
Antispam: Yes

SyntSMB( Model, Solver, dt, TransientSimulation )

Post 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
Post Reply