!*************************************************************************************!
! TITLE: CFVR4204 - GB50011: Spectrum definition
! SUBTITLE: Spectrum definition according to GB 50011
!
! DESCRIPTION: CivilFEM spectrum is defined and compared to the correct values.
!
! ELEMENT TYPE: NONE
! MODULES:
! UNITS: User
! KEYWORD1: GB 50011
!
!*************************************************************************************!

   FINISH
  ~CFCLEAR,,1
  NomFile='CFVR4204'
  /TITLE, %NomFile%, GB50011: Spectrum definition

! ----------------------------------------------------------------------
! Model definition
! ----------------------------------------------------------------------
! CivilFEM Setup: Code & Units
  ~UNITS,,LENG,CM
  ~UNITS,,TIME,S
  ~UNITS,,FORC,KP
  ~CODESEL,,,,,GB50011
/PREP7
! Preprocessor
! ----------------------------------------------------------------------

! Seismic spectrum definition
  ~DEFSPEC,ALL,0.20,RARE,I,1,0.04,0.7

  pi=3.141592654

! ----------------------------------------------------------------------
! DATA CHECK
! ----------------------------------------------------------------------
! Data comparision number
  NComp = 20
  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

*DO,II,1,20
 ! Obtained values
   ~CFGET,VALUE(20-II+1,2),SEISM,,SPEC,SDH,,II,1
*ENDDO

! Correct values
! ----------------------------------------------------------------------
  VALUE(01,1) = 0.715921875
  VALUE(02,1) = 0.7441875
  VALUE(03,1) = 0.772453125
  VALUE(04,1) = 0.80071875
  VALUE(05,1) = 0.828984375
  VALUE(06,1) = 0.85725
  VALUE(07,1) = 0.885515625
  VALUE(08,1) = 0.91378125
  VALUE(09,1) = 0.942046875
  VALUE(10,1) = 0.9703125
  VALUE(11,1) = 0.9703125
  VALUE(12,1) = 0.608244911
  VALUE(13,1) = 0.447173139
  VALUE(14,1) = 0.355374678
  VALUE(15,1) = 0.29580555
  VALUE(16,1) = 0.253906769
  VALUE(17,1) = 0.222768272
  VALUE(18,1) = 0.192487022
  VALUE(19,1) = 0.162205772
  VALUE(20,1) = 0.131924522

! Warning and error tolerances
! ----------------------------------------------------------------------
   *DO,II,1,20
     TOLER(II, 1)= 1E-5 $ TOLER(II, 2)= 1E-5
   *ENDDO

! ----------------------------------------------------------------------
! Results Comparison
!-----------------------------------------------------------------------
  COMPARA.MAC