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

xy112

高斯激光热源下的n层材料瞬态热传导方程
* }9 G8 T2 X) h  s' Y; l高斯, 激光, 热源, 方程, 热传导高斯, 激光, 热源, 方程, 热传导
出于需要， 写了个 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决定了其内部的热源分布。
  
6 z3 c& }4 T9 ]+ b( |) f本程序可已根据问题选择所需的材料层数，所有参数都可以改变。  
2 {* T  W' ]0 T- e8 i8 J& E- t  
! z- L- [  p* g6 d) D! u以下是命令流，如果有不合适的地方，欢迎大家提出、讨论！
  
/ P6 d$ ^7 j0 Y  U5 }$ ~: e" Q/ MFINISH  ! Make sure we are at BEGIN level   
/Clear
: b( b1 _) d( x0 l- ]( m. G/FILNAME,Media,0  
  
- ^. o8 p6 q6 g: {3 L$ ]1 L; I/PREP7
/PNUM,VOLU,1
  
! f6 c# i6 a: e; P, w$ m1 l$ D! In simulation, the X-Z is the media plane, while Y is the depth direction
9 |1 M' j/ t2 b; Q! A: B  
% o. w' p# o4 C5 r. G! i5 U8 g/ K5 l* ILayer=5                       !Total layer number
& u/ ~* D3 o, K  w" [  Cnm= 1e-9                    !length unit
Power=0.005               ! Laser incident power   W
1 K! M/ e& O$ w1 Wr=5e-7                        ! Laser spot radius      m
V=0                           ! Laser moving veolcity, along z direction,  m/s
2 M0 H5 e( I- h* t7 ]- S  
pw=Power*1000                 !For display Title
* L* V" s5 Z0 u- W1 Urh=r/1e-6                           !For display Title
% z7 t/ C6 ?: @. Z, D& t  k  
/TITLE,%Layer% layer media, Laser Power=%pw%mW, Heating radius=%rh%um, Velocity=%V%m/s  
  
*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
*DIM,Cp,ARRAY,Layer           ! Array input for each layer heat capacity
*DIM,Den,ARRAY,Layer          ! Array input for each layer density
7 d/ y5 Q5 x2 Q) d# ?*DIM,div,ARRAY,Layer          ! Array input for each layer thickness division
*DIM,dscale,ARRAY,Layer       ! Array input for each layer thickness division scale
: E; H2 H) B  p) E" \6 n: H; P. S*DIM,A0,ARRAY,Layer           ! Array input for each layer a0 coefficient for heat generation curve
4 N' N' W0 d. {1 ]*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
3 w" `, }) a0 F$ ]; u" u*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
8 [" t+ X3 y1 o& }  TEnd_Time=2e-8                 !Final time for simulation
5 f6 @1 O- U% `) K% @8 zStep=20                       !Time interval  
T0=50                         !Environment temperature   
% y: f! Q' a( _! r4 ^  
Length=4000*nm                !Length of the block in z direction
6 Y" k7 T9 s. k. M& ^Width=2000*nm                 !Width of the block in x direction
6 j$ l  A/ ]: n4 A3 k  
!Layer properties, from top to bottom  
  
  z; I! O3 f. |2 [0 a# Y!First layer thermal property
' ]) T5 @5 B( V2 u9 jNum=1                        
" i. F1 K4 p6 `9 E2 W# jk_L(Num)=1.4                  
k_N(Num)=1.4
Den(Num)=4000
4 N$ H  {& B! `- BCp(Num)=500
( b& r/ R& }; CThickness(Num)=50*nm
1 R! s1 }% ?$ Q2 W8 v' Jdiv(Num)=5
! a- W( K8 p% O; Y) `7 h, {/ E! p* L$ wdscale(Num)=1                 !last/first in -Y direction
8 X5 Z8 I0 X' v% y6 i, B( tA0(Num)=0
6 }2 S+ N+ c# @2 J- Y* E3 UA1(Num)=0
2 F5 N; u) q  D0 }& f3 p/ X  V! ]A2(Num)=0
/ U1 s; a4 z8 }0 O7 Q( \, lA3(Num)=0
8 _7 b7 U; {8 x! pA4(Num)=0
  i4 |6 V. q: Z, X( ~- z/ L  
* w# B# _( b& H$ J" Q  p9 @!Second layer thermal property
7 y( \6 ^2 {1 m2 ], S- e1 _# v& DNum=2
k_L(Num)=16
k_N(Num)=16
9 p5 ^) f7 L/ b1 ~2 oDen(Num)=3000
# |( I( N, ?. p+ P, _Cp(Num)=900
4 o  i5 n+ N& g! fThickness(Num)=20*nm
9 a* o, x7 z0 ~- j% adiv(Num)=2
dscale(Num)=1
  t, E8 W4 d0 d! l; VA0(Num)=807.929E+6
$ W& c  C+ c; I2 ~2 OA1(Num)=-43.896E+6
A2(Num)=934.529E+3
9 d; e/ J  f* hA3(Num)=-8.864E+3
A4(Num)=31.703
  
+ Y3 P* d+ {& O8 J, _+ T!Third layer thermal property
5 [3 v8 U& _2 J( W0 rNum=3
k_L(Num)=1.5
k_N(Num)=1.5
Den(Num)=4000
Cp(Num)=500
Thickness(Num)=100*nm
) j4 b9 [+ {! B" L1 Q$ r+ ?; R$ Fdiv(Num)=5
6 t4 m+ v2 l5 h. p  k2 I, p. Kdscale(Num)=1
A0(Num)=0
A1(Num)=0
A2(Num)=0
A3(Num)=0
" ?1 z% R" w/ s3 J2 ]A4(Num)=0
) G# E. z9 u, r( s  
- O; e% u" J" n/ x: N' X& g!Fourth layer thermal property
Num=4
k_L(Num)=100
k_N(Num)=100
; X- @* R" A8 ]! I" F/ f5 qDen(Num)=3900
Cp(Num)=700
Thickness(Num)=50*nm
! R% W3 H1 o$ Fdiv(Num)=2
dscale(Num)=2
* B* M# s/ j7 H; R5 q6 ^A0(Num)=1.627E+9
A1(Num)=-32.146E+6
' R6 K, V& i4 X" m. z, BA2(Num)=238.195E+3
A3(Num)=-784.173
& m  g3 k( g5 W0 A  P( U6 ]) qA4(Num)=967.716E-3
1 d7 c6 I% q+ m: g9 I; N  
4 ^! M0 g$ U/ F' p2 [- g$ c!Fifth layer thermal property
Num=5
k_L(Num)=0.2
" [! s( f  C+ n! _k_N(Num)=0.2
Den(Num)=3900
Cp(Num)=450
' r( H6 z' N( V9 L$ C7 ?6 W7 dThickness(Num)=100*nm
4 A( M& b# N: _) b: L, e- s) kdiv(Num)=5
dscale(Num)=2
) _  E- N! a, S+ u4 ^. |A0(Num)=0
A1(Num)=0
A2(Num)=0
, Q' O3 `1 I- P* m# ^1 NA3(Num)=0
* ^% [4 O/ K+ ~( k. LA4(Num)=0
  
! Define the thermal element
2 r" p- x* R! @6 S. u) \. oET,1,SOLID70  
KEYOPT,1,2,0  
* v: B  _6 |0 Z* s. |0 d0 @& fKEYOPT,1,4,0  
7 s* i* C. _) w6 aKEYOPT,1,7,0  
KEYOPT,1,8,0  
3 s4 T' b7 T- N, M  
!For post-precessing
# `  f/ ~4 D4 i/ JnPoint=200
display_r=1e-6
  
( m) D4 S, t; {+ g6 e  w! j; ]TSum1=0
TSum2=0
" ]* c& q. y$ A1 t*do,Num,1,Layer,1
- w$ u4 s8 q) Z9 |; f, n+ e9 K4 [  UIMP,Num,DENS, , ,Den(Num),                                        !Define density of material Num
2 t! g6 l# [7 x! ?+ E) M  UIMP,Num,KXX,KYY,KZZ,k_L(Num),k_N(Num),k_L(Num),                   !Define thermal conductivity of material Num
; t; `  R* x2 ^/ X1 o* J  UIMP,Num,C, , ,Cp(Num),                                            !Define heat capacity of material Num   
  TSum1=TSum1+Thickness(Num)  
  BLOCK,-Width/2,Width/2,-TSum1,-TSum2,-Length/2,Length/2,           !Draw the layer of Num  
% g$ a- {  Z' [+ J# g  TSum2=TSum2+Thickness(Num)
7 s1 I8 K. J; }3 U) C1 N' B1 b; ~- ]*Enddo
  
' P: D1 F; J3 I1 l$ P& \d_total=TSum1                 !Total thickness of the multilayer                                                  
/ i! [1 r5 B" v5 h! r! `9 P# a  
% E/ f2 v8 w3 n9 G$ f+ T9 G5 qVLSCALE,ALL, , , scale, scale, scale,  ,1,1,
  
( y+ Q! D, `; z- ]2 p4 {, I*IF,Layer,GT,1,THEN
" U3 Y& f- ~4 x' j( D9 m) d8 {   VGLUE,ALL
9 }* D# g$ C6 C0 b7 Q0 o. a+ @1 j*ENDIF
  
VLSCALE,ALL, , , 1/scale, 1/scale, 1/scale, ,1,1,
1 Q) y' n  K2 Y0 N  b% Z4 @  
TSum1=0
# p! m% R+ P( O2 aTSum2=0
  
Size=100*nm                             !Area mesh size, in X-Z plane
aesize,all,Size,  
% M9 I& _- v# P# s- o  
# e  G, Y7 c/ y' S+ ^6 c7 v4 f*do,Num,1,Layer,1
   TSum1=TSum1+THICKNESS(Num)
8 D; U' u/ R% `6 j9 d. [! x   lsel,all
4 d  e3 ~" c6 E. b7 u5 @' v3 R  
   !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
; h: K# H3 ^6 l2 @7 q0 F  o         lsel,r,loc,x,i*Width/2,i*Width/2
         lsel,r,loc,z,j*Length/2,j*Length/2
         *GET,Line_Num,LINE, 0, NUM,MIN,
5 D5 E, j  C0 m1 n  l- d5 e         *GET,KP1,LINE,Line_Num,KP,1,
         *GET,KP2,LINE,Line_Num,KP,2,
         *GET,KP1_Y,KP,KP1, LOC, Y,  
# M! g- \* B, ~: X         *GET,KP2_Y,KP,KP2, LOC, Y,
7 ~  F. @, ?0 r) Y         *IF,KP1_Y,GT,KP2_Y,THEN
7 N3 i# l; u2 X) @& w  o             lesize,all,THICKNESS(Num)/div(Num), , ,dscale(Num),1, , ,1,
         *ELSE
1 ]7 ?, }. N' j) L             lesize,all,THICKNESS(Num)/div(Num), , ,1/dscale(Num),1, , ,1,
         *ENDIF
5 a4 B# [/ d6 _/ n- I7 u      *ENDDO
   *ENDDO   
  
   asel,all
: Q: g, o' Q" F8 H" s' p  
1 E( c7 g( \) L+ o. A. v/ S9 k   esize,Size,0,
7 M- G: M4 e$ l" W/ S* Y  
   VSEL,S,LOC,Y,-TSum1,-TSum2
  
   TYPE,1
   MAT,Num
   mshkey,1  
   mshape,0,3d  
& p) Q5 J: U$ I, u, z! _. c( h; c   vmesh,ALL             !Mesh the layer Num
( |/ ~! }; O( L  
   allsel,all
, `" k( z2 G" |  
   TSum2=TSum2+Thickness(Num)
, J2 e* ~  A& ?( g; \& J2 v  
7 h9 M( @5 r( R*enddo
  
/REP,FAST  
4 Z5 c: \$ B! S: }) ]% h/psf,hflu,,2,0  
$ ~8 o" I8 n/ c/PNUM,LINE,1  
finish  
, X* A. N& U+ A& p" A/solu  
  
antype,trans             !Transient simulation
tunif,T0                 !Ambient temperature
TREF,T0
( F( O% n/ Y0 Q% P: P$ z. E4 Ntimint,on !time integration = on  
5 x" F- a- l& |" v& v: ]* bkbc,1  
OUTRES,ALL,ALL,          !Output all the result
  
/view, 1 ,1,1,1  
/ang, 1  
. @% e2 T6 t) X4 {  k% x9 l/rep,fast  
  a* N4 U' L6 d+ ~( A! |7 Z  
*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
8 T6 J& r7 c9 x5 `: l2 ]: e6 E1 J*enddo
  
! X3 ^$ k4 \  a+ C6 v! D& y8 r!Start simulation
! V- D: \! Q% b. E- z7 d*do,tm,1e-15,End_Time,End_Time/Step  
6 E: B6 d0 y+ ~$ k  
   time,tm  
   autots,0  
& a6 r! L2 K9 p1 q/ Y$ X   deltim, , , ,1  
   kbc,0  
7 k1 f* p9 K4 i8 X' K$ h: L3 U  x  
) m, F8 \2 E" k0 v, y   allsel,all  
" U9 S" A( D9 C2 ]1 Z& K   TSum1=0
) q5 }, Q! ^5 D( b   TSum2=0
# `+ J2 \0 I" }5 B# K% n/ A* h5 c   Num=0
& E( Y5 ]) i8 A+ t4 D$ P  
5 L+ M0 E4 b( s   *if,tm,LE,Pulse_Duration+1e-12,then  
% N! Q* E1 X  Z4 [1 x$ m# Y        
( G9 M+ D+ q8 P5 Y6 m8 y3 v      !Apply the heat generation in each layer  
) }& R2 S- [3 o" {& ?/ ]- w      Num=1
      TSum1=TSum1+THICKNESS(Num)
3 z9 s) m2 e5 N' y: f. S      nsel,s,LOC,y,-TSum1,-TSum2   
      BF,all,HGEN,%Layer1%
8 z/ U8 i* H- Q; O& o# D      TSum2=TSum2+Thickness(Num)
+ u+ A6 ~$ ?# p      
      Num=2
/ P" R% s* h$ Q4 l- B+ J) w! [; y      TSum1=TSum1+THICKNESS(Num)
3 u5 A! {6 K4 k      nsel,s,LOC,y,-TSum1,-TSum2   
% V+ V+ b. e) x5 R      BF,all,HGEN,%Layer2%
1 M% V: h& M3 m; ~" u) q      TSum2=TSum2+Thickness(Num)
   
      Num=3   
      TSum1=TSum1+THICKNESS(Num)
/ S' Z7 H! I2 O" ^  P3 r2 Z. n" D, W      nsel,s,LOC,y,-TSum1,-TSum2   
' W3 i6 q& K+ ^' o1 a      BF,all,HGEN,%Layer3%
      TSum2=TSum2+Thickness(Num)
   
: r7 m6 Y# R- }3 Q& B' l6 L      Num=4
, {" a8 `( Q% d5 [  c* n' w9 P& D      TSum1=TSum1+THICKNESS(Num)
      nsel,s,LOC,y,-TSum1,-TSum2   
" J! O; M# y# m. B' G; D      BF,all,HGEN,%Layer4%
      TSum2=TSum2+Thickness(Num)
   
6 q* j: n( S! X) ?, y5 S2 l      Num=5
      TSum1=TSum1+THICKNESS(Num)
      nsel,s,LOC,y,-TSum1,-TSum2   
- ^% W9 Z0 K) c# f      BF,all,HGEN,%Layer5%
      TSum2=TSum2+Thickness(Num)
' J+ l" a. C+ f; D  
4 y" Z% Y. `" b4 u! k   *else
5 F5 A, @5 M% P% v      allsel,all  
      BFDELE, ALL,HGEN
   *endif  
% d  |7 m/ e; W# V8 f0 {7 K  
   allsel,all
   asel,s,LOC,y,-1e-9,0
5 h( i3 _+ k7 |9 U: ~# `$ q!   DA,ALL,TEMP,T0  
   allsel,all
( a" j8 S( M! c! L& ^  
* O# m4 ^, D: ~8 A7 ?   /status,solu  
   solve  
*enddo  
  
finish  
/POST1
  
' w. n* F: o4 C" c!*ULIB, Post,mac
! t! y! N/ c+ y% V& ?!*USE,TTV    !Temperature decay at each layer thickness
!*USE,TRV,1  !Normal Temperature profile at each layer at given time  
6 Z2 S9 m# G$ r* ~!*USE,TRL,1  !Plane Temperature profile at each layer at given time  
6 B  d. W: f' T- R  
!Laser power absorption in each layer
3 z# ]4 L. 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)      
  
+ ^# x0 I( |7 C% T4 u* V; Z//UG.mac
9 n! }9 @6 i0 n- P7 Y& r- `  
$ A) V' I6 s: X!FileName=ARG1
!Power=ARG2
!r=ARG3
!V=ARG4
; x- Y# }+ A& r: I4 j4 t!A0=ARG5
!A1=ARG6
!A2=ARG7
!A3=ARG8
!A4=ARG9
3 Z) n$ o/ ~- T2 |  
6 h5 I" J8 B  g; H: i*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNC_C1
*DEL,_FNC_C2
% h! x, e& }2 u4 h*DEL,_FNC_C3
$ C! g7 c3 ^2 r3 _( K*DEL,_FNC_C4
6 w/ ~( i- G& V! L) f- r*DEL,_FNC_C5
*DEL,_FNC_C6
+ b0 I' b  u! n7 H' U6 o2 K*DEL,_FNC_C7
*DEL,_FNC_C8
3 Q8 q2 I: E0 s3 \1 K*SET,_FNCNAME,ARG1   
3 _- y% |6 S/ f! H2 q7 V4 O* o. O*DIM,_FNC_C1,,1  
  ^% f, ?' v  i% z- Z*DIM,_FNC_C2,,1  
$ d# H# x; j8 c# P/ J, \/ h  A*DIM,_FNC_C3,,1  
. w# k; n6 N/ P0 n1 r+ \) ]# d1 ], e*DIM,_FNC_C4,,1  
*DIM,_FNC_C5,,1  
% E) o! l5 g1 c*DIM,_FNC_C6,,1  
$ \9 C- o0 |+ d/ ^*DIM,_FNC_C7,,1  
*DIM,_FNC_C8,,1  
*SET,_FNC_C1(1),ARG2   
*SET,_FNC_C2(1),ARG3
+ R7 S: ?6 Q. q- k9 O" ~6 H*SET,_FNC_C3(1),ARG4
*SET,_FNC_C4(1),ARG5
: D: X! q5 C4 b1 v# T& Q: R- V6 S' p*SET,_FNC_C5(1),ARG6
  [( \% w# P# g*SET,_FNC_C6(1),ARG7
1 t/ a* {2 n8 u0 a  K*SET,_FNC_C7(1),ARG8
*SET,_FNC_C8(1),ARG9
# m& S8 W, }# r# v& o/ T- p9 v! /INPUT,HG.func
" k& N- I; `  U: X7 b*DIM,%_FNCNAME%,TABLE,10,61,1   
!   
" {7 Y( k# D/ T2 v0 `) P& ?! Begin 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)   
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
! e/ q; N3 V( r4 E- t*SET,%_FNCNAME%(2,0,1), 0.0  
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
& \9 q* B: X# l4 O& T- W0 l- r*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
4 L& ~2 m: J9 a/ O  Q*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
+ J7 i! @/ W8 q. U% W' @*SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)%
5 i8 h7 ?; d( f2 P0 g*SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)%
' b% V0 Q5 t* a& I: s*SET,%_FNCNAME%(8,0,1), %_FNC_C6(1)%
  s/ K' |; y$ B7 p! [*SET,%_FNCNAME%(9,0,1), %_FNC_C7(1)%
*SET,%_FNCNAME%(10,0,1), %_FNC_C8(1)%   
" T0 P( B3 C( a0 c% N4 Y& S*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
2 o! j- {% C. h; ?! ?*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  
*SET,%_FNCNAME%(0,8,1), 0.0, -4, 0, 1, -1, 2, -2
) u8 G! Y5 h4 O! S, x0 Y% K2 \*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 2, 0, 0, 2   
' n, @) B% I5 @% A8 E+ \*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
*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
: \+ f3 Z+ |4 F' H6 N/ x*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   
" Z; ^3 o9 s, |  W4 g/ t*SET,%_FNCNAME%(0,16,1), 0.0, -2, 0, 1, -4, 3, -1   
; v" E8 r* m5 J5 u*SET,%_FNCNAME%(0,17,1), 0.0, -1, 0, 2, 0, 0, 18
/ y$ V7 ^+ |7 b5 O5 ?( G8 Z*SET,%_FNCNAME%(0,18,1), 0.0, -4, 0, 1, 18, 17, -1   
2 R: n. S) p+ H+ w2 ?6 v: o0 |*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
6 w+ z. Q. x, H*SET,%_FNCNAME%(0,21,1), 0.0, -2, 0, 1, -3, 3, -1   
3 n: l; S! G3 l4 t4 O0 r6 O*SET,%_FNCNAME%(0,22,1), 0.0, -1, 0, 0, 0, 0, 0  
4 ~1 `- L4 q% N2 I& M2 ?5 x*SET,%_FNCNAME%(0,23,1), 0.0, -3, 0, 1, 0, 0, -1
  A9 `: F1 t# E" O( j1 q*SET,%_FNCNAME%(0,24,1), 0.0, -4, 0, 1, -1, 2, -3   
1 i4 q' V/ c; v, `" x*SET,%_FNCNAME%(0,25,1), 0.0, -1, 0, 1, -4, 3, 3
) O' M7 B# J; v. ~*SET,%_FNCNAME%(0,26,1), 0.0, -3, 0, 1e-9, 0, 0, -1  
) i: Y: o5 ~9 I) n3 `*SET,%_FNCNAME%(0,27,1), 0.0, -4, 0, 1, -1, 4, -3   
  D- T9 K. a  R! S0 u! g7 D*SET,%_FNCNAME%(0,28,1), 0.0, -1, 0, 1, 21, 3, -4   
% O' v) F3 J+ u  {, q- T6 c*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  
9 ^6 F8 u9 g+ z0 _*SET,%_FNCNAME%(0,31,1), 0.0, -4, 0, 1, 0, 0, -1
" H& X6 d" |6 @) X*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  
*SET,%_FNCNAME%(0,35,1), 0.0, -5, 0, 1, -1, 4, -4   
*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   
1 A8 [7 x  j6 D& {*SET,%_FNCNAME%(0,39,1), 0.0, -4, 0, 1, -3, 1, -1   
4 T8 H! ]9 h; b) h9 w: l( B2 p*SET,%_FNCNAME%(0,40,1), 0.0, -1, 0, 0, 0, 0, 0  
& g" i7 ]0 y! @7 X: z! Z6 |*SET,%_FNCNAME%(0,41,1), 0.0, -3, 0, 1, 0, 0, -1
  [! V. H5 g# ?: l$ t0 K6 B1 D*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 o, g8 M( x" O, z+ U6 I6 m*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   
! K' k# `, a  q3 k- Y4 G*SET,%_FNCNAME%(0,46,1), 0.0, -1, 0, 3, 0, 0, -5
4 w2 P; W5 [( q  g+ I*SET,%_FNCNAME%(0,47,1), 0.0, -3, 0, 1, -5, 17, -1   
" v/ f8 S* |0 O7 o*SET,%_FNCNAME%(0,48,1), 0.0, -1, 0, 1, 23, 3, -3   
; D6 t, Y4 Y5 _*SET,%_FNCNAME%(0,49,1), 0.0, -3, 0, 1, -4, 1, -1   
2 U, [' z9 |* o*SET,%_FNCNAME%(0,50,1), 0.0, -1, 0, 0, 0, 0, 0  
( v0 P8 a$ `/ d+ R' v& h*SET,%_FNCNAME%(0,51,1), 0.0, -4, 0, 1, 0, 0, -1
5 v$ D* a/ L# z) |0 p*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
5 C# g; r: l5 O3 [4 @1 B*SET,%_FNCNAME%(0,54,1), 0.0, -4, 0, 1e-9, 0, 0, -1  
1 _% V' x7 k+ q# F/ L*SET,%_FNCNAME%(0,55,1), 0.0, -5, 0, 1, -1, 4, -4   
  G! v! r: @! \, H5 t8 |*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   
0 i" P! N$ P) e  F- i) \% C*SET,%_FNCNAME%(0,58,1), 0.0, -1, 0, 1, 24, 3, -4   
. m! s- t( g* ?; p*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
2 l8 k3 w/ i3 b: R9 j& ]! End of equation: Power/r^2/{PI}*exp(-({X}^2+({Z}-V*{TIME})^2)/r^2)*(A0+A1*(-   
6 T# ?6 y2 V  o3 I! {Y}/1e-9)+A2*(-{Y}/1e-9)^2+A3*(-{Y}/1e-9)^3+ A4*(-{Y}/1e-9)^4)
) T: u- X. X4 q4 ?!-->