LoadFunctionLibrary.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
FUNCTION loadxdirection(Model, inindex, Tvalue) RESULT(load)

   USE defutils
   USE generalmodule
 
   IMPLICIT None

   TYPE(model_t) :: model
   INTEGER :: inindex, funci
   REAL(KIND=dp) :: tvalue, load, diff

   model%UDFUsed = .true.

   load =  model%NodeLoadsPass(1,model%ElmerToMyNodes(inindex))
   !WRITE(*,*) load

   IF (myverbosity > 3) WRITE(6,*)'Inside [LoadFunctionLibrary]...[LoadXDirection]'

   !Calculate diff for testing purposes
   IF (tvalue < 8.001) THEN
     diff = float(model%ElmerToMyNodes(inindex)-1)*3.0d0
     diff = diff*tvalue/8.0d0
     If (dabs(load - diff) > 1e-12) model%GetTestLoads = .false.
   ELSE
     diff = float(model%ElmerToMyNodes(inindex)-1)*3.0d0
     WRITE(*,*) 'diff = ',diff
     diff = diff + (tvalue- 8.0d0)/92.0d0
     WRITE(*,*) 'diff = ',diff
     If (dabs(load - diff) > 1e-12) model%GetTestLoads = .false.
   END IF

   IF (myverbosity > 3) THEN
     WRITE(6,*) 'inindex = ', inindex
     WRITE(*,*) 'Load 1 for node:', &
       model%NodeLoadsPass(1,model%ElmerToMyNodes(inindex))
     WRITE(*,*) 'Time = ', tvalue
     WRITE(*,*) 'diff = ',diff
   END IF

END FUNCTION loadxdirection

FUNCTION loadydirection(Model, inindex, Tvalue) RESULT(load)

   USE defutils
   USE generalmodule
 
   IMPLICIT None

   TYPE(model_t) :: model
   INTEGER :: inindex, funci
   REAL(KIND=dp) :: tvalue, load, diff

   ! M.S. Added
   !WRITE (*,*) 'THIS COMMENT IS ADDED BY MASOUD'
   ! M.S. Ended

   model%UDFUsed = .true.

   load = model%NodeLoadsPass(2,model%ElmerToMyNodes(inindex)) 

   IF (myverbosity > 3) WRITE(6,*)'Inside [LoadFunctionLibrary]...[LoadYDirection]'

   !Calculate diff for testing purposes
   IF (tvalue < 8.001) THEN
     diff = float(model%ElmerToMyNodes(inindex)-1)*3.0d0 + 1.0d0
     diff = diff*tvalue/8.0d0
     If (dabs(load - diff) > 1e-12) model%GetTestLoads = .false.
   ELSE
     diff = float(model%ElmerToMyNodes(inindex)-1)*3.0d0 + 1.0d0
     WRITE(*,*) 'diff = ',diff
     diff = diff + (tvalue-8.0d0)/92.0d0
     WRITE(*,*) 'diff = ',diff
     If (dabs(load - diff) > 1e-12) model%GetTestLoads = .false.
   END IF

   IF (myverbosity > 3) THEN
     WRITE(6,*) 'inindex = ', inindex
     WRITE(*,*) 'Load 2 for node:', &
       model%NodeLoadsPass(2,model%ElmerToMyNodes(inindex))
     WRITE(*,*) 'Time = ', tvalue
     WRITE(*,*) 'diff = ',diff
   END IF


END FUNCTION loadydirection

FUNCTION loadzdirection(Model, inindex, Tvalue) RESULT(load)

   USE defutils
   USE generalmodule
 
   IMPLICIT None

   TYPE(model_t) :: model
   INTEGER :: inindex, funci
   REAL(KIND=dp) :: tvalue, load, diff

   model%UDFUsed = .true.

   load =  model%NodeLoadsPass(3,model%ElmerToMyNodes(inindex))

   IF (myverbosity > 3) WRITE(6,*)'Inside [LoadFunctionLibrary]...[LoadZDirection]'

   !Calculate diff for testing purposes
   IF (tvalue < 8.001) THEN
     diff = float(model%ElmerToMyNodes(inindex)-1)*3.0d0 + 2.0d0
     diff = diff*tvalue/8.0d0
     If (dabs(load - diff) > 1e-12) model%GetTestLoads = .false.
   ELSE
     diff = float(model%ElmerToMyNodes(inindex)-1)*3.0d0 + 2.0d0
     WRITE(*,*) 'diff = ',diff
     diff = diff + (tvalue- 8.0d0)/92.0d0
     WRITE(*,*) 'diff = ',diff
     If (dabs(load - diff) > 1e-12) model%GetTestLoads = .false.
   END IF

   IF (myverbosity > 3) THEN
     WRITE(6,*) 'inindex = ', inindex
     WRITE(*,*) 'Load 3 for node:', &
       model%NodeLoadsPass(3,model%ElmerToMyNodes(inindex))
     WRITE(*,*) 'Time = ', tvalue
     WRITE(*,*) 'diff = ',diff
   END IF

END FUNCTION loadzdirection