diff options
Diffstat (limited to 'Dragon')
| -rwxr-xr-x | Dragon/data/salmacro.access | 1 | ||||
| -rw-r--r-- | Dragon/data/salmacro.x2m | 3 | ||||
| -rw-r--r-- | Dragon/data/salmacro_proc/TDCM69.c2m | 113 | ||||
| -rw-r--r-- | Dragon/data/salmacro_proc/irena_tiso_s30.dat | 2424 | ||||
| -rw-r--r-- | Dragon/src/EXCELP.f | 28 | ||||
| -rw-r--r-- | Dragon/src/MCGSIG.f | 4 | ||||
| -rw-r--r-- | Dragon/src/MUSACG.f90 | 10 | ||||
| -rw-r--r-- | Dragon/src/PIJWIJ.f | 18 | ||||
| -rw-r--r-- | Dragon/src/SALACG.f90 | 30 | ||||
| -rw-r--r-- | Dragon/src/SALTCG.f | 16 |
10 files changed, 2605 insertions, 42 deletions
diff --git a/Dragon/data/salmacro.access b/Dragon/data/salmacro.access index 21f437a..4fd6844 100755 --- a/Dragon/data/salmacro.access +++ b/Dragon/data/salmacro.access @@ -14,5 +14,6 @@ ln -s $1/data/salmacro_proc/lattice_hex_tiso_S30_7cells.dat . ln -s $1/data/salmacro_proc/lattice_hex_tspc_S30_7cells.dat . ln -s $1/data/salmacro_proc/lattice_hex_tiso_S30_7cells_macro.dat . ln -s $1/data/salmacro_proc/lattice_hex_tiso_S30_tdcm68.dat . +ln -s $1/data/salmacro_proc/irena_tiso_s30.dat . ls -l echo "salmacro access script terminated" diff --git a/Dragon/data/salmacro.x2m b/Dragon/data/salmacro.x2m index eeb678c..df4d85c 100644 --- a/Dragon/data/salmacro.x2m +++ b/Dragon/data/salmacro.x2m @@ -6,7 +6,7 @@ PROCEDURE TDCM02 TDCM03 TDCM07 TDCM08 TDCM31 TDCM32 TDCM33 TDCM43 TDCM44 TDCM45 TDCM46 TDCM47 TDCM48 TDCM49 TDCM50 TDCM51 TDCM52 TDCM53 TDCM54 TDCM55 TDCM56 TDCM57 TDCM58 TDCM59 TDCM60 TDCM61 TDCM62 TDCM63 TDCM64 TDCM65 TDCM66 TDCM67 - TDCM68 ; + TDCM68 TDCM69 ; * TDCM02 ; TDCM03 ; @@ -41,4 +41,5 @@ TDCM65 ; TDCM66 ; TDCM67 ; TDCM68 ; +TDCM69 ; QUIT "LIST" . diff --git a/Dragon/data/salmacro_proc/TDCM69.c2m b/Dragon/data/salmacro_proc/TDCM69.c2m new file mode 100644 index 0000000..590d461 --- /dev/null +++ b/Dragon/data/salmacro_proc/TDCM69.c2m @@ -0,0 +1,113 @@ +*------------------------------------------------------------------------------------------------- +* Name : TDCM69.x2m +* Author : B. Cui +* MODEL TO TEST GLOW GEOMETRIES WITH MACRO-GEOMETRIES +*------------------------------------------------------------------------------------------------- +* Define STRUCTURES and MODULES used +*------------------------------------------------------------------------------------------------- +SEQ_ASCII TDT_LAT_DAT :: FILE './irena_tiso_s30.dat' ; +SEQ_ASCII TDT_LAT_PS :: FILE './irena_tiso_s30.eps' ; +MODULE + G2S: MAC: SALT: ASM: FLU: EDI: END: ; +LINKED_LIST + MACRO4G_0 MACRO4G TRACK CP FLUX OUT ; +SEQ_BINARY + TRACK_f ; +PROCEDURE assertS ; +*---- +* Input Macrolib +*---- +MACRO4G_0 := MAC: :: + EDIT 1 NGRO 4 NMIX 9 NIFI 1 ANIS 2 + READ INPUT + MIX 1 (* NA *) + TOTAL 5.316E-01 9.352E-01 9.931E-01 1.5931E+00 + TRANC 2.880E-01 5.529E-01 4.936E-01 4.4340E-01 + SCAT 1 1 (*1->1*) 4.973E-01 + 2 2 (*2->2*) 6.570E-01 (*1->2*) 3.134E-02 + 4 4 (*4->3*) 1.163E-02 (*3->3*) 7.176E-01 (*2->3*) 2.528E-01 + (*1->3*) 2.304E-03 + 4 4 (*4->4*) 1.565E+00 (*3->4*) 2.707E-01 (*2->4*) 2.412E-02 + (*1->4*) 2.228E-04 + 1 1 (*1->1*) 2.735E-01 + 2 2 (*2->2*) 4.253E-01 (*1->2*) 1.410E-02 + 4 4 (*4->3*) 5.064E-03 (*3->3*) 4.286E-01 (*2->3*) 1.203E-01 + (*1->3*) 3.050E-04 + 4 4 (*4->4*) 4.384E-01 (*3->4*) 6.426E-02 (*2->4*) 3.373E-03 + (*1->4*) 1.245E-05 + MIX 2 (* CLAD *) + TOTAL 4.240E-01 1.0476E+00 2.7503E+00 1.1194E+01 + TRANC 0.000E+00 0.000E+00 0.000E+00 0.000E+00 + SCAT 1 1 (*1->1*) 3.738E-01 + 2 2 (*2->2*) 7.102E-01 (*1->2*) 2.439E-04 + 3 4 (*4->3*) 4.985E-03 (*3->3*) 8.956E-01 (*2->3*) 7.370E-04 + 2 4 (*4->4*) 8.996E-01 (*3->4*) 3.248E-03 + 1 1 (*1->1*) 0.000E+00 + 1 2 (*2->2*) -6.161E-02 + 1 3 (*3->3*) -2.326E-01 + 1 4 (*4->4*) 0.000E+00 + MIX 3 (* VOID *) + TOTAL 4.240E-02 1.0476E-01 2.7503E-01 1.1194E+00 + TRANC 0.000E+00 0.000E+00 0.000E+00 0.000E+00 + SCAT 1 1 (*1->1*) 3.738E-02 + 2 2 (*2->2*) 9.102E-02 (*1->2*) 2.439E-04 + 3 4 (*4->3*) 4.985E-04 (*3->3*) 9.956E-02 (*2->3*) 7.370E-04 + 2 4 (*4->4*) 9.996E-02 (*3->4*) 3.248E-04 + 1 1 (*1->1*) 0.000E+00 + 1 2 (*2->2*) 0.000E+00 + 1 3 (*3->3*) 0.000E+00 + 1 4 (*4->4*) 0.000E+00 + MIX 4 (* FUEL *) + TOTAL 4.097E-01 6.286E-01 4.594E-01 7.508E-01 + TRANC 4.841E-02 7.939E-03 8.080E-03 8.038E-03 + NUSIGF 1.931E-02 1.259E-02 2.129E-01 2.508E-01 + CHI 1.000E+00 0.000E+00 0.000E+00 0.000E+00 + SCAT 1 1 (*1->1*) 3.908E-01 + 2 2 (*2->2*) 4.407E-01 (*1->2*) 8.546E-04 + 3 4 (*4->3*) 5.318E-03 (*3->3*) 3.785E-01 (*2->3*) 8.057E-03 + 2 4 (*4->4*) 3.983E-01 (*3->4*) 1.203E-02 + 1 1 (*1->1*) 4.793E-02 + 2 2 (*2->2*) 5.419E-03 (*1->2*) -2.640E-04 + 3 4 (*4->3*) -7.338E-04 (*3->3*) 9.014E-03 (*2->3*) -2.453E-03 + 2 4 (*4->4*) 8.769E-03 (*3->4*) -2.572E-03 +; +MACRO4G := MAC: MACRO4G_0 :: + EDIT 1 NMIX 9 + MIX 5 4 + MIX 6 4 + MIX 7 4 + MIX 8 3 + MIX 9 2 +; +*----------------------------------------------------- +* Tracking of the imported geometry +*----------------------------------------------------- +TDT_LAT_PS := G2S: TDT_LAT_DAT :: DRAWMIX ; + +TRACK TRACK_f := SALT: TDT_LAT_DAT :: + TITL 'TDCM69 imported' + EDIT 6 + ANIS 1 + TISO 8 10.0 + REND LONG + IC EPSJ 1.0E-4 + ; +*------------------------------------------------------------------------------------------------- +* Flux calculation +*------------------------------------------------------------------------------------------------- +CP := ASM: MACRO4G TRACK TRACK_f :: PIJ ; +FLUX := FLU: MACRO4G CP TRACK TRACK_f :: TYPE K ; + +*------------------------------------------------------------------------------------------------- +* Editing +*------------------------------------------------------------------------------------------------- +OUT := EDI: MACRO4G TRACK FLUX :: + EDIT 4 + MERG COMP + COND + ; + +assertS FLUX :: 'K-EFFECTIVE' 1 1.007306 ; +ECHO "test TDCM69 completed" ; +END: ; +QUIT "LIST" . diff --git a/Dragon/data/salmacro_proc/irena_tiso_s30.dat b/Dragon/data/salmacro_proc/irena_tiso_s30.dat new file mode 100644 index 0000000..df17691 --- /dev/null +++ b/Dragon/data/salmacro_proc/irena_tiso_s30.dat @@ -0,0 +1,2424 @@ + + dat input file for DRAGON5 +------------------------------------------------------------ +* typge, nbfo, node, elem, macr, nreg, z, mac2 + 1, 12, 207, 481, 26, 207, 0, 1 +* index kindex + 0 0 1 +* eps eps0 + 1.000000E-05 1.000000E-05 +* flux region number per geometry region (mesh) + 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, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, + 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, + 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, + 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, + 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, + 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, + 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, + 157, 158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 168, + 169, 170, 171, 172, 173, 174, 175, 176, 177, 178, 179, 180, + 181, 182, 183, 184, 185, 186, 187, 188, 189, 190, 191, 192, + 193, 194, 195, 196, 197, 198, 199, 200, 201, 202, 203, 204, + 205, 206, 207 +* names of macros + MAC001 MAC002 MAC003 MAC004 + MAC005 MAC006 MAC007 MAC008 + MAC009 MAC010 MAC011 MAC012 + MAC013 MAC014 MAC015 MAC016 + MAC017 MAC018 MAC019 MAC020 + MAC021 MAC022 MAC023 MAC024 + MAC025 MAC026 +* macro order number per flux region + 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, + 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, + 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, + 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, + 7, 7, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, + 8, 8, 8, 8, 9, 9, 9, 9, 9, 9, 9, 9, + 10, 10, 10, 10, 10, 10, 10, 10, 11, 11, 11, 11, + 11, 11, 11, 11, 12, 12, 12, 12, 12, 12, 12, 12, + 13, 13, 13, 13, 13, 13, 13, 13, 14, 14, 14, 14, + 14, 14, 14, 14, 15, 15, 15, 15, 15, 15, 15, 15, + 16, 16, 16, 16, 16, 16, 16, 16, 17, 17, 17, 17, + 17, 17, 17, 17, 18, 18, 18, 18, 18, 18, 18, 18, + 19, 19, 19, 19, 19, 19, 19, 19, 20, 20, 20, 20, + 20, 20, 20, 20, 21, 21, 21, 21, 21, 21, 21, 21, + 22, 22, 22, 22, 22, 22, 22, 22, 23, 23, 23, 23, + 23, 23, 23, 23, 24, 24, 24, 24, 24, 24, 24, 24, + 25, 25, 25, 25, 25, 25, 25, 25, 26, 26, 26, 26, + 26, 26, 26 + elements +* ELEM 1 line segment + 1, 202, 0 +* + 7.4962660E+00, 4.3279712E+00, 2.4987553E+00, -4.3279712E+00 +* ELEM 2 line segment + 1, 0, 202 +* + 9.7352137E+00, 0.0000000E+00, 2.5980762E-01, 0.0000000E+00 +* ELEM 3 line segment + 1, 202, 201 +* + 9.7352137E+00, 0.0000000E+00, -2.4338034E+00, 4.2154712E+00 +* ELEM 4 line segment + 1, 0, 202 +* + 7.4962660E+00, 4.3279712E+00, -1.9485572E-01, -1.1250000E-01 +* ELEM 5 line segment + 1, 153, 121 +* + 8.0250000E+00, 0.0000000E+00, 0.0000000E+00, 3.0888239E-01 +* ELEM 6 line segment + 1, 161, 153 +* + 8.5600000E+00, 6.1776479E-01, -5.3500000E-01, -3.0888239E-01 +* ELEM 7 line segment + 1, 205, 153 +* + 9.0083977E+00, 3.5888229E-01, -4.4839763E-01, 2.5888250E-01 +* ELEM 8 line segment + 1, 201, 153 +* + 9.0950000E+00, 2.0888260E-01, -8.6602367E-02, 1.4999970E-01 +* ELEM 9 line segment + 1, 207, 153 +* + 9.0950000E+00, 0.0000000E+00, 0.0000000E+00, 2.0888260E-01 +* ELEM 10 line segment + 1, 0, 153 +* + 9.0450000E+00, 0.0000000E+00, 5.0000002E-02, 0.0000000E+00 +* ELEM 11 circular arc + 3, 154, 153 +* + 8.5600000E+00, 0.0000000E+00, 4.8500000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 12 line segment + 1, 0, 153 +* + 8.0250000E+00, 0.0000000E+00, 5.0000002E-02, 0.0000000E+00 +* ELEM 13 line segment + 1, 161, 121 +* + 8.0250000E+00, 3.0888239E-01, -5.3500000E-01, 3.0888240E-01 +* ELEM 14 line segment + 1, 0, 121 +* + 7.9750000E+00, 0.0000000E+00, 5.0000002E-02, 0.0000000E+00 +* ELEM 15 circular arc + 3, 122, 121 +* + 7.4900000E+00, 0.0000000E+00, 4.8500000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 16 line segment + 1, 0, 121 +* + 6.9550000E+00, 0.0000000E+00, 5.0000002E-02, 0.0000000E+00 +* ELEM 17 line segment + 1, 121, 89 +* + 6.9550000E+00, 0.0000000E+00, 0.0000000E+00, 3.0888239E-01 +* ELEM 18 line segment + 1, 129, 121 +* + 7.4900000E+00, 6.1776479E-01, -5.3500000E-01, -3.0888239E-01 +* ELEM 19 line segment + 1, 129, 89 +* + 6.9550000E+00, 3.0888239E-01, -5.3500000E-01, 3.0888240E-01 +* ELEM 20 line segment + 1, 0, 89 +* + 6.9050000E+00, 0.0000000E+00, 5.0000002E-02, 0.0000000E+00 +* ELEM 21 circular arc + 3, 90, 89 +* + 6.4200000E+00, 0.0000000E+00, 4.8500000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 22 line segment + 1, 0, 89 +* + 5.8850000E+00, 0.0000000E+00, 5.0000002E-02, 0.0000000E+00 +* ELEM 23 line segment + 1, 89, 65 +* + 5.8850000E+00, 0.0000000E+00, 0.0000000E+00, 3.0888239E-01 +* ELEM 24 line segment + 1, 97, 89 +* + 6.4200000E+00, 6.1776479E-01, -5.3500000E-01, -3.0888239E-01 +* ELEM 25 line segment + 1, 97, 65 +* + 5.8850000E+00, 3.0888239E-01, -5.3500000E-01, 3.0888240E-01 +* ELEM 26 line segment + 1, 0, 65 +* + 5.8350000E+00, 0.0000000E+00, 5.0000002E-02, 0.0000000E+00 +* ELEM 27 circular arc + 3, 66, 65 +* + 5.3500000E+00, 0.0000000E+00, 4.8500000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 28 line segment + 1, 0, 65 +* + 4.8150000E+00, 0.0000000E+00, 5.0000002E-02, 0.0000000E+00 +* ELEM 29 line segment + 1, 65, 41 +* + 4.8150000E+00, 0.0000000E+00, 0.0000000E+00, 3.0888239E-01 +* ELEM 30 line segment + 1, 73, 65 +* + 5.3500000E+00, 6.1776479E-01, -5.3500000E-01, -3.0888239E-01 +* ELEM 31 line segment + 1, 73, 41 +* + 4.8150000E+00, 3.0888239E-01, -5.3500000E-01, 3.0888240E-01 +* ELEM 32 line segment + 1, 0, 41 +* + 4.7650000E+00, 0.0000000E+00, 5.0000002E-02, 0.0000000E+00 +* ELEM 33 circular arc + 3, 42, 41 +* + 4.2800000E+00, 0.0000000E+00, 4.8500000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 34 line segment + 1, 0, 41 +* + 3.7450000E+00, 0.0000000E+00, 5.0000002E-02, 0.0000000E+00 +* ELEM 35 line segment + 1, 41, 25 +* + 3.7450000E+00, 0.0000000E+00, 0.0000000E+00, 3.0888239E-01 +* ELEM 36 line segment + 1, 49, 41 +* + 4.2800000E+00, 6.1776479E-01, -5.3500000E-01, -3.0888239E-01 +* ELEM 37 line segment + 1, 49, 25 +* + 3.7450000E+00, 3.0888239E-01, -5.3500000E-01, 3.0888240E-01 +* ELEM 38 line segment + 1, 0, 25 +* + 3.6950000E+00, 0.0000000E+00, 5.0000002E-02, 0.0000000E+00 +* ELEM 39 circular arc + 3, 26, 25 +* + 3.2100000E+00, 0.0000000E+00, 4.8500000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 40 line segment + 1, 0, 25 +* + 2.6750000E+00, 0.0000000E+00, 5.0000002E-02, 0.0000000E+00 +* ELEM 41 line segment + 1, 25, 9 +* + 2.6750000E+00, 0.0000000E+00, 0.0000000E+00, 3.0888239E-01 +* ELEM 42 line segment + 1, 33, 25 +* + 3.2100000E+00, 6.1776479E-01, -5.3500000E-01, -3.0888239E-01 +* ELEM 43 line segment + 1, 33, 9 +* + 2.6750000E+00, 3.0888239E-01, -5.3500000E-01, 3.0888240E-01 +* ELEM 44 line segment + 1, 0, 9 +* + 2.6250000E+00, 0.0000000E+00, 5.0000002E-02, 0.0000000E+00 +* ELEM 45 circular arc + 3, 10, 9 +* + 2.1400000E+00, 0.0000000E+00, 4.8500000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 46 line segment + 1, 0, 9 +* + 1.6050000E+00, 0.0000000E+00, 5.0000002E-02, 0.0000000E+00 +* ELEM 47 line segment + 1, 9, 1 +* + 1.6050000E+00, 0.0000000E+00, 0.0000000E+00, 3.0888239E-01 +* ELEM 48 line segment + 1, 17, 9 +* + 2.1400000E+00, 6.1776479E-01, -5.3500000E-01, -3.0888239E-01 +* ELEM 49 line segment + 1, 17, 1 +* + 1.6050000E+00, 3.0888239E-01, -5.3500000E-01, 3.0888240E-01 +* ELEM 50 line segment + 1, 0, 1 +* + 1.5550000E+00, 0.0000000E+00, 5.0000002E-02, 0.0000000E+00 +* ELEM 51 circular arc + 3, 2, 1 +* + 1.0700000E+00, 0.0000000E+00, 4.8500000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 52 line segment + 1, 0, 1 +* + 5.3500000E-01, 0.0000000E+00, 5.0000002E-02, 0.0000000E+00 +* ELEM 53 line segment + 1, 193, 1 +* + 5.3500000E-01, 3.0888240E-01, 0.0000000E+00, -3.0888240E-01 +* ELEM 54 line segment + 1, 0, 1 +* + 1.0700000E+00, 6.1776479E-01, -5.3500000E-01, -3.0888239E-01 +* ELEM 55 line segment + 1, 201, 207 +* + 9.2155985E+00, 0.0000000E+00, -1.2059842E-01, 2.0888260E-01 +* ELEM 56 line segment + 1, 0, 207 +* + 9.0950000E+00, 0.0000000E+00, 1.2059842E-01, 0.0000000E+00 +* ELEM 57 line segment + 1, 0, 193 +* + 4.8500000E-01, 0.0000000E+00, 5.0000002E-02, 0.0000000E+00 +* ELEM 58 circular arc + 3, 194, 193 +* + 0.0000000E+00, 0.0000000E+00, 4.8500000E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 59 line segment + 1, 0, 193 +* + 5.3500000E-01, 3.0888239E-01, -1.1497768E-01, -6.6382395E-02 +* ELEM 60 line segment + 1, 205, 161 +* + 8.5600000E+00, 6.1776479E-01, 0.0000000E+00, 5.1776499E-01 +* ELEM 61 line segment + 1, 161, 129 +* + 7.4900000E+00, 6.1776479E-01, 0.0000000E+00, 6.1776479E-01 +* ELEM 62 line segment + 1, 169, 161 +* + 8.0250000E+00, 1.5444120E+00, -5.3500000E-01, -3.0888239E-01 +* ELEM 63 line segment + 1, 206, 161 +* + 8.4733977E+00, 1.2855295E+00, -4.4839763E-01, 2.5888250E-01 +* ELEM 64 line segment + 1, 201, 161 +* + 8.5600000E+00, 1.1355298E+00, -8.6602367E-02, 1.4999970E-01 +* ELEM 65 circle + 2, 162, 161 +* + 8.0250000E+00, 9.2664718E-01, 4.8500000E-01, 0.0000000E+00 +* ELEM 66 line segment + 1, 137, 129 +* + 6.9550000E+00, 1.5444120E+00, -5.3500000E-01, -3.0888239E-01 +* ELEM 67 line segment + 1, 169, 129 +* + 7.4900000E+00, 1.2355296E+00, -5.3500000E-01, 3.0888239E-01 +* ELEM 68 line segment + 1, 129, 97 +* + 6.4200000E+00, 6.1776479E-01, 0.0000000E+00, 6.1776479E-01 +* ELEM 69 circle + 2, 130, 129 +* + 6.9550000E+00, 9.2664718E-01, 4.8500000E-01, 0.0000000E+00 +* ELEM 70 line segment + 1, 105, 97 +* + 5.8850000E+00, 1.5444120E+00, -5.3500000E-01, -3.0888239E-01 +* ELEM 71 line segment + 1, 137, 97 +* + 6.4200000E+00, 1.2355296E+00, -5.3500000E-01, 3.0888239E-01 +* ELEM 72 line segment + 1, 97, 73 +* + 5.3500000E+00, 6.1776479E-01, 0.0000000E+00, 6.1776479E-01 +* ELEM 73 circle + 2, 98, 97 +* + 5.8850000E+00, 9.2664718E-01, 4.8500000E-01, 0.0000000E+00 +* ELEM 74 line segment + 1, 81, 73 +* + 4.8150000E+00, 1.5444120E+00, -5.3500000E-01, -3.0888239E-01 +* ELEM 75 line segment + 1, 105, 73 +* + 5.3500000E+00, 1.2355296E+00, -5.3500000E-01, 3.0888240E-01 +* ELEM 76 line segment + 1, 73, 49 +* + 4.2800000E+00, 6.1776479E-01, 0.0000000E+00, 6.1776479E-01 +* ELEM 77 circle + 2, 74, 73 +* + 4.8150000E+00, 9.2664718E-01, 4.8500000E-01, 0.0000000E+00 +* ELEM 78 line segment + 1, 57, 49 +* + 3.7450000E+00, 1.5444120E+00, -5.3500000E-01, -3.0888239E-01 +* ELEM 79 line segment + 1, 81, 49 +* + 4.2800000E+00, 1.2355296E+00, -5.3500000E-01, 3.0888240E-01 +* ELEM 80 line segment + 1, 49, 33 +* + 3.2100000E+00, 6.1776479E-01, 0.0000000E+00, 6.1776479E-01 +* ELEM 81 circle + 2, 50, 49 +* + 3.7450000E+00, 9.2664718E-01, 4.8500000E-01, 0.0000000E+00 +* ELEM 82 line segment + 1, 0, 33 +* + 2.6750000E+00, 1.5444120E+00, -5.3500000E-01, -3.0888239E-01 +* ELEM 83 line segment + 1, 57, 33 +* + 3.2100000E+00, 1.2355296E+00, -5.3500000E-01, 3.0888239E-01 +* ELEM 84 line segment + 1, 33, 17 +* + 2.1400000E+00, 6.1776479E-01, 0.0000000E+00, 6.1776479E-01 +* ELEM 85 circle + 2, 34, 33 +* + 2.6750000E+00, 9.2664718E-01, 4.8500000E-01, 0.0000000E+00 +* ELEM 86 line segment + 1, 0, 17 +* + 1.1849777E+00, 6.8414718E-01, -1.1497768E-01, -6.6382395E-02 +* ELEM 87 circular arc + 3, 18, 17 +* + 1.6050000E+00, 9.2664718E-01, 4.8500000E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 88 circular arc + 3, 18, 17 +* + 1.6050000E+00, 9.2664718E-01, 4.8500000E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 89 line segment + 1, 0, 17 +* + 2.1400000E+00, 1.2355296E+00, -1.1497768E-01, -6.6382395E-02 +* ELEM 90 line segment + 1, 201, 205 +* + 9.0083977E+00, 3.5888229E-01, -4.4839763E-01, 7.7664749E-01 +* ELEM 91 line segment + 1, 206, 169 +* + 8.0250000E+00, 1.5444120E+00, 0.0000000E+00, 5.1776499E-01 +* ELEM 92 line segment + 1, 169, 137 +* + 6.9550000E+00, 1.5444120E+00, 0.0000000E+00, 6.1776479E-01 +* ELEM 93 line segment + 1, 177, 169 +* + 7.4900000E+00, 2.4710592E+00, -5.3500000E-01, -3.0888239E-01 +* ELEM 94 line segment + 1, 204, 169 +* + 7.9383977E+00, 2.2121767E+00, -4.4839763E-01, 2.5888250E-01 +* ELEM 95 line segment + 1, 201, 169 +* + 8.0250000E+00, 2.0621770E+00, -8.6602367E-02, 1.4999970E-01 +* ELEM 96 circle + 2, 170, 169 +* + 7.4900000E+00, 1.8532944E+00, 4.8500000E-01, 0.0000000E+00 +* ELEM 97 line segment + 1, 145, 137 +* + 6.4200000E+00, 2.4710592E+00, -5.3500000E-01, -3.0888239E-01 +* ELEM 98 line segment + 1, 177, 137 +* + 6.9550000E+00, 2.1621768E+00, -5.3500000E-01, 3.0888239E-01 +* ELEM 99 line segment + 1, 137, 105 +* + 5.8850000E+00, 1.5444120E+00, 0.0000000E+00, 6.1776479E-01 +* ELEM 100 circle + 2, 138, 137 +* + 6.4200000E+00, 1.8532944E+00, 4.8500000E-01, 0.0000000E+00 +* ELEM 101 line segment + 1, 113, 105 +* + 5.3500000E+00, 2.4710592E+00, -5.3500000E-01, -3.0888239E-01 +* ELEM 102 line segment + 1, 145, 105 +* + 5.8850000E+00, 2.1621768E+00, -5.3500000E-01, 3.0888239E-01 +* ELEM 103 line segment + 1, 105, 81 +* + 4.8150000E+00, 1.5444120E+00, 0.0000000E+00, 6.1776479E-01 +* ELEM 104 circle + 2, 106, 105 +* + 5.3500000E+00, 1.8532944E+00, 4.8500000E-01, 0.0000000E+00 +* ELEM 105 line segment + 1, 0, 81 +* + 4.2800000E+00, 2.4710592E+00, -5.3500000E-01, -3.0888239E-01 +* ELEM 106 line segment + 1, 113, 81 +* + 4.8150000E+00, 2.1621768E+00, -5.3500000E-01, 3.0888239E-01 +* ELEM 107 line segment + 1, 81, 57 +* + 3.7450000E+00, 1.5444120E+00, 0.0000000E+00, 6.1776479E-01 +* ELEM 108 circle + 2, 82, 81 +* + 4.2800000E+00, 1.8532944E+00, 4.8500000E-01, 0.0000000E+00 +* ELEM 109 line segment + 1, 0, 57 +* + 2.7899777E+00, 1.6107944E+00, -1.1497768E-01, -6.6382395E-02 +* ELEM 110 circular arc + 3, 58, 57 +* + 3.2100000E+00, 1.8532944E+00, 4.8500000E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 111 circular arc + 3, 58, 57 +* + 3.2100000E+00, 1.8532944E+00, 4.8500000E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 112 line segment + 1, 0, 57 +* + 3.7450000E+00, 2.1621768E+00, -1.1497768E-01, -6.6382395E-02 +* ELEM 113 line segment + 1, 201, 206 +* + 8.4733977E+00, 1.2855295E+00, -4.4839763E-01, 7.7664749E-01 +* ELEM 114 line segment + 1, 204, 177 +* + 7.4900000E+00, 2.4710592E+00, 0.0000000E+00, 5.1776499E-01 +* ELEM 115 line segment + 1, 177, 145 +* + 6.4200000E+00, 2.4710592E+00, 0.0000000E+00, 6.1776479E-01 +* ELEM 116 line segment + 1, 185, 177 +* + 6.9550000E+00, 3.3977063E+00, -5.3500000E-01, -3.0888239E-01 +* ELEM 117 line segment + 1, 203, 177 +* + 7.4033977E+00, 3.1388238E+00, -4.4839763E-01, 2.5888250E-01 +* ELEM 118 line segment + 1, 201, 177 +* + 7.4900000E+00, 2.9888242E+00, -8.6602367E-02, 1.4999970E-01 +* ELEM 119 circle + 2, 178, 177 +* + 6.9550000E+00, 2.7799416E+00, 4.8500000E-01, 0.0000000E+00 +* ELEM 120 line segment + 1, 0, 145 +* + 5.8850000E+00, 3.3977063E+00, -5.3500000E-01, -3.0888239E-01 +* ELEM 121 line segment + 1, 185, 145 +* + 6.4200000E+00, 3.0888239E+00, -5.3500000E-01, 3.0888239E-01 +* ELEM 122 line segment + 1, 145, 113 +* + 5.3500000E+00, 2.4710592E+00, 0.0000000E+00, 6.1776479E-01 +* ELEM 123 circle + 2, 146, 145 +* + 5.8850000E+00, 2.7799416E+00, 4.8500000E-01, 0.0000000E+00 +* ELEM 124 line segment + 1, 0, 113 +* + 4.3949777E+00, 2.5374416E+00, -1.1497768E-01, -6.6382395E-02 +* ELEM 125 circular arc + 3, 114, 113 +* + 4.8150000E+00, 2.7799416E+00, 4.8500000E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 126 circular arc + 3, 114, 113 +* + 4.8150000E+00, 2.7799416E+00, 4.8500000E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 127 line segment + 1, 0, 113 +* + 5.3500000E+00, 3.0888239E+00, -1.1497768E-01, -6.6382395E-02 +* ELEM 128 line segment + 1, 201, 204 +* + 7.9383977E+00, 2.2121767E+00, -4.4839763E-01, 7.7664749E-01 +* ELEM 129 line segment + 1, 0, 185 +* + 5.9999777E+00, 3.4640887E+00, -1.1497768E-01, -6.6382395E-02 +* ELEM 130 circular arc + 3, 186, 185 +* + 6.4200000E+00, 3.7065887E+00, 4.8500000E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 131 circular arc + 3, 186, 185 +* + 6.4200000E+00, 3.7065887E+00, 4.8500000E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 132 line segment + 1, 0, 185 +* + 6.9116988E+00, 3.9904712E+00, -7.1676497E-02, -4.1382445E-02 +* ELEM 133 line segment + 1, 201, 185 +* + 6.9550000E+00, 3.9154713E+00, -4.3301184E-02, 7.4999850E-02 +* ELEM 134 line segment + 1, 203, 185 +* + 6.9550000E+00, 3.3977063E+00, 0.0000000E+00, 5.1776499E-01 +* ELEM 135 line segment + 1, 201, 203 +* + 7.4033977E+00, 3.1388238E+00, -4.4839763E-01, 7.7664749E-01 +* ELEM 136 line segment + 1, 0, 2 +* + 1.5050000E+00, 0.0000000E+00, 5.0000000E-02, 0.0000000E+00 +* ELEM 137 circular arc + 3, 3, 2 +* + 1.0700000E+00, 0.0000000E+00, 4.3500000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 138 line segment + 1, 0, 2 +* + 5.8500000E-01, 0.0000000E+00, 5.0000000E-02, 0.0000000E+00 +* ELEM 139 line segment + 1, 0, 3 +* + 1.4925000E+00, 0.0000000E+00, 1.2500000E-02, 0.0000000E+00 +* ELEM 140 circular arc + 3, 4, 3 +* + 1.0700000E+00, 0.0000000E+00, 4.2250000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 141 line segment + 1, 0, 3 +* + 6.3500000E-01, 0.0000000E+00, 1.2500000E-02, 0.0000000E+00 +* ELEM 142 line segment + 1, 0, 4 +* + 1.4825360E+00, 0.0000000E+00, 9.9640000E-03, 0.0000000E+00 +* ELEM 143 circular arc + 3, 5, 4 +* + 1.0700000E+00, 0.0000000E+00, 4.1253600E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 144 line segment + 1, 0, 4 +* + 6.4750000E-01, 0.0000000E+00, 9.9640000E-03, 0.0000000E+00 +* ELEM 145 line segment + 1, 0, 5 +* + 1.4510840E+00, 0.0000000E+00, 3.1452000E-02, 0.0000000E+00 +* ELEM 146 circular arc + 3, 6, 5 +* + 1.0700000E+00, 0.0000000E+00, 3.8108400E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 147 line segment + 1, 0, 5 +* + 6.5746400E-01, 0.0000000E+00, 3.1452000E-02, 0.0000000E+00 +* ELEM 148 line segment + 1, 0, 6 +* + 1.3787120E+00, 0.0000000E+00, 7.2372000E-02, 0.0000000E+00 +* ELEM 149 circular arc + 3, 7, 6 +* + 1.0700000E+00, 0.0000000E+00, 3.0871200E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 150 line segment + 1, 0, 6 +* + 6.8891600E-01, 0.0000000E+00, 7.2372000E-02, 0.0000000E+00 +* ELEM 151 line segment + 1, 0, 7 +* + 1.1800000E+00, 0.0000000E+00, 1.9871200E-01, 0.0000000E+00 +* ELEM 152 circular arc + 3, 8, 7 +* + 1.0700000E+00, 0.0000000E+00, 1.1000000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 153 line segment + 1, 0, 7 +* + 7.6128800E-01, 0.0000000E+00, 1.9871200E-01, 0.0000000E+00 +* ELEM 154 line segment + 1, 0, 8 +* + 9.6000000E-01, 0.0000000E+00, 2.2000000E-01, 0.0000000E+00 +* ELEM 155 line segment + 1, 0, 10 +* + 2.5750000E+00, 0.0000000E+00, 5.0000000E-02, 0.0000000E+00 +* ELEM 156 circular arc + 3, 11, 10 +* + 2.1400000E+00, 0.0000000E+00, 4.3500000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 157 line segment + 1, 0, 10 +* + 1.6550000E+00, 0.0000000E+00, 5.0000000E-02, 0.0000000E+00 +* ELEM 158 line segment + 1, 0, 11 +* + 2.5625000E+00, 0.0000000E+00, 1.2500000E-02, 0.0000000E+00 +* ELEM 159 circular arc + 3, 12, 11 +* + 2.1400000E+00, 0.0000000E+00, 4.2250000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 160 line segment + 1, 0, 11 +* + 1.7050000E+00, 0.0000000E+00, 1.2500000E-02, 0.0000000E+00 +* ELEM 161 line segment + 1, 0, 12 +* + 2.5525360E+00, 0.0000000E+00, 9.9640000E-03, 0.0000000E+00 +* ELEM 162 circular arc + 3, 13, 12 +* + 2.1400000E+00, 0.0000000E+00, 4.1253600E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 163 line segment + 1, 0, 12 +* + 1.7175000E+00, 0.0000000E+00, 9.9640000E-03, 0.0000000E+00 +* ELEM 164 line segment + 1, 0, 13 +* + 2.5210840E+00, 0.0000000E+00, 3.1452000E-02, 0.0000000E+00 +* ELEM 165 circular arc + 3, 14, 13 +* + 2.1400000E+00, 0.0000000E+00, 3.8108400E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 166 line segment + 1, 0, 13 +* + 1.7274640E+00, 0.0000000E+00, 3.1452000E-02, 0.0000000E+00 +* ELEM 167 line segment + 1, 0, 14 +* + 2.4487120E+00, 0.0000000E+00, 7.2372000E-02, 0.0000000E+00 +* ELEM 168 circular arc + 3, 15, 14 +* + 2.1400000E+00, 0.0000000E+00, 3.0871200E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 169 line segment + 1, 0, 14 +* + 1.7589160E+00, 0.0000000E+00, 7.2372000E-02, 0.0000000E+00 +* ELEM 170 line segment + 1, 0, 15 +* + 2.2500000E+00, 0.0000000E+00, 1.9871200E-01, 0.0000000E+00 +* ELEM 171 circular arc + 3, 16, 15 +* + 2.1400000E+00, 0.0000000E+00, 1.1000000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 172 line segment + 1, 0, 15 +* + 1.8312880E+00, 0.0000000E+00, 1.9871200E-01, 0.0000000E+00 +* ELEM 173 line segment + 1, 0, 16 +* + 2.0300000E+00, 0.0000000E+00, 2.2000000E-01, 0.0000000E+00 +* ELEM 174 line segment + 1, 0, 18 +* + 1.2282790E+00, 7.0914718E-01, -4.3301270E-02, -2.5000000E-02 +* ELEM 175 circular arc + 3, 19, 18 +* + 1.6050000E+00, 9.2664718E-01, 4.3500000E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 176 circular arc + 3, 19, 18 +* + 1.6050000E+00, 9.2664718E-01, 4.3500000E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 177 line segment + 1, 0, 18 +* + 2.0250223E+00, 1.1691472E+00, -4.3301270E-02, -2.5000000E-02 +* ELEM 178 line segment + 1, 0, 19 +* + 1.2391043E+00, 7.1539718E-01, -1.0825318E-02, -6.2500000E-03 +* ELEM 179 circular arc + 3, 20, 19 +* + 1.6050000E+00, 9.2664718E-01, 4.2250000E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 180 circular arc + 3, 20, 19 +* + 1.6050000E+00, 9.2664718E-01, 4.2250000E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 181 line segment + 1, 0, 19 +* + 1.9817211E+00, 1.1441472E+00, -1.0825318E-02, -6.2500000E-03 +* ELEM 182 line segment + 1, 0, 20 +* + 1.2477333E+00, 7.2037918E-01, -8.6290771E-03, -4.9820000E-03 +* ELEM 183 circular arc + 3, 21, 20 +* + 1.6050000E+00, 9.2664718E-01, 4.1253600E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 184 circular arc + 3, 21, 20 +* + 1.6050000E+00, 9.2664718E-01, 4.1253600E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 185 line segment + 1, 0, 20 +* + 1.9708957E+00, 1.1378972E+00, -8.6290771E-03, -4.9820000E-03 +* ELEM 186 line segment + 1, 0, 21 +* + 1.2749716E+00, 7.3610518E-01, -2.7238231E-02, -1.5726000E-02 +* ELEM 187 circular arc + 3, 22, 21 +* + 1.6050000E+00, 9.2664718E-01, 3.8108400E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 188 circular arc + 3, 22, 21 +* + 1.6050000E+00, 9.2664718E-01, 3.8108400E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 189 line segment + 1, 0, 21 +* + 1.9622667E+00, 1.1329152E+00, -2.7238231E-02, -1.5726000E-02 +* ELEM 190 line segment + 1, 0, 22 +* + 1.3376476E+00, 7.7229118E-01, -6.2675991E-02, -3.6186000E-02 +* ELEM 191 circular arc + 3, 23, 22 +* + 1.6050000E+00, 9.2664718E-01, 3.0871200E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 192 circular arc + 3, 23, 22 +* + 1.6050000E+00, 9.2664718E-01, 3.0871200E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 193 line segment + 1, 0, 22 +* + 1.9350284E+00, 1.1171892E+00, -6.2675991E-02, -3.6186000E-02 +* ELEM 194 line segment + 1, 0, 23 +* + 1.5097372E+00, 8.7164718E-01, -1.7208964E-01, -9.9356000E-02 +* ELEM 195 circular arc + 3, 24, 23 +* + 1.6050000E+00, 9.2664718E-01, 1.1000000E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 196 circular arc + 3, 24, 23 +* + 1.6050000E+00, 9.2664718E-01, 1.1000000E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 197 line segment + 1, 0, 23 +* + 1.8723524E+00, 1.0810032E+00, -1.7208964E-01, -9.9356000E-02 +* ELEM 198 line segment + 1, 0, 24 +* + 1.7002628E+00, 9.8164718E-01, -1.9052559E-01, -1.1000000E-01 +* ELEM 199 line segment + 1, 0, 26 +* + 3.6450000E+00, 0.0000000E+00, 5.0000000E-02, 0.0000000E+00 +* ELEM 200 circular arc + 3, 27, 26 +* + 3.2100000E+00, 0.0000000E+00, 4.3500000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 201 line segment + 1, 0, 26 +* + 2.7250000E+00, 0.0000000E+00, 5.0000000E-02, 0.0000000E+00 +* ELEM 202 line segment + 1, 0, 27 +* + 3.6325000E+00, 0.0000000E+00, 1.2500000E-02, 0.0000000E+00 +* ELEM 203 circular arc + 3, 28, 27 +* + 3.2100000E+00, 0.0000000E+00, 4.2250000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 204 line segment + 1, 0, 27 +* + 2.7750000E+00, 0.0000000E+00, 1.2500000E-02, 0.0000000E+00 +* ELEM 205 line segment + 1, 0, 28 +* + 3.6225360E+00, 0.0000000E+00, 9.9640000E-03, 0.0000000E+00 +* ELEM 206 circular arc + 3, 29, 28 +* + 3.2100000E+00, 0.0000000E+00, 4.1253600E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 207 line segment + 1, 0, 28 +* + 2.7875000E+00, 0.0000000E+00, 9.9640000E-03, 0.0000000E+00 +* ELEM 208 line segment + 1, 0, 29 +* + 3.5910840E+00, 0.0000000E+00, 3.1452000E-02, 0.0000000E+00 +* ELEM 209 circular arc + 3, 30, 29 +* + 3.2100000E+00, 0.0000000E+00, 3.8108400E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 210 line segment + 1, 0, 29 +* + 2.7974640E+00, 0.0000000E+00, 3.1452000E-02, 0.0000000E+00 +* ELEM 211 line segment + 1, 0, 30 +* + 3.5187120E+00, 0.0000000E+00, 7.2372000E-02, 0.0000000E+00 +* ELEM 212 circular arc + 3, 31, 30 +* + 3.2100000E+00, 0.0000000E+00, 3.0871200E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 213 line segment + 1, 0, 30 +* + 2.8289160E+00, 0.0000000E+00, 7.2372000E-02, 0.0000000E+00 +* ELEM 214 line segment + 1, 0, 31 +* + 3.3200000E+00, 0.0000000E+00, 1.9871200E-01, 0.0000000E+00 +* ELEM 215 circular arc + 3, 32, 31 +* + 3.2100000E+00, 0.0000000E+00, 1.1000000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 216 line segment + 1, 0, 31 +* + 2.9012880E+00, 0.0000000E+00, 1.9871200E-01, 0.0000000E+00 +* ELEM 217 line segment + 1, 0, 32 +* + 3.1000000E+00, 0.0000000E+00, 2.2000000E-01, 0.0000000E+00 +* ELEM 218 circle + 2, 35, 34 +* + 2.6750000E+00, 9.2664718E-01, 4.3500000E-01, 0.0000000E+00 +* ELEM 219 circle + 2, 36, 35 +* + 2.6750000E+00, 9.2664718E-01, 4.2250000E-01, 0.0000000E+00 +* ELEM 220 circle + 2, 37, 36 +* + 2.6750000E+00, 9.2664718E-01, 4.1253600E-01, 0.0000000E+00 +* ELEM 221 circle + 2, 38, 37 +* + 2.6750000E+00, 9.2664718E-01, 3.8108400E-01, 0.0000000E+00 +* ELEM 222 circle + 2, 39, 38 +* + 2.6750000E+00, 9.2664718E-01, 3.0871200E-01, 0.0000000E+00 +* ELEM 223 circle + 2, 40, 39 +* + 2.6750000E+00, 9.2664718E-01, 1.1000000E-01, 0.0000000E+00 +* ELEM 224 line segment + 1, 0, 42 +* + 4.7150000E+00, 0.0000000E+00, 5.0000000E-02, 0.0000000E+00 +* ELEM 225 circular arc + 3, 43, 42 +* + 4.2800000E+00, 0.0000000E+00, 4.3500000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 226 line segment + 1, 0, 42 +* + 3.7950000E+00, 0.0000000E+00, 5.0000000E-02, 0.0000000E+00 +* ELEM 227 line segment + 1, 0, 43 +* + 4.7025000E+00, 0.0000000E+00, 1.2500000E-02, 0.0000000E+00 +* ELEM 228 circular arc + 3, 44, 43 +* + 4.2800000E+00, 0.0000000E+00, 4.2250000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 229 line segment + 1, 0, 43 +* + 3.8450000E+00, 0.0000000E+00, 1.2500000E-02, 0.0000000E+00 +* ELEM 230 line segment + 1, 0, 44 +* + 4.6925360E+00, 0.0000000E+00, 9.9640000E-03, 0.0000000E+00 +* ELEM 231 circular arc + 3, 45, 44 +* + 4.2800000E+00, 0.0000000E+00, 4.1253600E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 232 line segment + 1, 0, 44 +* + 3.8575000E+00, 0.0000000E+00, 9.9640000E-03, 0.0000000E+00 +* ELEM 233 line segment + 1, 0, 45 +* + 4.6610840E+00, 0.0000000E+00, 3.1452000E-02, 0.0000000E+00 +* ELEM 234 circular arc + 3, 46, 45 +* + 4.2800000E+00, 0.0000000E+00, 3.8108400E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 235 line segment + 1, 0, 45 +* + 3.8674640E+00, 0.0000000E+00, 3.1452000E-02, 0.0000000E+00 +* ELEM 236 line segment + 1, 0, 46 +* + 4.5887120E+00, 0.0000000E+00, 7.2372000E-02, 0.0000000E+00 +* ELEM 237 circular arc + 3, 47, 46 +* + 4.2800000E+00, 0.0000000E+00, 3.0871200E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 238 line segment + 1, 0, 46 +* + 3.8989160E+00, 0.0000000E+00, 7.2372000E-02, 0.0000000E+00 +* ELEM 239 line segment + 1, 0, 47 +* + 4.3900000E+00, 0.0000000E+00, 1.9871200E-01, 0.0000000E+00 +* ELEM 240 circular arc + 3, 48, 47 +* + 4.2800000E+00, 0.0000000E+00, 1.1000000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 241 line segment + 1, 0, 47 +* + 3.9712880E+00, 0.0000000E+00, 1.9871200E-01, 0.0000000E+00 +* ELEM 242 line segment + 1, 0, 48 +* + 4.1700000E+00, 0.0000000E+00, 2.2000000E-01, 0.0000000E+00 +* ELEM 243 circle + 2, 51, 50 +* + 3.7450000E+00, 9.2664718E-01, 4.3500000E-01, 0.0000000E+00 +* ELEM 244 circle + 2, 52, 51 +* + 3.7450000E+00, 9.2664718E-01, 4.2250000E-01, 0.0000000E+00 +* ELEM 245 circle + 2, 53, 52 +* + 3.7450000E+00, 9.2664718E-01, 4.1253600E-01, 0.0000000E+00 +* ELEM 246 circle + 2, 54, 53 +* + 3.7450000E+00, 9.2664718E-01, 3.8108400E-01, 0.0000000E+00 +* ELEM 247 circle + 2, 55, 54 +* + 3.7450000E+00, 9.2664718E-01, 3.0871200E-01, 0.0000000E+00 +* ELEM 248 circle + 2, 56, 55 +* + 3.7450000E+00, 9.2664718E-01, 1.1000000E-01, 0.0000000E+00 +* ELEM 249 line segment + 1, 0, 58 +* + 2.8332790E+00, 1.6357944E+00, -4.3301270E-02, -2.5000000E-02 +* ELEM 250 circular arc + 3, 59, 58 +* + 3.2100000E+00, 1.8532944E+00, 4.3500000E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 251 circular arc + 3, 59, 58 +* + 3.2100000E+00, 1.8532944E+00, 4.3500000E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 252 line segment + 1, 0, 58 +* + 3.6300223E+00, 2.0957944E+00, -4.3301270E-02, -2.5000000E-02 +* ELEM 253 line segment + 1, 0, 59 +* + 2.8441043E+00, 1.6420444E+00, -1.0825318E-02, -6.2500000E-03 +* ELEM 254 circular arc + 3, 60, 59 +* + 3.2100000E+00, 1.8532944E+00, 4.2250000E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 255 circular arc + 3, 60, 59 +* + 3.2100000E+00, 1.8532944E+00, 4.2250000E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 256 line segment + 1, 0, 59 +* + 3.5867211E+00, 2.0707944E+00, -1.0825318E-02, -6.2500000E-03 +* ELEM 257 line segment + 1, 0, 60 +* + 2.8527334E+00, 1.6470264E+00, -8.6290771E-03, -4.9820000E-03 +* ELEM 258 circular arc + 3, 61, 60 +* + 3.2100000E+00, 1.8532944E+00, 4.1253600E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 259 circular arc + 3, 61, 60 +* + 3.2100000E+00, 1.8532944E+00, 4.1253600E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 260 line segment + 1, 0, 60 +* + 3.5758957E+00, 2.0645444E+00, -8.6290771E-03, -4.9820000E-03 +* ELEM 261 line segment + 1, 0, 61 +* + 2.8799716E+00, 1.6627524E+00, -2.7238231E-02, -1.5726000E-02 +* ELEM 262 circular arc + 3, 62, 61 +* + 3.2100000E+00, 1.8532944E+00, 3.8108400E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 263 circular arc + 3, 62, 61 +* + 3.2100000E+00, 1.8532944E+00, 3.8108400E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 264 line segment + 1, 0, 61 +* + 3.5672667E+00, 2.0595624E+00, -2.7238231E-02, -1.5726000E-02 +* ELEM 265 line segment + 1, 0, 62 +* + 2.9426476E+00, 1.6989384E+00, -6.2675991E-02, -3.6186000E-02 +* ELEM 266 circular arc + 3, 63, 62 +* + 3.2100000E+00, 1.8532944E+00, 3.0871200E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 267 circular arc + 3, 63, 62 +* + 3.2100000E+00, 1.8532944E+00, 3.0871200E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 268 line segment + 1, 0, 62 +* + 3.5400284E+00, 2.0438364E+00, -6.2675991E-02, -3.6186000E-02 +* ELEM 269 line segment + 1, 0, 63 +* + 3.1147372E+00, 1.7982944E+00, -1.7208964E-01, -9.9356000E-02 +* ELEM 270 circular arc + 3, 64, 63 +* + 3.2100000E+00, 1.8532944E+00, 1.1000000E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 271 circular arc + 3, 64, 63 +* + 3.2100000E+00, 1.8532944E+00, 1.1000000E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 272 line segment + 1, 0, 63 +* + 3.4773524E+00, 2.0076504E+00, -1.7208964E-01, -9.9356000E-02 +* ELEM 273 line segment + 1, 0, 64 +* + 3.3052628E+00, 1.9082944E+00, -1.9052559E-01, -1.1000000E-01 +* ELEM 274 line segment + 1, 0, 66 +* + 5.7850000E+00, 0.0000000E+00, 5.0000000E-02, 0.0000000E+00 +* ELEM 275 circular arc + 3, 67, 66 +* + 5.3500000E+00, 0.0000000E+00, 4.3500000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 276 line segment + 1, 0, 66 +* + 4.8650000E+00, 0.0000000E+00, 5.0000000E-02, 0.0000000E+00 +* ELEM 277 line segment + 1, 0, 67 +* + 5.7725000E+00, 0.0000000E+00, 1.2500000E-02, 0.0000000E+00 +* ELEM 278 circular arc + 3, 68, 67 +* + 5.3500000E+00, 0.0000000E+00, 4.2250000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 279 line segment + 1, 0, 67 +* + 4.9150000E+00, 0.0000000E+00, 1.2500000E-02, 0.0000000E+00 +* ELEM 280 line segment + 1, 0, 68 +* + 5.7625360E+00, 0.0000000E+00, 9.9640000E-03, 0.0000000E+00 +* ELEM 281 circular arc + 3, 69, 68 +* + 5.3500000E+00, 0.0000000E+00, 4.1253600E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 282 line segment + 1, 0, 68 +* + 4.9275000E+00, 0.0000000E+00, 9.9640000E-03, 0.0000000E+00 +* ELEM 283 line segment + 1, 0, 69 +* + 5.7310840E+00, 0.0000000E+00, 3.1452000E-02, 0.0000000E+00 +* ELEM 284 circular arc + 3, 70, 69 +* + 5.3500000E+00, 0.0000000E+00, 3.8108400E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 285 line segment + 1, 0, 69 +* + 4.9374640E+00, 0.0000000E+00, 3.1452000E-02, 0.0000000E+00 +* ELEM 286 line segment + 1, 0, 70 +* + 5.6587120E+00, 0.0000000E+00, 7.2372000E-02, 0.0000000E+00 +* ELEM 287 circular arc + 3, 71, 70 +* + 5.3500000E+00, 0.0000000E+00, 3.0871200E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 288 line segment + 1, 0, 70 +* + 4.9689160E+00, 0.0000000E+00, 7.2372000E-02, 0.0000000E+00 +* ELEM 289 line segment + 1, 0, 71 +* + 5.4600000E+00, 0.0000000E+00, 1.9871200E-01, 0.0000000E+00 +* ELEM 290 circular arc + 3, 72, 71 +* + 5.3500000E+00, 0.0000000E+00, 1.1000000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 291 line segment + 1, 0, 71 +* + 5.0412880E+00, 0.0000000E+00, 1.9871200E-01, 0.0000000E+00 +* ELEM 292 line segment + 1, 0, 72 +* + 5.2400000E+00, 0.0000000E+00, 2.2000000E-01, 0.0000000E+00 +* ELEM 293 circle + 2, 75, 74 +* + 4.8150000E+00, 9.2664718E-01, 4.3500000E-01, 0.0000000E+00 +* ELEM 294 circle + 2, 76, 75 +* + 4.8150000E+00, 9.2664718E-01, 4.2250000E-01, 0.0000000E+00 +* ELEM 295 circle + 2, 77, 76 +* + 4.8150000E+00, 9.2664718E-01, 4.1253600E-01, 0.0000000E+00 +* ELEM 296 circle + 2, 78, 77 +* + 4.8150000E+00, 9.2664718E-01, 3.8108400E-01, 0.0000000E+00 +* ELEM 297 circle + 2, 79, 78 +* + 4.8150000E+00, 9.2664718E-01, 3.0871200E-01, 0.0000000E+00 +* ELEM 298 circle + 2, 80, 79 +* + 4.8150000E+00, 9.2664718E-01, 1.1000000E-01, 0.0000000E+00 +* ELEM 299 circle + 2, 83, 82 +* + 4.2800000E+00, 1.8532944E+00, 4.3500000E-01, 0.0000000E+00 +* ELEM 300 circle + 2, 84, 83 +* + 4.2800000E+00, 1.8532944E+00, 4.2250000E-01, 0.0000000E+00 +* ELEM 301 circle + 2, 85, 84 +* + 4.2800000E+00, 1.8532944E+00, 4.1253600E-01, 0.0000000E+00 +* ELEM 302 circle + 2, 86, 85 +* + 4.2800000E+00, 1.8532944E+00, 3.8108400E-01, 0.0000000E+00 +* ELEM 303 circle + 2, 87, 86 +* + 4.2800000E+00, 1.8532944E+00, 3.0871200E-01, 0.0000000E+00 +* ELEM 304 circle + 2, 88, 87 +* + 4.2800000E+00, 1.8532944E+00, 1.1000000E-01, 0.0000000E+00 +* ELEM 305 line segment + 1, 0, 90 +* + 6.8550000E+00, 0.0000000E+00, 5.0000000E-02, 0.0000000E+00 +* ELEM 306 circular arc + 3, 91, 90 +* + 6.4200000E+00, 0.0000000E+00, 4.3500000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 307 line segment + 1, 0, 90 +* + 5.9350000E+00, 0.0000000E+00, 5.0000000E-02, 0.0000000E+00 +* ELEM 308 line segment + 1, 0, 91 +* + 6.8425000E+00, 0.0000000E+00, 1.2500000E-02, 0.0000000E+00 +* ELEM 309 circular arc + 3, 92, 91 +* + 6.4200000E+00, 0.0000000E+00, 4.2250000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 310 line segment + 1, 0, 91 +* + 5.9850000E+00, 0.0000000E+00, 1.2500000E-02, 0.0000000E+00 +* ELEM 311 line segment + 1, 0, 92 +* + 6.8325360E+00, 0.0000000E+00, 9.9640000E-03, 0.0000000E+00 +* ELEM 312 circular arc + 3, 93, 92 +* + 6.4200000E+00, 0.0000000E+00, 4.1253600E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 313 line segment + 1, 0, 92 +* + 5.9975000E+00, 0.0000000E+00, 9.9640000E-03, 0.0000000E+00 +* ELEM 314 line segment + 1, 0, 93 +* + 6.8010840E+00, 0.0000000E+00, 3.1452000E-02, 0.0000000E+00 +* ELEM 315 circular arc + 3, 94, 93 +* + 6.4200000E+00, 0.0000000E+00, 3.8108400E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 316 line segment + 1, 0, 93 +* + 6.0074640E+00, 0.0000000E+00, 3.1452000E-02, 0.0000000E+00 +* ELEM 317 line segment + 1, 0, 94 +* + 6.7287120E+00, 0.0000000E+00, 7.2372000E-02, 0.0000000E+00 +* ELEM 318 circular arc + 3, 95, 94 +* + 6.4200000E+00, 0.0000000E+00, 3.0871200E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 319 line segment + 1, 0, 94 +* + 6.0389160E+00, 0.0000000E+00, 7.2372000E-02, 0.0000000E+00 +* ELEM 320 line segment + 1, 0, 95 +* + 6.5300000E+00, 0.0000000E+00, 1.9871200E-01, 0.0000000E+00 +* ELEM 321 circular arc + 3, 96, 95 +* + 6.4200000E+00, 0.0000000E+00, 1.1000000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 322 line segment + 1, 0, 95 +* + 6.1112880E+00, 0.0000000E+00, 1.9871200E-01, 0.0000000E+00 +* ELEM 323 line segment + 1, 0, 96 +* + 6.3100000E+00, 0.0000000E+00, 2.2000000E-01, 0.0000000E+00 +* ELEM 324 circle + 2, 99, 98 +* + 5.8850000E+00, 9.2664718E-01, 4.3500000E-01, 0.0000000E+00 +* ELEM 325 circle + 2, 100, 99 +* + 5.8850000E+00, 9.2664718E-01, 4.2250000E-01, 0.0000000E+00 +* ELEM 326 circle + 2, 101, 100 +* + 5.8850000E+00, 9.2664718E-01, 4.1253600E-01, 0.0000000E+00 +* ELEM 327 circle + 2, 102, 101 +* + 5.8850000E+00, 9.2664718E-01, 3.8108400E-01, 0.0000000E+00 +* ELEM 328 circle + 2, 103, 102 +* + 5.8850000E+00, 9.2664718E-01, 3.0871200E-01, 0.0000000E+00 +* ELEM 329 circle + 2, 104, 103 +* + 5.8850000E+00, 9.2664718E-01, 1.1000000E-01, 0.0000000E+00 +* ELEM 330 circle + 2, 107, 106 +* + 5.3500000E+00, 1.8532944E+00, 4.3500000E-01, 0.0000000E+00 +* ELEM 331 circle + 2, 108, 107 +* + 5.3500000E+00, 1.8532944E+00, 4.2250000E-01, 0.0000000E+00 +* ELEM 332 circle + 2, 109, 108 +* + 5.3500000E+00, 1.8532944E+00, 4.1253600E-01, 0.0000000E+00 +* ELEM 333 circle + 2, 110, 109 +* + 5.3500000E+00, 1.8532944E+00, 3.8108400E-01, 0.0000000E+00 +* ELEM 334 circle + 2, 111, 110 +* + 5.3500000E+00, 1.8532944E+00, 3.0871200E-01, 0.0000000E+00 +* ELEM 335 circle + 2, 112, 111 +* + 5.3500000E+00, 1.8532944E+00, 1.1000000E-01, 0.0000000E+00 +* ELEM 336 line segment + 1, 0, 114 +* + 4.4382790E+00, 2.5624416E+00, -4.3301270E-02, -2.5000000E-02 +* ELEM 337 circular arc + 3, 115, 114 +* + 4.8150000E+00, 2.7799416E+00, 4.3500000E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 338 circular arc + 3, 115, 114 +* + 4.8150000E+00, 2.7799416E+00, 4.3500000E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 339 line segment + 1, 0, 114 +* + 5.2350223E+00, 3.0224416E+00, -4.3301270E-02, -2.5000000E-02 +* ELEM 340 line segment + 1, 0, 115 +* + 4.4491043E+00, 2.5686916E+00, -1.0825318E-02, -6.2500000E-03 +* ELEM 341 circular arc + 3, 116, 115 +* + 4.8150000E+00, 2.7799416E+00, 4.2250000E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 342 circular arc + 3, 116, 115 +* + 4.8150000E+00, 2.7799416E+00, 4.2250000E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 343 line segment + 1, 0, 115 +* + 5.1917211E+00, 2.9974416E+00, -1.0825318E-02, -6.2500000E-03 +* ELEM 344 line segment + 1, 0, 116 +* + 4.4577334E+00, 2.5736736E+00, -8.6290771E-03, -4.9820000E-03 +* ELEM 345 circular arc + 3, 117, 116 +* + 4.8150000E+00, 2.7799416E+00, 4.1253600E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 346 circular arc + 3, 117, 116 +* + 4.8150000E+00, 2.7799416E+00, 4.1253600E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 347 line segment + 1, 0, 116 +* + 5.1808957E+00, 2.9911916E+00, -8.6290771E-03, -4.9820000E-03 +* ELEM 348 line segment + 1, 0, 117 +* + 4.4849716E+00, 2.5893996E+00, -2.7238231E-02, -1.5726000E-02 +* ELEM 349 circular arc + 3, 118, 117 +* + 4.8150000E+00, 2.7799416E+00, 3.8108400E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 350 circular arc + 3, 118, 117 +* + 4.8150000E+00, 2.7799416E+00, 3.8108400E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 351 line segment + 1, 0, 117 +* + 5.1722667E+00, 2.9862096E+00, -2.7238231E-02, -1.5726000E-02 +* ELEM 352 line segment + 1, 0, 118 +* + 4.5476476E+00, 2.6255856E+00, -6.2675991E-02, -3.6186000E-02 +* ELEM 353 circular arc + 3, 119, 118 +* + 4.8150000E+00, 2.7799416E+00, 3.0871200E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 354 circular arc + 3, 119, 118 +* + 4.8150000E+00, 2.7799416E+00, 3.0871200E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 355 line segment + 1, 0, 118 +* + 5.1450284E+00, 2.9704836E+00, -6.2675991E-02, -3.6186000E-02 +* ELEM 356 line segment + 1, 0, 119 +* + 4.7197372E+00, 2.7249416E+00, -1.7208964E-01, -9.9356000E-02 +* ELEM 357 circular arc + 3, 120, 119 +* + 4.8150000E+00, 2.7799416E+00, 1.1000000E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 358 circular arc + 3, 120, 119 +* + 4.8150000E+00, 2.7799416E+00, 1.1000000E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 359 line segment + 1, 0, 119 +* + 5.0823524E+00, 2.9342976E+00, -1.7208964E-01, -9.9356000E-02 +* ELEM 360 line segment + 1, 0, 120 +* + 4.9102628E+00, 2.8349416E+00, -1.9052559E-01, -1.1000000E-01 +* ELEM 361 line segment + 1, 0, 122 +* + 7.9250000E+00, 0.0000000E+00, 5.0000000E-02, 0.0000000E+00 +* ELEM 362 circular arc + 3, 123, 122 +* + 7.4900000E+00, 0.0000000E+00, 4.3500000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 363 line segment + 1, 0, 122 +* + 7.0050000E+00, 0.0000000E+00, 5.0000000E-02, 0.0000000E+00 +* ELEM 364 line segment + 1, 0, 123 +* + 7.9125000E+00, 0.0000000E+00, 1.2500000E-02, 0.0000000E+00 +* ELEM 365 circular arc + 3, 124, 123 +* + 7.4900000E+00, 0.0000000E+00, 4.2250000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 366 line segment + 1, 0, 123 +* + 7.0550000E+00, 0.0000000E+00, 1.2500000E-02, 0.0000000E+00 +* ELEM 367 line segment + 1, 0, 124 +* + 7.9025360E+00, 0.0000000E+00, 9.9640000E-03, 0.0000000E+00 +* ELEM 368 circular arc + 3, 125, 124 +* + 7.4900000E+00, 0.0000000E+00, 4.1253600E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 369 line segment + 1, 0, 124 +* + 7.0675000E+00, 0.0000000E+00, 9.9640000E-03, 0.0000000E+00 +* ELEM 370 line segment + 1, 0, 125 +* + 7.8710840E+00, 0.0000000E+00, 3.1452000E-02, 0.0000000E+00 +* ELEM 371 circular arc + 3, 126, 125 +* + 7.4900000E+00, 0.0000000E+00, 3.8108400E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 372 line segment + 1, 0, 125 +* + 7.0774640E+00, 0.0000000E+00, 3.1452000E-02, 0.0000000E+00 +* ELEM 373 line segment + 1, 0, 126 +* + 7.7987120E+00, 0.0000000E+00, 7.2372000E-02, 0.0000000E+00 +* ELEM 374 circular arc + 3, 127, 126 +* + 7.4900000E+00, 0.0000000E+00, 3.0871200E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 375 line segment + 1, 0, 126 +* + 7.1089160E+00, 0.0000000E+00, 7.2372000E-02, 0.0000000E+00 +* ELEM 376 line segment + 1, 0, 127 +* + 7.6000000E+00, 0.0000000E+00, 1.9871200E-01, 0.0000000E+00 +* ELEM 377 circular arc + 3, 128, 127 +* + 7.4900000E+00, 0.0000000E+00, 1.1000000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 378 line segment + 1, 0, 127 +* + 7.1812880E+00, 0.0000000E+00, 1.9871200E-01, 0.0000000E+00 +* ELEM 379 line segment + 1, 0, 128 +* + 7.3800000E+00, 0.0000000E+00, 2.2000000E-01, 0.0000000E+00 +* ELEM 380 circle + 2, 131, 130 +* + 6.9550000E+00, 9.2664718E-01, 4.3500000E-01, 0.0000000E+00 +* ELEM 381 circle + 2, 132, 131 +* + 6.9550000E+00, 9.2664718E-01, 4.2250000E-01, 0.0000000E+00 +* ELEM 382 circle + 2, 133, 132 +* + 6.9550000E+00, 9.2664718E-01, 4.1253600E-01, 0.0000000E+00 +* ELEM 383 circle + 2, 134, 133 +* + 6.9550000E+00, 9.2664718E-01, 3.8108400E-01, 0.0000000E+00 +* ELEM 384 circle + 2, 135, 134 +* + 6.9550000E+00, 9.2664718E-01, 3.0871200E-01, 0.0000000E+00 +* ELEM 385 circle + 2, 136, 135 +* + 6.9550000E+00, 9.2664718E-01, 1.1000000E-01, 0.0000000E+00 +* ELEM 386 circle + 2, 139, 138 +* + 6.4200000E+00, 1.8532944E+00, 4.3500000E-01, 0.0000000E+00 +* ELEM 387 circle + 2, 140, 139 +* + 6.4200000E+00, 1.8532944E+00, 4.2250000E-01, 0.0000000E+00 +* ELEM 388 circle + 2, 141, 140 +* + 6.4200000E+00, 1.8532944E+00, 4.1253600E-01, 0.0000000E+00 +* ELEM 389 circle + 2, 142, 141 +* + 6.4200000E+00, 1.8532944E+00, 3.8108400E-01, 0.0000000E+00 +* ELEM 390 circle + 2, 143, 142 +* + 6.4200000E+00, 1.8532944E+00, 3.0871200E-01, 0.0000000E+00 +* ELEM 391 circle + 2, 144, 143 +* + 6.4200000E+00, 1.8532944E+00, 1.1000000E-01, 0.0000000E+00 +* ELEM 392 circle + 2, 147, 146 +* + 5.8850000E+00, 2.7799416E+00, 4.3500000E-01, 0.0000000E+00 +* ELEM 393 circle + 2, 148, 147 +* + 5.8850000E+00, 2.7799416E+00, 4.2250000E-01, 0.0000000E+00 +* ELEM 394 circle + 2, 149, 148 +* + 5.8850000E+00, 2.7799416E+00, 4.1253600E-01, 0.0000000E+00 +* ELEM 395 circle + 2, 150, 149 +* + 5.8850000E+00, 2.7799416E+00, 3.8108400E-01, 0.0000000E+00 +* ELEM 396 circle + 2, 151, 150 +* + 5.8850000E+00, 2.7799416E+00, 3.0871200E-01, 0.0000000E+00 +* ELEM 397 circle + 2, 152, 151 +* + 5.8850000E+00, 2.7799416E+00, 1.1000000E-01, 0.0000000E+00 +* ELEM 398 line segment + 1, 0, 154 +* + 8.9950000E+00, 0.0000000E+00, 5.0000000E-02, 0.0000000E+00 +* ELEM 399 circular arc + 3, 155, 154 +* + 8.5600000E+00, 0.0000000E+00, 4.3500000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 400 line segment + 1, 0, 154 +* + 8.0750000E+00, 0.0000000E+00, 5.0000000E-02, 0.0000000E+00 +* ELEM 401 line segment + 1, 0, 155 +* + 8.9825000E+00, 0.0000000E+00, 1.2500000E-02, 0.0000000E+00 +* ELEM 402 circular arc + 3, 156, 155 +* + 8.5600000E+00, 0.0000000E+00, 4.2250000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 403 line segment + 1, 0, 155 +* + 8.1250000E+00, 0.0000000E+00, 1.2500000E-02, 0.0000000E+00 +* ELEM 404 line segment + 1, 0, 156 +* + 8.9725360E+00, 0.0000000E+00, 9.9640000E-03, 0.0000000E+00 +* ELEM 405 circular arc + 3, 157, 156 +* + 8.5600000E+00, 0.0000000E+00, 4.1253600E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 406 line segment + 1, 0, 156 +* + 8.1375000E+00, 0.0000000E+00, 9.9640000E-03, 0.0000000E+00 +* ELEM 407 line segment + 1, 0, 157 +* + 8.9410840E+00, 0.0000000E+00, 3.1452000E-02, 0.0000000E+00 +* ELEM 408 circular arc + 3, 158, 157 +* + 8.5600000E+00, 0.0000000E+00, 3.8108400E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 409 line segment + 1, 0, 157 +* + 8.1474640E+00, 0.0000000E+00, 3.1452000E-02, 0.0000000E+00 +* ELEM 410 line segment + 1, 0, 158 +* + 8.8687120E+00, 0.0000000E+00, 7.2372000E-02, 0.0000000E+00 +* ELEM 411 circular arc + 3, 159, 158 +* + 8.5600000E+00, 0.0000000E+00, 3.0871200E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 412 line segment + 1, 0, 158 +* + 8.1789160E+00, 0.0000000E+00, 7.2372000E-02, 0.0000000E+00 +* ELEM 413 line segment + 1, 0, 159 +* + 8.6700000E+00, 0.0000000E+00, 1.9871200E-01, 0.0000000E+00 +* ELEM 414 circular arc + 3, 160, 159 +* + 8.5600000E+00, 0.0000000E+00, 1.1000000E-01, 0.0000000E+00, 1.8000000E+02 +* ELEM 415 line segment + 1, 0, 159 +* + 8.2512880E+00, 0.0000000E+00, 1.9871200E-01, 0.0000000E+00 +* ELEM 416 line segment + 1, 0, 160 +* + 8.4500000E+00, 0.0000000E+00, 2.2000000E-01, 0.0000000E+00 +* ELEM 417 circle + 2, 163, 162 +* + 8.0250000E+00, 9.2664718E-01, 4.3500000E-01, 0.0000000E+00 +* ELEM 418 circle + 2, 164, 163 +* + 8.0250000E+00, 9.2664718E-01, 4.2250000E-01, 0.0000000E+00 +* ELEM 419 circle + 2, 165, 164 +* + 8.0250000E+00, 9.2664718E-01, 4.1253600E-01, 0.0000000E+00 +* ELEM 420 circle + 2, 166, 165 +* + 8.0250000E+00, 9.2664718E-01, 3.8108400E-01, 0.0000000E+00 +* ELEM 421 circle + 2, 167, 166 +* + 8.0250000E+00, 9.2664718E-01, 3.0871200E-01, 0.0000000E+00 +* ELEM 422 circle + 2, 168, 167 +* + 8.0250000E+00, 9.2664718E-01, 1.1000000E-01, 0.0000000E+00 +* ELEM 423 circle + 2, 171, 170 +* + 7.4900000E+00, 1.8532944E+00, 4.3500000E-01, 0.0000000E+00 +* ELEM 424 circle + 2, 172, 171 +* + 7.4900000E+00, 1.8532944E+00, 4.2250000E-01, 0.0000000E+00 +* ELEM 425 circle + 2, 173, 172 +* + 7.4900000E+00, 1.8532944E+00, 4.1253600E-01, 0.0000000E+00 +* ELEM 426 circle + 2, 174, 173 +* + 7.4900000E+00, 1.8532944E+00, 3.8108400E-01, 0.0000000E+00 +* ELEM 427 circle + 2, 175, 174 +* + 7.4900000E+00, 1.8532944E+00, 3.0871200E-01, 0.0000000E+00 +* ELEM 428 circle + 2, 176, 175 +* + 7.4900000E+00, 1.8532944E+00, 1.1000000E-01, 0.0000000E+00 +* ELEM 429 circle + 2, 179, 178 +* + 6.9550000E+00, 2.7799416E+00, 4.3500000E-01, 0.0000000E+00 +* ELEM 430 circle + 2, 180, 179 +* + 6.9550000E+00, 2.7799416E+00, 4.2250000E-01, 0.0000000E+00 +* ELEM 431 circle + 2, 181, 180 +* + 6.9550000E+00, 2.7799416E+00, 4.1253600E-01, 0.0000000E+00 +* ELEM 432 circle + 2, 182, 181 +* + 6.9550000E+00, 2.7799416E+00, 3.8108400E-01, 0.0000000E+00 +* ELEM 433 circle + 2, 183, 182 +* + 6.9550000E+00, 2.7799416E+00, 3.0871200E-01, 0.0000000E+00 +* ELEM 434 circle + 2, 184, 183 +* + 6.9550000E+00, 2.7799416E+00, 1.1000000E-01, 0.0000000E+00 +* ELEM 435 line segment + 1, 0, 186 +* + 6.0432790E+00, 3.4890887E+00, -4.3301270E-02, -2.5000000E-02 +* ELEM 436 circular arc + 3, 187, 186 +* + 6.4200000E+00, 3.7065887E+00, 4.3500000E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 437 circular arc + 3, 187, 186 +* + 6.4200000E+00, 3.7065887E+00, 4.3500000E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 438 line segment + 1, 0, 186 +* + 6.8400223E+00, 3.9490887E+00, -4.3301270E-02, -2.5000000E-02 +* ELEM 439 line segment + 1, 0, 187 +* + 6.0541043E+00, 3.4953387E+00, -1.0825318E-02, -6.2500000E-03 +* ELEM 440 circular arc + 3, 188, 187 +* + 6.4200000E+00, 3.7065887E+00, 4.2250000E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 441 circular arc + 3, 188, 187 +* + 6.4200000E+00, 3.7065887E+00, 4.2250000E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 442 line segment + 1, 0, 187 +* + 6.7967211E+00, 3.9240887E+00, -1.0825318E-02, -6.2500000E-03 +* ELEM 443 line segment + 1, 0, 188 +* + 6.0627334E+00, 3.5003207E+00, -8.6290771E-03, -4.9820000E-03 +* ELEM 444 circular arc + 3, 189, 188 +* + 6.4200000E+00, 3.7065887E+00, 4.1253600E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 445 circular arc + 3, 189, 188 +* + 6.4200000E+00, 3.7065887E+00, 4.1253600E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 446 line segment + 1, 0, 188 +* + 6.7858958E+00, 3.9178387E+00, -8.6290771E-03, -4.9820000E-03 +* ELEM 447 line segment + 1, 0, 189 +* + 6.0899716E+00, 3.5160467E+00, -2.7238231E-02, -1.5726000E-02 +* ELEM 448 circular arc + 3, 190, 189 +* + 6.4200000E+00, 3.7065887E+00, 3.8108400E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 449 circular arc + 3, 190, 189 +* + 6.4200000E+00, 3.7065887E+00, 3.8108400E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 450 line segment + 1, 0, 189 +* + 6.7772667E+00, 3.9128567E+00, -2.7238231E-02, -1.5726000E-02 +* ELEM 451 line segment + 1, 0, 190 +* + 6.1526476E+00, 3.5522327E+00, -6.2675991E-02, -3.6186000E-02 +* ELEM 452 circular arc + 3, 191, 190 +* + 6.4200000E+00, 3.7065887E+00, 3.0871200E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 453 circular arc + 3, 191, 190 +* + 6.4200000E+00, 3.7065887E+00, 3.0871200E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 454 line segment + 1, 0, 190 +* + 6.7500284E+00, 3.8971307E+00, -6.2675991E-02, -3.6186000E-02 +* ELEM 455 line segment + 1, 0, 191 +* + 6.3247372E+00, 3.6515887E+00, -1.7208964E-01, -9.9356000E-02 +* ELEM 456 circular arc + 3, 192, 191 +* + 6.4200000E+00, 3.7065887E+00, 1.1000000E-01, 2.1000000E+02, 1.5000000E+02 +* ELEM 457 circular arc + 3, 192, 191 +* + 6.4200000E+00, 3.7065887E+00, 1.1000000E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 458 line segment + 1, 0, 191 +* + 6.6873525E+00, 3.8609447E+00, -1.7208964E-01, -9.9356000E-02 +* ELEM 459 line segment + 1, 0, 192 +* + 6.5152628E+00, 3.7615887E+00, -1.9052559E-01, -1.1000000E-01 +* ELEM 460 line segment + 1, 0, 194 +* + 4.3500000E-01, 0.0000000E+00, 5.0000000E-02, 0.0000000E+00 +* ELEM 461 circular arc + 3, 195, 194 +* + 0.0000000E+00, 0.0000000E+00, 4.3500000E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 462 line segment + 1, 0, 194 +* + 4.2002232E-01, 2.4250000E-01, -4.3301270E-02, -2.5000000E-02 +* ELEM 463 line segment + 1, 0, 195 +* + 4.2250000E-01, 0.0000000E+00, 1.2500000E-02, 0.0000000E+00 +* ELEM 464 circular arc + 3, 196, 195 +* + 0.0000000E+00, 0.0000000E+00, 4.2250000E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 465 line segment + 1, 0, 195 +* + 3.7672105E-01, 2.1750000E-01, -1.0825318E-02, -6.2500000E-03 +* ELEM 466 line segment + 1, 0, 196 +* + 4.1253600E-01, 0.0000000E+00, 9.9640000E-03, 0.0000000E+00 +* ELEM 467 circular arc + 3, 197, 196 +* + 0.0000000E+00, 0.0000000E+00, 4.1253600E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 468 line segment + 1, 0, 196 +* + 3.6589573E-01, 2.1125000E-01, -8.6290771E-03, -4.9820000E-03 +* ELEM 469 line segment + 1, 0, 197 +* + 3.8108400E-01, 0.0000000E+00, 3.1452000E-02, 0.0000000E+00 +* ELEM 470 circular arc + 3, 198, 197 +* + 0.0000000E+00, 0.0000000E+00, 3.8108400E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 471 line segment + 1, 0, 197 +* + 3.5726666E-01, 2.0626800E-01, -2.7238231E-02, -1.5726000E-02 +* ELEM 472 line segment + 1, 0, 198 +* + 3.0871200E-01, 0.0000000E+00, 7.2372000E-02, 0.0000000E+00 +* ELEM 473 circular arc + 3, 199, 198 +* + 0.0000000E+00, 0.0000000E+00, 3.0871200E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 474 line segment + 1, 0, 198 +* + 3.3002842E-01, 1.9054200E-01, -6.2675991E-02, -3.6186000E-02 +* ELEM 475 line segment + 1, 0, 199 +* + 1.1000000E-01, 0.0000000E+00, 1.9871200E-01, 0.0000000E+00 +* ELEM 476 circular arc + 3, 200, 199 +* + 0.0000000E+00, 0.0000000E+00, 1.1000000E-01, 0.0000000E+00, 3.0000000E+01 +* ELEM 477 line segment + 1, 0, 199 +* + 2.6735243E-01, 1.5435600E-01, -1.7208964E-01, -9.9356000E-02 +* ELEM 478 line segment + 1, 0, 200 +* + 0.0000000E+00, 0.0000000E+00, 1.1000000E-01, 0.0000000E+00 +* ELEM 479 line segment + 1, 0, 200 +* + 9.5262794E-02, 5.5000000E-02, -9.5262794E-02, -5.5000000E-02 +* ELEM 480 line segment + 1, 0, 201 +* + 7.3014103E+00, 4.2154712E+00, -3.8971143E-01, -2.2500000E-01 +* ELEM 481 line segment + 1, 0, 201 +* + 9.2155985E+00, 0.0000000E+00, 5.1961524E-01, 0.0000000E+00 +* boundaries conditions: defaul nbbcda allsur + 0, 3, 0 +* albedo + 0.0 +* type number of elements + 4, 74 +* elements +479 +477 +474 +471 +468 +465 +462 +59 +54 +86 +174 +178 +182 +186 +190 +194 +198 +197 +193 +189 +185 +181 +177 +89 +82 +109 +249 +253 +257 +261 +265 +269 +273 +272 +268 +264 +260 +256 +252 +112 +105 +124 +336 +340 +344 +348 +352 +356 +360 +359 +355 +351 +347 +343 +339 +127 +120 +129 +435 +439 +443 +447 +451 +455 +459 +458 +454 +450 +446 +442 +438 +132 +480 +4 +* tx, ty, angle + 0.0000000E+00 0.0000000E+00 3.0000000E+01 +* type number of elements + 4, 131 +* elements +478 +475 +472 +469 +466 +463 +460 +57 +52 +138 +141 +144 +147 +150 +153 +154 +151 +148 +145 +142 +139 +136 +50 +46 +157 +160 +163 +166 +169 +172 +173 +170 +167 +164 +161 +158 +155 +44 +40 +201 +204 +207 +210 +213 +216 +217 +214 +211 +208 +205 +202 +199 +38 +34 +226 +229 +232 +235 +238 +241 +242 +239 +236 +233 +230 +227 +224 +32 +28 +276 +279 +282 +285 +288 +291 +292 +289 +286 +283 +280 +277 +274 +26 +22 +307 +310 +313 +316 +319 +322 +323 +320 +317 +314 +311 +308 +305 +20 +16 +363 +366 +369 +372 +375 +378 +379 +376 +373 +370 +367 +364 +361 +14 +12 +400 +403 +406 +409 +412 +415 +416 +413 +410 +407 +404 +401 +398 +10 +56 +481 +2 +* tx, ty, angle + 0.0000000E+00 0.0000000E+00 0.0000000E+00 +* type number of elements + 4, 1 +* elements +1 +* tx, ty, angle + 9.9950213E+00 0.0000000E+00 1.2000000E+02 +# 1 - MIX1_Na +# 2 - MIX6_CLAD +# 3 - MIX9_He +# 4 - MIX5_MOX +# 5 - MIX4_MOX +# 6 - MIX3_MOX +# 7 - MIX2_MOX +# 8 - MIX8_He +# 9 - MIX7_EM10 +* medium number per region + 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 1 + 2 + 3 + 4 + 5 + 6 + 7 + 8 + 9 + 1 + 1 + 1 + 1 + 1 + 1 +------------------------------------------------------------ + this is the last line generated by conversion diff --git a/Dragon/src/EXCELP.f b/Dragon/src/EXCELP.f index 77a9fab..9781c59 100644 --- a/Dragon/src/EXCELP.f +++ b/Dragon/src/EXCELP.f @@ -85,18 +85,18 @@ *---- * LOCAL VARIABLES *---- - INTEGER IOUT, ICPALL, ICPEND, MXGAUS, NSTATE + INTEGER IOUT, ICPALL, ICPEND, MXGAUS, NSTATE, MAXCDA PARAMETER ( IOUT=6, ICPALL=4, ICPEND=3, MXGAUS=64, - > NSTATE=40 ) + > NSTATE=40, MAXCDA=30 ) CHARACTER NAMSBR*6 PARAMETER ( NAMSBR='EXCELP') INTEGER MKI1, MKI2, MKI3, MKI4, MKI5 PARAMETER (MKI1=600,MKI2=600,MKI3=600,MKI4=600,MKI5=600) - INTEGER ISTATE(NSTATE),ICODE(6) + INTEGER ISTATE(NSTATE),ICODE(MAXCDA) INTEGER NPROB,ISBG,KSBG,ITYPBC - REAL ALBEDO(6),EXTKOP(NSTATE),CUTOF,RCUTOF,ASCRP, - > YGSS,XGSS(MXGAUS),WGSS(MXGAUS),WGSSX(MXGAUS), - > ALBG(6) + REAL ALBEDO(MAXCDA),EXTKOP(NSTATE),CUTOF,RCUTOF, + > ASCRP,YGSS,XGSS(MXGAUS),WGSS(MXGAUS), + > WGSSX(MXGAUS),ALBG(MAXCDA) LOGICAL SWVOID, LPIJK CHARACTER CTRKT*4, COMENT*80 DOUBLE PRECISION DANG0,DASCRP @@ -104,7 +104,8 @@ INTEGER JJ,MSYM,IL,NALLOC,ITRAK,IANG,IC,IPRT,ISPEC, > IUN,KSPEC,LOPT,MXSEG,NALBG,NANGL,NCOMNT,NCOR, > NCORT,NDIM,NGSS,NREG2,NSCRP,NTRK,NUNKNO,JGSS, - > JUN,IFMT,MXSUB,ISA,IBATCH,IL1 ,III,IND,I,J + > JUN,IFMT,MXSUB,ISA,IBATCH,IL1,III,IND,I,J, + > ILONG,ITYLCM *---- * Variables for NXT: inline tracking *---- @@ -153,7 +154,7 @@ 1 + MIN(I+NSOUT+1,J+NSOUT+1) *---- * RECOVER EXCELL SPECIFIC TRACKING INFORMATION. -* ALBEDO: SURFACE ALBEDOS (REAL(6)) +* ALBEDO: SURFACE ALBEDOS (REAL(MAXCDA)) * KSPEC : KIND OF PIJ INTEGRATION (0:ISOTROPE,1:SPECULAR) * CUTOF : MFP CUTOFF FOR SPECULAR INTEGRATION *---- @@ -162,6 +163,9 @@ KSPEC=ISTATE(10) CALL LCMGET(IPTRK,'EXCELTRACKOP',EXTKOP) CUTOF=EXTKOP(1) + CALL LCMLEN(IPTRK,'ICODE',ILONG,ITYLCM) + IF(ILONG.GT.MAXCDA) CALL XABORT('EXCELP: MAXCDA OVERFLOW(1).') + ICODE(:MAXCDA)=0 CALL LCMGET(IPTRK,'ICODE',ICODE) CALL LCMGET(IPTRK,'ALBEDO',ALBG) * @@ -269,16 +273,16 @@ SWVOID= .FALSE. DO ISBG=1,NSBG IF(NPSYS(ISBG).NE.0) THEN - DO ISA=1,6 - ALBEDO(ISA)=ALBG(ISA) - ENDDO + ALBEDO(:MAXCDA)=ALBG(:MAXCDA) IF(NALBP .GT. 0) THEN - DO ISA=1,6 + DO ISA=1,MAXCDA IF(ICODE(ISA).GT.0) ALBEDO(ISA)=ALBP(ICODE(ISA),ISBG) ENDDO ENDIF DO IUN= -NSOUT, -1 SIGT00(IUN,ISBG)= 0.0 + IF(-MATALB(IUN).GT.MAXCDA) CALL XABORT('EXCELP: MAXCDA OV' + > //'ERFLOW(2).') SIGTAL(IUN,ISBG)= ALBEDO(-MATALB(IUN)) SWNZBC= SWNZBC.OR.(SIGTAL(IUN,ISBG).NE.0.0) ENDDO diff --git a/Dragon/src/MCGSIG.f b/Dragon/src/MCGSIG.f index 9509841..3589cd6 100644 --- a/Dragon/src/MCGSIG.f +++ b/Dragon/src/MCGSIG.f @@ -41,12 +41,14 @@ * LOCAL VARIABLES *--- TYPE(C_PTR) JPSYS - INTEGER I,II,ISA,ICODE(6) + INTEGER I,II,ISA,ICODE(6),ILONG,ITYLCM REAL ALBG(6),ALBEDO(6) REAL, ALLOCATABLE, DIMENSION(:) :: ALBP *--- * RECOVER ALBEDO INFORMATION FROM TRACKING *--- + CALL LCMLEN(IPTRK,'ICODE',ILONG,ITYLCM) + IF(ILONG.GT.6) CALL XABORT('MCGSIG: ALBEDO OVERFLOW.') CALL LCMGET(IPTRK,'ICODE',ICODE) CALL LCMGET(IPTRK,'ALBEDO',ALBG) * diff --git a/Dragon/src/MUSACG.f90 b/Dragon/src/MUSACG.f90 index 9ceaeb8..3d560e4 100644 --- a/Dragon/src/MUSACG.f90 +++ b/Dragon/src/MUSACG.f90 @@ -54,12 +54,14 @@ SUBROUTINE MUSACG(ITRACK,IFTRK,IPRINT,IMACRO,NBSLIN,RCUTOF,GG,LGINF,NBNODE_MACRO INTEGER, PARAMETER :: NSTATE=40 INTEGER, PARAMETER :: FOUT=6 INTEGER, PARAMETER :: NDIM=2 ! NUMBER OF DIMENSIONS + INTEGER, PARAMETER :: MAXCDA=30 ! MAXIMUM NUMBER OF PERIMETERS INTEGER ELEM, OK, TYPE REAL(PDB) :: X1,X2,Y1,Y2,DET1,DET2 REAL(PDB) :: DGMESHX(2),DGMESHY(2) LOGICAL :: LTEST INTEGER, DIMENSION(NSTATE) :: I_STATE,IEDIMG CHARACTER(LEN=72) :: TEXT72 + CHARACTER(LEN=131) :: HSMG !---- ! Allocatable arrays !---- @@ -384,7 +386,11 @@ SUBROUTINE MUSACG(ITRACK,IFTRK,IPRINT,IMACRO,NBSLIN,RCUTOF,GG,LGINF,NBNODE_MACRO WRITE(FOUT,'(3HEND)') WRITE(FOUT,'(5H--cut,70(1H-),I5)') IMACRO ENDIF - IF(GG_MAC%NBBCDA.GT.6) CALL XABORT('MUSACG: The unfolded geometry has more than 6 perimeters') + IF(GG_MAC%NBBCDA.GT.MAXCDA) THEN + WRITE(HSMG,'(33HMUSACG: The unfolded geometry has,I3,14H perimeters (>,I3,2H).)') & + & GG_MAC%NBBCDA,MAXCDA + CALL XABORT(HSMG) + ENDIF !**** !* compute node perimeters for the macro ALLOCATE (GG_MAC%PPERIM_NODE(GG_MAC%NB_NODE+1),STAT=OK) @@ -649,7 +655,7 @@ SUBROUTINE MUSACG(ITRACK,IFTRK,IPRINT,IMACRO,NBSLIN,RCUTOF,GG,LGINF,NBNODE_MACRO MATALB(-I)=-1 GALBED(1)=REAL(GG_MAC%ALBEDO) ELSE - IF(INDEX.GT.6) CALL XABORT('MUSACG: SDIRE overflow.') + IF(INDEX.GT.MAXCDA) CALL XABORT('MUSACG: INDEX overflow.') IF(INDEX > GG_MAC%NALBG) THEN CALL XABORT('MUSACG: Albedo array overflow(2).') ENDIF diff --git a/Dragon/src/PIJWIJ.f b/Dragon/src/PIJWIJ.f index c4d9e11..bea36a2 100644 --- a/Dragon/src/PIJWIJ.f +++ b/Dragon/src/PIJWIJ.f @@ -88,13 +88,14 @@ *---- * LOCAL VARIABLES *---- - INTEGER IOUT, ICPALL, ICPEND, MXGAUS, NSTATE + INTEGER IOUT, ICPALL, ICPEND, MXGAUS, NSTATE, MAXCDA PARAMETER ( IOUT=6, ICPALL=4, ICPEND=3, MXGAUS=64, - > NSTATE=40 ) + > NSTATE=40, MAXCDA=30 ) CHARACTER NAMSBR*6 PARAMETER ( NAMSBR='PIJWIJ') - INTEGER ILONG,ITYPE,NPROB,ISBG,ISTATE(NSTATE),ICODE(6) - REAL FACT,ALBEDO(6),ALBG(6) + INTEGER ILONG,ITYPE,NPROB,ISBG,ISTATE(NSTATE), + > ICODE(MAXCDA) + REAL FACT,ALBEDO(MAXCDA),ALBG(MAXCDA) LOGICAL LSKIP,SWNZBC,SWVOID * INTEGER MSYM,IU,IL,ISOUT,IIN,I,J,IBM,IOP,INDPIJ,IJKS, @@ -125,6 +126,9 @@ ISTATE(:NSTATE)=0 CALL LCMGET(IPTRK,'STATE-VECTOR',ISTATE) KSPEC=ISTATE(10) + CALL LCMLEN(IPTRK,'ICODE',ILONG,ITYPE) + IF(ILONG.GT.MAXCDA) CALL XABORT('PIJWIJ: MAXCDA OVERFLOW(1).') + ICODE(:MAXCDA)=0 CALL LCMGET(IPTRK,'ICODE',ICODE) CALL LCMGET(IPTRK,'ALBEDO',ALBG) *---- @@ -135,11 +139,9 @@ SWVOID= .FALSE. DO ISBG=1,NSBG IF(NPSYS(ISBG).NE.0) THEN - DO ISA=1,6 - ALBEDO(ISA)=ALBG(ISA) - ENDDO + ALBEDO(:MAXCDA)=ALBG(:MAXCDA) IF(NALBP .GT. 0) THEN - DO ISA=1,6 + DO ISA=1,MAXCDA IF(ICODE(ISA).GT.0) ALBEDO(ISA)=ALBP(ICODE(ISA),ISBG) ENDDO ENDIF diff --git a/Dragon/src/SALACG.f90 b/Dragon/src/SALACG.f90 index 0dbc708..679446c 100644 --- a/Dragon/src/SALACG.f90 +++ b/Dragon/src/SALACG.f90 @@ -45,13 +45,14 @@ SUBROUTINE SALACG(FGEO ,ITRACK, RCUTOF, IPRINT, GG) ! Local variables !---- INTEGER, PARAMETER :: NSTATE=40 - INTEGER, PARAMETER :: NDIM=2 ! NUMBER OF DIMENSIONS - INTEGER, PARAMETER :: NALBG=6 ! NUMBER OF ALBEDOS + INTEGER, PARAMETER :: NDIM=2 ! NUMBER OF DIMENSIONS + INTEGER, PARAMETER :: MAXCDA=30 ! MAXIMUM NUMBER OF ALBEDOS LOGICAL LGINF INTEGER, DIMENSION(NSTATE) :: I_STATE,IEDIMG - INTEGER OK,I,J,NREG,ELEM,NFREG,LEAK,NSOUT,ICODE(NALBG),INDEX,MMAX - REAL GALBED(NALBG) + INTEGER NALBG,OK,I,J,NREG,ELEM,NFREG,LEAK,NSOUT,ICODE(MAXCDA),INDEX,MMAX + REAL GALBED(MAXCDA) CHARACTER(LEN=72) TEXT72 + CHARACTER(LEN=131) HSMG REAL(PDB) :: DGMESHX(2),DGMESHY(2) !---- ! Allocatable arrays @@ -234,10 +235,16 @@ SUBROUTINE SALACG(FGEO ,ITRACK, RCUTOF, IPRINT, GG) VOLSUR(1:NREG)=VOLUME(:NREG) DEALLOCATE(VOLUME) ! boundary conditions structures - ICODE(:NALBG)=(/ (-I,I=1,NALBG) /) - GALBED(:NALBG)=REAL(GG%ALBEDO) + ICODE(:MAXCDA)=(/ (-I,I=1,MAXCDA) /) + GALBED(:MAXCDA)=REAL(GG%ALBEDO) + IF(GG%NALBG.GT.MAXCDA) THEN + WRITE(HSMG,'(24HSALACG: The geometry has,I3,17H albedo values (>,I3,2H).)') & + & GG%NALBG,MAXCDA + CALL XABORT(HSMG) + ENDIF IF(ISPEC == 0) THEN - IF(GG%NALBG > 6) CALL XABORT('SALACG: Albedo array overflow(1).') + NALBG=GG%NALBG + IF(TYPGEO.EQ.0) NALBG=6 DO I=1,NSOUT KEYMRG(-I)=-I VOLSUR(-I)=GG%SURF2(I) @@ -247,10 +254,8 @@ SUBROUTINE SALACG(FGEO ,ITRACK, RCUTOF, IPRINT, GG) MATALB(-I)=-1 GALBED(1)=REAL(GG%ALBEDO) ELSE - IF(INDEX.GT.6) CALL XABORT('SALACG: SDIRE overflow.') - IF(INDEX > GG%NALBG) THEN - CALL XABORT('SALACG: Albedo array overflow(2).') - ENDIF + IF(INDEX.GT.MAXCDA) CALL XABORT('SALACG: INDEX overflow(1).') + IF(INDEX.GT.GG%NALBG) CALL XABORT('SALACG: INDEX overflow(2).') MATALB(-I)=-INDEX IF(SIZE(GG%BCDATA) > 0) THEN GALBED(INDEX)=REAL(GG%BCDATA(6,INDEX)) @@ -260,13 +265,14 @@ SUBROUTINE SALACG(FGEO ,ITRACK, RCUTOF, IPRINT, GG) ENDIF ENDDO ELSE + NALBG=6 DO I=1,NSOUT VOLSUR(-I)=0.0 KEYMRG(-I)=-I MATALB(-I)=-1 ENDDO GALBED(:NALBG)=1.0 - ENDIF + ENDIF MATALB(0)=0 KEYMRG(0)=0 VOLSUR(0)=0._PDB diff --git a/Dragon/src/SALTCG.f b/Dragon/src/SALTCG.f index 123a288..dcf496f 100644 --- a/Dragon/src/SALTCG.f +++ b/Dragon/src/SALTCG.f @@ -57,10 +57,10 @@ PARAMETER (IOUT=6,NAMSBR='SALTCG') INTEGER NSTATE PARAMETER (NSTATE=40) - INTEGER NMAX0 + INTEGER NMAX0,MAXCDA DOUBLE PRECISION PI,DZERO,DONE,DTWO,DSUM PARAMETER (PI=3.14159265358979, DZERO=0.0D0,DONE=1.0D0, - > DTWO=2.0D0,NMAX0=100000) + > DTWO=2.0D0,NMAX0=100000,MAXCDA=30) *---- * Functions *---- @@ -68,8 +68,8 @@ *---- * Local variables *---- - INTEGER ISTATE(NSTATE),IEDIMG(NSTATE),ICODE(6) - REAL RSTATT(NSTATE),ALBEDO(6) + INTEGER ISTATE(NSTATE),IEDIMG(NSTATE),ICODE(MAXCDA) + REAL RSTATT(NSTATE),ALBEDO(MAXCDA) INTEGER RENO,LTRK,AZMOAQ,ISYMM,POLQUA,POLOAQ,AZMQUA, > AZMNBA,OK DOUBLE PRECISION DENUSR,RCUTOF,DENLIN,SPACLN,WEIGHT @@ -81,7 +81,7 @@ > MAXSUB,MAXSGL,NBDR,ILONG,ITYLCM,IPER(3) INTEGER JJ,KK,NCOR,NQUAD,NANGL,NBANGL,LINMAX DOUBLE PRECISION DQUAD(4),ABSC(3,2),RCIRC,SIDEH,ANGLE - CHARACTER CTRK*4,COMENT*80 + CHARACTER CTRK*4,COMENT*80,HSMG*131 INTEGER IFMT,NEREG,NESUR *---- * Allocatable arrays @@ -132,7 +132,11 @@ * Get main tracking records *---- CALL LCMLEN(IPTRK,'ICODE ',ILONG,ITYLCM) - IF(ILONG.GT.6) CALL XABORT('SALTCG: ALBEDO OVERFLOW.') + IF(ILONG.GT.MAXCDA) THEN + WRITE(HSMG,'(24HSALTCG: The geometry has,I3,15H albedo values , + 1 2H(>,I3,2H).)') ILONG,MAXCDA + CALL XABORT(HSMG) + ENDIF CALL LCMGET(IPTRK,'ICODE ',ICODE ) CALL LCMGET(IPTRK,'ALBEDO ',ALBEDO) CALL LCMSIX(IPTRK,'NXTRecords ',1) |