Télécharger la plaquette Cast3M


fr en

Télécharger test_iwprd3D_sol.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : test_iwpr3d_sol.dgibi
  2. *
  3. **********************************************************
  4. *  Test pour la loi de comportement IWPR3D_SOL
  5. *  On considere un essai numerique cyclique en imposant 
  6. *  un tensuer de deformation donne:
  7. *
  8. *  eps = gam1[(ex^ez)+(ez^ex)]
  9. *  avec gam1=1e-4
  10. *
  11. *  L'etat de contraite initiale est:
  12. *  sig0 = sig0[(ex^ex)+(ey^ey)+(ez^ez)]
  13. *  avec sig0 = -3.0e4 Pa
  14. *  les parametres de la loi IWPR3D_SOL sont:
  15. *     n1 = 0.1 k1 = 0.75 
  16. *  et les seuil de chaque criteres sont:
  17. *     tau1 = 200.      m1 = (sqrt(3)*tau1)/sig0 
  18. *     tau2 = 400.      m2 = (sqrt(3)*tau2)/sig0 
  19. *     tau3 = 600.      m3 = (sqrt(3)*tau3)/sig0 
  20. *     tau4 = 800.      m4 = (sqrt(3)*tau4)/sig0 
  21. *     tau5 = 900.      m5 = (sqrt(3)*tau5)/sig0 
  22. *     tau6 = 1000.     m6 = (sqrt(3)*tau6)/sig0 
  23. *     tau7 = 1200.     m7 = (sqrt(3)*tau7)/sig0 
  24. *     tau8 = 1400.     m8 = (sqrt(3)*tau8)/sig0 
  25. *     tau9 = 1600.     m9 = (sqrt(3)*tau9)/sig0 
  26. *
  27. *  et les caracteristiques associées à chaque critere
  28. *     H1 = 4000. D1=0.1 eta1=0.142 eta1=0.142
  29. *     H2 = 3500. D2=0.1 eta2=0.142 eta2=0.142
  30. *     H3 = 3000. D3=0.1 eta3=0.142 eta3=0.142
  31. *     H4 = 2800. D4=0.1 eta4=0.142 eta4=0.142
  32. *     H5 = 2600. D5=0.1 eta5=0.142 eta5=0.142
  33. *     H6 = 2200. D6=0.1 eta6=0.142 eta6=0.142
  34. *     H7 = 1800. D7=0.1 eta7=0.142 eta7=0.142
  35. *     H8 = 1000. D8=0.1 eta8=0.142 eta8=0.142
  36. *     H9 = 800.  D9=0.1 eta9=0.142 eta9=0.142
  37. *
  38. *  Develloppé par Alberto FRAU /DEN/DANS/DM2S/SEMT/EMSI
  39. *
  40. **********************************************************
  41. *
  42. OPTI DIME 3 MODE TRID ELEM CUB8;
  43. OPTI TRAC PSC;
  44. OPTI EPTR 12;
  45. *
  46. VSSOL1 = 200.;
  47. NU1 = 0.3;
  48. RHO1 = 2000.;
  49. G1 = (VSSOL1*VSSOL1)*(RHO1);
  50. E1 = ((G1)*((2.0)*(1.0 + NU1)));
  51. *
  52. DS1 = 1.0;
  53. *
  54. DT1 = 0.01;
  55. *
  56. D_IMP = 1.E-4;
  57. *
  58. S_XX_INI = -3.E4;
  59. S_YY_INI = -3.E4;
  60. S_ZZ_INI = -3.E4;
  61. *
  62. VAR_GRAP = 0;
  63. *
  64. *********************************************************************
  65. *   Maillage
  66. *********************************************************************
  67. *
  68. P0 = 0. 0. 0.;
  69. P1 = P0 PLUS (DS1 0. 0.);
  70. P2 = P0 PLUS (DS1 DS1 0.);
  71. P3 = P0 PLUS (0. DS1 0.);
  72. *
  73. L1 = D 1 P0 P1;
  74. L2 = D 1 P1 P2;
  75. L3 = D 1 P2 P3;
  76. L4 = D 1 P3 P0;
  77. *
  78. SUR1 = DALL L1 L2 L3 L4;
  79. *
  80. VOL1 = SUR1 VOLU 1 TRANS (0. 0. DS1);
  81. *
  82. FLX_1 = L1 TRANS 1 (0. 0. DS1);
  83. FLX_2 = L3 TRANS 1 (0. 0. DS1);
  84. FLY_1 = L4 TRANS 1 (0. 0. DS1);
  85. FLY_2 = L2 TRANS 1 (0. 0. DS1);
  86. FB_IN = SUR1;
  87. FB_SU = SUR1 PLUS (0. 0. DS1);
  88. *
  89. VOL1 = VOL1 COUL VERT;
  90. FLX_1 = FLX_1 COUL ROSE;
  91. FLX_2 = FLX_2 COUL ROSE;
  92. FLY_1 = FLY_1 COUL BLEU;
  93. FLY_2 = FLY_2 COUL BLEU;
  94. FB_IN = FB_IN COUL ROUG;
  95. FB_SU = FB_SU COUL ROUG;
  96. *
  97. ELIM 0.001 (VOL1 ET FLX_1 ET FLX_2 ET FLY_1 ET FLY_2 
  98.                                    ET FB_IN ET FB_SU); 
  99. *
  100. SI (VAR_GRAP EGA 1);
  101.   TRAC (VOL1);
  102.   TRAC (VOL1 ET FLX_1 ET FLX_2);
  103.   TRAC (VOL1 ET FLX_1 ET FLX_2 ET FLY_1 ET FLY_2);
  104.   TRAC (VOL1 ET FLX_1 ET FLX_2 ET FLY_1 ET FLY_2
  105.              ET FB_IN ET FB_SU);
  106. FINSI;
  107. *
  108. *********************************************************************
  109. *   Modele et Materiau
  110. *********************************************************************
  111. *
  112. *
  113. MOD1 = MODE VOL1 MECANIQUE ELASTIQUE PLASTIQUE IWPR3D_SOL CUB8;
  114. *
  115. XM1 = (((3.0)**(0.5))*(200.0))/(30000.0);
  116. XM2 = (((3.0)**(0.5))*(400.0))/(30000.0);
  117. XM3 = (((3.0)**(0.5))*(600.0))/(30000.0);
  118. XM4 = (((3.0)**(0.5))*(800.0))/(30000.0);
  119. XM5 = (((3.0)**(0.5))*(900.0))/(30000.0);
  120. XM6 = (((3.0)**(0.5))*(1000.0))/(30000.0);
  121. XM7 = (((3.0)**(0.5))*(1200.0))/(30000.0);
  122. XM8 = (((3.0)**(0.5))*(1400.0))/(30000.0);
  123. XM9 = (((3.0)**(0.5))*(1600.0))/(30000.0);
  124. *
  125. MAT1 = MATE MOD1 'YOUNG' E1  'NU' NU1 'RHO' RHO1;
  126. MAT1 = MAT1 ET 
  127. (MATE MOD1 'E0' ((E1)/((30000)**(0.1))) 'ALP0' 1.0 'C' 0. 'N1' 0.1 
  128.     'PREF' -1.0 'K1' 0.75);
  129. MAT1 = MAT1 ET 
  130. (MATE MOD1 'M1' XM1 'M2' XM2 'M3' XM3 'M4' XM4 'M5' XM5 
  131.            'M6' XM6 'M7' XM7 'M8' XM8 'M9' XM9 'M10' 19.0);
  132. MAT1 = MAT1 ET 
  133. (MATE MOD1 'HC1' 4000.  'HC2' 3500.   'HC3' 3000.   'HC4' 2800.
  134.            'HC5' 2600.  'HC6' 2200.   'HC7' 1800.   'HC8' 1000.
  135.            'HC9' 800.);
  136. MAT1 = MAT1 ET 
  137. (MATE MOD1 'HE1' 4000.  'HE2' 3500. 'HE3' 3000. 'HE4' 2800.
  138.            'HE5' 2600.  'HE6' 2200. 'HE7' 1800. 'HE8' 1000.
  139.            'HE9' 800.);
  140. MAT1 = MAT1 ET 
  141. (MATE MOD1 'DA1' 0.1 'DA2' 0.1 'DA3' 0.1 'DA4' 0.1 'DA5' 0.1
  142.            'DA6' 0.1 'DA7' 0.1 'DA8' 0.1 'DA9' 0.1);
  143. MAT1 = MAT1 ET 
  144. (MATE MOD1   'E_C1' 0.142 'E_C2' 0.142 'E_C3' 0.142 'E_C4' 0.142  
  145. 'E_C5' 0.142 'E_C6' 0.142 'E_C7' 0.142 'E_C8' 0.142 'E_C9' 0.142);
  146. MAT1 = MAT1 ET 
  147. (MATE MOD1   'E_E1' 0.142 'E_E2' 0.142 'E_E3' 0.142 'E_E4' 0.142 
  148. 'E_E5' 0.142 'E_E6' 0.142 'E_E7' 0.142 'E_E8' 0.142 'E_E9' 0.142);
  149. *
  150. *********************************************************************
  151. *   Condictions initiales
  152. *********************************************************************
  153. *
  154. SIG_INI = MANU 'CHML' MOD1 'SMXX' S_XX_INI 
  155.                              'SMYY' S_YY_INI 
  156.                              'SMZZ' S_ZZ_INI
  157.                              'SMXY' 0.0 'SMXZ' 0.0 'SMYZ' 0.0;
  158. F_EX_INI = BSIGM MOD1 SIG_INI;
  159. *
  160. *********************************************************************
  161. *   Chargement
  162. *********************************************************************
  163. *
  164. BL1 = BLOQUER 'DEPLA' FB_IN;
  165. BL2 = BLOQUER 'UY' FB_SU;
  166. BL3 = BLOQUER 'UZ' FB_SU;
  167. BL4 = BLOQUER 'UX' FB_SU;
  168. *
  169. DP1 = DEPI BL4 1.0;
  170. *
  171. L_T1 = PROG 0. PAS DT1 13.;
  172. L_E1 = IPOL L_T1 (PROG 0. 1. 2.  3. 4. 5. 6.  7. 8. 9. 10. 11. 12. 13.)
  173.                  (PROG 0. 1. 0. -1. 0. 1. 0. -1. 0. 1.  0. -1.  0.  1.);
  174. L_E2 = PROG (DIME L_T1) * 1.0;
  175. *
  176. EV1 = EVOL MANU L_T1 ((L_E1)*(D_IMP));
  177. EV2 = EVOL MANU L_T1 ((L_E2));
  178. *
  179. CHA1 = CHAR 'DIMP' DP1 EV1;
  180. CHA2 = CHAR 'FORCE' F_EX_INI EV2;
  181. *
  182. *********************************************************************
  183. *   Calcul PASAPAS
  184. *********************************************************************
  185. *
  186. CH_INI = MANU CHPO VOL1 3 'UX' 0. 'UY' 0. 'UZ' 0.;
  187. *
  188. LLT3 = PROG 0. PAS      DT1      13.;
  189. LLT4 = PROG 0. PAS ((DT1)/(1.0)) 13.;
  190. *
  191. TAB1                           = TABLE;
  192. TAB1.'MODELE'                  = MOD1;
  193. TAB1.'CARACTERISTIQUES'        = MAT1;
  194. TAB1.'CHARGEMENT'              = CHA1 ET CHA2;
  195. TAB1.'TEMPS_CALCULES'          = LLT3;
  196. TAB1.'TEMPS_SAUVES'            = LLT4;
  197. TAB1.'CONVERGENCE_FORCEE'      = FAUX;
  198. TAB1.'BLOCAGES_MECANIQUES'     = BL1 ET BL2 ET BL3 ET BL4;
  199. TAB1.'MES_SAUVEGARDES'         = TABLE;
  200. TAB1.'MES_SAUVEGARDES'.'DEFIN' = VRAI;
  201. TAB1.'MOVA'                    = 'GAM1';
  202. TAB1.'DEPLACEMENTS'            = TABLE;
  203. TAB1.'DEPLACEMENTS'. 0         = CH_INI;
  204. TAB1.'CONTRAINTES'             = TABLE;
  205. TAB1.'CONTRAINTES'. 0          = SIG_INI;
  206. TAB1.'PRECISION'               = 1E-4;
  207. TAB1.'HYPOTHESE_DEFORMATIONS'  = 'LINEAIRE';
  208. *
  209. PASAPAS TAB1;
  210. *
  211. *********************************************************************
  212. *   Posttraitement
  213. *********************************************************************
  214. *
  215. L_T1 = PROG;
  216. L_EXX = PROG;
  217. L_EYY = PROG;
  218. L_EZZ = PROG;
  219. L_EXY = PROG;
  220. L_EXZ = PROG;
  221. L_EYZ = PROG;
  222. L_SXX = PROG;
  223. L_SYY = PROG;
  224. L_SZZ = PROG;
  225. L_SXY = PROG;
  226. L_SXZ = PROG;
  227. L_SYZ = PROG;
  228. *
  229. I = 0;
  230. NN1 = DIME (TAB1.DEPLACEMENTS);
  231. *
  232. REPE IND1 (NN1);
  233. *
  234.   L_T1 = L_T1 ET (PROG (TAB1.TEMPS. I));
  235. *
  236.   DEP1 = TAB1.DEPLACEMENTS. I;
  237.   EPS1 = EPSI DEP1 MOD1 'LINE';
  238.   SIG1 = TAB1.CONTRAINTES. I;
  239. *
  240.   L_EXX = L_EXX ET (PROG (EXTR EPS1 'EPXX' 1 1 1));
  241.   L_EYY = L_EYY ET (PROG (EXTR EPS1 'EPYY' 1 1 1));
  242.   L_EZZ = L_EZZ ET (PROG (EXTR EPS1 'EPZZ' 1 1 1));
  243.   L_EXY = L_EXY ET (PROG (EXTR EPS1 'GAXY' 1 1 1));
  244.   L_EXZ = L_EXZ ET (PROG (EXTR EPS1 'GAXZ' 1 1 1));
  245.   L_EYZ = L_EYZ ET (PROG (EXTR EPS1 'GAYZ' 1 1 1));
  246. *
  247.   L_SXX = L_SXX ET (PROG (EXTR SIG1 'SMXX' 1 1 1));
  248.   L_SYY = L_SYY ET (PROG (EXTR SIG1 'SMYY' 1 1 1));
  249.   L_SZZ = L_SZZ ET (PROG (EXTR SIG1 'SMZZ' 1 1 1));
  250.   L_SXY = L_SXY ET (PROG (EXTR SIG1 'SMXY' 1 1 1));
  251.   L_SXZ = L_SXZ ET (PROG (EXTR SIG1 'SMXZ' 1 1 1));
  252.   L_SYZ = L_SYZ ET (PROG (EXTR SIG1 'SMYZ' 1 1 1));
  253. *
  254.   I = I + 1
  255. FIN IND1;
  256. *
  257. EV_EXX = EVOL (BLEU) MANU 'Temps [s]' L_T1 'Exx []' L_EXX;
  258. EV_EYY = EVOL (BLEU) MANU 'Temps [s]' L_T1 'Eyy []' L_EYY;
  259. EV_EZZ = EVOL (BLEU) MANU 'Temps [s]' L_T1 'Ezz []' L_EZZ;
  260. EV_EXY = EVOL (BLEU) MANU 'Temps [s]' L_T1 'Exy []' L_EXY;
  261. EV_EXZ = EVOL (BLEU) MANU 'Temps [s]' L_T1 'Exz []' L_EXZ;
  262. EV_EYZ = EVOL (BLEU) MANU 'Temps [s]' L_T1 'Eyz []' L_EYZ;
  263. *
  264. EV_SXX = EVOL (ROUG) MANU 'Temps [s]' L_T1 'Sxx [Pa]' L_SXX;
  265. EV_SYY = EVOL (ROUG) MANU 'Temps [s]' L_T1 'Syy [Pa]' L_SYY;
  266. EV_SZZ = EVOL (ROUG) MANU 'Temps [s]' L_T1 'Szz [Pa]' L_SZZ;
  267. EV_SXY = EVOL (ROUG) MANU 'Temps [s]' L_T1 'Sxy [Pa]' L_SXY;
  268. EV_SXZ = EVOL (ROUG) MANU 'Temps [s]' L_T1 'Sxz [Pa]' L_SXZ;
  269. EV_SYZ = EVOL (ROUG) MANU 'Temps [s]' L_T1 'Syz [Pa]' L_SYZ;
  270. *
  271. EV_1 = EVOL (ROSE) MANU 'Exz []' L_EXZ 'Sxz [Pa]' L_SXZ;
  272. EV_2 = EVOL (ROSE) MANU 'Exz []' L_EXZ 'p [Pa]' 
  273.                                    ((L_SXX + L_SYY + L_SZZ)/(3.0));
  274. *
  275. SI (VAR_GRAP EGA 1);
  276.   DESS EV_EXX 'TITR' 'Exx';
  277.   DESS EV_EYY 'TITR' 'Eyy';
  278.   DESS EV_EZZ 'TITR' 'Ezz';
  279.   DESS EV_EXY 'TITR' 'Exy';
  280.   DESS EV_EXZ 'TITR' 'Exz';
  281.   DESS EV_EYZ 'TITR' 'Eyz';
  282. *
  283.   DESS EV_SXX 'TITR' 'Sxx';
  284.   DESS EV_SYY 'TITR' 'Syy';
  285.   DESS EV_SZZ 'TITR' 'Szz';
  286.   DESS EV_SXY 'TITR' 'Sxy';
  287.   DESS EV_SXZ 'TITR' 'Sxz';
  288.   DESS EV_SYZ 'TITR' 'Syz';
  289. *
  290.   DESS EV_1 'TITR' 'Exz - Sxz';
  291.   DESS EV_2 'TITR' 'Exz - p';
  292. *
  293. FINSI;
  294. *
  295. *********************************************************************
  296. *   Contruction de la solution de reference
  297. *********************************************************************
  298. *
  299. L_SXZT = (PROG
  300.       0.00000000e+00      8.00000000e+01      1.60000000e+02
  301.       2.12839070e+02      2.60203100e+02      3.05291070e+02
  302.       3.52371660e+02      3.96740470e+02      4.40781330e+02
  303.       4.84086590e+02      5.27096150e+02      5.68406540e+02
  304.       6.08937740e+02      6.49230570e+02      6.88908370e+02);
  305. L_SXZT = L_SXZT ET (PROG
  306.       7.28652200e+02      7.67707310e+02      8.05779040e+02
  307.       8.43343830e+02      8.80292300e+02      9.16612880e+02
  308.       9.50543070e+02      9.84382610e+02      1.01808480e+03
  309.       1.05141200e+03      1.08465590e+03      1.11536850e+03
  310.       1.14532130e+03      1.17519080e+03      1.20491190e+03);
  311. L_SXZT = L_SXZT ET (PROG
  312.       1.23433030e+03      1.26367880e+03      1.28733980e+03
  313.       1.30749550e+03      1.32755070e+03      1.34750650e+03
  314.       1.36735080e+03      1.38697860e+03      1.40652480e+03
  315.       1.42600100e+03      1.44539400e+03      1.46449510e+03
  316.       1.48114470e+03      1.49774810e+03      1.51431050e+03);
  317. L_SXZT = L_SXZT ET (PROG
  318.       1.53083650e+03      1.54733000e+03      1.56379460e+03
  319.       1.58023360e+03      1.59664980e+03      1.61304570e+03
  320.       1.62942370e+03      1.64578580e+03      1.66213370e+03
  321.       1.67846920e+03      1.69479360e+03      1.71110830e+03
  322.       1.72741440e+03      1.74371300e+03      1.76000510e+03);
  323. L_SXZT = L_SXZT ET (PROG
  324.       1.77629150e+03      1.79257290e+03      1.80885010e+03
  325.       1.82512380e+03      1.84139450e+03      1.85766270e+03
  326.       1.87392890e+03      1.89019360e+03      1.90645720e+03
  327.       1.92272000e+03      1.93898240e+03      1.95524480e+03
  328.       1.97150730e+03      1.98777040e+03      2.00403420e+03);
  329. L_SXZT = L_SXZT ET (PROG
  330.       2.02029900e+03      2.03656490e+03      2.05283240e+03
  331.       2.06910140e+03      2.08537220e+03      2.10164500e+03
  332.       2.11791990e+03      2.13419710e+03      2.15047680e+03
  333.       2.16675900e+03      2.18304400e+03      2.19933170e+03
  334.       2.21562250e+03      2.23191620e+03      2.24821320e+03);
  335. L_SXZT = L_SXZT ET (PROG
  336.       2.26451340e+03      2.28081710e+03      2.29712420e+03
  337.       2.31343480e+03      2.32974920e+03      2.34606720e+03
  338.       2.36238910e+03      2.37871500e+03      2.39504470e+03
  339.       2.41137860e+03      2.42771660e+03      2.34778880e+03
  340.       2.26786090e+03      2.18793310e+03      2.10800520e+03);
  341. L_SXZT = L_SXZT ET (PROG
  342.       2.04283480e+03      1.99707210e+03      1.95141180e+03
  343.       1.90572340e+03      1.85994860e+03      1.81407540e+03
  344.       1.76811290e+03      1.72334380e+03      1.67878060e+03
  345.       1.63569170e+03      1.59266830e+03      1.54966200e+03
  346.       1.50663980e+03      1.46357950e+03      1.42046400e+03);
  347. L_SXZT = L_SXZT ET (PROG
  348.       1.37722510e+03      1.33653190e+03      1.29600770e+03
  349.       1.25607350e+03      1.21616880e+03      1.17627200e+03
  350.       1.13636690e+03      1.09644030e+03      1.05647860e+03
  351.       1.01643240e+03      9.78694220e+02      9.39895240e+02
  352.       9.02095010e+02      8.64372590e+02      8.25892880e+02);
  353. L_SXZT = L_SXZT ET (PROG
  354.       7.90553610e+02      7.53565560e+02      7.17343480e+02
  355.       6.81188090e+02      6.44171820e+02      6.10099840e+02
  356.       5.75797600e+02      5.42109140e+02      5.08434210e+02
  357.       4.74767080e+02      4.41103010e+02      4.07437980e+02
  358.       3.73768270e+02      3.40089410e+02      3.06390860e+02);
  359. L_SXZT = L_SXZT ET (PROG
  360.       2.73862150e+02      2.41176290e+02      2.10389400e+02
  361.       1.80620380e+02      1.50863640e+02      1.21115570e+02
  362.       9.13731080e+01      6.16336650e+01      3.18949020e+01
  363.       2.15442210e+00     -2.75911010e+01     -5.73508650e+01
  364.      -8.59856250e+01     -1.14697550e+02     -1.39112180e+02);
  365. L_SXZT = L_SXZT ET (PROG
  366.      -1.58847760e+02     -1.78571410e+02     -1.98285320e+02
  367.      -2.17991280e+02     -2.37690740e+02     -2.57384890e+02
  368.      -2.77074670e+02     -2.96760870e+02     -3.16444120e+02
  369.      -3.36125010e+02     -3.55804130e+02     -3.75482230e+02
  370.      -3.95160660e+02     -4.14843900e+02     -4.33736680e+02);
  371. L_SXZT = L_SXZT ET (PROG
  372.      -4.52617800e+02     -4.71540760e+02     -4.90675650e+02
  373.      -5.06641090e+02     -5.22597350e+02     -5.38546040e+02
  374.      -5.54488440e+02     -5.70425560e+02     -5.86358220e+02
  375.      -6.02287050e+02     -6.18212540e+02     -6.34135100e+02
  376.      -6.50055000e+02     -6.65972480e+02     -6.81887700e+02);
  377. L_SXZT = L_SXZT ET (PROG
  378.      -6.97800790e+02     -7.13711830e+02     -7.29620890e+02
  379.      -7.45528010e+02     -7.61433210e+02     -7.77336520e+02
  380.      -7.93237930e+02     -8.09137470e+02     -8.25035130e+02
  381.      -8.40930920e+02     -8.56824840e+02     -8.72716900e+02
  382.      -8.88607110e+02     -9.04495490e+02     -9.20382050e+02);
  383. L_SXZT = L_SXZT ET (PROG
  384.      -9.36266820e+02     -9.52149830e+02     -9.68031110e+02
  385.      -9.83910680e+02     -9.99788610e+02     -1.01566490e+03
  386.      -1.03153970e+03     -1.04741290e+03     -1.06328470e+03
  387.      -1.07915510e+03     -1.09502410e+03     -1.11089190e+03
  388.      -1.12675850e+03     -1.14262390e+03     -1.15848830e+03);
  389. L_SXZT = L_SXZT ET (PROG
  390.      -1.17435170e+03     -1.19021430e+03     -1.20607590e+03
  391.      -1.22193690e+03     -1.23779720e+03     -1.25365690e+03
  392.      -1.26951610e+03     -1.28537500e+03     -1.30123350e+03
  393.      -1.31709170e+03     -1.33294990e+03     -1.34880790e+03
  394.      -1.36466600e+03     -1.38052430e+03     -1.39638270e+03);
  395. L_SXZT = L_SXZT ET (PROG
  396.      -1.41224150e+03     -1.42810060e+03     -1.44396020e+03
  397.      -1.45982040e+03     -1.47568130e+03     -1.49154290e+03
  398.      -1.50740530e+03     -1.52326870e+03     -1.53913310e+03
  399.      -1.55499850e+03     -1.57086520e+03     -1.58673320e+03
  400.      -1.60260250e+03     -1.61847320e+03     -1.63434550e+03);
  401. L_SXZT = L_SXZT ET (PROG
  402.      -1.65021940e+03     -1.66609500e+03     -1.68197240e+03
  403.      -1.69785160e+03     -1.71373280e+03     -1.72961600e+03
  404.      -1.74550130e+03     -1.76138880e+03     -1.77727860e+03
  405.      -1.79317070e+03     -1.80906530e+03     -1.82496230e+03
  406.      -1.84086190e+03     -1.85676410e+03     -1.87266910e+03);
  407. L_SXZT = L_SXZT ET (PROG
  408.      -1.88857690e+03     -1.90448750e+03     -1.92040100e+03
  409.      -1.93631760e+03     -1.95223720e+03     -1.96816000e+03
  410.      -1.98408600e+03     -2.00001530e+03     -2.01594790e+03
  411.      -2.03188390e+03     -2.04782340e+03     -2.06376640e+03
  412.      -2.07971310e+03     -2.09566340e+03     -2.11161740e+03);
  413. L_SXZT = L_SXZT ET (PROG
  414.      -2.12757530e+03     -2.14353690e+03     -2.15950250e+03
  415.      -2.17547210e+03     -2.19144560e+03     -2.20742330e+03
  416.      -2.22340510e+03     -2.23939110e+03     -2.25538130e+03
  417.      -2.27137580e+03     -2.28737470e+03     -2.30337800e+03
  418.      -2.31938580e+03     -2.33539800e+03     -2.35141480e+03);
  419. L_SXZT = L_SXZT ET (PROG
  420.      -2.36743630e+03     -2.28770520e+03     -2.20797410e+03
  421.      -2.12824300e+03     -2.04851200e+03     -1.98715270e+03
  422.      -1.94203080e+03     -1.89699580e+03     -1.85191550e+03
  423.      -1.80673460e+03     -1.76144450e+03     -1.71604710e+03
  424.      -1.67222980e+03     -1.62870880e+03     -1.58623760e+03);
  425. L_SXZT = L_SXZT ET (PROG
  426.      -1.54383720e+03     -1.50145170e+03     -1.45904280e+03
  427.      -1.41658420e+03     -1.37404960e+03     -1.33258950e+03
  428.      -1.29085750e+03     -1.25150480e+03     -1.21218240e+03
  429.      -1.17286970e+03     -1.13355120e+03     -1.09421490e+03
  430.      -1.05485040e+03     -1.01544110e+03     -9.77155120e+02);
  431. L_SXZT = L_SXZT ET (PROG
  432.      -9.38238380e+02     -9.00744510e+02     -8.63320810e+02
  433.      -8.25867240e+02     -7.89527150e+02     -7.52675570e+02
  434.      -7.16763710e+02     -6.80914190e+02     -6.45012980e+02
  435.      -6.10394570e+02     -5.75611400e+02     -5.42467990e+02
  436.      -5.09337000e+02     -4.76212450e+02     -4.43089440e+02);
  437. L_SXZT = L_SXZT ET (PROG
  438.      -4.09963820e+02     -3.76831860e+02     -3.43689320e+02
  439.      -3.10528210e+02     -2.77286810e+02     -2.46141910e+02
  440.      -2.15038800e+02     -1.85769100e+02     -1.56509860e+02
  441.      -1.27257250e+02     -9.80080750e+01     -6.87596320e+01
  442.      -3.95094940e+01     -1.02552060e+01      1.90065120e+01);
  443. L_SXZT = L_SXZT ET (PROG
  444.       4.82831510e+01      7.66897660e+01      1.05112840e+02
  445.       1.31754240e+02      1.51123650e+02      1.70483010e+02
  446.       1.89834660e+02      2.09180480e+02      2.28522030e+02
  447.       2.47860530e+02      2.67196970e+02      2.86532150e+02
  448.       3.05866750e+02      3.25201370e+02      3.44536640e+02);
  449. L_SXZT = L_SXZT ET (PROG
  450.       3.63873390e+02      3.83213100e+02      4.02559970e+02
  451.       4.21212020e+02      4.39844700e+02      4.58502200e+02
  452.       4.77271770e+02      4.93693170e+02      5.09384420e+02
  453.       5.25069860e+02      5.40750690e+02      5.56427860e+02
  454.       5.72102140e+02      5.87774110e+02      6.03444230e+02);
  455. L_SXZT = L_SXZT ET (PROG
  456.       6.19112860e+02      6.34780260e+02      6.50446620e+02
  457.       6.66112100e+02      6.81776800e+02      6.97440800e+02
  458.       7.13104130e+02      7.28766830e+02      7.44428910e+02
  459.       7.60090380e+02      7.75751250e+02      7.91411500e+02
  460.       8.07071140e+02      8.22730160e+02      8.38388560e+02);
  461. L_SXZT = L_SXZT ET (PROG
  462.       8.54046330e+02      8.69703490e+02      8.85360040e+02
  463.       9.01015980e+02      9.16671330e+02      9.32326120e+02
  464.       9.47980360e+02      9.63634080e+02      9.79287300e+02
  465.       9.94940070e+02      1.01059240e+03      1.02624440e+03
  466.       1.04189600e+03      1.05754740e+03      1.07319850e+03);
  467. L_SXZT = L_SXZT ET (PROG
  468.       1.08884950e+03      1.10450030e+03      1.12015110e+03
  469.       1.13580180e+03      1.15145270e+03      1.16710360e+03
  470.       1.18275470e+03      1.19840610e+03      1.21405780e+03
  471.       1.22970990e+03      1.24536250e+03      1.26101560e+03
  472.       1.27666940e+03      1.29232380e+03      1.30797890e+03);
  473. L_SXZT = L_SXZT ET (PROG
  474.       1.32363500e+03      1.33929190e+03      1.35494980e+03
  475.       1.37060870e+03      1.38626880e+03      1.40193020e+03
  476.       1.41759280e+03      1.43325670e+03      1.44892210e+03
  477.       1.46458900e+03      1.48025760e+03      1.49592770e+03
  478.       1.51159960e+03      1.52727340e+03      1.54294900e+03);
  479. L_SXZT = L_SXZT ET (PROG
  480.       1.55862660e+03      1.57430620e+03      1.58998790e+03
  481.       1.60567180e+03      1.62135800e+03      1.63704640e+03
  482.       1.65273730e+03      1.66843070e+03      1.68412650e+03
  483.       1.69982500e+03      1.71552620e+03      1.73123010e+03
  484.       1.74693680e+03      1.76264640e+03      1.77835890e+03);
  485. L_SXZT = L_SXZT ET (PROG
  486.       1.79407440e+03      1.80979300e+03      1.82551470e+03
  487.       1.84123970e+03      1.85696790e+03      1.87269940e+03
  488.       1.88843430e+03      1.90417260e+03      1.91991450e+03
  489.       1.93565990e+03      1.95140900e+03      1.96716170e+03
  490.       1.98291820e+03      1.99867850e+03      2.01444260e+03);
  491. L_SXZT = L_SXZT ET (PROG
  492.       2.03021070e+03      2.04598270e+03      2.06175870e+03
  493.       2.07753880e+03      2.09332310e+03      2.10911150e+03
  494.       2.12490410e+03      2.14070110e+03      2.15650230e+03
  495.       2.17230800e+03      2.18811810e+03      2.20393260e+03
  496.       2.21975170e+03      2.23557540e+03      2.25140370e+03);
  497. L_SXZT = L_SXZT ET (PROG
  498.       2.26723670e+03      2.28307440e+03      2.29891680e+03
  499.       2.31476410e+03      2.33061620e+03      2.34647320e+03
  500.       2.26686710e+03      2.18726110e+03      2.10765510e+03
  501.       2.02804910e+03      1.96912920e+03      1.92442230e+03
  502.       1.87979080e+03      1.83510320e+03      1.79030690e+03);
  503. L_SXZT = L_SXZT ET (PROG
  504.       1.74539470e+03      1.70035570e+03      1.65720720e+03
  505.       1.61436420e+03      1.57229270e+03      1.53029560e+03
  506.       1.48831170e+03      1.44629900e+03      1.40422740e+03
  507.       1.36205970e+03      1.32127700e+03      1.28050720e+03
  508.       1.24154620e+03      1.20261670e+03      1.16369580e+03);
  509. L_SXZT = L_SXZT ET (PROG
  510.       1.12476620e+03      1.08581420e+03      1.04682690e+03
  511.       1.00777620e+03      9.70239200e+02      9.31996960e+02
  512.       8.95022090e+02      8.58123340e+02      8.21076800e+02
  513.       7.86004300e+02      7.49958440e+02      7.14576920e+02
  514.       6.79260390e+02      6.43705940e+02      6.10360700e+02);
  515. L_SXZT = L_SXZT ET (PROG
  516.       5.76674250e+02      5.43878600e+02      5.11094190e+02
  517.       4.78314260e+02      4.45533230e+02      4.12746320e+02
  518.       3.79948990e+02      3.47135530e+02      3.14292570e+02
  519.       2.82532720e+02      2.50665450e+02      2.20351380e+02
  520.       1.91388230e+02      1.62433480e+02      1.33482850e+02);
  521. L_SXZT = L_SXZT ET (PROG
  522.       1.04532750e+02      7.55800860e+01      4.66219880e+01
  523.       1.76553310e+01     -1.13246870e+01     -4.03327340e+01
  524.      -6.81639330e+01     -9.60781200e+01     -1.19788730e+02
  525.      -1.38951110e+02     -1.58105570e+02     -1.77254600e+02
  526.      -1.96400220e+02     -2.15544050e+02     -2.34687400e+02);
  527. L_SXZT = L_SXZT ET (PROG
  528.      -2.53831310e+02     -2.72976650e+02     -2.92124160e+02
  529.      -3.11274530e+02     -3.30428550e+02     -3.49587380e+02
  530.      -3.68753360e+02     -3.87935860e+02     -4.06257040e+02
  531.      -4.24571520e+02     -4.42940410e+02     -4.61432070e+02
  532.      -4.77005150e+02     -4.92573330e+02     -5.08138020e+02);
  533. L_SXZT = L_SXZT ET (PROG
  534.      -5.23700300e+02     -5.39261030e+02     -5.54820900e+02
  535.      -5.70380420e+02     -5.85939980e+02     -6.01499890e+02
  536.      -6.17060360e+02     -6.32621560e+02     -6.48183590e+02
  537.      -6.63746530e+02     -6.79310410e+02     -6.94875260e+02
  538.      -7.10441080e+02     -7.26007860e+02     -7.41575600e+02);
  539. L_SXZT = L_SXZT ET (PROG
  540.      -7.57144260e+02     -7.72713840e+02     -7.88284300e+02
  541.      -8.03855620e+02     -8.19427790e+02     -8.35000780e+02
  542.      -8.50574580e+02     -8.66149180e+02     -8.81724580e+02
  543.      -8.97300770e+02     -9.12877750e+02     -9.28455530e+02
  544.      -9.44034110e+02     -9.59613510e+02     -9.75193760e+02);
  545. L_SXZT = L_SXZT ET (PROG
  546.      -9.90774850e+02     -1.00635680e+03     -1.02193970e+03
  547.      -1.03752360e+03     -1.05310840e+03     -1.06869420e+03
  548.      -1.08428100e+03     -1.09986900e+03     -1.11545810e+03
  549.      -1.13104840e+03     -1.14663980e+03     -1.16223260e+03
  550.      -1.17782670e+03     -1.19342220e+03     -1.20901910e+03);
  551. L_SXZT = L_SXZT ET (PROG
  552.      -1.22461740e+03     -1.24021730e+03     -1.25581880e+03
  553.      -1.27142200e+03     -1.28702680e+03     -1.30263350e+03
  554.      -1.31824190e+03     -1.33385220e+03     -1.34946450e+03
  555.      -1.36507880e+03     -1.38069510e+03     -1.39631350e+03
  556.      -1.41193420e+03     -1.42755700e+03     -1.44318220e+03);
  557. L_SXZT = L_SXZT ET (PROG
  558.      -1.45880970e+03     -1.47443970e+03     -1.49007210e+03
  559.      -1.50570710e+03     -1.52134460e+03     -1.53698480e+03
  560.      -1.55262780e+03     -1.56827350e+03     -1.58392200e+03
  561.      -1.59957350e+03     -1.61522780e+03     -1.63088520e+03
  562.      -1.64654570e+03     -1.66220930e+03     -1.67787600e+03);
  563. L_SXZT = L_SXZT ET (PROG
  564.      -1.69354600e+03     -1.70921930e+03     -1.72489590e+03
  565.      -1.74057590e+03     -1.75625940e+03     -1.77194640e+03
  566.      -1.78763700e+03     -1.80333110e+03     -1.81902900e+03
  567.      -1.83473050e+03     -1.85043580e+03     -1.86614500e+03
  568.      -1.88185800e+03     -1.89757490e+03     -1.91329580e+03);
  569. L_SXZT = L_SXZT ET (PROG
  570.      -1.92902080e+03     -1.94474980e+03     -1.96048290e+03
  571.      -1.97622020e+03     -1.99196170e+03     -2.00770750e+03
  572.      -2.02345750e+03     -2.03921200e+03     -2.05497080e+03
  573.      -2.07073410e+03     -2.08650180e+03     -2.10227410e+03
  574.      -2.11805100e+03     -2.13383240e+03     -2.14961860e+03);
  575. L_SXZT = L_SXZT ET (PROG
  576.      -2.16540940e+03     -2.18120500e+03     -2.19700540e+03
  577.      -2.21281060e+03     -2.22862060e+03     -2.24443560e+03
  578.      -2.26025550e+03     -2.27608040e+03     -2.29191030e+03
  579.      -2.30774530e+03     -2.32358530e+03     -2.24402140e+03
  580.      -2.16445750e+03     -2.08489350e+03     -2.00532960e+03);
  581. L_SXZT = L_SXZT ET (PROG
  582.      -1.94724250e+03     -1.90268340e+03     -1.85819590e+03
  583.      -1.81364970e+03     -1.76899320e+03     -1.72421920e+03
  584.      -1.67930860e+03     -1.63641850e+03     -1.59387320e+03
  585.      -1.55194360e+03     -1.51008990e+03     -1.46824920e+03
  586.      -1.42637790e+03     -1.38444480e+03     -1.34240720e+03);
  587. L_SXZT = L_SXZT ET (PROG
  588.      -1.30189260e+03     -1.26146370e+03     -1.22264030e+03
  589.      -1.18384860e+03     -1.14506490e+03     -1.10627130e+03
  590.      -1.06745340e+03     -1.02859700e+03     -9.89668210e+02
  591.      -9.52430140e+02     -9.14374280e+02     -8.77582230e+02
  592.      -8.40867130e+02     -8.03888640e+02     -7.69333970e+02);
  593. L_SXZT = L_SXZT ET (PROG
  594.      -7.33442020e+02     -6.98245910e+02     -6.63113670e+02
  595.      -6.27509130e+02     -5.94579260e+02     -5.61227250e+02
  596.      -5.28539390e+02     -4.95861270e+02     -4.63185870e+02
  597.      -4.30507380e+02     -3.97820770e+02     -3.65121180e+02
  598.      -3.32402260e+02     -2.99647780e+02     -2.68105190e+02);
  599. L_SXZT = L_SXZT ET (PROG
  600.      -2.36402690e+02     -2.06501080e+02     -1.77617570e+02
  601.      -1.48740500e+02     -1.19865460e+02     -9.09887150e+01
  602.      -6.21070360e+01     -3.32173970e+01     -4.31641630e+00
  603.       2.46013210e+01      5.35544060e+01      8.11996960e+01
  604.       1.08966280e+02      1.31829330e+02      1.50957680e+02);
  605. L_SXZT = L_SXZT ET (PROG
  606.       1.70080070e+02      1.89199010e+02      2.08316540e+02
  607.       2.27434320e+02      2.46553680e+02      2.65675690e+02
  608.       2.84801230e+02      3.03931080e+02      3.23065960e+02
  609.       3.42206720e+02      3.61354610e+02      3.80512310e+02
  610.       3.99692140e+02      4.17947610e+02      4.36204660e+02);
  611. L_SXZT = L_SXZT ET (PROG
  612.       4.54530790e+02      4.72728290e+02      4.88315540e+02
  613.       5.03900400e+02      5.19484110e+02      5.35067680e+02
  614.       5.50651890e+02      5.66237330e+02      5.81824450e+02
  615.       5.97413600e+02      6.13005030e+02      6.28598900e+02
  616.       6.44195350e+02      6.59794440e+02      6.75396210e+02);
  617. L_SXZT = L_SXZT ET (PROG
  618.       6.91000670e+02      7.06607820e+02      7.22217630e+02
  619.       7.37830080e+02      7.53445120e+02      7.69062710e+02
  620.       7.84682800e+02      8.00305360e+02      8.15930350e+02
  621.       8.31557710e+02      8.47187410e+02      8.62819430e+02
  622.       8.78453720e+02      8.94090260e+02      9.09729030e+02);
  623. L_SXZT = L_SXZT ET (PROG
  624.       9.25370020e+02      9.41013200e+02      9.56658580e+02
  625.       9.72306140e+02      9.87955880e+02      1.00360780e+03
  626.       1.01926190e+03      1.03491820e+03      1.05057680e+03
  627.       1.06623750e+03      1.08190050e+03      1.09756580e+03
  628.       1.11323340e+03      1.12890330e+03      1.14457550e+03);
  629. L_SXZT = L_SXZT ET (PROG
  630.       1.16025010e+03      1.17592710e+03      1.19160650e+03
  631.       1.20728840e+03      1.22297280e+03      1.23865970e+03
  632.       1.25434930e+03      1.27004140e+03      1.28573620e+03
  633.       1.30143370e+03      1.31713400e+03      1.33283700e+03
  634.       1.34854280e+03      1.36425160e+03      1.37996320e+03);
  635. L_SXZT = L_SXZT ET (PROG
  636.       1.39567780e+03      1.41139540e+03      1.42711600e+03
  637.       1.44283970e+03      1.45856660e+03      1.47429660e+03
  638.       1.49002990e+03      1.50576640e+03      1.52150620e+03
  639.       1.53724940e+03      1.55299600e+03      1.56874610e+03
  640.       1.58449970e+03      1.60025670e+03      1.61601740e+03);
  641. L_SXZT = L_SXZT ET (PROG
  642.       1.63178170e+03      1.64754970e+03      1.66332140e+03
  643.       1.67909690e+03      1.69487610e+03      1.71065930e+03
  644.       1.72644630e+03      1.74223720e+03      1.75803220e+03
  645.       1.77383110e+03      1.78963410e+03      1.80544130e+03
  646.       1.82125250e+03      1.83706800e+03      1.85288770e+03);
  647. L_SXZT = L_SXZT ET (PROG
  648.       1.86871170e+03      1.88453990e+03      1.90037260e+03
  649.       1.91620960e+03      1.93205110e+03      1.94789700e+03
  650.       1.96374740e+03      1.97960240e+03      1.99546200e+03
  651.       2.01132620e+03      2.02719510e+03      2.04306860e+03
  652.       2.05894690e+03      2.07483000e+03      2.09071790e+03);
  653. L_SXZT = L_SXZT ET (PROG
  654.       2.10661070e+03      2.12250830e+03      2.13841080e+03
  655.       2.15431830e+03      2.17023080e+03      2.18614830e+03
  656.       2.20207080e+03      2.21799850e+03      2.23393120e+03
  657.       2.24986910e+03      2.26581220e+03      2.28176060e+03
  658.       2.29771410e+03      2.31367300e+03      2.32963720e+03);
  659. L_SXZT = L_SXZT ET (PROG
  660.       2.34560670e+03      2.26600390e+03      2.18640120e+03
  661.       2.10679840e+03      2.02719570e+03      1.96836960e+03
  662.       1.92369750e+03      1.87910140e+03      1.83445220e+03
  663.       1.78969740e+03      1.74482860e+03      1.69983100e+03
  664.       1.65678260e+03      1.61403870e+03      1.57199880e+03);
  665. L_SXZT = L_SXZT ET (PROG
  666.       1.53003390e+03      1.48808340e+03      1.44610490e+03
  667.       1.40406830e+03      1.36193570e+03      1.32124930e+03
  668.       1.28054770e+03      1.24161310e+03      1.20270940e+03
  669.       1.16381380e+03      1.12490870e+03      1.08598040e+03
  670.       1.04701570e+03      1.00798600e+03      9.70539260e+02);
  671. L_SXZT = L_SXZT ET (PROG
  672.       9.32441050e+02      8.95467630e+02      8.58569690e+02
  673.       8.21507750e+02      7.86565740e+02      7.50602190e+02
  674.       7.15206550e+02      6.79874470e+02      6.44277870e+02
  675.       6.11031600e+02      5.77417810e+02      5.44603120e+02
  676.       5.11796710e+02      4.78991810e+02      4.46182840e+02);
  677. L_SXZT = L_SXZT ET (PROG
  678.       4.13364990e+02      3.80533700e+02      3.47683300e+02
  679.       3.14800650e+02      2.83037460e+02      2.51147750e+02
  680.       2.20775160e+02      1.91752230e+02      1.62734090e+02
  681.       1.33716480e+02      1.04695840e+02      7.56690770e+01
  682.       4.66333630e+01      1.75856180e+01     -1.14787850e+01);
  683. L_SXZT = L_SXZT ET (PROG
  684.      -4.05736840e+01     -6.84719410e+01     -9.64657050e+01
  685.      -1.20267270e+02     -1.39526190e+02     -1.58780670e+02
  686.      -1.78033180e+02     -1.97285680e+02     -2.16539800e+02
  687.      -2.35796840e+02     -2.55057840e+02     -2.74323680e+02
  688.      -2.93595080e+02     -3.12872730e+02     -3.32157410e+02);
  689. L_SXZT = L_SXZT ET (PROG
  690.      -3.51450180e+02     -3.70753160e+02     -3.90074380e+02
  691.      -4.08535530e+02     -4.26997830e+02     -4.45519570e+02
  692.      -4.64186830e+02     -4.79917950e+02     -4.95648440e+02
  693.      -5.11379550e+02     -5.27112260e+02     -5.42847320e+02
  694.      -5.58585300e+02     -5.74326650e+02     -5.90071690e+02);
  695. L_SXZT = L_SXZT ET (PROG
  696.      -6.05820640e+02     -6.21573650e+02     -6.37330830e+02
  697.      -6.53092220e+02     -6.68857850e+02     -6.84627700e+02
  698.      -7.00401750e+02     -7.16179950e+02     -7.31962240e+02
  699.      -7.47748580e+02     -7.63538890e+02     -7.79333110e+02
  700.      -7.95131170e+02     -8.10933020e+02     -8.26738580e+02);
  701. L_SXZT = L_SXZT ET (PROG
  702.      -8.42547800e+02     -8.58360630e+02     -8.74177020e+02
  703.      -8.89996920e+02     -9.05820280e+02     -9.21647080e+02
  704.      -9.37477270e+02     -9.53310830e+02     -9.69147740e+02
  705.      -9.84987980e+02     -1.00083150e+03     -1.01667840e+03
  706.      -1.03252850e+03     -1.04838190e+03     -1.06423860e+03);
  707. L_SXZT = L_SXZT ET (PROG
  708.      -1.08009860e+03     -1.09596190e+03     -1.11182840e+03
  709.      -1.12769830e+03     -1.14357150e+03     -1.15944790e+03
  710.      -1.17532780e+03     -1.19121100e+03     -1.20709760e+03
  711.      -1.22298750e+03     -1.23888090e+03     -1.25477770e+03
  712.      -1.27067800e+03     -1.28658180e+03     -1.30248910e+03);
  713. L_SXZT = L_SXZT ET (PROG
  714.      -1.31840000e+03     -1.33431440e+03     -1.35023240e+03
  715.      -1.36615400e+03     -1.38207940e+03     -1.39800840e+03
  716.      -1.41394110e+03     -1.42987760e+03     -1.44581780e+03
  717.      -1.46176190e+03     -1.47770980e+03     -1.49366160e+03
  718.      -1.50961740e+03     -1.52557700e+03     -1.54154070e+03);
  719. L_SXZT = L_SXZT ET (PROG
  720.      -1.55750840e+03     -1.57348010e+03     -1.58945590e+03
  721.      -1.60543580e+03     -1.62141990e+03     -1.63740820e+03
  722.      -1.65340060e+03     -1.66939740e+03     -1.68539840e+03
  723.      -1.70140370e+03     -1.71741340e+03     -1.73342750e+03
  724.      -1.74944600e+03     -1.76546890e+03     -1.78149630e+03);
  725. L_SXZT = L_SXZT ET (PROG
  726.      -1.79752830e+03     -1.81356480e+03     -1.82960580e+03
  727.      -1.84565150e+03     -1.86170190e+03     -1.87775690e+03
  728.      -1.89381660e+03     -1.90988110e+03     -1.92595030e+03
  729.      -1.94202430e+03     -1.95810320e+03     -1.97418700e+03
  730.      -1.99027560e+03     -2.00636920e+03     -2.02246770e+03);
  731. L_SXZT = L_SXZT ET (PROG
  732.      -2.03857120e+03     -2.05467980e+03     -2.07079330e+03
  733.      -2.08691200e+03     -2.10303580e+03     -2.11916470e+03
  734.      -2.13529870e+03     -2.15143800e+03     -2.16758240e+03
  735.      -2.18373210e+03     -2.19988710e+03     -2.21604740e+03
  736.      -2.23221300e+03     -2.24838400e+03     -2.26456040e+03);
  737. L_SXZT = L_SXZT ET (PROG
  738.      -2.28074220e+03     -2.29692940e+03     -2.31312200e+03
  739.      -2.32932020e+03     -2.34552390e+03     -2.36173310e+03
  740.      -2.28202410e+03     -2.20231510e+03     -2.12260620e+03
  741.      -2.04289720e+03     -1.98202240e+03     -1.93701810e+03
  742.      -1.89209980e+03     -1.84713990e+03     -1.80208420e+03);
  743. L_SXZT = L_SXZT ET (PROG
  744.      -1.75692230e+03     -1.71164920e+03     -1.66810010e+03
  745.      -1.62473930e+03     -1.58237790e+03     -1.54008890e+03
  746.      -1.49781650e+03     -1.45552160e+03     -1.41317730e+03
  747.      -1.37075480e+03     -1.32954940e+03     -1.28811140e+03
  748.      -1.24885820e+03     -1.20963430e+03     -1.17041880e+03);
  749. L_SXZT = L_SXZT ET (PROG
  750.      -1.13119570e+03     -1.09195280e+03     -1.05267870e+03
  751.      -1.01335530e+03     -9.75310730e+02     -9.36713540e+02
  752.      -8.99287670e+02     -8.61931690e+02     -8.24528580e+02
  753.      -7.88531980e+02     -7.51927490e+02     -7.16057670e+02
  754.      -6.80248230e+02     -6.44362450e+02     -6.10066870e+02);
  755. L_SXZT = L_SXZT ET (PROG
  756.      -5.75710520e+02     -5.42577660e+02     -5.09451500e+02
  757.      -4.76325960e+02     -4.43195980e+02     -4.10057300e+02
  758.      -3.76906030e+02     -3.43737790e+02     -3.10544340e+02
  759.      -2.78450010e+02     -2.45976840e+02     -2.14661020e+02
  760.      -1.85322180e+02     -1.55986210e+02     -1.26649390e+02);
  761. L_SXZT = L_SXZT ET (PROG
  762.      -9.73085700e+01     -6.79610740e+01     -3.86045080e+01
  763.      -9.23644540e+00      2.01463280e+01      4.95512180e+01
  764.       7.80079940e+01      1.06509910e+02      1.32712670e+02
  765.       1.52234580e+02      1.71753460e+02      1.91271570e+02
  766.       2.10790740e+02      2.30312490e+02      2.49838010e+02);
  767. L_SXZT = L_SXZT ET (PROG
  768.       2.69368290e+02      2.88904130e+02      3.08446190e+02
  769.       3.27995070e+02      3.47551360e+02      3.67115820e+02
  770.       3.86689760e+02      4.06276810e+02      4.25144500e+02
  771.       4.44008510e+02      4.62911230e+02      4.81940030e+02
  772.       4.98703860e+02      5.14683880e+02      5.30666070e+02);
  773. L_SXZT = L_SXZT ET (PROG
  774.       5.46651380e+02      5.62640560e+02      5.78634160e+02
  775.       5.94632610e+02      6.10636200e+02      6.26645160e+02
  776.       6.42659610e+02      6.58679640e+02      6.74705280e+02
  777.       6.90736530e+02      7.06773370e+02      7.22815740e+02
  778.       7.38863580e+02      7.54916830e+02      7.70975400e+02);
  779. L_SXZT = L_SXZT ET (PROG
  780.       7.87039220e+02      8.03108190e+02      8.19182230e+02
  781.       8.35261260e+02      8.51345210e+02      8.67433990e+02
  782.       8.83527530e+02      8.99625760e+02      9.15728630e+02
  783.       9.31836060e+02      9.47948020e+02      9.64064430e+02
  784.       9.80185270e+02      9.96310490e+02      1.01244010e+03);
  785. L_SXZT = L_SXZT ET (PROG
  786.       1.02857390e+03      1.04471210e+03      1.06085450e+03
  787.       1.07700120e+03      1.09315200e+03      1.10930710e+03
  788.       1.12546640e+03      1.14162990e+03      1.15779760e+03
  789.       1.17396940e+03      1.19014540e+03      1.20632560e+03
  790.       1.22251000e+03      1.23869860e+03      1.25489140e+03);
  791. L_SXZT = L_SXZT ET (PROG
  792.       1.27108840e+03      1.28728960e+03      1.30349500e+03
  793.       1.31970470e+03      1.33591860e+03      1.35213690e+03
  794.       1.36835940e+03      1.38458620e+03      1.40081730e+03
  795.       1.41705280e+03      1.43329260e+03      1.44953690e+03
  796.       1.46578550e+03      1.48203860e+03      1.49829610e+03);
  797. L_SXZT = L_SXZT ET (PROG
  798.       1.51455800e+03      1.53082450e+03      1.54709550e+03
  799.       1.56337100e+03      1.57965110e+03      1.59593580e+03
  800.       1.61222510e+03      1.62851900e+03      1.64481760e+03
  801.       1.66112090e+03      1.67742880e+03      1.69374150e+03
  802.       1.71005900e+03      1.72638120e+03      1.74270830e+03);
  803. L_SXZT = L_SXZT ET (PROG
  804.       1.75904020e+03      1.77537690e+03      1.79171850e+03
  805.       1.80806510e+03      1.82441650e+03      1.84077290e+03
  806.       1.85713430e+03      1.87350070e+03      1.88987210e+03
  807.       1.90624860e+03      1.92263020e+03      1.93901680e+03
  808.       1.95540860e+03      1.97180560e+03      1.98820770e+03);
  809. L_SXZT = L_SXZT ET (PROG
  810.       2.00461500e+03      2.02102750e+03      2.03744530e+03
  811.       2.05386830e+03      2.07029670e+03      2.08673030e+03
  812.       2.10316930e+03      2.11961370e+03      2.13606340e+03
  813.       2.15251860e+03      2.16897910e+03      2.18544520e+03
  814.       2.20191670e+03      2.21839370e+03      2.23487620e+03);
  815. L_SXZT = L_SXZT ET (PROG
  816.       2.25136430e+03      2.26785790e+03      2.28435710e+03
  817.       2.30086190e+03      2.31737240e+03      2.33388850e+03
  818.       2.35041020e+03      2.36693770e+03      2.38347080e+03
  819.       2.40000970e+03      2.41655430e+03);
  820. *
  821. *********************************************************************
  822. *   Test
  823. *********************************************************************
  824. *
  825. EV_TEST = EVOL (VERT) MANU 'Exz []' L_EXZ 'Sxz [Pa]' L_SXZT;
  826. SI (VAR_GRAP EGA 1);
  827.   DESS (EV_1 ET EV_TEST);
  828. FINSI;
  829. VAL_TEST = MAXI (ABS(L_SXZ - L_SXZT));
  830. VAL_TEST = VAL_TEST/(MAXI (ABS(L_SXZ)));
  831. *
  832. SI (VAL_TEST > 1.E-5);
  833.   ERRE 5;
  834. FINSI;
  835. *
  836. FIN;
  837.  
  838.  
  839.  
  840.  
  841.  
  842.  
  843.  
  844.  
  845.  
  846.  
  847.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales