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

xy112

高斯激光热源下的n层材料瞬态热传导方程
% _' _! s5 t+ A# Y& V, ?  M高斯, 激光, 热源, 方程, 热传导高斯, 激光, 热源, 方程, 热传导
出于需要， 写了个 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).
  
* k; C# Q, S" @1 N. T9 Y2 Y高斯分布体热源在程序中由HG.mac加载。不同层的a0,a1,a2,a3,a4决定了其内部的热源分布。
7 x9 G' U" ~% F: D3 y4 o  
$ m  E0 ?3 h* ^/ S2 M7 c本程序可已根据问题选择所需的材料层数，所有参数都可以改变。  
  
以下是命令流，如果有不合适的地方，欢迎大家提出、讨论！
  
FINISH  ! Make sure we are at BEGIN level   
7 z, @, P8 l( p2 n3 O: ~& p/Clear
/FILNAME,Media,0  
  
' H/ ^7 e; f+ b; t* S# A/PREP7
; o/ y3 \6 u% v/ W- E3 c, S# `/PNUM,VOLU,1
  
6 _4 B) c) C* ]! In simulation, the X-Z is the media plane, while Y is the depth direction
, D/ [6 V* T! G0 U  h- E  
Layer=5                       !Total layer number
: i, }. y" a, e+ a0 Rnm= 1e-9                    !length unit
Power=0.005               ! Laser incident power   W
6 X6 o, D& B5 X: f/ r* Nr=5e-7                        ! Laser spot radius      m
V=0                           ! Laser moving veolcity, along z direction,  m/s
  
, W' ^5 u! c( Q3 j3 Tpw=Power*1000                 !For display Title
6 p$ z7 t% g+ M9 L8 E: Qrh=r/1e-6                           !For display Title
  
5 _# V6 s& s# V0 T- L9 i/TITLE,%Layer% layer media, Laser Power=%pw%mW, Heating radius=%rh%um, Velocity=%V%m/s  
4 ]) ^3 o' A8 ^' S& k) h  l  
, B. W- |5 Q# F  t5 o) C2 w*DIM,Thickness,ARRAY,Layer    ! Array input for each layer thickness
*DIM,k_L,ARRAY,Layer          ! Array input for each layer lateral thermal conductivity  
/ ]0 \/ o9 T* X" L6 m7 m% r*DIM,k_N,ARRAY,Layer          ! Array input for each layer normal thermal conductivity
*DIM,Cp,ARRAY,Layer           ! Array input for each layer heat capacity
: i% K* x5 m9 K6 d*DIM,Den,ARRAY,Layer          ! Array input for each layer density
*DIM,div,ARRAY,Layer          ! Array input for each layer thickness division
* a' S5 y6 |, s4 X. g) v6 v*DIM,dscale,ARRAY,Layer       ! Array input for each layer thickness division scale
*DIM,A0,ARRAY,Layer           ! Array input for each layer a0 coefficient for heat generation curve
5 M' ~# e0 {- a( ~1 x8 R' ^*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
% u9 N- S5 r7 O: c% i2 w8 g*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
0 Q2 b  O! n3 I( Y5 F, EEnd_Time=2e-8                 !Final time for simulation
4 k  e" O" s5 i4 YStep=20                       !Time interval  
# y7 ^/ e% \7 u' L, [& _T0=50                         !Environment temperature   
$ W1 F3 J' M# z, x- Z  
* }; b8 U2 v/ s' ^! n( ]9 FLength=4000*nm                !Length of the block in z direction
Width=2000*nm                 !Width of the block in x direction
  B4 Y: Q1 x% j  L, m  
!Layer properties, from top to bottom  
  
!First layer thermal property
Num=1                        
8 J1 `9 S  }% O, P2 _/ _1 Jk_L(Num)=1.4                  
k_N(Num)=1.4
7 A* j7 W# Y1 f$ d/ K; ]  }Den(Num)=4000
Cp(Num)=500
2 _. ^2 D; }' m0 @; YThickness(Num)=50*nm
$ ]. V/ m' I9 y2 E& b: wdiv(Num)=5
& b1 Y6 B9 c; Q& Rdscale(Num)=1                 !last/first in -Y direction
A0(Num)=0
2 _9 q6 t6 F, W0 iA1(Num)=0
A2(Num)=0
A3(Num)=0
A4(Num)=0
. `8 o1 E/ ]* E( \/ c, ]  
& C4 l! Y- D( [4 u: u0 z!Second layer thermal property
- Y3 ^! m0 d* D# A3 T! O( A7 Y1 P% INum=2
* A/ Y, D+ J+ N- Sk_L(Num)=16
: c" P3 v1 i+ i! S/ uk_N(Num)=16
Den(Num)=3000
; h$ {2 q/ P1 u/ p4 jCp(Num)=900
Thickness(Num)=20*nm
div(Num)=2
3 s; E1 ]. q& ^. [# H8 T# M2 v& ddscale(Num)=1
) Q2 D6 V; x3 \. F8 |! q: ^6 @A0(Num)=807.929E+6
A1(Num)=-43.896E+6
8 W! S: T" ?: |; q& B- {A2(Num)=934.529E+3
A3(Num)=-8.864E+3
( Y" V+ `& N6 u0 v4 lA4(Num)=31.703
( v0 B( z4 n1 ~6 p  o  v+ `9 d) w+ p  
!Third layer thermal property
9 `; E) R3 u8 M5 s! MNum=3
k_L(Num)=1.5
k_N(Num)=1.5
7 {1 ]/ H+ P2 M& ^% I$ MDen(Num)=4000
Cp(Num)=500
  _( ^" b! Q( P* j* J/ j& kThickness(Num)=100*nm
div(Num)=5
dscale(Num)=1
A0(Num)=0
A1(Num)=0
A2(Num)=0
- t$ q( M2 V8 o; b+ vA3(Num)=0
A4(Num)=0
  
!Fourth layer thermal property
Num=4
' Q' M4 Y6 S# p' P3 tk_L(Num)=100
k_N(Num)=100
Den(Num)=3900
: o. W& ~% t& gCp(Num)=700
. X9 [* k! U* A1 M9 kThickness(Num)=50*nm
div(Num)=2
) v1 a! _0 @4 v' H# j) a$ e# `# Tdscale(Num)=2
A0(Num)=1.627E+9
A1(Num)=-32.146E+6
A2(Num)=238.195E+3
A3(Num)=-784.173
5 @$ Y* W4 o2 V; v5 eA4(Num)=967.716E-3
  
' ]/ @: ?5 S: \" i* X) m!Fifth layer thermal property
Num=5
/ ^+ Y  H# g& ^1 vk_L(Num)=0.2
( A- i1 ?' v# W/ I$ [3 P1 ?k_N(Num)=0.2
Den(Num)=3900
- \, m' M- l& S  R( aCp(Num)=450
Thickness(Num)=100*nm
4 J+ {+ g; ~" ~  C, fdiv(Num)=5
dscale(Num)=2
A0(Num)=0
A1(Num)=0
A2(Num)=0
A3(Num)=0
A4(Num)=0
4 @+ t: o6 C" i: _5 H5 E, `  
! Define the thermal element
7 ?3 u1 ^- N" p. |5 X# l0 o: jET,1,SOLID70  
KEYOPT,1,2,0  
* K# l$ M8 `, E0 }, _KEYOPT,1,4,0  
KEYOPT,1,7,0  
4 w! p- M* n5 T: r, }) \# nKEYOPT,1,8,0  
  
!For post-precessing
# c. R. [$ u5 K# ]+ B: L4 \nPoint=200
display_r=1e-6
  
9 N6 {! ?! K; J. lTSum1=0
TSum2=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
  UIMP,Num,C, , ,Cp(Num),                                            !Define heat capacity of material Num   
( @# o/ G" }4 O/ E2 b7 Z! ]  TSum1=TSum1+Thickness(Num)  
  BLOCK,-Width/2,Width/2,-TSum1,-TSum2,-Length/2,Length/2,           !Draw the layer of Num  
  TSum2=TSum2+Thickness(Num)
1 k) ]5 V4 y) z* Z$ T! J*Enddo
  
d_total=TSum1                 !Total thickness of the multilayer                                                  
: f9 _' @+ H  |8 S  
VLSCALE,ALL, , , scale, scale, scale,  ,1,1,
8 }9 p- v8 r+ d% x" Q6 A, @  
*IF,Layer,GT,1,THEN
   VGLUE,ALL
*ENDIF
  
: M8 d% k! f6 B, H+ FVLSCALE,ALL, , , 1/scale, 1/scale, 1/scale, ,1,1,
+ a2 e3 a& n& v, L$ d$ V) [: z  
: x. F  R6 n% i6 ETSum1=0
TSum2=0
  
  H" ?4 |: ]0 f( I8 T7 \Size=100*nm                             !Area mesh size, in X-Z plane
& a' |  C6 v- T; p% ]7 m- qaesize,all,Size,  
, A( E5 m0 Q6 z3 @. `$ c5 O& n  
( p9 _8 d  {  Y' |' U*do,Num,1,Layer,1
9 o2 V% n3 }  Z% i# a   TSum1=TSum1+THICKNESS(Num)
5 N; [. V5 e( s* j   lsel,all
0 C$ O% }0 m. k. c1 k  
   !Determine the mesh scale and divison for each layer at Y direction
1 }! {& {( J  b   *do,i,-1,1,2
      *do, j,-1,1,2
* x5 M0 E, V, u3 {  s         lsel,s,loc,y,-TSum1,-TSum2
         lsel,r,loc,x,i*Width/2,i*Width/2
1 u( v" ^9 k! z9 t2 ~& P- {. G         lsel,r,loc,z,j*Length/2,j*Length/2
         *GET,Line_Num,LINE, 0, NUM,MIN,
/ O2 Q( ?; K2 D+ u3 M6 v         *GET,KP1,LINE,Line_Num,KP,1,
         *GET,KP2,LINE,Line_Num,KP,2,
, D9 g1 c  ~- R         *GET,KP1_Y,KP,KP1, LOC, Y,  
         *GET,KP2_Y,KP,KP2, LOC, Y,
9 }& ?4 @) S  O" x0 G* ]8 H         *IF,KP1_Y,GT,KP2_Y,THEN
             lesize,all,THICKNESS(Num)/div(Num), , ,dscale(Num),1, , ,1,
         *ELSE
: _& \6 z1 Q$ b. d             lesize,all,THICKNESS(Num)/div(Num), , ,1/dscale(Num),1, , ,1,
# |1 q, J2 h" k% k: a- i4 C+ C         *ENDIF
9 N. v  i, B; a$ u      *ENDDO
, t6 V/ v' a; Y0 i/ i   *ENDDO   
' {, t5 q' `9 v! m0 a$ \  
0 r" g- ^# g* z   asel,all
* T. p4 {/ j. p5 m0 ]! t: P  
4 d" o$ l( L  I   esize,Size,0,
  E& G0 F! Z  h  
+ V# w& m( y- N  y# b% O- [   VSEL,S,LOC,Y,-TSum1,-TSum2
  
   TYPE,1
" z: c2 n+ G3 P) n" C# \2 s) B   MAT,Num
   mshkey,1  
  v' \% X2 h- C( O! A+ L$ F   mshape,0,3d  
4 [9 \1 B4 {; w+ l0 v: ~   vmesh,ALL             !Mesh the layer Num
) ~0 {& _, t9 t& ^  
   allsel,all
  @' m: i3 J+ O0 [. m5 \; a% D9 X  
& n- P2 z" t# s$ R   TSum2=TSum2+Thickness(Num)
  
*enddo
  
/REP,FAST  
( `% }+ i) \7 R, o7 s$ K/psf,hflu,,2,0  
/PNUM,LINE,1  
2 w" g1 U$ d# P' pfinish  
' r9 U* K% p& S7 N* J/solu  
  
antype,trans             !Transient simulation
tunif,T0                 !Ambient temperature
! f2 z% B1 U7 iTREF,T0
timint,on !time integration = on  
kbc,1  
OUTRES,ALL,ALL,          !Output all the result
  
/view, 1 ,1,1,1  
/ang, 1  
/ d% m% c" U5 m! N2 ]6 ?/rep,fast  
  
*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
) J2 M- `+ S6 b. `*enddo
0 F4 T0 F$ ~8 J( _) q  
!Start simulation
( b+ b; K: t* C*do,tm,1e-15,End_Time,End_Time/Step  
- [( f9 C/ t" `6 R4 F5 z  
   time,tm  
( M2 f0 L1 N) s* w$ O- z   autots,0  
   deltim, , , ,1  
, [. n* C! ?+ B% k   kbc,0  
: I9 D7 ~2 _/ B  E  
4 R/ n9 P- m1 s! r+ N   allsel,all  
- P. S; y3 l) z) P$ c   TSum1=0
+ h; Q$ [1 ?  u# U. \3 `* B% L   TSum2=0
   Num=0
, p8 F' ~! V: r9 [3 i  
( \6 t1 S9 R8 E% V; L1 E$ X6 C+ D   *if,tm,LE,Pulse_Duration+1e-12,then  
        
      !Apply the heat generation in each layer  
      Num=1
      TSum1=TSum1+THICKNESS(Num)
6 h+ c7 X6 M. D      nsel,s,LOC,y,-TSum1,-TSum2   
      BF,all,HGEN,%Layer1%
      TSum2=TSum2+Thickness(Num)
      
      Num=2
4 \/ f# l7 I: ?; ^      TSum1=TSum1+THICKNESS(Num)
' T& |3 S+ E0 [! f      nsel,s,LOC,y,-TSum1,-TSum2   
; Q: f2 v4 ?, \8 q      BF,all,HGEN,%Layer2%
# P* s5 B2 l0 ^1 ?* r  Z      TSum2=TSum2+Thickness(Num)
8 E( }+ A" U3 j, Q5 g   
      Num=3   
7 l% b& e/ m5 N  z3 X' z4 u      TSum1=TSum1+THICKNESS(Num)
      nsel,s,LOC,y,-TSum1,-TSum2   
      BF,all,HGEN,%Layer3%
6 Y, M' W$ V  u$ B      TSum2=TSum2+Thickness(Num)
   
$ D+ l! A8 M* e& o8 r- R7 X- j      Num=4
      TSum1=TSum1+THICKNESS(Num)
      nsel,s,LOC,y,-TSum1,-TSum2   
      BF,all,HGEN,%Layer4%
. f% j5 w5 p# _) m9 o      TSum2=TSum2+Thickness(Num)
' u# |1 h( {% C4 o   
      Num=5
$ A# _6 y( }% R      TSum1=TSum1+THICKNESS(Num)
      nsel,s,LOC,y,-TSum1,-TSum2   
      BF,all,HGEN,%Layer5%
8 n; B( d0 D% b( j; `, a      TSum2=TSum2+Thickness(Num)
  
   *else
      allsel,all  
, @0 n4 c# I4 e8 {      BFDELE, ALL,HGEN
6 Z$ F, B2 A  j   *endif  
5 l" G2 m/ y: b  
3 s9 a9 W, U# K2 e6 Q3 _0 ]) r/ h3 |   allsel,all
! p! W& p9 B/ R/ [7 u   asel,s,LOC,y,-1e-9,0
+ f) o7 Q) T7 H" H!   DA,ALL,TEMP,T0  
6 G4 k. R6 P' f0 N7 {8 c   allsel,all
  
   /status,solu  
   solve  
7 A  R" d3 t6 K  y* d*enddo  
  
finish  
) H& y; [0 b, k, P/POST1
  
!*ULIB, Post,mac
) C! Z/ g! e$ U- w, R7 b, M!*USE,TTV    !Temperature decay at each layer thickness
- U* L/ Y) s/ z+ w1 g0 J!*USE,TRV,1  !Normal Temperature profile at each layer at given time  
" r3 H' `% Q0 i  z!*USE,TRL,1  !Plane Temperature profile at each layer at given time  
  
! p3 _" j: M1 B/ o) c! p!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)      
! V8 X0 H3 h# `, N/ ]  {$ M; o  
! F8 d; }  V3 g% _+ r//UG.mac
' f8 S' u, z' t$ V5 Y7 l  
!FileName=ARG1
' l; D2 R4 F  n4 T!Power=ARG2
!r=ARG3
!V=ARG4
5 h" W, b% @9 }!A0=ARG5
!A1=ARG6
!A2=ARG7
!A3=ARG8
!A4=ARG9
' P. j% u+ r6 @4 L  
6 Y5 a# r/ O4 P' H; @- N) F1 ?*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNC_C1
6 F  X+ C8 w: I4 z7 X. `; L/ Z*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNC_C4
% b7 X4 v$ K# B*DEL,_FNC_C5
% ?& t+ R* O. v5 d*DEL,_FNC_C6
  A# D/ q0 n( p4 I*DEL,_FNC_C7
*DEL,_FNC_C8
$ U1 O. p! V) U0 g6 g& i*SET,_FNCNAME,ARG1   
" ~+ X, {' v. d' I. I*DIM,_FNC_C1,,1  
*DIM,_FNC_C2,,1  
*DIM,_FNC_C3,,1  
*DIM,_FNC_C4,,1  
- ^7 u9 [) A+ Z' s9 `0 f*DIM,_FNC_C5,,1  
9 K4 J+ |5 B6 L0 P2 r# v*DIM,_FNC_C6,,1  
+ ]! Y, t( _" `*DIM,_FNC_C7,,1  
*DIM,_FNC_C8,,1  
% Z9 H; T$ R4 o0 a. K*SET,_FNC_C1(1),ARG2   
- V5 j" L3 q+ X3 k/ |) J0 W  y*SET,_FNC_C2(1),ARG3
7 v4 I9 R# O, s$ {+ a. Z2 {*SET,_FNC_C3(1),ARG4
*SET,_FNC_C4(1),ARG5
  J" x8 M+ M1 z*SET,_FNC_C5(1),ARG6
$ e5 [# F2 F' W0 F0 Q2 r' K*SET,_FNC_C6(1),ARG7
* C8 C+ z7 s& ~+ J/ Q' ?6 E*SET,_FNC_C7(1),ARG8
*SET,_FNC_C8(1),ARG9
! /INPUT,HG.func
*DIM,%_FNCNAME%,TABLE,10,61,1   
!   
! Begin of equation: Power/r^2/{PI}*exp(-({X}^2+({Z}-V*{TIME})^2)/r^2)*(A0+A1*   
) K# ]" W  P/ K3 V' z! (-{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   
*SET,%_FNCNAME%(2,0,1), 0.0  
4 L9 f& m# D0 v# a*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
- p: F5 S$ A5 o% k! a*SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)%
- M2 L) I4 G+ z8 _& t*SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)%
*SET,%_FNCNAME%(8,0,1), %_FNC_C6(1)%
0 @. O' }0 M( {' D. V3 ~1 x*SET,%_FNCNAME%(9,0,1), %_FNC_C7(1)%
*SET,%_FNCNAME%(10,0,1), %_FNC_C8(1)%   
# o4 @% G  [! F*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   
) T- W3 {2 [, z*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   
+ s& T) C: r1 F( `3 @. @1 n. B2 n*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   
/ b1 S7 J4 Q+ W# T*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
8 K3 z' e3 ?4 B# B*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 2, 0, 0, 2   
4 t  S! x6 b4 u* a: o*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
" O3 N! d0 q0 Y7 ^*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   
  H0 b2 A& J( G2 {$ k*SET,%_FNCNAME%(0,15,1), 0.0, -1, 0, 1, -2, 1, -6   
/ P7 q" H6 D# t) Y4 h: N  X*SET,%_FNCNAME%(0,16,1), 0.0, -2, 0, 1, -4, 3, -1   
" M2 n7 ~8 r( H9 S6 N' j*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   
/ p( R' j- ~' @1 D9 S*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   
6 X# z% r( o# Q8 P7 j, C, P*SET,%_FNCNAME%(0,22,1), 0.0, -1, 0, 0, 0, 0, 0  
0 ?6 Q$ [$ \2 K3 _) Y*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
' R# Z& ?2 v1 Y9 N- c*SET,%_FNCNAME%(0,26,1), 0.0, -3, 0, 1e-9, 0, 0, -1  
*SET,%_FNCNAME%(0,27,1), 0.0, -4, 0, 1, -1, 4, -3   
. i  Y" d: ~( 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   
*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
. H8 T; A% ?" O+ F3 R*SET,%_FNCNAME%(0,32,1), 0.0, -5, 0, 1, -1, 2, -4   
7 `# E& g2 Z& v( O& b*SET,%_FNCNAME%(0,33,1), 0.0, -1, 0, 1, -5, 3, 3
" D+ R( Z# H% r6 E6 h*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   
* \9 ~% y, w- l  A*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   
1 B4 p7 L8 F) F  t*SET,%_FNCNAME%(0,38,1), 0.0, -1, 0, 1, 22, 3, -4   
3 F' d' v3 @: r) n) n2 m# o*SET,%_FNCNAME%(0,39,1), 0.0, -4, 0, 1, -3, 1, -1   
# i& x3 `- q8 t' z+ T*SET,%_FNCNAME%(0,40,1), 0.0, -1, 0, 0, 0, 0, 0  
# |+ i/ @" a9 G7 A  S6 Y1 R& L* c! t*SET,%_FNCNAME%(0,41,1), 0.0, -3, 0, 1, 0, 0, -1
. E* R; f$ _, C0 x, J0 G. U*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
*SET,%_FNCNAME%(0,44,1), 0.0, -3, 0, 1e-9, 0, 0, -1  
4 T$ U8 c, @8 x. M*SET,%_FNCNAME%(0,45,1), 0.0, -5, 0, 1, -1, 4, -3   
. Z, `; I7 ]: I/ Q1 o& W  I$ p*SET,%_FNCNAME%(0,46,1), 0.0, -1, 0, 3, 0, 0, -5
/ h- n4 n: L8 e, y' `3 X( }: Y*SET,%_FNCNAME%(0,47,1), 0.0, -3, 0, 1, -5, 17, -1   
$ j* g* e' a8 j1 [* l& Z( e* V" g8 v* O*SET,%_FNCNAME%(0,48,1), 0.0, -1, 0, 1, 23, 3, -3   
3 ]& n; ?  i8 R& b*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  
; u% V8 n0 x( g. A*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   
1 r! o: \  w5 P5 K7 R1 v*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
*SET,%_FNCNAME%(0,57,1), 0.0, -4, 0, 1, -5, 17, -1   
0 ]) F8 W8 }* Q7 [- }6 G*SET,%_FNCNAME%(0,58,1), 0.0, -1, 0, 1, 24, 3, -4   
  d/ q8 U/ }- D9 E* y9 w*SET,%_FNCNAME%(0,59,1), 0.0, -4, 0, 1, -3, 1, -1   
1 u8 _) V; x7 O# D6 r8 m- a*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
0 G5 R3 \( Z; L) ~! End of equation: Power/r^2/{PI}*exp(-({X}^2+({Z}-V*{TIME})^2)/r^2)*(A0+A1*(-   
( u6 J! q9 X3 _; z! {Y}/1e-9)+A2*(-{Y}/1e-9)^2+A3*(-{Y}/1e-9)^3+ A4*(-{Y}/1e-9)^4)
!-->