!*******************************************************************************
! TITLE:  CFVR2958 - Shell shear check according to SP63133 (compression & tension)
! SUBTITLE: Check Shear plus axial force of a reinforced concrete shell according to SP63133 
!
! ELEMENT TYPE: Shell_181
! MODULES:
! UNITS: User
! KEYWORD1: Shell
! KEYWORD2: Shell Vertex
! KEYWORD3: Reinforced Concrete
! KEYWORD3: Shear check
! KEYWORD4: SP63133
!*******************************************************************************
   FINISH
  ~CFCLEAR,,1
   NomFile='CFVR2958'
  /TITLE, %NomFile%, Shell shear check with SP63133 (compression & tension)

! ------------------------------------------------------------------------------
! Model definition and solve
! ------------------------------------------------------------------------------

! Parameters
   A = 3
   B = 6
   T = 0.2         ! Thickness
   R = 0.03        ! Cover 
   AST = 3.0708e-3
   ASB = 6.296e-3
   EsIZ = 1
   S =0.1
   ASW = 0.00022831/S

   
! Code
~CODESEL,,SP63133,

! Units
 ~UNITS,,FORC,KN 
  
/PREP7
! Material
~CFMP,1,LIB,CONCRETE,SP63133,B25
~CFMP,2,LIB,REINF,SP63133,A400  
~CFMP,1,SP_C    ,GAMb    ,,18.5/14.5
~CFMP,1,SP_C    ,GAMbt   ,,1.55/1.05
  
! Element types
ET,1,181
KEYOPT,1,3,2
KEYOPT,1,8,2
   
! Shell Vertex 
~SHLRNF,1,T,1,2,R,AST,ASB,AST,ASB,0,0.000,45.00
~SHLSHR,1,2,,,,,,,ASW,ASW 

! Shell properties
~BMSHPRO,1,SHELL,1,1,1,1, 181,,, ,Shell 1   
  
! Geometry
RECTNG,0,A,0,B
  
! Mesh
AESIZE,1,EsIZ,  
AMESH,1
  
FINISH
/SOLU 
 
! BCs
NSEL,S,LOC,X,0
NSEL,A,LOC,X,A
D,ALL,UX  
NSEL,S,LOC,Y,0
D,ALL, , , , , ,UX,UY,UZ,ROTZ, ,   
NSEL,S,LOC,Y,B
D,ALL, , , , , ,UX,UZ,ROTZ, , ,
ALLSEL
  
! CASE 1: COMPRESSION 1
Q = 60
H = 100

! Vertical Pressure
!SFE,ALL,1,PRES, ,-Q

! Compression pressure to forces
!NSEL,S,LOC,Y,B
!*GET,NNODE,NODE,,COUNT
!F,ALL,FY,-H*A/NNODE

ALLSEL
SOLVE		! CASE 1

FDELE,ALL,ALL
SFEDELE,ALL,ALL,ALL

! CASE 2: COMPRESSION 2
Q = 60
H = 2000

! Vertical Pressure
!SFE,ALL,1,PRES, ,-Q

! Compression pressure to forces
! NSEL,S,LOC,Y,B
! *GET,NNODE,NODE,,COUNT
! F,ALL,FY,-H*A/NNODE

ALLSEL
SOLVE		! CASE 2

FDELE,ALL,ALL
SFEDELE,ALL,ALL,ALL

! CASE 3:  COMPRESSION
Q = 60
H = 2320

! Vertical Pressure
! SFE,ALL,1,PRES, ,-Q

! Compression pressure to forces
! NSEL,S,LOC,Y,B
! *GET,NNODE,NODE,,COUNT
! F,ALL,FY,-H*A/NNODE

ALLSEL
SOLVE		! CASE 3

FDELE,ALL,ALL
SFEDELE,ALL,ALL,ALL

! CASE 4:  TENSION
Q = 60
H = 200

! Vertical Pressure
! SFE,ALL,1,PRES, ,-Q

! Compression pressure to forces
! NSEL,S,LOC,Y,B
! *GET,NNODE,NODE,,COUNT
! F,ALL,FY,H*A/NNODE

ALLSEL
SOLVE		! CASE 4

FDELE,ALL,ALL
SFEDELE,ALL,ALL,ALL


FINISH
/POST1

! CASE 1
~CFSET,,1

ESEL,S, , ,     15  ! Lateral endside element
~FMREAD,CFVR2958,FM 
~CHKCON,SHELL,SHEAR

! CASE 2
~CFSET,,2

ESEL,S, , ,     16  ! Lateral endside element
~FMREAD,CFVR2958,FM 
~CHKCON,SHELL,SHEAR

! CASE 3
~CFSET,,3

ESEL,S, , ,     17  ! Lateral endside element
~FMREAD,CFVR2958,FM 
~CHKCON,SHELL,SHEAR

! CASE 4
~CFSET,,4

ESEL,S, , ,     18  ! Lateral endside element
~FMREAD,CFVR2958,FM 
~CHKCON,SHELL,SHEAR
 
!--------------------------------------------------------------------------------------
! DATA CHECK
!--------------------------------------------------------------------------------------
! Data comparison number
  NComp = 68
  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

! Labels
!--------------------------------------------------------------------------------------
! CASE 1 CHECK
  LABEL(1 ) = 'C1_Q1_X'
  LABEL(2 ) = 'C1_CRTQ1_X'
  LABEL(3 ) = 'C1_QB_X'
  LABEL(4 ) = 'C1_QSW_X'
  LABEL(5 ) = 'C1_Q2_X'
  LABEL(6 ) = 'C1_CRTQ2_X'
  LABEL(7 ) = 'C1_CRT_TOTX'
  LABEL(8 ) = 'C1_PHI_N_X'
  LABEL(9 ) = 'C1_Q1_Y'
  LABEL(10) = 'C1_CRTQ1_Y'
  LABEL(11) = 'C1_QB_Y'
  LABEL(12) = 'C1_QSW_Y'
  LABEL(13) = 'C1_Q2_Y'
  LABEL(14) = 'C1_CRTQ2_Y'
  LABEL(15) = 'C1_CRT_TOTY'
  LABEL(16) = 'C1_PHI_N_Y'
  LABEL(17) = 'C1_CRT_TOT'
 
  
! CASE 2 CHECK
  LABEL(18) = 'C2_Q1_X'
  LABEL(19) = 'C2_CRTQ1_X'
  LABEL(20) = 'C2_QB_X'
  LABEL(21) = 'C2_QSW_X'
  LABEL(22) = 'C2_Q2_X'
  LABEL(23) = 'C2_CRTQ2_X'
  LABEL(24) = 'C2_CRT_TOTX'
  LABEL(25) = 'C2_PHI_N_X'
  LABEL(26) = 'C2_Q1_Y'
  LABEL(27) = 'C2_CRTQ1_Y'
  LABEL(28) = 'C2_QB_Y'
  LABEL(29) = 'C2_QSW_Y'
  LABEL(30) = 'C2_Q2_Y'
  LABEL(31) = 'C2_CRTQ2_Y'
  LABEL(32) = 'C2_CRT_TOTY'
  LABEL(33) = 'C2_PHI_N_Y'
  LABEL(34) = 'C2_CRT_TOT'
  
! CASE 3 check
  LABEL(35) = 'C3_Q1_X'
  LABEL(36) = 'C3_CRTQ1_X'
  LABEL(37) = 'C3_QB_X'
  LABEL(38) = 'C3_QSW_X'
  LABEL(39) = 'C3_Q2_X'
  LABEL(40) = 'C3_CRTQ2_X'
  LABEL(41) = 'C3_CRT_TOTX'
  LABEL(42) = 'C3_PHI_N_X'
  LABEL(43) = 'C3_Q1_Y'
  LABEL(44) = 'C3_CRTQ1_Y'
  LABEL(45) = 'C3_QB_Y'
  LABEL(46) = 'C3_QSW_Y'
  LABEL(47) = 'C3_Q2_Y'
  LABEL(48) = 'C3_CRTQ2_Y'
  LABEL(49) = 'C3_CRT_TOTY'
  LABEL(50) = 'C3_PHI_N_Y'
  LABEL(51) = 'C3_CRT_TOT'
  
! CASE 3 check
  LABEL(52) = 'C4_Q1_X'
  LABEL(53) = 'C4_CRTQ1_X'
  LABEL(54) = 'C4_QB_X'
  LABEL(55) = 'C4_QSW_X'
  LABEL(56) = 'C4_Q2_X'
  LABEL(57) = 'C4_CRTQ2_X'
  LABEL(58) = 'C4_CRT_TOTX'
  LABEL(59) = 'C4_PHI_N_X'
  LABEL(60) = 'C4_Q1_Y'
  LABEL(61) = 'C4_CRTQ1_Y'
  LABEL(62) = 'C4_QB_Y'
  LABEL(63) = 'C4_QSW_Y'
  LABEL(64) = 'C4_Q2_Y'
  LABEL(65) = 'C4_CRTQ2_Y'
  LABEL(66) = 'C4_CRT_TOTY'
  LABEL(67) = 'C4_PHI_N_Y'
  LABEL(68) = 'C4_CRT_TOT'
  

! Correct values
! CASE 1 CHECK
  VALUE(1 ,1)= 744.677
  VALUE(2 ,1)= 204.4/VALUE(1 ,1)
  VALUE(3 ,1)= 89.250
  VALUE(4 ,1)= 108
  VALUE(5 ,1)= 198.63
  VALUE(6 ,1)= 204.4/VALUE(5 ,1)
  VALUE(7 ,1)= MAX(VALUE(2 ,1),VALUE(6 ,1)) ! max CRTQ1_X,CRTQ2_X
  VALUE(8 ,1)= 1.007
  VALUE(9 ,1)= 764.64
  VALUE(10,1)= 394.5/VALUE(9 ,1)
  VALUE(11,1)= 89.250
  VALUE(12,1)= 108
  VALUE(13,1)= 203.96
  VALUE(14,1)= 394.5/VALUE(13,1)
  VALUE(15,1)= MAX(VALUE(10 ,1),VALUE(14 ,1)) ! max CRTQ1_Y,CRTQ2_Y
  VALUE(16,1)= 1.034
  VALUE(17,1)= VALUE(7 ,1) + VALUE(15,1)
  
! CASE 2 CHECK
  VALUE(18,1)= 841.551
  VALUE(19,1)= 204.4/VALUE(18 ,1)
  VALUE(20,1)= 89.250
  VALUE(21,1)= 108
  VALUE(22,1)= 224.472
  VALUE(23,1)= 204.4/VALUE(22 ,1)
  VALUE(24,1)= MAX(VALUE(19 ,1),VALUE(23 ,1)) ! max CRTQ1_X,CRTQ2_X
  VALUE(25,1)= 1.138
  VALUE(26,1)= 924.375
  VALUE(27,1)= 394.5/VALUE(26 ,1)
  VALUE(28,1)= 89.250
  VALUE(29,1)= 108
  VALUE(30,1)= 246.564
  VALUE(31,1)= 394.5/VALUE(30,1)
  VALUE(32,1)= MAX(VALUE(27 ,1),VALUE(31 ,1)) ! max CRTQ1_Y,CRTQ2_Y
  VALUE(33,1)= 1.25
  VALUE(34,1)= VALUE(24 ,1)+VALUE(32,1)
  
! CASE 3
  VALUE(35,1)= 857.8
  VALUE(36,1)= 204.4/VALUE(35 ,1)
  VALUE(37,1)= 89.250
  VALUE(38,1)= 108
  VALUE(39,1)= 228.8
  VALUE(40,1)= 204.4/VALUE(39 ,1)
  VALUE(41,1)= MAX(VALUE(36 ,1),VALUE(40 ,1)) ! max CRTQ1_X,CRTQ2_X
  VALUE(42,1)= 1.16
  VALUE(43,1)= 739.5
  VALUE(44,1)= 394.5/VALUE(43 ,1)
  VALUE(45,1)= 89.250
  VALUE(46,1)= 108
  VALUE(47,1)= 197.25    ! 
  VALUE(48,1)= 394.5/VALUE(47,1)
  VALUE(49,1)= MAX(VALUE(44 ,1),VALUE(48 ,1)) ! max CRTQ1_Y,CRTQ2_Y
  VALUE(50,1)= 1.0
  VALUE(51,1)= VALUE(41 ,1)+VALUE(49,1)
  
! CASE 4
  VALUE(52,1)= 369.750
  VALUE(53,1)= 204.4/VALUE(52 ,1)
  VALUE(54,1)= 89.250
  VALUE(55,1)= 108
  VALUE(56,1)= 98.625
  VALUE(57,1)= 204.4/VALUE(56 ,1)
  VALUE(58,1)= MAX(VALUE(53 ,1),VALUE(57 ,1)) ! max CRTQ1_X,CRTQ2_X
  VALUE(59,1)= 0.5
  VALUE(60,1)= 387
  VALUE(61,1)= 394.5/VALUE(60 ,1)
  VALUE(62,1)= 89.250
  VALUE(63,1)= 108
  VALUE(64,1)= 103.360    ! 
  VALUE(65,1)= 394.5/VALUE(64,1)
  VALUE(66,1)= MAX(VALUE(61 ,1),VALUE(65 ,1)) ! max CRTQ1_Y,CRTQ2_Y
  VALUE(67,1)= 0.524
  VALUE(68,1)= VALUE(58 ,1)+VALUE(66,1)
  
 
! Comparison
 ~CFSET,1
 ~CFGET,VALUE(1 ,2),ELEMENT,15,RESULT,Q1_X,I
 ~CFGET,VALUE(2 ,2),ELEMENT,15,RESULT,CRTQ1_X,I
 ~CFGET,VALUE(3 ,2),ELEMENT,15,RESULT,QB_X,I
 ~CFGET,VALUE(4 ,2),ELEMENT,15,RESULT,QSW_X,I
 ~CFGET,VALUE(5 ,2),ELEMENT,15,RESULT,Q2_X,I
 ~CFGET,VALUE(6 ,2),ELEMENT,15,RESULT,CRTQ2_X,I
 ~CFGET,VALUE(7 ,2),ELEMENT,15,RESULT,CRT_TOTX,I
 ~CFGET,VALUE(8 ,2),ELEMENT,15,RESULT,PHI_N_X,I
 ~CFGET,VALUE(9 ,2),ELEMENT,15,RESULT,Q1_Y,I
 ~CFGET,VALUE(10,2),ELEMENT,15,RESULT,CRTQ1_Y,I
 ~CFGET,VALUE(11,2),ELEMENT,15,RESULT,QB_Y,I
 ~CFGET,VALUE(12,2),ELEMENT,15,RESULT,QSW_Y,I
 ~CFGET,VALUE(13,2),ELEMENT,15,RESULT,Q2_Y,I
 ~CFGET,VALUE(14,2),ELEMENT,15,RESULT,CRTQ2_Y,I
 ~CFGET,VALUE(15,2),ELEMENT,15,RESULT,CRT_TOTY,I
 ~CFGET,VALUE(16,2),ELEMENT,15,RESULT,PHI_N_Y,I
 ~CFGET,VALUE(17,2),ELEMENT,15,RESULT,CRT_TOT,I
 
 
 ~CFSET,2
 ~CFGET,VALUE(18,2),ELEMENT,16,RESULT,Q1_X,I
 ~CFGET,VALUE(19,2),ELEMENT,16,RESULT,CRTQ1_X,I
 ~CFGET,VALUE(20,2),ELEMENT,16,RESULT,QB_X,I
 ~CFGET,VALUE(21,2),ELEMENT,16,RESULT,QSW_X,I
 ~CFGET,VALUE(22,2),ELEMENT,16,RESULT,Q2_X,I
 ~CFGET,VALUE(23,2),ELEMENT,16,RESULT,CRTQ2_X,I
 ~CFGET,VALUE(24,2),ELEMENT,16,RESULT,CRT_TOTX,I
 ~CFGET,VALUE(25,2),ELEMENT,16,RESULT,PHI_N_X,I
 ~CFGET,VALUE(26,2),ELEMENT,16,RESULT,Q1_Y,I
 ~CFGET,VALUE(27,2),ELEMENT,16,RESULT,CRTQ1_Y,I
 ~CFGET,VALUE(28,2),ELEMENT,16,RESULT,QB_Y,I
 ~CFGET,VALUE(29,2),ELEMENT,16,RESULT,QSW_Y,I
 ~CFGET,VALUE(30,2),ELEMENT,16,RESULT,Q2_Y,I
 ~CFGET,VALUE(31,2),ELEMENT,16,RESULT,CRTQ2_Y,I
 ~CFGET,VALUE(32,2),ELEMENT,16,RESULT,CRT_TOTY,I
 ~CFGET,VALUE(33,2),ELEMENT,16,RESULT,PHI_N_Y,I
 ~CFGET,VALUE(34,2),ELEMENT,16,RESULT,CRT_TOT,I
 
 ~CFSET,3
 ~CFGET,VALUE(35,2),ELEMENT,17,RESULT,Q1_X,I
 ~CFGET,VALUE(36,2),ELEMENT,17,RESULT,CRTQ1_X,I
 ~CFGET,VALUE(37,2),ELEMENT,17,RESULT,QB_X,I
 ~CFGET,VALUE(38,2),ELEMENT,17,RESULT,QSW_X,I
 ~CFGET,VALUE(39,2),ELEMENT,17,RESULT,Q2_X,I
 ~CFGET,VALUE(40,2),ELEMENT,17,RESULT,CRTQ2_X,I
 ~CFGET,VALUE(41,2),ELEMENT,17,RESULT,CRT_TOTX,I
 ~CFGET,VALUE(42,2),ELEMENT,17,RESULT,PHI_N_X,I
 ~CFGET,VALUE(43,2),ELEMENT,17,RESULT,Q1_Y,I
 ~CFGET,VALUE(44,2),ELEMENT,17,RESULT,CRTQ1_Y,I
 ~CFGET,VALUE(45,2),ELEMENT,17,RESULT,QB_Y,I
 ~CFGET,VALUE(46,2),ELEMENT,17,RESULT,QSW_Y,I
 ~CFGET,VALUE(47,2),ELEMENT,17,RESULT,Q2_Y,I
 ~CFGET,VALUE(48,2),ELEMENT,17,RESULT,CRTQ2_Y,I
 ~CFGET,VALUE(49,2),ELEMENT,17,RESULT,CRT_TOTY,I
 ~CFGET,VALUE(50,2),ELEMENT,17,RESULT,PHI_N_Y,I
 ~CFGET,VALUE(51,2),ELEMENT,17,RESULT,CRT_TOT,I
 
 ~CFSET,4
 ~CFGET,VALUE(52,2),ELEMENT,18,RESULT,Q1_X,I
 ~CFGET,VALUE(53,2),ELEMENT,18,RESULT,CRTQ1_X,I
 ~CFGET,VALUE(54,2),ELEMENT,18,RESULT,QB_X,I
 ~CFGET,VALUE(55,2),ELEMENT,18,RESULT,QSW_X,I
 ~CFGET,VALUE(56,2),ELEMENT,18,RESULT,Q2_X,I
 ~CFGET,VALUE(57,2),ELEMENT,18,RESULT,CRTQ2_X,I
 ~CFGET,VALUE(58,2),ELEMENT,18,RESULT,CRT_TOTX,I
 ~CFGET,VALUE(59,2),ELEMENT,18,RESULT,PHI_N_X,I
 ~CFGET,VALUE(60,2),ELEMENT,18,RESULT,Q1_Y,I
 ~CFGET,VALUE(61,2),ELEMENT,18,RESULT,CRTQ1_Y,I
 ~CFGET,VALUE(62,2),ELEMENT,18,RESULT,QB_Y,I
 ~CFGET,VALUE(63,2),ELEMENT,18,RESULT,QSW_Y,I
 ~CFGET,VALUE(64,2),ELEMENT,18,RESULT,Q2_Y,I
 ~CFGET,VALUE(65,2),ELEMENT,18,RESULT,CRTQ2_Y,I
 ~CFGET,VALUE(66,2),ELEMENT,18,RESULT,CRT_TOTY,I
 ~CFGET,VALUE(67,2),ELEMENT,18,RESULT,PHI_N_Y,I
 ~CFGET,VALUE(68,2),ELEMENT,18,RESULT,CRT_TOT,I
 

! Warning and error tolerances
! CASE 1 CHECK
TOLER(1 ,1)= VALUE(1 ,1)*2E-2 $ TOLER( 1,2)=VALUE(1 ,1)*2E-2 
TOLER(2 ,1)= VALUE(2 ,1)*1E-2 $ TOLER( 2,2)=VALUE(2 ,1)*1E-2 
TOLER(3 ,1)= VALUE(3 ,1)*2E-2 $ TOLER( 3,2)=VALUE(3 ,1)*2E-2
TOLER(4 ,1)= VALUE(4 ,1)*1E-2 $ TOLER( 4,2)=VALUE(4 ,1)*1E-2
TOLER(5 ,1)= VALUE(5 ,1)*1E-2 $ TOLER( 5,2)=VALUE(5 ,1)*1E-2
TOLER(6 ,1)= VALUE(6 ,1)*1E-2 $ TOLER( 6,2)=VALUE(6 ,1)*1E-2
TOLER(7 ,1)= VALUE(7 ,1)*1E-2 $ TOLER( 7,2)=VALUE(7 ,1)*1E-2
TOLER(8 ,1)= VALUE(8 ,1)*1E-2 $ TOLER( 8,2)=VALUE(8 ,1)*1E-2 
TOLER(9 ,1)= VALUE(9 ,1)*2E-2 $ TOLER( 9,2)=VALUE(9 ,1)*2E-2 
TOLER(10,1)= VALUE(10,1)*1E-2 $ TOLER(10,2)=VALUE(10,1)*1E-2
TOLER(11,1)= VALUE(11,1)*2E-2 $ TOLER(11,2)=VALUE(11,1)*2E-2
TOLER(12,1)= VALUE(12,1)*1E-2 $ TOLER(12,2)=VALUE(12,1)*1E-2
TOLER(13,1)= VALUE(13,1)*1E-2 $ TOLER(13,2)=VALUE(13,1)*1E-2
TOLER(14,1)= VALUE(14,1)*1E-2 $ TOLER(14,2)=VALUE(14,1)*1E-2
TOLER(15,1)= VALUE(15,1)*1E-2 $ TOLER(15,2)=VALUE(15,1)*1E-2
TOLER(16,1)= VALUE(16,1)*1E-2 $ TOLER(16,2)=VALUE(16,1)*1E-2 
TOLER(17,1)= VALUE(17,1)*1E-2 $ TOLER(17,2)=VALUE(17,1)*1E-2 


! CASE 2 CHECK
TOLER(18,1)= VALUE(18,1)*2E-2 $ TOLER(18,2)=VALUE(18,1)*2E-2 
TOLER(19,1)= VALUE(19,1)*1E-2 $ TOLER(19,2)=VALUE(19,1)*1E-2 
TOLER(20,1)= VALUE(20,1)*2E-2 $ TOLER(20,2)=VALUE(20,1)*2E-2
TOLER(21,1)= VALUE(21,1)*1E-2 $ TOLER(21,2)=VALUE(21,1)*1E-2
TOLER(22,1)= VALUE(22,1)*1E-2 $ TOLER(22,2)=VALUE(22,1)*1E-2
TOLER(23,1)= VALUE(23,1)*1E-2 $ TOLER(23,2)=VALUE(23,1)*1E-2
TOLER(24,1)= VALUE(24,1)*1E-2 $ TOLER(24,2)=VALUE(24,1)*1E-2
TOLER(25,1)= VALUE(25,1)*1E-2 $ TOLER(25,2)=VALUE(25,1)*1E-2 
TOLER(26,1)= VALUE(26,1)*2E-2 $ TOLER(26,2)=VALUE(26,1)*2E-2 
TOLER(27,1)= VALUE(27,1)*1E-2 $ TOLER(27,2)=VALUE(27,1)*1E-2
TOLER(28,1)= VALUE(28,1)*2E-2 $ TOLER(28,2)=VALUE(28,1)*2E-2
TOLER(29,1)= VALUE(29,1)*1E-2 $ TOLER(29,2)=VALUE(29,1)*1E-2
TOLER(30,1)= VALUE(30,1)*1E-2 $ TOLER(30,2)=VALUE(30,1)*1E-2
TOLER(31,1)= VALUE(31,1)*1E-2 $ TOLER(31,2)=VALUE(31,1)*1E-2
TOLER(32,1)= VALUE(32,1)*1E-2 $ TOLER(32,2)=VALUE(32,1)*1E-2
TOLER(33,1)= VALUE(33,1)*1E-2 $ TOLER(33,2)=VALUE(33,1)*1E-2 
TOLER(34,1)= VALUE(34,1)*1E-2 $ TOLER(34,2)=VALUE(34,1)*1E-2 


! CASE 3 CHECK
TOLER(35,1)= VALUE(35,1)*2E-2 $ TOLER(35,2)=VALUE(35,1)*2E-2 
TOLER(36,1)= VALUE(36,1)*1E-2 $ TOLER(36,2)=VALUE(36,1)*1E-2 
TOLER(37,1)= VALUE(37,1)*2E-2 $ TOLER(37,2)=VALUE(37,1)*2E-2
TOLER(38,1)= VALUE(38,1)*1E-2 $ TOLER(38,2)=VALUE(38,1)*1E-2
TOLER(39,1)= VALUE(39,1)*1E-2 $ TOLER(39,2)=VALUE(39,1)*1E-2
TOLER(40,1)= VALUE(40,1)*1E-2 $ TOLER(40,2)=VALUE(40,1)*1E-2
TOLER(41,1)= VALUE(41,1)*1E-2 $ TOLER(41,2)=VALUE(41,1)*1E-2
TOLER(42,1)= VALUE(42,1)*1E-2 $ TOLER(42,2)=VALUE(42,1)*1E-2 
TOLER(43,1)= VALUE(43,1)*2E-2 $ TOLER(43,2)=VALUE(43,1)*2E-2 
TOLER(44,1)= VALUE(44,1)*1E-2 $ TOLER(44,2)=VALUE(44,1)*1E-2
TOLER(45,1)= VALUE(45,1)*2E-2 $ TOLER(45,2)=VALUE(45,1)*2E-2
TOLER(46,1)= VALUE(46,1)*1E-2 $ TOLER(46,2)=VALUE(46,1)*1E-2
TOLER(47,1)= VALUE(47,1)*1E-2 $ TOLER(47,2)=VALUE(47,1)*1E-2
TOLER(48,1)= VALUE(48,1)*1E-2 $ TOLER(48,2)=VALUE(48,1)*1E-2
TOLER(49,1)= VALUE(49,1)*1E-2 $ TOLER(49,2)=VALUE(49,1)*1E-2
TOLER(50,1)= VALUE(50,1)*1E-2 $ TOLER(50,2)=VALUE(50,1)*1E-2 
TOLER(51,1)= VALUE(51,1)*1E-2 $ TOLER(51,2)=VALUE(51,1)*1E-2 

! CASE 4 CHECK
TOLER(52,1)= VALUE(52,1)*2E-2 $ TOLER(52,2)=VALUE(52,1)*2E-2 
TOLER(53,1)= VALUE(53,1)*1E-2 $ TOLER(53,2)=VALUE(53,1)*1E-2 
TOLER(54,1)= VALUE(54,1)*2E-2 $ TOLER(54,2)=VALUE(54,1)*2E-2
TOLER(55,1)= VALUE(55,1)*1E-2 $ TOLER(55,2)=VALUE(55,1)*1E-2
TOLER(56,1)= VALUE(56,1)*1E-2 $ TOLER(56,2)=VALUE(56,1)*1E-2
TOLER(57,1)= VALUE(57,1)*1E-2 $ TOLER(57,2)=VALUE(57,1)*1E-2
TOLER(58,1)= VALUE(58,1)*1E-2 $ TOLER(58,2)=VALUE(58,1)*1E-2
TOLER(59,1)= VALUE(59,1)*1E-2 $ TOLER(59,2)=VALUE(59,1)*1E-2 
TOLER(60,1)= VALUE(60,1)*2E-2 $ TOLER(60,2)=VALUE(60,1)*2E-2 
TOLER(61,1)= VALUE(61,1)*1E-2 $ TOLER(61,2)=VALUE(61,1)*1E-2
TOLER(62,1)= VALUE(62,1)*2E-2 $ TOLER(62,2)=VALUE(62,1)*2E-2
TOLER(63,1)= VALUE(63,1)*1E-2 $ TOLER(63,2)=VALUE(63,1)*1E-2
TOLER(64,1)= VALUE(64,1)*1E-2 $ TOLER(64,2)=VALUE(64,1)*1E-2
TOLER(65,1)= VALUE(65,1)*1E-2 $ TOLER(65,2)=VALUE(65,1)*1E-2
TOLER(66,1)= VALUE(66,1)*1E-2 $ TOLER(66,2)=VALUE(66,1)*1E-2
TOLER(67,1)= VALUE(67,1)*1E-2 $ TOLER(67,2)=VALUE(67,1)*1E-2 
TOLER(68,1)= VALUE(68,1)*1E-2 $ TOLER(68,2)=VALUE(68,1)*1E-2 


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