summaryrefslogtreecommitdiff
path: root/Dragon/data/VanDerGucht_proc/svdg12.c2m
blob: c6d5f73052823b01d102398ae445472d71f9f889 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
XSM_FILE XSM1 XSM2 XSM3 XSM4 ;
PROCEDURE svdg12_auto svdg12_shiba svdg12_uss svdg12_rse ;
MODULE VDG: DELETE: ABORT: END: ;
SEQ_ASCII SECLIB_XSM :: FILE './SECLIB_XSM' ;
XSM_FILE DUMMYDA ;
REAL zaver REFVALUE DELTA ;
*
* import the cross section library (autolib)
DUMMYDA := SECLIB_XSM ;
*
* TEST AUTO: SELF-SHIELDING MODULE.
XSM1 := svdg12_auto ;
*
* TEST SHI: SELF-SHIELDING MODULE.
XSM2 := svdg12_shiba ;
*
* TEST USS/PTSL: SELF-SHIELDING MODULE.
XSM3 := svdg12_uss ;
*
* TEST USS/RSE: SELF-SHIELDING MODULE.
XSM4 := svdg12_rse ;
*
* COMPARE THE TWO RUNS.
VDG: XSM1 XSM2 :: EDIT 2 PICK AVER >>zaver<< ;
ECHO "SHI: averaged error=" zaver "%" ;

EVALUATE REFVALUE := 1.722 ;
EVALUATE DELTA := zaver REFVALUE - REFVALUE / ABS ;
IF DELTA 5.0E-2 < THEN
  ECHO "TEST SUCCESSFUL; DELTA=" DELTA ;
ELSE
  ECHO "------------" ;
  ECHO "TEST FAILURE" ;
  ECHO "------------" ;
  ECHO "REFERENCE=" REFVALUE " CALCULATED=" zaver ;
  ABORT: ;
ENDIF ;
*
* COMPARE THE TWO RUNS.
VDG: XSM1 XSM3 :: EDIT 2 PICK AVER >>zaver<< ;
ECHO "USS/PTSL: averaged error=" zaver "%" ;

EVALUATE REFVALUE := 7.616 ;
EVALUATE DELTA := zaver REFVALUE - REFVALUE / ABS ;
IF DELTA 5.0E-2 < THEN
  ECHO "TEST SUCCESSFUL; DELTA=" DELTA ;
ELSE
  ECHO "------------" ;
  ECHO "TEST FAILURE" ;
  ECHO "------------" ;
  ECHO "REFERENCE=" REFVALUE " CALCULATED=" zaver ;
  ABORT: ;
ENDIF ;
*
* COMPARE THE TWO RUNS.
VDG: XSM1 XSM4 :: EDIT 2 PICK AVER >>zaver<< ;
ECHO "USS/RSE: averaged error=" zaver "%" ;

EVALUATE REFVALUE := 3.54031E-1 ;
EVALUATE DELTA := zaver REFVALUE - REFVALUE / ABS ;
IF DELTA 5.0E-2 < THEN
  ECHO "TEST SUCCESSFUL; DELTA=" DELTA ;
ELSE
  ECHO "------------" ;
  ECHO "TEST FAILURE" ;
  ECHO "------------" ;
  ECHO "REFERENCE=" REFVALUE " CALCULATED=" zaver ;
  ABORT: ;
ENDIF ;
*
XSM1 XSM2 XSM3 XSM4 := DELETE: XSM1 XSM2 XSM3 XSM4 ;
ECHO "test svdg12 completed" ;
END: ;
QUIT "LIST" .