高斯激光热源下的n层材料瞬态热传导方程

xy112

高斯激光热源下的n层材料瞬态热传导方程
高斯, 激光, 热源, 方程, 热传导高斯, 激光, 热源, 方程, 热传导
出于需要， 写了个 n层材料的瞬态热传导方程， 由表面激光加热。每层的热源由其光学性质N= n + ik决定. 通过解Maxwell方程获得激光在材料中的强度分布， 由多项式 a0+a1*y+a2*y^2+a3*y^3+a4*y^4拟和获得。
P=Power/r^2/{PI}*exp(-({X}^2+({Z}-V*{TIME})^2)/r^2)*( A0+A1*(-{Y}/1e-9)+A2*(-{Y}/1e-9)^2+A3*(-{Y}/1e-9)^3+ A4*(-{Y}/1e-9)^4).
  
高斯分布体热源在程序中由HG.mac加载。不同层的a0,a1,a2,a3,a4决定了其内部的热源分布。
  
本程序可已根据问题选择所需的材料层数，所有参数都可以改变。  
  
) T: U% w0 p3 n: t  Y以下是命令流，如果有不合适的地方，欢迎大家提出、讨论！
  
FINISH  ! Make sure we are at BEGIN level   
! t: y, k9 p' p/ A7 q  d: t/Clear
/FILNAME,Media,0  
  
/PREP7
; s$ u/ j3 a: Z8 s# H8 I/PNUM,VOLU,1
  
! In simulation, the X-Z is the media plane, while Y is the depth direction
* x. B" q, @1 ~$ p3 \4 Y- G3 C  
4 X! R( F7 ]; g8 P% xLayer=5                       !Total layer number
nm= 1e-9                    !length unit
" x5 f9 p+ V$ i- KPower=0.005               ! Laser incident power   W
r=5e-7                        ! Laser spot radius      m
2 k/ w3 B5 t3 C, VV=0                           ! Laser moving veolcity, along z direction,  m/s
  
pw=Power*1000                 !For display Title
rh=r/1e-6                           !For display Title
  
/TITLE,%Layer% layer media, Laser Power=%pw%mW, Heating radius=%rh%um, Velocity=%V%m/s  
2 e% k' ~$ g/ w  b& i2 {/ N$ m6 J  
6 P' f& H4 @1 F*DIM,Thickness,ARRAY,Layer    ! Array input for each layer thickness
*DIM,k_L,ARRAY,Layer          ! Array input for each layer lateral thermal conductivity  
. @; q- y7 Q  f* n*DIM,k_N,ARRAY,Layer          ! Array input for each layer normal thermal conductivity
  `1 e4 x5 S2 @3 ?# i% R' S*DIM,Cp,ARRAY,Layer           ! Array input for each layer heat capacity
2 `6 j/ n  N$ M9 }2 V7 H$ ]6 c) L*DIM,Den,ARRAY,Layer          ! Array input for each layer density
. ^) _$ v! e1 R- g! d*DIM,div,ARRAY,Layer          ! Array input for each layer thickness division
*DIM,dscale,ARRAY,Layer       ! Array input for each layer thickness division scale
- |' T* X# r: Z  A; h' L*DIM,A0,ARRAY,Layer           ! Array input for each layer a0 coefficient for heat generation curve
*DIM,A1,ARRAY,Layer           ! Array input for each layer a1 coefficient for heat generation curve
*DIM,A2,ARRAY,Layer           ! Array input for each layer a2 coefficient for heat generation curve
*DIM,A3,ARRAY,Layer           ! Array input for each layer a3 coefficient for heat generation curve
- F) B0 L. J" Z  p" k8 v) M9 {*DIM,A4,ARRAY,Layer           ! Array input for each layer a4 coefficient for heat generation curve
% |! g) D6 L9 a2 q' s" Y4 {6 ]3 {  
scale=10000                   !for vovlap operation
Pulse_Duration=1e-8           !Laser pulse duration
End_Time=2e-8                 !Final time for simulation
. G) k  G  r* }! h1 XStep=20                       !Time interval  
T0=50                         !Environment temperature   
+ e3 T+ r& r* M) w& ]1 X7 d8 Q6 |+ o  
Length=4000*nm                !Length of the block in z direction
2 @/ c2 D( O% H9 x4 s# @2 B3 bWidth=2000*nm                 !Width of the block in x direction
  
!Layer properties, from top to bottom  
  
!First layer thermal property
( ~% U3 ]4 Y5 L( b) Y8 h: `Num=1                        
7 o$ l" h, B  z% Nk_L(Num)=1.4                  
k_N(Num)=1.4
Den(Num)=4000
* p$ ]* C: M* N! `7 bCp(Num)=500
Thickness(Num)=50*nm
div(Num)=5
dscale(Num)=1                 !last/first in -Y direction
# C. z* J8 j" E9 o, ~) O' LA0(Num)=0
5 A6 Y( G7 @, j% z: T; yA1(Num)=0
A2(Num)=0
A3(Num)=0
, W" y, X3 R& e2 q+ B0 D5 i6 C+ g( mA4(Num)=0
& b0 r5 a9 o  l: i% T% w  
!Second layer thermal property
# E& S/ N2 D% d+ u/ I1 V: ~Num=2
$ |2 A1 a! l8 h" v+ sk_L(Num)=16
5 o- b. ?9 i/ w: ck_N(Num)=16
1 V6 [* k( N9 W6 r) t" s2 YDen(Num)=3000
Cp(Num)=900
Thickness(Num)=20*nm
5 M" o% S- |  t1 R3 ^/ h, p& Cdiv(Num)=2
: v/ {, K! _) r% ?dscale(Num)=1
A0(Num)=807.929E+6
A1(Num)=-43.896E+6
/ X5 @% g6 r6 u9 y8 \A2(Num)=934.529E+3
3 ~3 |9 H3 N5 ^A3(Num)=-8.864E+3
A4(Num)=31.703
  
!Third layer thermal property
Num=3
k_L(Num)=1.5
k_N(Num)=1.5
& t( z# m% o5 O% _8 q1 P; n. IDen(Num)=4000
8 {0 n2 n' N% j8 G5 k9 ]- b9 ICp(Num)=500
- @4 g5 k1 H0 D1 u( y* ~& [/ a* fThickness(Num)=100*nm
div(Num)=5
dscale(Num)=1
A0(Num)=0
: Z; f, p0 M" x% {. A* [! s6 XA1(Num)=0
A2(Num)=0
A3(Num)=0
A4(Num)=0
# ^) p# P. I0 P7 d  
" Z( N8 f3 _' v1 p" \5 u! c- x!Fourth layer thermal property
$ Z6 n6 Y) l, [5 uNum=4
' B+ b5 {/ W3 b* \2 p" t1 tk_L(Num)=100
+ }; V, ^8 U- U. ]k_N(Num)=100
3 y$ w, v" B' H! qDen(Num)=3900
) k1 s( n$ I' m! l9 v8 L/ G, z( CCp(Num)=700
Thickness(Num)=50*nm
div(Num)=2
dscale(Num)=2
* T% P3 q$ M0 `8 W% ^$ A" Y, fA0(Num)=1.627E+9
A1(Num)=-32.146E+6
A2(Num)=238.195E+3
A3(Num)=-784.173
A4(Num)=967.716E-3
4 U* Y$ G, x* V& [% C6 F  
!Fifth layer thermal property
1 F5 H1 F% s; M$ X/ W, P$ \Num=5
) J, t- Z+ }3 C' ^- f% R  Ck_L(Num)=0.2
% Q. ~/ q. [. c$ [k_N(Num)=0.2
Den(Num)=3900
! E5 k! t5 O6 e3 s  N, C; ICp(Num)=450
6 V# ^; T! G: Y& pThickness(Num)=100*nm
% v5 }: q, W* |" c4 {: y  l$ rdiv(Num)=5
dscale(Num)=2
A0(Num)=0
A1(Num)=0
A2(Num)=0
A3(Num)=0
; _3 K0 i9 F2 j% d1 JA4(Num)=0
( G8 M/ o5 J: J- M  
! Define the thermal element
7 N# _+ e: U: d* D' Q3 w7 J0 ZET,1,SOLID70  
% p* d+ {' u. c3 P( m* F7 nKEYOPT,1,2,0  
& h2 l9 J5 p9 k) qKEYOPT,1,4,0  
' v+ v& o5 v! fKEYOPT,1,7,0  
, b$ ?" P' Q! M" S$ f0 sKEYOPT,1,8,0  
  
!For post-precessing
nPoint=200
" X6 y8 s! S- R% A/ O: Z9 bdisplay_r=1e-6
' @% r" ]& R) Z7 o  
3 `  D* Y% b: s& rTSum1=0
+ U$ b. ]& ?# o2 N8 {' oTSum2=0
*do,Num,1,Layer,1
  UIMP,Num,DENS, , ,Den(Num),                                        !Define density of material Num
  UIMP,Num,KXX,KYY,KZZ,k_L(Num),k_N(Num),k_L(Num),                   !Define thermal conductivity of material Num
7 V7 x4 @  Q4 K( T  UIMP,Num,C, , ,Cp(Num),                                            !Define heat capacity of material Num   
  TSum1=TSum1+Thickness(Num)  
; n) d, ?% k/ F4 c2 q  BLOCK,-Width/2,Width/2,-TSum1,-TSum2,-Length/2,Length/2,           !Draw the layer of Num  
# X8 N: r  |1 y  TSum2=TSum2+Thickness(Num)
: {5 y; l0 B- ~) l*Enddo
  
0 T; u# p* G) I2 Jd_total=TSum1                 !Total thickness of the multilayer                                                  
  
VLSCALE,ALL, , , scale, scale, scale,  ,1,1,
5 b, E5 d* g8 i+ t" i  
*IF,Layer,GT,1,THEN
   VGLUE,ALL
*ENDIF
  
VLSCALE,ALL, , , 1/scale, 1/scale, 1/scale, ,1,1,
: M- h0 O. D+ b9 h" g) e  
TSum1=0
, k8 e) V, W6 s6 S$ yTSum2=0
  
Size=100*nm                             !Area mesh size, in X-Z plane
aesize,all,Size,  
- \2 g* [& e- X) G/ Q0 r  
*do,Num,1,Layer,1
% @6 i/ `8 f0 M9 N   TSum1=TSum1+THICKNESS(Num)
9 V, B- \1 y/ @2 r   lsel,all
. I. U, P/ @6 I: R$ B' @  
: Q; E+ w1 g0 O  W4 j8 N1 y- f   !Determine the mesh scale and divison for each layer at Y direction
' R( T' F+ m$ W2 T% G( X5 G   *do,i,-1,1,2
      *do, j,-1,1,2
/ O1 n, M+ {5 W5 R! ^         lsel,s,loc,y,-TSum1,-TSum2
) h$ r; H  U' s! N6 B         lsel,r,loc,x,i*Width/2,i*Width/2
; \& Q; Y5 N- [3 e, y         lsel,r,loc,z,j*Length/2,j*Length/2
         *GET,Line_Num,LINE, 0, NUM,MIN,
; Z9 |$ N! |5 W; T/ L) Q( {3 _1 J$ q         *GET,KP1,LINE,Line_Num,KP,1,
         *GET,KP2,LINE,Line_Num,KP,2,
         *GET,KP1_Y,KP,KP1, LOC, Y,  
5 _2 {3 J1 ^" N         *GET,KP2_Y,KP,KP2, LOC, Y,
. M! T% c) k# T9 W1 k6 w         *IF,KP1_Y,GT,KP2_Y,THEN
             lesize,all,THICKNESS(Num)/div(Num), , ,dscale(Num),1, , ,1,
% B/ O5 y& D, [: z8 w/ x9 N* j         *ELSE
             lesize,all,THICKNESS(Num)/div(Num), , ,1/dscale(Num),1, , ,1,
         *ENDIF
      *ENDDO
0 t& _" U8 H  R2 y( h) E/ H   *ENDDO   
# {, D+ G8 i! r& }( I) @  
   asel,all
4 ]' i. J7 f6 @( D5 ?5 [  
   esize,Size,0,
  
   VSEL,S,LOC,Y,-TSum1,-TSum2
  
6 M0 w, Z' r" ~# Z/ K. a; B   TYPE,1
* w8 F9 N7 ]( W, y5 w9 [   MAT,Num
   mshkey,1  
   mshape,0,3d  
% m, J- m. g; F4 A) d1 l8 p   vmesh,ALL             !Mesh the layer Num
. q0 z) w4 b9 \& ?" b  
( y5 W1 h4 d1 K) s' P( N; o   allsel,all
  
   TSum2=TSum2+Thickness(Num)
  
*enddo
  
+ H$ d- s- _, m0 ]4 y7 X/REP,FAST  
+ f! H1 L; q$ S. G/psf,hflu,,2,0  
9 x1 J% y! P$ l8 R  P/PNUM,LINE,1  
finish  
/solu  
  
0 e& U% P/ K5 x* E+ w% P! Pantype,trans             !Transient simulation
. [' `5 A1 j0 utunif,T0                 !Ambient temperature
" C; y; v8 ?) W% eTREF,T0
timint,on !time integration = on  
4 t! E7 o9 z# j% i3 a+ ]* Okbc,1  
7 x# x/ l1 ^7 oOUTRES,ALL,ALL,          !Output all the result
8 u7 j6 ]* ?8 M  
/view, 1 ,1,1,1  
/ang, 1  
9 J+ ^. }6 k4 e; A3 w/rep,fast  
1 }) G, C+ I* E+ J+ R# g2 |1 l  
% s( K, o5 G. }% s- |/ j; T- C, l*do,Num,1,Layer,1
6 f  H6 P; _9 D1 r   *USE,HG.mac,'Layer%Num%',Power,r,V,A0(Num),A1(Num),A2(Num),A3(Num),A4(Num)     !Load heat generation for each layer
*enddo
  
!Start simulation
*do,tm,1e-15,End_Time,End_Time/Step  
2 x4 p/ O/ s2 ], P  G! |! N) A# c  
   time,tm  
   autots,0  
   deltim, , , ,1  
4 P3 ^; H* T) `8 }9 L   kbc,0  
- N: H, O* V5 l% ~2 U( ]  
   allsel,all  
* o- Y) J( g% Q0 Z; M   TSum1=0
   TSum2=0
   Num=0
  
% T% @8 u/ ?5 r7 x& L( I  J9 \   *if,tm,LE,Pulse_Duration+1e-12,then  
        
      !Apply the heat generation in each layer  
1 ?% b! u0 _& M3 O; X" V4 _- e+ S# @' ~      Num=1
      TSum1=TSum1+THICKNESS(Num)
4 S. n% G2 ]: w  r      nsel,s,LOC,y,-TSum1,-TSum2   
5 b+ R) u' l! X4 S      BF,all,HGEN,%Layer1%
      TSum2=TSum2+Thickness(Num)
      
1 C8 I. x; c! G0 f      Num=2
/ N" ]0 a7 |0 r- a! o; M5 L      TSum1=TSum1+THICKNESS(Num)
+ L4 ?/ r4 U9 n0 ], ]7 N/ O; z      nsel,s,LOC,y,-TSum1,-TSum2   
      BF,all,HGEN,%Layer2%
      TSum2=TSum2+Thickness(Num)
& e7 D5 @6 {5 a8 e   
+ m) P- r9 ~+ T$ r. B3 _      Num=3   
      TSum1=TSum1+THICKNESS(Num)
5 h2 o. X* z6 Z4 ~  r      nsel,s,LOC,y,-TSum1,-TSum2   
      BF,all,HGEN,%Layer3%
8 t! q' N5 g+ z! R0 E4 U      TSum2=TSum2+Thickness(Num)
   
      Num=4
# E$ q/ X  [1 W! h2 }      TSum1=TSum1+THICKNESS(Num)
/ @, G: u+ f: M, o  q6 {) F      nsel,s,LOC,y,-TSum1,-TSum2   
      BF,all,HGEN,%Layer4%
4 e. o' m- _7 y4 D9 a$ G& N      TSum2=TSum2+Thickness(Num)
: v% w- Y% [3 G  }7 m+ ^0 K   
      Num=5
      TSum1=TSum1+THICKNESS(Num)
" |+ H1 W$ Q2 a/ g      nsel,s,LOC,y,-TSum1,-TSum2   
5 d( L8 e; F2 r2 N) V5 i8 e2 N" q0 {9 V      BF,all,HGEN,%Layer5%
! V( Q6 V( Y' g  ]& J  T- `5 i      TSum2=TSum2+Thickness(Num)
  
   *else
7 T/ S' Z+ h% J( n  X$ a3 _      allsel,all  
, v1 p: b9 U, a4 ]) z3 X/ ]8 }      BFDELE, ALL,HGEN
   *endif  
  
5 L& F+ U% n0 O9 a* P, l   allsel,all
   asel,s,LOC,y,-1e-9,0
7 j; o/ H6 N( B, h/ @$ m!   DA,ALL,TEMP,T0  
   allsel,all
  
1 W; Q9 H% ]- B9 {4 ]% ?$ |) `; U2 {   /status,solu  
   solve  
+ x; B, A( N" D& }*enddo  
# k/ o+ j+ k% _0 x! m9 ]  
finish  
( V! N- ^8 R& Y  y  `, Z/POST1
8 }2 t+ x: M0 c* a5 \4 A% j% Z+ c0 e  
!*ULIB, Post,mac
8 t# ~$ A6 ]8 [: v% [' @- i# l!*USE,TTV    !Temperature decay at each layer thickness
8 W  R4 ]' w: |  b& a- R" |!*USE,TRV,1  !Normal Temperature profile at each layer at given time  
, T- f: U* p/ C* |* h9 b!*USE,TRL,1  !Plane Temperature profile at each layer at given time  
; l. j4 m) ?4 R7 Z. {; a  
( w; _1 p4 A9 w1 m' y8 ?!Laser power absorption in each layer
5 Q/ `! s8 z. }/ r4 v2 H, V!Power/r^2/{PI}*exp(-({X}^2+({Z}-V*{TIME})^2)/r^2)*( A0+A1*(-{Y}/1e-9)+A2*(-{Y}/1e-9)^2+A3*(-{Y}/1e-9)^3+ A4*(-{Y}/1e-9)^4)      
0 k  m3 Z" D, t2 j2 t# {+ d2 q  
//UG.mac
  V% m! x; m8 \# L2 h  
) P* ~2 \- q1 H* U!FileName=ARG1
!Power=ARG2
!r=ARG3
!V=ARG4
!A0=ARG5
8 g7 ]! q9 G  t2 [* @, H$ U* ~9 T!A1=ARG6
!A2=ARG7
!A3=ARG8
0 h. H2 R: J" L. y0 N  L6 S+ P!A4=ARG9
* B) @0 Q% Y! N; K& z. a  
& F: O9 c. G: J; K* E*DEL,_FNCNAME   
" C  {% P5 _/ V6 i5 n*DEL,_FNCMTID   
*DEL,_FNC_C1
*DEL,_FNC_C2
8 I/ ^# ~5 g  m8 Q  `*DEL,_FNC_C3
*DEL,_FNC_C4
*DEL,_FNC_C5
*DEL,_FNC_C6
*DEL,_FNC_C7
*DEL,_FNC_C8
*SET,_FNCNAME,ARG1   
0 G6 N" Y+ S; Y*DIM,_FNC_C1,,1  
0 _9 k4 s: I7 ~& |9 L4 Y*DIM,_FNC_C2,,1  
7 f/ |6 m! L, m$ e*DIM,_FNC_C3,,1  
2 `# |9 @4 A4 q/ W+ A, a, i*DIM,_FNC_C4,,1  
*DIM,_FNC_C5,,1  
. Z$ ~% p. l% G. L*DIM,_FNC_C6,,1  
9 {8 i6 A+ n/ `*DIM,_FNC_C7,,1  
*DIM,_FNC_C8,,1  
*SET,_FNC_C1(1),ARG2   
*SET,_FNC_C2(1),ARG3
5 Y4 g2 w5 i* ~& k) d8 Q*SET,_FNC_C3(1),ARG4
  x1 p5 S3 ]* g# I, F4 D*SET,_FNC_C4(1),ARG5
8 d; Y' }  \  k: E  K2 F0 W*SET,_FNC_C5(1),ARG6
' x+ W7 D# l" `: @9 V: c; v*SET,_FNC_C6(1),ARG7
*SET,_FNC_C7(1),ARG8
*SET,_FNC_C8(1),ARG9
! /INPUT,HG.func
*DIM,%_FNCNAME%,TABLE,10,61,1   
# {/ e% i, w4 k: _9 X2 Q!   
! Begin of equation: Power/r^2/{PI}*exp(-({X}^2+({Z}-V*{TIME})^2)/r^2)*(A0+A1*   
+ E' r( q5 ?, {9 C8 C  n! (-{Y}/1e-9)+A2*(-{Y}/1e-9)^2+A3*(-{Y}/1e-9)^3+ A4*(-{Y}/1e-9)^4)   
* _0 s: m% n- b! Z; y) B" x( u*SET,%_FNCNAME%(0,0,1), 0.0, -999   
  V* U) \" q. c*SET,%_FNCNAME%(2,0,1), 0.0  
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
6 U9 ~! m% c# g6 y*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
*SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)%
*SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)%
*SET,%_FNCNAME%(8,0,1), %_FNC_C6(1)%
8 M2 E- y5 P: n& e: d5 P" ?5 V5 x*SET,%_FNCNAME%(9,0,1), %_FNC_C7(1)%
, g/ l3 ~/ [' @*SET,%_FNCNAME%(10,0,1), %_FNC_C8(1)%   
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 2, 0, 0, 18  
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 18, 17, -1   
*SET,%_FNCNAME%(0,3,1),   0, -1, 0, 1, 17, 4, -2
8 m3 Z% \, h* m" p9 ~% ]. ]*SET,%_FNCNAME%(0,4,1), 0.0, -2, 0, 3.14159265358979310, 0, 0, -1   
*SET,%_FNCNAME%(0,5,1), 0.0, -3, 0, 1, -1, 4, -2
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 0, 0, 0, 0   
*SET,%_FNCNAME%(0,7,1), 0.0, -2, 0, 1, 0, 0, -1  
; Z7 {; X/ v; `9 R4 C4 Y% z*SET,%_FNCNAME%(0,8,1), 0.0, -4, 0, 1, -1, 2, -2
% P& N6 l' ]. F0 q*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 2, 0, 0, 2   
# U) |; ^; X$ e! `6 F2 _: p5 G*SET,%_FNCNAME%(0,10,1), 0.0, -2, 0, 1, 2, 17, -1   
*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 1, 19, 3, 1
0 v$ O/ f2 g3 ]% E*SET,%_FNCNAME%(0,12,1), 0.0, -5, 0, 1, 4, 2, -1
*SET,%_FNCNAME%(0,13,1), 0.0, -1, 0, 2, 0, 0, -5
*SET,%_FNCNAME%(0,14,1), 0.0, -6, 0, 1, -5, 17, -1   
*SET,%_FNCNAME%(0,15,1), 0.0, -1, 0, 1, -2, 1, -6   
*SET,%_FNCNAME%(0,16,1), 0.0, -2, 0, 1, -4, 3, -1   
! x! B! g( w( ~6 V2 l; N, E*SET,%_FNCNAME%(0,17,1), 0.0, -1, 0, 2, 0, 0, 18
*SET,%_FNCNAME%(0,18,1), 0.0, -4, 0, 1, 18, 17, -1   
*SET,%_FNCNAME%(0,19,1), 0.0, -1, 0, 1, -2, 4, -4   
*SET,%_FNCNAME%(0,20,1), 0.0, -1, 7, 1, -1, 0, 0
: A8 Y0 a" @9 _* W- W8 ]* H*SET,%_FNCNAME%(0,21,1), 0.0, -2, 0, 1, -3, 3, -1   
*SET,%_FNCNAME%(0,22,1), 0.0, -1, 0, 0, 0, 0, 0  
3 E# z9 R; `$ U" O  L% ]9 k*SET,%_FNCNAME%(0,23,1), 0.0, -3, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,24,1), 0.0, -4, 0, 1, -1, 2, -3   
*SET,%_FNCNAME%(0,25,1), 0.0, -1, 0, 1, -4, 3, 3
7 V8 w) `( |* d: X- }0 `8 g2 X+ y*SET,%_FNCNAME%(0,26,1), 0.0, -3, 0, 1e-9, 0, 0, -1  
$ [" l* f$ Z& M- B9 M6 c*SET,%_FNCNAME%(0,27,1), 0.0, -4, 0, 1, -1, 4, -3   
* _* k/ X4 R6 C2 P8 k5 g* Y% {*SET,%_FNCNAME%(0,28,1), 0.0, -1, 0, 1, 21, 3, -4   
*SET,%_FNCNAME%(0,29,1), 0.0, -3, 0, 1, 20, 1, -1   
: }4 x6 }" h0 y; Y+ n$ {*SET,%_FNCNAME%(0,30,1), 0.0, -1, 0, 0, 0, 0, 0  
*SET,%_FNCNAME%(0,31,1), 0.0, -4, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,32,1), 0.0, -5, 0, 1, -1, 2, -4   
*SET,%_FNCNAME%(0,33,1), 0.0, -1, 0, 1, -5, 3, 3
1 c( P) b  r& J$ d0 m7 B# n: W*SET,%_FNCNAME%(0,34,1), 0.0, -4, 0, 1e-9, 0, 0, -1  
*SET,%_FNCNAME%(0,35,1), 0.0, -5, 0, 1, -1, 4, -4   
) e8 d" T" I' t/ S6 Q& R5 O9 K* r( |*SET,%_FNCNAME%(0,36,1), 0.0, -1, 0, 2, 0, 0, -5
0 D. r$ g! g& k$ L2 G1 R7 m; d  G*SET,%_FNCNAME%(0,37,1), 0.0, -4, 0, 1, -5, 17, -1   
; T( h3 C9 J' k*SET,%_FNCNAME%(0,38,1), 0.0, -1, 0, 1, 22, 3, -4   
*SET,%_FNCNAME%(0,39,1), 0.0, -4, 0, 1, -3, 1, -1   
  s' F9 i: \9 E2 P% b. q*SET,%_FNCNAME%(0,40,1), 0.0, -1, 0, 0, 0, 0, 0  
4 H& Z: v2 U) E8 G9 d% A4 y*SET,%_FNCNAME%(0,41,1), 0.0, -3, 0, 1, 0, 0, -1
: l- }. y3 D. J0 [& E0 q! O7 S*SET,%_FNCNAME%(0,42,1), 0.0, -5, 0, 1, -1, 2, -3   
*SET,%_FNCNAME%(0,43,1), 0.0, -1, 0, 1, -5, 3, 3
4 C( z( k) s, @2 z8 L*SET,%_FNCNAME%(0,44,1), 0.0, -3, 0, 1e-9, 0, 0, -1  
7 x- H8 u' M1 |& Q& m5 a3 Y% H( T1 k*SET,%_FNCNAME%(0,45,1), 0.0, -5, 0, 1, -1, 4, -3   
*SET,%_FNCNAME%(0,46,1), 0.0, -1, 0, 3, 0, 0, -5
0 h/ y( @# H1 f# J" y5 w*SET,%_FNCNAME%(0,47,1), 0.0, -3, 0, 1, -5, 17, -1   
*SET,%_FNCNAME%(0,48,1), 0.0, -1, 0, 1, 23, 3, -3   
1 s$ l+ c6 @$ r7 r/ @*SET,%_FNCNAME%(0,49,1), 0.0, -3, 0, 1, -4, 1, -1   
2 }! _( y# ~6 H0 K. s& r- \4 k5 V# J*SET,%_FNCNAME%(0,50,1), 0.0, -1, 0, 0, 0, 0, 0  
*SET,%_FNCNAME%(0,51,1), 0.0, -4, 0, 1, 0, 0, -1
9 I9 A$ |6 L/ @" i5 M9 E" W* ^( q*SET,%_FNCNAME%(0,52,1), 0.0, -5, 0, 1, -1, 2, -4   
% x& E3 U8 e+ H" j*SET,%_FNCNAME%(0,53,1), 0.0, -1, 0, 1, -5, 3, 3
*SET,%_FNCNAME%(0,54,1), 0.0, -4, 0, 1e-9, 0, 0, -1  
*SET,%_FNCNAME%(0,55,1), 0.0, -5, 0, 1, -1, 4, -4   
- r4 u8 w7 z2 I  A2 ~4 u0 j( ^*SET,%_FNCNAME%(0,56,1), 0.0, -1, 0, 4, 0, 0, -5
*SET,%_FNCNAME%(0,57,1), 0.0, -4, 0, 1, -5, 17, -1   
*SET,%_FNCNAME%(0,58,1), 0.0, -1, 0, 1, 24, 3, -4   
*SET,%_FNCNAME%(0,59,1), 0.0, -4, 0, 1, -3, 1, -1   
*SET,%_FNCNAME%(0,60,1), 0.0, -1, 0, 1, -2, 3, -4   
*SET,%_FNCNAME%(0,61,1), 0.0, 99, 0, 1, -1, 0, 0
! End of equation: Power/r^2/{PI}*exp(-({X}^2+({Z}-V*{TIME})^2)/r^2)*(A0+A1*(-   
! {Y}/1e-9)+A2*(-{Y}/1e-9)^2+A3*(-{Y}/1e-9)^3+ A4*(-{Y}/1e-9)^4)
1 f! v" Z* Z0 G( N1 w0 B4 b7 `!-->