!***********************************************************************!
! TITLE: CFVR3903 - Testing fluxes 4 node 2D laplace element
! SUBTITLE:
!
! DESCRIPTION: 	Test fluxes in nodes of a PLANE42 2D model. The model consist of
! DESCRIPTION:  both soil triangle and square.
! DESCRIPTION: 	
! DESCRIPTION:  For this purpose it calculates and applies a ground pressure
! DESCRIPTION:  on both plane elements.
!
! ELEMENT TYPE: PLANE42
! MODULES:
! UNITS: SI
! KEYWORD1: Material Models
! KEYWORD2: Materials
!
!***********************************************************************!
  FINISH
  ~CFCLEAR,,1
  NomFile='CFVR3903'
  /TITLE, %NomFile%, Testing fluxes 4 node 2D laplace element

  ~SEEPAGE,2D

/prep7

~cfmp,1,lib,soil,,sc
~cfmp,1,soil,kx,,10
~cfmp,1,soil,ky,,10

!      X   Y   Z
n,10, 10, 10,  0
n,20, 20, 10,  0
n,30, 20, 20,  0
n,40, 10, 20,  0
n,50, 50, 10,  0
n,60, 60, 10,  0
n,70, 55, 20,  0

e,10,20,30,40
e,50,60,70,70

/SOLU
d,10,head,100
d,20,head,100
d,30,head,10
d,40,head,10
d,70,head,100
d,60,head,10
d,50,head,10

~LPSOLVE
/post1
file,file,seep
set,1
/dscale,,off
plnsol,head

etable,fx1,smisc,1
etable,fx2,smisc,2
etable,fx3,smisc,3
etable,fx4,smisc,4

!--------------------------------------------------------------------------------------
! DATA CHECK
!--------------------------------------------------------------------------------------
! Data comparison number
  NComp    =  8
  NComp_ch =  0

! Marix dim.
  *DIM,LABEL,CHAR,Ncomp,1
  *DIM,LABEL_CH,CHAR,Ncomp_ch,1
  *DIM,VALUE,,Ncomp,3
  *DIM,VALUE_CH,CHAR,Ncomp_ch,3
  *DIM,TOLER,,Ncomp,2

! Labels
  LABEL(1,1) = 'E1F1'
  LABEL(2,1) = 'E1F2'
  LABEL(3,1) = 'E1F3'
  LABEL(4,1) = 'E1F4'
  LABEL(5,1) = 'E2F1'
  LABEL(6,1) = 'E2F2'
  LABEL(7,1) = 'E2F3'
  LABEL(8,1) = 'E2F4'

! Correct values
  VALUE(1,1)= -900.0
  VALUE(2,1)= 0.0
  VALUE(3,1)= 900.0
  VALUE(4,1)= 0.0
  VALUE(5,1)= 900.0
  VALUE(6,1)= -450.0
  VALUE(7,1)= 0.0
  VALUE(8,1)= -450.0

! Obtained values
  *get,VALUE(1,2),element,1,etable,fx1
  *get,VALUE(2,2),element,1,etable,fx2
  *get,VALUE(3,2),element,1,etable,fx3
  *get,VALUE(4,2),element,1,etable,fx4
  *get,VALUE(5,2),element,2,etable,fx1
  *get,VALUE(6,2),element,2,etable,fx2
  *get,VALUE(7,2),element,2,etable,fx3
  *get,VALUE(8,2),element,2,etable,fx4

! Warning and error tolerances
  *DO,II,1,8
    TOLER(II, 1)= 1E-9
    TOLER(II, 2)= 1E-9
  *ENDDO

!--------------------------------------------------------------------------------------
! Results comparison
!--------------------------------------------------------------------------------------
  COMPARA.MAC