!*****************************************************************************************************!
! TITLE: CFVR5112 - Test Defined Sections with holes
! SUBTITLE: Example to test bridge sections with holes
!
! DESCRIPTION: Several bridge sections are defined in this example
! DESCRIPTION:
! DESCRIPTION: The aim of the example is to verify the geometry of the sections.
!
! ELEMENT TYPE: NONE
! MODULES: NLBR
! UNITS: SI
! KEYWORD1: Bridges
!
!*************************************************************************************!
  FINISH
  ~CFCLEAR,,1
  ~CFACTIV,NLBR,Y
  NomFile='CFVR5112'
  /TITLE, %NomFile%, Test Defined Sections with Holes

! ---------------------------------------------------------------------------------
! Initial data
! ---------------------------------------------------------------------------------
/PREP7

~BRSSLAB,1,RS,1,10 ,5
~BRSMDF,1,NAME,,,Sec1
~BRSSLAB,2,TS,2,20 ,50 ,10 ,3
~BRSMDF,2,NAME,,,Sec2
~BRSSLAB,3,TF,3,100 ,120 ,200 ,20 ,8 ,4
~BRSMDF,3,NAME,,,Sec3
~BRSSLAB,4,PS,4,10 ,15 ,20 ,30 ,10 ,3 ,2 ,2
~BRSMDF,4,NAME,,,Sec4
~BRSSLAB,5,PAL,5,30 ,39 ,48 ,59 ,25 ,10 ,5 ,3
~BRSSLAB,5,PAR,5,10 ,15 ,25 ,38 ,15 ,5 ,5 ,3
~BRSMDF,5,NAME,,,Sec5
~BRHL,1,C,-2.5,2.5,1 ,,3 ,-3
~BRSMDF,1,NAME,,,Sec1
~BRSMDF,1,NAME,,,Sec1
~BRHL,4,C,-3 ,10 ,1 ,,1 ,0
~BRHL,4,C,-3 ,7 ,1 ,,1 ,0
~BRHLMDF,4,2,C,-3 ,7 ,1.2,
~BRHLMDF,4,2,C,-3 ,7 ,1.4,
~BRHLMDF,4,2,C,-3 ,7 ,1.5,
~BRHL,4,C,-3 ,3 ,2 ,,1 ,0
~BRHLMDF,4,3,C,-3 ,3 ,2.2,
~BRHL,4,C,-3 ,-3 ,2.2,,1 ,0
~BRHL,4,C,-3 ,-7 ,1.5,,1 ,0
~BRHL,4,C,-3 ,-10 ,1 ,,1 ,0
~BRSMDF,4,NAME,,,Sec4
~BRSMDF,4,NAME,,,Sec4
~BRHL,3,C,-10 ,25 ,5 ,,1 ,0
~BRHLDEL,3,1
~BRHL,3,C,-10 ,30 ,2.5,,5 ,-6
~BRSMDF,3,NAME,,,Sec3
~BRSMDF,3,NAME,,,Sec3
~BRHLMDF,3,5,C,-10 ,-6 ,3 ,
~BRHLMDF,3,5,C,-10 ,-6 ,4 ,
~BRHL,3,C,-8 ,-14 ,2 ,,2 ,-8.5
~BRHLMDF,3,7,C,-8.5,-29.5,4 ,
~BRSMDF,3,NAME,,,Sec3
~BRHL,2,C,-4 ,10 ,3 ,,4 ,-6.1
~BRSMDF,2,NAME,,,Sec2
~BRSMDF,4,NAME,,,Sec4
~BRHL,5,C,-5 ,20 ,3 ,,1 ,0
~BRHLMDF,5,1,C,-10 ,30 ,6 ,
~BRHL,5,C,-4 ,10 ,2.5,,1 ,0
~BRHLMDF,5,2,C,-8 ,10 ,4 ,
~BRHLMDF,5,2,C,-8 ,10 ,5 ,
~BRHLMDF,5,2,C,-10 ,10 ,5 ,
~BRHLMDF,5,2,C,-12 ,10 ,5 ,
~BRHL,5,C,-9 ,-2 ,4 ,,1 ,0
~BRHL,5,C,-4 ,-20 ,1 ,,1 ,0
~BRSMDF,5,NAME,,,Sec5
~BRMVDL,5,2,39
~BRBC,5,ADD,2
~BRMVDL,5,5,-15
~BRMVDL,5,6,-25
~BRMVDL,5,1,48
~BRBC,5,ADD,1
~BRBC,5,ADD,6
~BRBC,5,ADD,7
~BRSMDF,5,NAME,,,Sec5
~BRBC,1,ADD,3
~BRBC,1,ADD,4
~BRMVDL,1,2,4.4
~BRSMDF,1,NAME,,,Sec1
~BRMVDL,2,7,-18.433
~BRBC,2,ADD,7
~BRMVDL,2,1,18
~BRBC,2,ADD,1
~BRSMDF,2,NAME,,,Sec2
~BRMVDL,3,9,-45
~BRMVDL,3,2,45
~BRBC,3,ADD,2
~BRBC,3,ADD,9
~BRBC,3,ADD,6
~BRSMDF,3,NAME,,,Sec3
~BRBC,4,ADD,1
~BRBC,4,ADD,2
~BRBC,4,ADD,5
~BRBC,4,ADD,6
~BRBC,4,ADD,8
~BRBC,4,ADD,9
~BRSMDF,4,NAME,,,Sec4

!--------------------------------------------------------------------------------------
! DATA CHECK
!--------------------------------------------------------------------------------------
! Data comparison number

  NComp = 199
  NComp_ch = 0

! Matrix 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
  *DIM,Nodes,,Ncomp,1
  *DIM,COORDX,,Ncomp,1

  j = 0
  VALUE(1,1) = 3,4,7,6,4
  *do,ii,1,5
  j=j+1
  LABEL(j)  = 'NHOLES%ii%'  $ ~CFGET,value( j,2),BRIDGE,,SECTION,NHOLES,,ii
  *enddo
  VALUE(6,1) = 6,7,10,9,7
  *do,ii,1,5
  j=j+1
  LABEL(j)  = 'NDIV%ii%'  $ ~CFGET,value( j,2),BRIDGE,,SECTION,NDIV,,ii
  *enddo


VALUE(11,1)= 1.0000,1.0000,1.0000,3.0000,3.0000,3.0000,3.0000,2.5000,2.5000,2.5000
VALUE(21,1)= 2.5000,4.0000,2.0000,4.0000,1.0000,1.5000,2.2000,2.2000,1.5000,1.0000
VALUE(31,1)= 6.0000,5.0000,4.0000,1.0000


  *do,ii,1,5
  *do,JJ,1,VALUE(II,1)
  j=j+1
  LABEL(j)  = 'R%ii%H%jj%'  $ ~CFGET,value( j,2),BRIDGE,,SECTION,HOLES,RADIO,ii,jj
  *enddo
  *enddo

VALUE(35,1)=-2.5000,-2.5000,-2.5000,-4.0000,-4.0000,-4.0000,-4.0000,-10.000,-10.000
VALUE(44,1)=-10.000,-10.000,-10.000,-8.0000,-8.5000,-3.0000,-3.0000,-3.0000,-3.0000
VALUE(53,1)=-3.0000,-3.0000,-10.000,-12.000,-9.0000,-4.0000

  *do,ii,1,5
  *do,JJ,1,VALUE(II,1)
  j=j+1
  LABEL(j)  = 'CY%ii%H%jj%'  $ ~CFGET,value( j,2),BRIDGE,,SECTION,HOLES,YCENT,ii,jj
  *enddo
  *enddo

VALUE(59,1)= 2.5000,-0.500,-3.500,10.000,3.9000,-2.200,-8.300,30.000,24.000,18.000
VALUE(69,1)= 12.000,-6.000,-14.00,-29.50,10.000,7.0000,3.0000,-3.000,-7.000,-10.00
VALUE(79,1)= 30.000,10.000,-2.000,-20.00

  *do,ii,1,5
  *do,JJ,1,VALUE(II,1)
  j=j+1
  LABEL(j)  = 'CZ%ii%H%jj%'  $ ~CFGET,value( j,2),BRIDGE,,SECTION,HOLES,ZCENT,ii,jj
  *enddo
  *enddo

VALUE(83,1)  = 6,8,10,12,14,16,8,10,12,14,16,18,20,10,12,14,16,18
VALUE(101,1) =20,22,24,26,28,12,14,16,18,20,22,24,26,28,11,12,14
VALUE(118,1) =16,17,18,20

  *do,ii,1,5
  *do,JJ,1,VALUE(II+5,1)
  j=j+1
  LABEL(j)  = 'DT%ii%D%jj%'  $ ~CFGET,value( j,2),BRIDGE,,SECTION,DIV,TOP,ii,jj
  *enddo
  *enddo

VALUE(122,1) =5,7,9,11,13,15,7,9,11,13,15,17,19,9,11,13,15,17
VALUE(140,1) =19,21,23,25,27,11,13,15,17,19,21,23,25,27,3
VALUE(155,1) = 4,13,15,7,8,19

  *do,ii,1,5
  *do,JJ,1,VALUE(II+5,1)
  j=j+1
  LABEL(j)  = 'DD%ii%D%jj%'  $ ~CFGET,value( j,2),BRIDGE,,SECTION,DIV,BOT,ii,jj
  *enddo
  *enddo

VALUE(163,1) = 1,1
VALUE(167,1) = 1
VALUE(173,1) = 1
VALUE(175,1) = 1
VALUE(179,1) = 1
VALUE(182,1) = 1
VALUE(184,1) = 1,1
VALUE(188,1) = 1,1
VALUE(191,1) = 1,1,1,1
VALUE(198,1) = 1,1

 *do,ii,1,5
  *do,JJ,1,VALUE(II+5,1)
   j=j+1
   LABEL(j)  = 'A%ii%D%jj%'  $ ~CFGET,value( j,2),BRIDGE,,SECTION,DIV,SUPPORT,ii,jj
  *enddo
 *enddo

 *DO,II,1,NComp
   TOLER(II,1)  = 1E-3
   TOLER(II,2)  = 1E-3
 *ENDDO

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