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

xy112

高斯激光热源下的n层材料瞬态热传导方程
高斯, 激光, 热源, 方程, 热传导高斯, 激光, 热源, 方程, 热传导
- H" l7 A0 q2 m, E出于需要， 写了个 n层材料的瞬态热传导方程， 由表面激光加热。每层的热源由其光学性质N= n + ik决定. 通过解Maxwell方程获得激光在材料中的强度分布， 由多项式 a0+a1*y+a2*y^2+a3*y^3+a4*y^4拟和获得。
5 M( e, ]7 p: V) FP=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).
  
/ I% N$ i! c- [' x高斯分布体热源在程序中由HG.mac加载。不同层的a0,a1,a2,a3,a4决定了其内部的热源分布。
  
) H0 _% n9 `' F' g% ?本程序可已根据问题选择所需的材料层数，所有参数都可以改变。  
4 O5 Y0 x6 n# j7 f, G9 Q5 P8 Q  
以下是命令流，如果有不合适的地方，欢迎大家提出、讨论！
. L2 C6 S& T0 o* U! I! l3 C/ }  
FINISH  ! Make sure we are at BEGIN level   
0 f$ n* Q: |" o; _/Clear
+ B8 F  z* C9 `7 T/FILNAME,Media,0  
3 k2 y, [" c1 z  
/PREP7
6 @; N, O9 Y9 U1 N% {/PNUM,VOLU,1
, R9 S9 |4 E4 Z  
! In simulation, the X-Z is the media plane, while Y is the depth direction
  
, |5 ]! A' Q# {, Y- p  [* WLayer=5                       !Total layer number
+ c- C. Z- m% E- ^. ynm= 1e-9                    !length unit
Power=0.005               ! Laser incident power   W
r=5e-7                        ! Laser spot radius      m
2 Y3 R1 e8 A( uV=0                           ! Laser moving veolcity, along z direction,  m/s
4 E3 ?- E- v, q  ^. a  
3 O. ~7 i; s6 _4 `  epw=Power*1000                 !For display Title
rh=r/1e-6                           !For display Title
( d, l5 g2 [' h' a5 `5 k# c  
# W- P) f0 w- B* r0 F6 U, y7 p/TITLE,%Layer% layer media, Laser Power=%pw%mW, Heating radius=%rh%um, Velocity=%V%m/s  
/ t) O: K) B$ R# M* y/ I' O! n  
*DIM,Thickness,ARRAY,Layer    ! Array input for each layer thickness
*DIM,k_L,ARRAY,Layer          ! Array input for each layer lateral thermal conductivity  
*DIM,k_N,ARRAY,Layer          ! Array input for each layer normal thermal conductivity
6 q2 }5 I8 z* ?( F& A; |*DIM,Cp,ARRAY,Layer           ! Array input for each layer heat capacity
*DIM,Den,ARRAY,Layer          ! Array input for each layer density
*DIM,div,ARRAY,Layer          ! Array input for each layer thickness division
*DIM,dscale,ARRAY,Layer       ! Array input for each layer thickness division scale
/ M/ H9 X5 n) y: @& S*DIM,A0,ARRAY,Layer           ! Array input for each layer a0 coefficient for heat generation curve
: n5 d# |. e+ B/ @*DIM,A1,ARRAY,Layer           ! Array input for each layer a1 coefficient for heat generation curve
  h* Z+ G, l7 T*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
*DIM,A4,ARRAY,Layer           ! Array input for each layer a4 coefficient for heat generation curve
  
scale=10000                   !for vovlap operation
Pulse_Duration=1e-8           !Laser pulse duration
End_Time=2e-8                 !Final time for simulation
/ b$ F( D$ e: I" ~# GStep=20                       !Time interval  
1 \4 E) Q  C* ]' |6 B0 ?2 _4 j5 H1 H3 TT0=50                         !Environment temperature   
- `' }& Y' F# A/ w% h  
7 C: P0 y' B8 C; E+ U; `7 eLength=4000*nm                !Length of the block in z direction
8 l1 ?1 y; a! w: P: L% wWidth=2000*nm                 !Width of the block in x direction
  
!Layer properties, from top to bottom  
  
!First layer thermal property
' i5 c! y! R* ^6 tNum=1                        
k_L(Num)=1.4                  
k_N(Num)=1.4
0 G+ i! M1 V  ~0 P8 C+ ]6 i: sDen(Num)=4000
# s7 l$ c" C# T! t8 T; m, \- HCp(Num)=500
2 s" \6 a2 B9 S6 u: Q9 u6 ^Thickness(Num)=50*nm
div(Num)=5
, F2 f( r6 b) |- q1 k/ ~; ~dscale(Num)=1                 !last/first in -Y direction
A0(Num)=0
1 ^) s- T. E8 `- v; J& p. `A1(Num)=0
( t1 F# A2 c1 X8 `" v+ RA2(Num)=0
3 A! p6 J, t. oA3(Num)=0
A4(Num)=0
  
; e! J3 @5 l% [$ H!Second layer thermal property
, Z4 L$ z5 J/ n% xNum=2
( t" z+ M- r+ Wk_L(Num)=16
k_N(Num)=16
7 \- v% z3 n3 S2 }  }6 lDen(Num)=3000
, T2 n- s* I' h2 R* l4 |Cp(Num)=900
) X8 y3 Q; a! t" ZThickness(Num)=20*nm
! h" e# e; |- {0 b0 x. }7 qdiv(Num)=2
8 U0 V% a* M1 ^) i. W3 S/ ^# mdscale(Num)=1
A0(Num)=807.929E+6
1 I# T5 H  C; ~A1(Num)=-43.896E+6
# h3 [- q) w% y/ IA2(Num)=934.529E+3
3 ~0 t6 G* ?. q" sA3(Num)=-8.864E+3
7 A4 X# i0 N6 p; GA4(Num)=31.703
& @# Z3 t  t; h+ J' S/ K7 c  
!Third layer thermal property
# `4 g; R2 Q5 H* ^, G) ]Num=3
k_L(Num)=1.5
k_N(Num)=1.5
Den(Num)=4000
Cp(Num)=500
: h0 ~7 t) v$ j& q/ [3 jThickness(Num)=100*nm
div(Num)=5
dscale(Num)=1
0 F: G0 ]; X( t2 K. KA0(Num)=0
! u5 K8 Z% R0 u6 \A1(Num)=0
A2(Num)=0
3 l9 X( {$ y# J9 j2 \5 q$ Y+ ?A3(Num)=0
2 w3 ~% O8 O5 y; l$ C1 t+ wA4(Num)=0
( c, P. f& Z$ c$ U! K' u! _8 u  
# z! p' G8 u  b- M9 Q' ?* T, w!Fourth layer thermal property
Num=4
k_L(Num)=100
0 @  }. M/ F: ^& w1 A; U: jk_N(Num)=100
Den(Num)=3900
% p5 X/ p6 j0 M1 G: t6 BCp(Num)=700
; C3 T* c: g& K6 `# qThickness(Num)=50*nm
div(Num)=2
7 l. n# U% l9 L( Jdscale(Num)=2
" i& j6 D# e" ]# V+ I6 rA0(Num)=1.627E+9
A1(Num)=-32.146E+6
7 u3 @, J) k8 W( S: C5 OA2(Num)=238.195E+3
A3(Num)=-784.173
3 y& T( o3 L5 r; P" g& }A4(Num)=967.716E-3
0 g  ]6 i. S' I* G  
9 \9 O. {5 M4 e9 `# z) v!Fifth layer thermal property
Num=5
k_L(Num)=0.2
k_N(Num)=0.2
' q( B: [8 m: q$ s% dDen(Num)=3900
6 e1 L# f2 I% PCp(Num)=450
) B) V5 c, n  K: i+ \& ZThickness(Num)=100*nm
# X& X7 @( ~& qdiv(Num)=5
dscale(Num)=2
5 M: Q0 B7 H9 X( b& P9 ZA0(Num)=0
A1(Num)=0
A2(Num)=0
/ e2 {5 n1 d* r% l4 N7 Z/ N# GA3(Num)=0
A4(Num)=0
  
! Define the thermal element
' k9 @: F* b/ C) M: g5 IET,1,SOLID70  
5 w" ?4 Y; O1 t/ i- a* oKEYOPT,1,2,0  
+ }) a( K) q" {3 H$ N) ^9 G4 \KEYOPT,1,4,0  
6 n8 ?) e) B. IKEYOPT,1,7,0  
KEYOPT,1,8,0  
2 E0 [5 L, k# i. n0 `9 w' f  
!For post-precessing
! N$ P) C6 j. C: O+ }/ [nPoint=200
3 S6 \) J2 {6 ~* |; N3 cdisplay_r=1e-6
1 o6 e& J- {4 O. R; |  
TSum1=0
% P( @$ K# }: A" M7 u6 m4 B% `0 tTSum2=0
*do,Num,1,Layer,1
3 r% ?8 \5 ]6 B  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
9 \' K+ u6 }- _1 h6 k& a  UIMP,Num,C, , ,Cp(Num),                                            !Define heat capacity of material Num   
" h. ^1 S8 @5 M+ d  TSum1=TSum1+Thickness(Num)  
/ c( ~% l, |7 }/ w5 p* N+ e0 g' V+ k  BLOCK,-Width/2,Width/2,-TSum1,-TSum2,-Length/2,Length/2,           !Draw the layer of Num  
  TSum2=TSum2+Thickness(Num)
/ E. c& P7 |( a*Enddo
  
d_total=TSum1                 !Total thickness of the multilayer                                                  
& k7 I# m2 B8 b, f  
VLSCALE,ALL, , , scale, scale, scale,  ,1,1,
  
5 Q2 e* o: M4 M! B*IF,Layer,GT,1,THEN
% X% @) C8 b. Y: W2 u# o   VGLUE,ALL
*ENDIF
  
; B! Z4 h3 g2 FVLSCALE,ALL, , , 1/scale, 1/scale, 1/scale, ,1,1,
. O5 D6 V" P- l- f  
- f. a' u! L( dTSum1=0
TSum2=0
* @2 E1 f/ P4 i$ A- L, d  
% r0 @% Z1 w6 w. s! dSize=100*nm                             !Area mesh size, in X-Z plane
, n) w  H, {2 e; L( Saesize,all,Size,  
! M3 ?% [: v7 |8 ^3 R' D  
6 i& s  j$ K& X3 m6 ^5 p& {/ _% G*do,Num,1,Layer,1
) Q4 _5 C3 J! @   TSum1=TSum1+THICKNESS(Num)
. i% B% `, e; r- A) `   lsel,all
  
3 x4 q' X" L' l3 t8 U! |" Y, K   !Determine the mesh scale and divison for each layer at Y direction
   *do,i,-1,1,2
- T& p: F  M. f0 L0 v8 c: G4 V- L      *do, j,-1,1,2
! N0 j0 b4 L& W5 B) d         lsel,s,loc,y,-TSum1,-TSum2
4 Z9 u& e! O" ?1 M7 j         lsel,r,loc,x,i*Width/2,i*Width/2
7 i3 M8 f# ~9 C- C         lsel,r,loc,z,j*Length/2,j*Length/2
         *GET,Line_Num,LINE, 0, NUM,MIN,
         *GET,KP1,LINE,Line_Num,KP,1,
5 \* b# c% I4 K1 K  r, w5 `         *GET,KP2,LINE,Line_Num,KP,2,
         *GET,KP1_Y,KP,KP1, LOC, Y,  
/ L1 m% i% \$ H, A* E         *GET,KP2_Y,KP,KP2, LOC, Y,
         *IF,KP1_Y,GT,KP2_Y,THEN
             lesize,all,THICKNESS(Num)/div(Num), , ,dscale(Num),1, , ,1,
         *ELSE
3 o" }- t6 o  v/ w             lesize,all,THICKNESS(Num)/div(Num), , ,1/dscale(Num),1, , ,1,
         *ENDIF
      *ENDDO
   *ENDDO   
0 u; z+ f7 E# W: j7 H/ h* k7 M  
   asel,all
  
   esize,Size,0,
+ d5 r* j' g) y0 J& g9 J2 g  
   VSEL,S,LOC,Y,-TSum1,-TSum2
: B2 F, F6 h/ O* X3 s  
9 G- M1 T3 j- i; N5 C   TYPE,1
   MAT,Num
   mshkey,1  
  ~* d" C: O9 S, h  l+ y5 `   mshape,0,3d  
" f, _1 h7 D7 R6 O   vmesh,ALL             !Mesh the layer Num
  
   allsel,all
3 \& T/ U- |' s6 q* S  
' A: a2 y' w9 R   TSum2=TSum2+Thickness(Num)
  
5 ]6 I' y- u3 O) y; B  W*enddo
0 [  |% U; ]: p8 m  
6 \+ J. X! Y3 Q) S1 R, R* N/REP,FAST  
. F; F9 W1 c* }) w* l4 C/psf,hflu,,2,0  
/PNUM,LINE,1  
% o6 ?1 l! ?" i7 C$ u; K) ifinish  
/solu  
, u  h) p9 G2 H  
antype,trans             !Transient simulation
" o0 c" r6 x9 D( Jtunif,T0                 !Ambient temperature
TREF,T0
9 W8 Y: ]; ^5 ytimint,on !time integration = on  
kbc,1  
# X; D* y7 R/ _8 @* J% O3 G+ UOUTRES,ALL,ALL,          !Output all the result
  
, h: C' c# k8 v0 L3 x" L2 P/view, 1 ,1,1,1  
/ang, 1  
  M# |# z& Q: P3 V/ N8 \5 T/rep,fast  
5 J0 S% @# g3 j8 u% A  
*do,Num,1,Layer,1
8 b/ R7 l! q0 w! x/ W: H   *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
& f6 g- M( n. o3 z  
/ Z$ v1 [/ b! N5 D2 p; e' X% r!Start simulation
3 @# O4 y; L5 d# u7 I% S: R*do,tm,1e-15,End_Time,End_Time/Step  
  
) N1 a; F2 @# L! e) L* c+ a   time,tm  
   autots,0  
   deltim, , , ,1  
   kbc,0  
  
1 e0 N& l3 [8 ^, I# d7 M  j+ Z8 X: x   allsel,all  
! C' J! g( K1 [" o+ D" _   TSum1=0
   TSum2=0
2 T( g5 I( a+ s1 S6 ^  l* C   Num=0
  
   *if,tm,LE,Pulse_Duration+1e-12,then  
$ x# y% P1 T0 X& Y& N! h( I% N        
  ~1 E% g9 t. T- C      !Apply the heat generation in each layer  
& B9 `6 Z$ Z  D# S5 [! q& D% d6 v      Num=1
      TSum1=TSum1+THICKNESS(Num)
  `, [" k) K4 i+ ?) l      nsel,s,LOC,y,-TSum1,-TSum2   
- L# x8 Y& t2 }8 [8 c8 h      BF,all,HGEN,%Layer1%
) s1 B, q+ L7 z: f6 _+ m+ X      TSum2=TSum2+Thickness(Num)
      
2 z7 ?+ ~6 J. m+ k      Num=2
. w. K) p3 T! _0 N      TSum1=TSum1+THICKNESS(Num)
  J' j' U$ u2 j- |      nsel,s,LOC,y,-TSum1,-TSum2   
      BF,all,HGEN,%Layer2%
' P) U: L+ j% j6 j; J2 w      TSum2=TSum2+Thickness(Num)
   
: ~( i2 W" K& R2 m6 k3 T      Num=3   
      TSum1=TSum1+THICKNESS(Num)
2 }( ?$ t# O. d, C      nsel,s,LOC,y,-TSum1,-TSum2   
! ]/ f. c  J4 c1 X2 _8 K      BF,all,HGEN,%Layer3%
      TSum2=TSum2+Thickness(Num)
   
      Num=4
9 d& w/ C/ B6 H( C' P6 y0 H! e# V      TSum1=TSum1+THICKNESS(Num)
6 Z4 h; _- w! R4 k4 r/ O5 W      nsel,s,LOC,y,-TSum1,-TSum2   
3 Q5 C# i# e, |+ x: x1 ]      BF,all,HGEN,%Layer4%
      TSum2=TSum2+Thickness(Num)
! M" ~. W' E3 L* g* H- ~$ ^* }   
      Num=5
      TSum1=TSum1+THICKNESS(Num)
; L: n& W6 `3 X4 \7 o9 S8 ^      nsel,s,LOC,y,-TSum1,-TSum2   
      BF,all,HGEN,%Layer5%
3 |. n  ~& U% H4 @      TSum2=TSum2+Thickness(Num)
, e* Y! Q/ c0 D( O/ b6 o# Q  
! r0 o& V& O8 {* k/ F   *else
. [4 m6 c; n' }9 a3 a# Z- A; ^      allsel,all  
- a) n0 z' a8 d& c% C9 m) u      BFDELE, ALL,HGEN
   *endif  
  
" r3 D  Z0 @. @   allsel,all
9 i! x& ~: A; P( }$ Q  k0 Y   asel,s,LOC,y,-1e-9,0
!   DA,ALL,TEMP,T0  
   allsel,all
  y9 y/ H( ?3 Z' w$ [8 Q  
2 b  B  O, E+ v- k   /status,solu  
7 [, s* b7 T/ _$ D+ \; D   solve  
*enddo  
$ G  ^# z, f, W, Z% z  
finish  
/POST1
  
# V6 E/ l' v. T4 Q* x!*ULIB, Post,mac
!*USE,TTV    !Temperature decay at each layer thickness
!*USE,TRV,1  !Normal Temperature profile at each layer at given time  
$ C8 ?9 l$ L' u) J!*USE,TRL,1  !Plane Temperature profile at each layer at given time  
  
!Laser power absorption in each layer
!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)      
  
& ^  k0 B6 L6 z  s) R5 ?" Q8 e//UG.mac
8 A8 B+ @- H9 k0 e7 N  
3 _$ X  Y) r% S' X2 D2 ]% a" r9 Z( U!FileName=ARG1
!Power=ARG2
!r=ARG3
!V=ARG4
2 e: F$ b# W3 _! z, I1 K: T!A0=ARG5
  Y- W" M% {2 D) u2 A!A1=ARG6
!A2=ARG7
!A3=ARG8
' Y# D' Z9 J/ u4 L  T. j- _!A4=ARG9
  
9 ?( S* F( Y8 Z' ?; p4 I2 c6 Y4 v  K*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNC_C1
*DEL,_FNC_C2
7 \3 c1 e4 v8 H9 A*DEL,_FNC_C3
*DEL,_FNC_C4
*DEL,_FNC_C5
; v% R4 X( l9 P( ~*DEL,_FNC_C6
6 h/ Y- e. q7 M4 v; f*DEL,_FNC_C7
+ N: j7 Z# i8 B1 E$ C+ B*DEL,_FNC_C8
*SET,_FNCNAME,ARG1   
/ ^+ f- I; @: W1 {* C' y) e; V*DIM,_FNC_C1,,1  
: Y! z1 q  A  D*DIM,_FNC_C2,,1  
*DIM,_FNC_C3,,1  
* t& s- W0 X; \) P4 @1 W*DIM,_FNC_C4,,1  
*DIM,_FNC_C5,,1  
*DIM,_FNC_C6,,1  
6 _+ D4 Y8 z5 p7 X" x6 y9 g*DIM,_FNC_C7,,1  
*DIM,_FNC_C8,,1  
$ t1 J' h, C) a/ T! I0 V& a# R*SET,_FNC_C1(1),ARG2   
*SET,_FNC_C2(1),ARG3
*SET,_FNC_C3(1),ARG4
( D8 E3 |/ q" ]. U* L*SET,_FNC_C4(1),ARG5
/ u3 g  _2 T- _% S0 m! O/ V*SET,_FNC_C5(1),ARG6
( L$ |( l! Y0 d" w; |! _( s*SET,_FNC_C6(1),ARG7
*SET,_FNC_C7(1),ARG8
) L, c% V  }( P*SET,_FNC_C8(1),ARG9
* B$ y" _& l1 ]2 w# v/ P! /INPUT,HG.func
*DIM,%_FNCNAME%,TABLE,10,61,1   
6 J0 s) e$ O: {% d3 _7 M# m!   
! Begin of equation: Power/r^2/{PI}*exp(-({X}^2+({Z}-V*{TIME})^2)/r^2)*(A0+A1*   
4 |9 J/ r% {( X& X! (-{Y}/1e-9)+A2*(-{Y}/1e-9)^2+A3*(-{Y}/1e-9)^3+ A4*(-{Y}/1e-9)^4)   
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
6 T  C2 M7 K" b*SET,%_FNCNAME%(2,0,1), 0.0  
9 D& H. l7 q( E*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
! o* [7 `' k: B0 k% D% o/ q*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
* h2 T# Z9 i  _) f0 M; @! l*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
* @+ N( }2 d. @; x. m*SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)%
  w# j- t7 M  y$ s& h" O8 O*SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)%
/ v) [% J- @: Q' d: W& O*SET,%_FNCNAME%(8,0,1), %_FNC_C6(1)%
. l7 F3 p) b9 [5 F- Q*SET,%_FNCNAME%(9,0,1), %_FNC_C7(1)%
*SET,%_FNCNAME%(10,0,1), %_FNC_C8(1)%   
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 2, 0, 0, 18  
' p1 W% [& c$ Q" @, A*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 18, 17, -1   
( g) r( ?) ?  W) d) c*SET,%_FNCNAME%(0,3,1),   0, -1, 0, 1, 17, 4, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -2, 0, 3.14159265358979310, 0, 0, -1   
2 y+ ?# R. c& P8 k( J*SET,%_FNCNAME%(0,5,1), 0.0, -3, 0, 1, -1, 4, -2
3 q' L, @+ v' [$ S9 ^*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 0, 0, 0, 0   
6 s$ J6 ]7 V) K0 i- ?*SET,%_FNCNAME%(0,7,1), 0.0, -2, 0, 1, 0, 0, -1  
*SET,%_FNCNAME%(0,8,1), 0.0, -4, 0, 1, -1, 2, -2
*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 2, 0, 0, 2   
$ v( f  @2 [2 f: j$ Q7 Q- k6 y*SET,%_FNCNAME%(0,10,1), 0.0, -2, 0, 1, 2, 17, -1   
/ u$ o1 L" B4 W6 j*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 1, 19, 3, 1
' h  D8 Q# y! _" k% J; \, k+ V*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   
4 r, a  g! H2 _& _! [1 Z*SET,%_FNCNAME%(0,16,1), 0.0, -2, 0, 1, -4, 3, -1   
*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   
7 G9 F+ A4 j1 \% Q*SET,%_FNCNAME%(0,20,1), 0.0, -1, 7, 1, -1, 0, 0
" x) R3 ^7 {& m5 k*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 O3 f" M! O+ V  U( @; W) Z3 `*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   
+ O' ]& u) B2 D+ W/ z5 \" b% m*SET,%_FNCNAME%(0,25,1), 0.0, -1, 0, 1, -4, 3, 3
1 d: w1 g2 \4 k. k7 a# @*SET,%_FNCNAME%(0,26,1), 0.0, -3, 0, 1e-9, 0, 0, -1  
  v/ I+ o4 `. V*SET,%_FNCNAME%(0,27,1), 0.0, -4, 0, 1, -1, 4, -3   
*SET,%_FNCNAME%(0,28,1), 0.0, -1, 0, 1, 21, 3, -4   
: v1 A: }- V$ r( O8 O: q*SET,%_FNCNAME%(0,29,1), 0.0, -3, 0, 1, 20, 1, -1   
*SET,%_FNCNAME%(0,30,1), 0.0, -1, 0, 0, 0, 0, 0  
+ V7 P* L# f( \0 X: E7 j*SET,%_FNCNAME%(0,31,1), 0.0, -4, 0, 1, 0, 0, -1
9 M$ y3 Y( e) T  x9 u+ v*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
; \* f2 @) z+ Z' O3 ^+ v*SET,%_FNCNAME%(0,34,1), 0.0, -4, 0, 1e-9, 0, 0, -1  
# [! V3 w7 }# c8 _" X*SET,%_FNCNAME%(0,35,1), 0.0, -5, 0, 1, -1, 4, -4   
. a, j2 @. D8 I! V8 n& Q*SET,%_FNCNAME%(0,36,1), 0.0, -1, 0, 2, 0, 0, -5
*SET,%_FNCNAME%(0,37,1), 0.0, -4, 0, 1, -5, 17, -1   
*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   
*SET,%_FNCNAME%(0,40,1), 0.0, -1, 0, 0, 0, 0, 0  
- m4 L: N2 R% w. |: {9 A, M( s*SET,%_FNCNAME%(0,41,1), 0.0, -3, 0, 1, 0, 0, -1
$ u4 S5 _5 v1 m6 ~9 x1 _: Q) K*SET,%_FNCNAME%(0,42,1), 0.0, -5, 0, 1, -1, 2, -3   
. J( K9 T& D! S& ~6 D+ J0 U*SET,%_FNCNAME%(0,43,1), 0.0, -1, 0, 1, -5, 3, 3
/ I8 l. b8 d1 |  h. z8 H*SET,%_FNCNAME%(0,44,1), 0.0, -3, 0, 1e-9, 0, 0, -1  
*SET,%_FNCNAME%(0,45,1), 0.0, -5, 0, 1, -1, 4, -3   
' u' M- m* S; p& T8 W*SET,%_FNCNAME%(0,46,1), 0.0, -1, 0, 3, 0, 0, -5
*SET,%_FNCNAME%(0,47,1), 0.0, -3, 0, 1, -5, 17, -1   
- F+ @+ a& C% ?0 d1 a# E9 c% Y+ E*SET,%_FNCNAME%(0,48,1), 0.0, -1, 0, 1, 23, 3, -3   
*SET,%_FNCNAME%(0,49,1), 0.0, -3, 0, 1, -4, 1, -1   
* O! Z+ T; s4 z: c7 \*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
*SET,%_FNCNAME%(0,52,1), 0.0, -5, 0, 1, -1, 2, -4   
*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   
* K  S5 d8 U8 ?& L*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   
' R0 k8 w, l) b2 j*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*(-   
$ m# M) z* o# I# Y+ l$ E! @& y* s, Q5 ]! {Y}/1e-9)+A2*(-{Y}/1e-9)^2+A3*(-{Y}/1e-9)^3+ A4*(-{Y}/1e-9)^4)
; J2 ~. U  `- `) b! V!-->