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

xy112

高斯激光热源下的n层材料瞬态热传导方程
高斯, 激光, 热源, 方程, 热传导高斯, 激光, 热源, 方程, 热传导
出于需要， 写了个 n层材料的瞬态热传导方程， 由表面激光加热。每层的热源由其光学性质N= n + ik决定. 通过解Maxwell方程获得激光在材料中的强度分布， 由多项式 a0+a1*y+a2*y^2+a3*y^3+a4*y^4拟和获得。
3 f  x. `( D# W% W; aP=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决定了其内部的热源分布。
5 K* [, c/ k. O5 C  X  
本程序可已根据问题选择所需的材料层数，所有参数都可以改变。  
9 x( j+ f% _1 S9 l* W. U  
$ c0 z8 A* e3 v- u5 J; Y1 k9 G以下是命令流，如果有不合适的地方，欢迎大家提出、讨论！
  
FINISH  ! Make sure we are at BEGIN level   
! ?) f# \& L+ h# N( M/ k/Clear
/FILNAME,Media,0  
. X, Y9 H" |) f2 ~  p$ y* X  
9 f  A; \/ n5 b2 {0 {4 W& w& S/PREP7
! g$ C. ]% q& a/PNUM,VOLU,1
  
! In simulation, the X-Z is the media plane, while Y is the depth direction
  
Layer=5                       !Total layer number
nm= 1e-9                    !length unit
Power=0.005               ! Laser incident power   W
r=5e-7                        ! Laser spot radius      m
V=0                           ! Laser moving veolcity, along z direction,  m/s
  
7 [9 U' l9 O, Y  K; O5 fpw=Power*1000                 !For display Title
+ T3 k* W3 T  T) Y. \/ J' c  yrh=r/1e-6                           !For display Title
; @$ P8 B1 f7 F/ }  
/TITLE,%Layer% layer media, Laser Power=%pw%mW, Heating radius=%rh%um, Velocity=%V%m/s  
) ]5 D' K  f3 q- e) ~  
2 r$ G7 L/ i# B# V% ^*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
! u- D1 ]7 @% X1 o1 @: P2 V*DIM,Cp,ARRAY,Layer           ! Array input for each layer heat capacity
*DIM,Den,ARRAY,Layer          ! Array input for each layer density
" s0 w# ^0 z6 L6 o, h; R0 E) \+ H*DIM,div,ARRAY,Layer          ! Array input for each layer thickness division
* h; R5 h1 |0 j*DIM,dscale,ARRAY,Layer       ! Array input for each layer thickness division scale
3 z6 o- Y0 M' T3 E3 v$ r- p2 g*DIM,A0,ARRAY,Layer           ! Array input for each layer a0 coefficient for heat generation curve
/ ?! Q( Z; A& L: M*DIM,A1,ARRAY,Layer           ! Array input for each layer a1 coefficient for heat generation curve
9 M8 R+ x$ O: n! g) X*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
4 _$ ^( x* ^. P( P+ g0 b*DIM,A4,ARRAY,Layer           ! Array input for each layer a4 coefficient for heat generation curve
" r% @/ Z5 J; h- `5 h9 Y9 r$ [  
scale=10000                   !for vovlap operation
: `8 ~* ^( o6 r, K0 `# V( fPulse_Duration=1e-8           !Laser pulse duration
End_Time=2e-8                 !Final time for simulation
Step=20                       !Time interval  
T0=50                         !Environment temperature   
$ p3 G8 N. I2 v0 c  
Length=4000*nm                !Length of the block in z direction
$ F7 `$ C5 s, v9 ]0 b$ M, wWidth=2000*nm                 !Width of the block in x direction
8 Z6 f7 ?4 @1 ?$ ~" ^2 N# X. Z  
!Layer properties, from top to bottom  
  
!First layer thermal property
* u# _, Z3 j0 g( [9 S  }Num=1                        
k_L(Num)=1.4                  
! C( c3 F! ?# p# W2 p' tk_N(Num)=1.4
Den(Num)=4000
) R, C5 \5 p1 K# l. [Cp(Num)=500
( z! I' H( D* M/ }+ L( l4 BThickness(Num)=50*nm
div(Num)=5
; r. A' b4 |8 l+ D4 ~3 {: kdscale(Num)=1                 !last/first in -Y direction
A0(Num)=0
7 a  D# ^' \' }% jA1(Num)=0
) x+ S4 b6 {" B& R1 C. i' f7 ]2 J0 FA2(Num)=0
A3(Num)=0
A4(Num)=0
  
!Second layer thermal property
Num=2
k_L(Num)=16
k_N(Num)=16
Den(Num)=3000
Cp(Num)=900
Thickness(Num)=20*nm
div(Num)=2
. r; a: ^" ^6 g" Ldscale(Num)=1
A0(Num)=807.929E+6
4 A. U# y+ h# C8 P% P9 J' `% V, a' Z9 qA1(Num)=-43.896E+6
A2(Num)=934.529E+3
A3(Num)=-8.864E+3
* v' [3 S/ \! n, R* a6 o& F8 }) vA4(Num)=31.703
- h& H6 R: R8 a  
!Third layer thermal property
. u% B6 |2 A% o; F, c1 Z( mNum=3
; T. g7 Q& s6 y6 ?7 Hk_L(Num)=1.5
. Z$ C0 g% Y$ B6 q5 xk_N(Num)=1.5
Den(Num)=4000
Cp(Num)=500
Thickness(Num)=100*nm
. D" g! ]$ T' n( I) Ydiv(Num)=5
dscale(Num)=1
5 i. W. X. \/ p7 TA0(Num)=0
4 Z+ t) A- ~& V" ZA1(Num)=0
A2(Num)=0
A3(Num)=0
A4(Num)=0
, I( X6 h: L/ e0 O" N. X% c1 Q  
!Fourth layer thermal property
- M( I0 _3 T) h- c  fNum=4
k_L(Num)=100
k_N(Num)=100
Den(Num)=3900
Cp(Num)=700
: j: \3 ~- p7 d# F. d" O. a6 Y" hThickness(Num)=50*nm
div(Num)=2
4 m1 L# @; S+ ~8 w  p( Ddscale(Num)=2
A0(Num)=1.627E+9
* q4 I. S3 ~7 l1 g, L. CA1(Num)=-32.146E+6
A2(Num)=238.195E+3
A3(Num)=-784.173
1 B9 u' M! w# b! W( U) vA4(Num)=967.716E-3
  
) M: N+ S- t0 ~' w6 }% X+ C!Fifth layer thermal property
Num=5
: B8 J0 G5 D0 z1 G2 D1 nk_L(Num)=0.2
. |" @9 B1 m) n) M0 L# n" i3 Ck_N(Num)=0.2
Den(Num)=3900
Cp(Num)=450
Thickness(Num)=100*nm
div(Num)=5
dscale(Num)=2
1 d/ `% R- u5 H; i; z: HA0(Num)=0
A1(Num)=0
A2(Num)=0
, a& {! k& _2 K1 }  t1 u7 h( ?A3(Num)=0
A4(Num)=0
- h; g; _6 ^; C' C5 G  X  
! Define the thermal element
ET,1,SOLID70  
KEYOPT,1,2,0  
KEYOPT,1,4,0  
4 Q& r2 ]6 h- v9 S# e2 RKEYOPT,1,7,0  
KEYOPT,1,8,0  
( W/ v0 b  Y; @6 j1 C) v  
!For post-precessing
nPoint=200
display_r=1e-6
2 a- S3 X, F: C) }  
TSum1=0
TSum2=0
5 p. j+ @8 c1 r" k3 y, u: Z*do,Num,1,Layer,1
) |& }& v, `/ a$ @/ I  UIMP,Num,DENS, , ,Den(Num),                                        !Define density of material Num
3 g- I9 R( ~2 F6 p  UIMP,Num,KXX,KYY,KZZ,k_L(Num),k_N(Num),k_L(Num),                   !Define thermal conductivity of material Num
  UIMP,Num,C, , ,Cp(Num),                                            !Define heat capacity of material Num   
  TSum1=TSum1+Thickness(Num)  
( t2 C' l) _& F. l  BLOCK,-Width/2,Width/2,-TSum1,-TSum2,-Length/2,Length/2,           !Draw the layer of Num  
' q5 L6 k$ W0 I& e5 t. c  TSum2=TSum2+Thickness(Num)
*Enddo
, w1 [) m4 [: x9 a: s! n8 W  
d_total=TSum1                 !Total thickness of the multilayer                                                  
# ~7 t' a8 u: T4 n  
VLSCALE,ALL, , , scale, scale, scale,  ,1,1,
2 _, G1 e0 |- A8 I2 Y5 F  
*IF,Layer,GT,1,THEN
% u# N# u5 J  O7 y! s   VGLUE,ALL
*ENDIF
  
1 ~' n* R5 U) l  IVLSCALE,ALL, , , 1/scale, 1/scale, 1/scale, ,1,1,
# M' K0 p. [( J% F/ ^$ i, s: \  
TSum1=0
! |, g7 l9 u! D* MTSum2=0
  
0 _7 G; E( b' g% S$ VSize=100*nm                             !Area mesh size, in X-Z plane
aesize,all,Size,  
  L5 `( x; {* P- e& j; O  
*do,Num,1,Layer,1
   TSum1=TSum1+THICKNESS(Num)
$ F2 u- n  c$ L% l$ `   lsel,all
  
% g4 C) e0 y  n+ u   !Determine the mesh scale and divison for each layer at Y direction
   *do,i,-1,1,2
      *do, j,-1,1,2
         lsel,s,loc,y,-TSum1,-TSum2
* H0 \4 \$ ?" l% {8 S. z% k5 V         lsel,r,loc,x,i*Width/2,i*Width/2
8 F2 N8 y) _# d5 i         lsel,r,loc,z,j*Length/2,j*Length/2
0 V7 ^  ^2 O6 X! K' V1 l) T         *GET,Line_Num,LINE, 0, NUM,MIN,
         *GET,KP1,LINE,Line_Num,KP,1,
         *GET,KP2,LINE,Line_Num,KP,2,
9 [# o3 b" ?$ ~0 h% U         *GET,KP1_Y,KP,KP1, LOC, Y,  
         *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
             lesize,all,THICKNESS(Num)/div(Num), , ,1/dscale(Num),1, , ,1,
         *ENDIF
: |6 ~+ \8 i% b( E      *ENDDO
   *ENDDO   
) N* {0 a6 v0 m! M, F6 S  
% w; `6 b& C( A& [9 ?+ k   asel,all
% u; \. Q; p" L, Q2 l2 o% W0 g* M  
5 O7 t+ F1 Q! @0 G# h   esize,Size,0,
3 b4 e: W$ f# A# Z- K0 {  
1 Q0 r- @; U' b/ }& E: f   VSEL,S,LOC,Y,-TSum1,-TSum2
  
   TYPE,1
   MAT,Num
% F" g5 N  a2 b   mshkey,1  
0 k- H' w% {+ g   mshape,0,3d  
! f- I0 M% D% C* b) k' {( F   vmesh,ALL             !Mesh the layer Num
; E/ V% I2 }" x' \  
; m! k. G' t' l' i$ E6 _   allsel,all
) v9 d! C( ~) [( i  
   TSum2=TSum2+Thickness(Num)
: z) X, z  o8 r5 t3 W  Z  
+ s4 I7 b& E' P, m! x: x! X6 k*enddo
5 B  M: Q( Y& \% e  
/REP,FAST  
+ u% j" `; _0 j% J$ J$ a/psf,hflu,,2,0  
/PNUM,LINE,1  
finish  
/solu  
3 N8 Q0 U2 [8 |5 F* P  
antype,trans             !Transient simulation
  {5 P( f* H/ N' E+ e' e( c  ?6 `tunif,T0                 !Ambient temperature
; t/ n0 U0 p4 p0 n5 c& MTREF,T0
timint,on !time integration = on  
kbc,1  
OUTRES,ALL,ALL,          !Output all the result
  
/view, 1 ,1,1,1  
/ang, 1  
/rep,fast  
  L' O" n* a; S- D6 }6 M6 r; n  
*do,Num,1,Layer,1
   *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
1 q8 W3 z' b8 C  G5 ]  
!Start simulation
( X2 A" {' i1 D; i*do,tm,1e-15,End_Time,End_Time/Step  
  
/ T0 p* q7 v' D; l5 W2 h% r5 {   time,tm  
! }$ c3 l" D& |- J: n) U7 d4 F   autots,0  
   deltim, , , ,1  
# B3 `5 W& M8 {: f0 S! @   kbc,0  
  
   allsel,all  
   TSum1=0
   TSum2=0
   Num=0
  
- S6 H/ Z3 V: \* e0 k- y3 X$ [   *if,tm,LE,Pulse_Duration+1e-12,then  
# e0 x& s6 c3 D' J$ a        
; Q  z5 q/ q. E9 a      !Apply the heat generation in each layer  
      Num=1
      TSum1=TSum1+THICKNESS(Num)
) ]: J# o1 f4 T      nsel,s,LOC,y,-TSum1,-TSum2   
4 C9 x7 p( A8 e      BF,all,HGEN,%Layer1%
      TSum2=TSum2+Thickness(Num)
/ q. u( K0 @& p      
2 U0 M# x* Z# A' l3 P( z      Num=2
* D" W( z3 Q, b: w6 z. M      TSum1=TSum1+THICKNESS(Num)
      nsel,s,LOC,y,-TSum1,-TSum2   
2 p) {2 P. `# K+ ~$ d/ _      BF,all,HGEN,%Layer2%
      TSum2=TSum2+Thickness(Num)
   
9 z, q4 G; ~* B3 a& r* X/ A& b      Num=3   
7 s6 p8 S# F& ?2 A* U      TSum1=TSum1+THICKNESS(Num)
      nsel,s,LOC,y,-TSum1,-TSum2   
; n+ X4 F  I6 ^4 b; r6 b( }2 }0 B      BF,all,HGEN,%Layer3%
, {, @1 d. g; j( ~& F      TSum2=TSum2+Thickness(Num)
   
& {6 l0 B- l/ P9 L9 H2 V7 J" t      Num=4
" y+ i3 P/ y5 i4 J2 \      TSum1=TSum1+THICKNESS(Num)
      nsel,s,LOC,y,-TSum1,-TSum2   
+ _: L: ?# L3 }: [$ X      BF,all,HGEN,%Layer4%
3 u' @  x6 T8 U      TSum2=TSum2+Thickness(Num)
   
      Num=5
      TSum1=TSum1+THICKNESS(Num)
0 {8 S) L: a$ e. m9 a, o+ x      nsel,s,LOC,y,-TSum1,-TSum2   
      BF,all,HGEN,%Layer5%
      TSum2=TSum2+Thickness(Num)
: Q. V( G% |! S7 m- q1 b8 M' g# j  
) A) z1 S9 E" o4 _. M   *else
0 R4 g1 k  g& z! \      allsel,all  
      BFDELE, ALL,HGEN
   *endif  
  
   allsel,all
: R0 p4 b# l1 g; P& z- w   asel,s,LOC,y,-1e-9,0
!   DA,ALL,TEMP,T0  
- g5 A4 v! ?. v: N$ r$ s   allsel,all
+ u1 V) v! I6 E+ V  
   /status,solu  
   solve  
( \5 [& E2 K+ m( h*enddo  
3 N' v' x) t! Y  ^+ W  
' t5 w9 P. f5 rfinish  
/POST1
# Y7 L* q7 `1 T1 v& c$ G  
; ]. z! a2 y# u: i5 O$ v% H# s!*ULIB, Post,mac
5 V% A6 u; m+ R$ M% Q- W!*USE,TTV    !Temperature decay at each layer thickness
!*USE,TRV,1  !Normal Temperature profile at each layer at given time  
!*USE,TRL,1  !Plane Temperature profile at each layer at given time  
2 @  R) o9 Y, s$ w! Z7 V9 x! o5 v  
!Laser power absorption in each layer
  u2 h" b. m# S  r!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)      
# n5 r8 }5 e" _. x  e- d  
//UG.mac
  
!FileName=ARG1
7 A( w! R/ f$ ?% n% Y0 c! v!Power=ARG2
!r=ARG3
5 F) c; w$ k" F4 h1 }!V=ARG4
/ |: @0 \- N$ L  C/ f5 V# M!A0=ARG5
!A1=ARG6
!A2=ARG7
! f* L# F  D( x4 X!A3=ARG8
!A4=ARG9
  
/ I, i2 _! ~& M" b! Q6 J" n*DEL,_FNCNAME   
; d! `! e# L$ ]" D+ m*DEL,_FNCMTID   
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
- ^. ?' c9 T, ~- y2 x*DEL,_FNC_C4
*DEL,_FNC_C5
*DEL,_FNC_C6
7 h! \3 D5 K% d9 {* [" d*DEL,_FNC_C7
*DEL,_FNC_C8
*SET,_FNCNAME,ARG1   
*DIM,_FNC_C1,,1  
*DIM,_FNC_C2,,1  
*DIM,_FNC_C3,,1  
*DIM,_FNC_C4,,1  
*DIM,_FNC_C5,,1  
+ n( m5 L2 i; }+ Y8 u*DIM,_FNC_C6,,1  
: w8 l4 s9 S7 Y# U*DIM,_FNC_C7,,1  
9 i* r- Q' A7 K- V: _" ?*DIM,_FNC_C8,,1  
*SET,_FNC_C1(1),ARG2   
*SET,_FNC_C2(1),ARG3
6 Z- m; o. [1 t) ]  v*SET,_FNC_C3(1),ARG4
# B( e% }& m' z4 L, }*SET,_FNC_C4(1),ARG5
*SET,_FNC_C5(1),ARG6
" J  o" @7 F  O4 O  ^*SET,_FNC_C6(1),ARG7
1 v- ^6 n) c) g7 T1 ]) g7 ]*SET,_FNC_C7(1),ARG8
9 }- b- @2 X( U5 D/ B! [" p3 u/ Y, `*SET,_FNC_C8(1),ARG9
! /INPUT,HG.func
3 M: [# m) b2 _" w5 F, E& L*DIM,%_FNCNAME%,TABLE,10,61,1   
1 l) h8 ~9 D4 G& \3 S!   
5 S; ]- h3 Y- Y" y6 C# B! Begin of equation: Power/r^2/{PI}*exp(-({X}^2+({Z}-V*{TIME})^2)/r^2)*(A0+A1*   
+ M. z1 K( R9 B; q3 p  I! (-{Y}/1e-9)+A2*(-{Y}/1e-9)^2+A3*(-{Y}/1e-9)^3+ A4*(-{Y}/1e-9)^4)   
/ j3 O: b2 B- Q0 {5 M4 W, _3 i, D0 r*SET,%_FNCNAME%(0,0,1), 0.0, -999   
2 x; U! P$ X% h: T! a, U# @) |  W*SET,%_FNCNAME%(2,0,1), 0.0  
; u" }+ F( _3 O0 B: O( V*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
8 w- e9 ~0 N: u; r*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
( y* `+ G5 }; R* i4 r3 V*SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)%
*SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)%
3 U$ @- |, M( H! Y+ t- |0 s8 M*SET,%_FNCNAME%(8,0,1), %_FNC_C6(1)%
*SET,%_FNCNAME%(9,0,1), %_FNC_C7(1)%
*SET,%_FNCNAME%(10,0,1), %_FNC_C8(1)%   
$ T  F5 P2 Y9 E- N2 t*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 2, 0, 0, 18  
6 e1 _  s5 J2 Y; I*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 18, 17, -1   
( s% j, X" d4 n" R6 ?; c. j0 ^% _*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   
1 Z2 M" _: Q5 k% \*SET,%_FNCNAME%(0,5,1), 0.0, -3, 0, 1, -1, 4, -2
2 O$ [) N- E# L/ q# I# \* N' U*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  
*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   
*SET,%_FNCNAME%(0,10,1), 0.0, -2, 0, 1, 2, 17, -1   
; C: L  }9 Z3 j*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 1, 19, 3, 1
*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   
5 X+ C: `/ G* @7 D1 N1 b* h  q9 x*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   
6 D6 M$ m3 |" X- G" I$ x$ d*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
*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  
9 j+ b! u2 P. H( T: f$ a5 E*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
6 ?8 N" W) u' M+ v8 K+ y% O*SET,%_FNCNAME%(0,26,1), 0.0, -3, 0, 1e-9, 0, 0, -1  
/ @* M" b9 _2 c9 [*SET,%_FNCNAME%(0,27,1), 0.0, -4, 0, 1, -1, 4, -3   
- M! ]% g7 h! e( p0 S% k) u*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   
0 ~1 n& c6 q8 M+ e*SET,%_FNCNAME%(0,30,1), 0.0, -1, 0, 0, 0, 0, 0  
4 Q" C6 b4 x7 ]*SET,%_FNCNAME%(0,31,1), 0.0, -4, 0, 1, 0, 0, -1
! P" f# M5 j) `4 q*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
*SET,%_FNCNAME%(0,34,1), 0.0, -4, 0, 1e-9, 0, 0, -1  
$ s$ `7 t2 e' ~& I  d6 M) y' V; T*SET,%_FNCNAME%(0,35,1), 0.0, -5, 0, 1, -1, 4, -4   
7 X6 W9 e# z( o% ?5 T1 f*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   
0 }& ?: Y- a1 K*SET,%_FNCNAME%(0,39,1), 0.0, -4, 0, 1, -3, 1, -1   
5 H6 _- v, n' p*SET,%_FNCNAME%(0,40,1), 0.0, -1, 0, 0, 0, 0, 0  
*SET,%_FNCNAME%(0,41,1), 0.0, -3, 0, 1, 0, 0, -1
*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
( ~, S+ q7 K; p0 J*SET,%_FNCNAME%(0,44,1), 0.0, -3, 0, 1e-9, 0, 0, -1  
5 f; x& J$ q# f5 h; J1 m  Z*SET,%_FNCNAME%(0,45,1), 0.0, -5, 0, 1, -1, 4, -3   
9 y3 u$ y/ k/ ?0 _' H, T*SET,%_FNCNAME%(0,46,1), 0.0, -1, 0, 3, 0, 0, -5
0 e7 j: ^) a; E7 ~4 E  F*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   
' n% T, @  [# Q: r% T% m+ p*SET,%_FNCNAME%(0,49,1), 0.0, -3, 0, 1, -4, 1, -1   
*SET,%_FNCNAME%(0,50,1), 0.0, -1, 0, 0, 0, 0, 0  
; O( S# H) G8 a: u6 P& `*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   
*SET,%_FNCNAME%(0,56,1), 0.0, -1, 0, 4, 0, 0, -5
3 \5 b0 K% [* D! k: e0 ]*SET,%_FNCNAME%(0,57,1), 0.0, -4, 0, 1, -5, 17, -1   
; o; e9 T& m7 i: n" z7 q*SET,%_FNCNAME%(0,58,1), 0.0, -1, 0, 1, 24, 3, -4   
( `: s2 ?. W; u. D" T& Z1 g+ P  S3 n*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)
!-->