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

   FINISH
  ~CFCLEAR,,1
  NomFile='CFVR4208'
  /TITLE, %NomFile%, GB50011-2010: Spectrum definition

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

! Seismic spectrum definition
  ~DEFSPEC,ALL,0.05,RARE,I0,1,0.04,0.7 

! ----------------------------------------------------------------------
! 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.143344444
  VALUE(02,1) = 0.160688889
  VALUE(03,1) = 0.178033333
  VALUE(04,1) = 0.195377778
  VALUE(05,1) = 0.212722222
  VALUE(06,1) = 0.230066667
  VALUE(07,1) = 0.247411111
  VALUE(08,1) = 0.264755556
  VALUE(09,1) = 0.2821
  VALUE(10,1) = 0.299444444
  VALUE(11,1) = 0.299444444
  VALUE(12,1) = 0.187302716
  VALUE(13,1) = 0.137506354
  VALUE(14,1) = 0.109162017
  VALUE(15,1) = 0.090786757
  VALUE(16,1) = 0.077872366
  VALUE(17,1) = 0.06828092
  VALUE(18,1) = 0.058063749
  VALUE(19,1) = 0.047846577
  VALUE(20,1) = 0.037629405
  
! Warning and error tolerances
! ----------------------------------------------------------------------
   *DO,II,1,20
     TOLER(II, 1)= 1E-5 $ TOLER(II, 2)= 1E-5
   *ENDDO

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