sr!com.femlab.server.ModelFileHeaderD%LtagstLcom/femlab/util/FlStringList;Ltypesq~LvrsntLcom/femlab/util/FlVersion;xpwsrcom.femlab.util.FlVersion%/B = IbuildImajorLdatetLjava/lang/String;Lextq~Lnameq~Lrcsq~L reactionExtq~L reactionNameq~L scriptExtq~L scriptNameq~xpwtCOMSOL Script 1.1tat#COMSOL Reaction Engineering Lab 1.3tt COMSOL 3.3q~wt $Name: $t$Date: 2007/02/02 19:05:58 $xur[Ljava.lang.String;V{Gxpgt modelinfotxfemtguitg229tg227tg228tg222tg221tg220tg226tg225tg224tg223tg202tg203tg204tg216tg217tg218tg219tg211tg210tg213tg212tg215tg250tg214tg206tg205tg208tg207tg209tg249tg248tg247tg246tg245tg244tg243tg242tg241tg240tg238tg239tg235tg234tg237tg236tg231tg230tg233tg232tfem6q~q~q~q~q~q~q~q~q~q~q~q~q~q~ q~!q~"q~#q~$q~%q~&q~'q~)q~(q~*q~+q~,q~-q~.q~/q~0q~1q~2q~3q~4q~5q~6q~7q~8q~9q~:q~;q~q~?q~@q~Aq~Bq~Ct mfileinfouq~gq~t femstructt guistructtdrawq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Itgeomq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Iq~Exsrcom.femlab.api.client.ModelInfo^%Ldescrq~LdocURLq~[imaget[Bxpwur[BTxpOPNG  IHDRp IDATxY,Ivv"r_3=!)Y$Je-`mI ?A |0؀a昆LC~a[H%d#L]{\"'DvSSV~uavi@yvZ,|#4PX٣FfiE<:+ad.p#|(6Kx^ 7^DB7A#zSzӿW!a+_ݶbd4=jzᅙGxp/#cXC6[s|K3@]K2j," Z׷X,`:I> ZM2%- VM)YEk D)b\آHÔ`Ak:܈s8.~a>i/ɊmR^ӣd6ŧʙe n}҂`$xE(>"PyLҿWOTILOLB>Kѓ^ ڒ =. _oqnhNOz̟ƫ&CWB͒<(Mn%ɚٗ@fH!MUb2k0\ʞ0DL3v ؠ `3vBvٸA_ yrLA(?EP1(; dC$Uc7#' ٩1fta8*4aNX^SVľWNXp[&O |=hBOkRI v#6PUtF M wk&N;Jn&WE2ɵ u._V6 -ySzH(ҋmAy毯%Z/jg]ڜ+KL,UlVʚ.pX' {iN]ゕ?۳)%<|,NM[Ak3OY6eG: 8ހb7O.QM/kVFCp_R&σR$A \L Jh l^R:'- 0PUab!u1O"9Agc:$/;eeCdW{JF1romeUQ EWҴn. 1˞ +\C\=\v?f\Sm2n.:JoV6x@v1-tFP}+ԳI}/)_ECEf+?b`ePHd"lOQq%x8ġʐN)Ba]j213{P@&l%= f:W Eҫ4RF#(c,VY7dAnÌѧ8h%%8@;=<L!zDCiHqW`^wxyDچYoY~j܏bݷv*3p{k<@PR4qcyBن{^ncVA-ݣ3kOۋC{~ɬēr6 a1G!W|['~1AXBBetJMyIzeXJEXap Xy%Q u$ $TM}\efhV%B]Y\y$fשαUCY/x^YfIWTyt+zEL0[/SxU췤]UJVy~j*T̾|VR"A'Hp}M)p3v ˌK*Xz+63ڹޮR"()OYS_p͇yN:RPYB$p)3V q )[DB CHVqd b4wϟMIus-\18bR%Sj"ş(Վ{$r* /r3Oߓ"MCrYsM%7Ú{/8#H~Z/MRXnBl< RҤGX;kA@&HTs z&Ty^ JRzlʄD~"+5~20It(!"ah2BӍhhEMBLC{SNP3#4,#4O c ;Efkte-N$Sі[]A +M\N$e>LHZT_dI/80Tކ .c"ei FiI)H|"8h!TЯEutM1I∄j O)IYz]¡J賸޼F"oG/wsA0%Q( Z}w Hf c`'=qjwV%;>;L/|voSE܊c\&  3muhE y'iU#mk֩08orVlCF!+R_+fM8$"p`U} ﯭzͻ\5k4TyJF6fYňI6Yp % T"LD o k|ln!ω\\QGS2!(Q`,BA.s進_4egP'sʧk\u)q>`nO52݋߫]dQ ZSUTΛ;+ȖP?Fa,KD[&rnW@$ĬkG糃R2 0ڈ}C;*4Zy+֚nyfm4:: o"ߠj;Iv~Ea[WIK dŋ! 'eKj 5;Hw-xe+1fٰ@5D 9$^i_c`P L=_w :ߦN^#{_3b eWx eLr)ȥUb'wl3˚}e]׿˿ ۽nxh1>}*s}@(_)+ jq #V{ ɤ4d%fѥ;̔{W1W\B! LABGd.)8wh9.GY]545WُIhͰ&Kk gRc jo=b I=f/ƦXgLE}N2JDW(L(598> (;}-9(!8DPmD NHP'ͻ 6=VWqNt,b>do,i q(;xmo@ EQ|+__[A Pd\S`hA?sbYqrq,:t)jݺ>B݇N/~ՍkTp C<Vj, {#0sF4z}(jŷR6BjOGO#zE\^n8rP>D: wDž˲h9j戥dފ {3aԩt=s_Sgip=ocȢ\:Rq^7q9AO=Pj&!l! `KG3A}:|cHF sNZH>3J5%$fjVX.b&CX1U.N8(5BWC>Pj 8ϧ'q%0pTqy)#JMZ/gXP`04U_M3I@&;*-K꒕@z aA.bUs<аˎǬa!f}@y]Fb쌼#LGys+vS\7O LXg-D}C2v(z:L!t,;S%a=4 @C7Ԓ(R>}f[q]2N Cf^I8Oӯxt2mշ>h39K^gl.] %g%E35% )Q[` dv$mwsŊt__#vYMI)B(q\+XqP1"+*cH)7 u_>·kLT $ nPxr.qAS`PRR"" F"ℕ `;@/$tBQN.o%"ئgIXq&c^EN| <:ߍNV, ⵴tޙU{ZG]YpiQ jE?ma'eh];Ԁ'C!2es)e>7P.]X 6ed1oh_X"Ɗ8Sf`c1%}s$+&uUOLd%EDZTPq2.M'xZ&,NkAݝ~ҏa V#(8G4wm/7WKH0%zBq`A?]y66efe\+rkPO^::2VqT(Ե#bWd, e^9/AqEY&Z-Y098#9Q<߳:cOgskP%fs0v: UxqHG55__g/ܒ-e:wd,.j c,R"bՊ4#iLѫys _ItX,HX='uCO K,uOn 6c|7ߕ `Aa3e=I zeSbSJQK:^uwzNs(FԳĭR^fP"%hV415(;6T\1ʵvY\4+g)6m3=:=!);פhL`7sNkRu?- 0T@9%Pw^әa_x5|+IݼX418 ٓٔ2Vܼ%R2%'08L w6~ܝS|GHҘim6 vXY56¹xblIJJt!ow0~sYMO/S6cT֣eTq[{xFØgt[DaO=st:w19(oRj!MJQ]7ȭ)ktV {!u2 W@ܷA7=R߲+eA)!Qyh[cqkW$f2 2gġ$ÆH7_F9quCqw, ' uO4>%6eb7PޠF})i! xY. sTp57utตS…x| d= l};zGeFbgzkɯ$[}PV@z1ͳ^RΤ][( k J 'ou61Ov&1ɧSp#V(< .*(᫡)U@RP.f &:rЦ"Vu) fI\;dHLbTT#X|οQ oM|AF4k&LL 4FGT£BTp## 7qyyʟ֦̆ J$=JpH1{h'Epf剀/y@4;_jSq2a[stIu3 R$Ց)OTf_^#>oQHv [dP,=1 e7≮tn}G?I$zPp !o*&-i ,#^ŵ ǎhT#ܰn!@'w $IM2:tW) ?pp%(,}c K>v9TϺ6_ 2{iC:T"Ye\%+8s`j5S26@>Mh70O dl_z#E$(~@+)ryՕO~φ(2=vYbaŻ"p IDAT<x7qӋAvdHx LQ`ؐ@e*v9dFŽ&h\ `s$HP}|V/Jyo;J& ,g S` ]8Q}~^j.)53ñ`T2`LbB@`z9gmeABf5‡mX|K_eR~̃R~S LzIg=Áؚfu/.z8nS-S'k%x<RD|J!uaف3$z>(a6}k_;מ?l%.buAmw3 mԦAx2Ј5iΎ_;?RsrQsR4ҏx7g+6 L .#_cDitSCJڬ@YSe(,&/jYtkvw4 2 |d7qYcq i[3{gPޥ {]/Ɣ>$؃NH}&J@gƔoR@Q ;D J$*CA#Џp~~;$i=d Q:9uLq; <30m`Räq #pXY9]cd1i5TYVj``ipLYBDC !}KO[) ? C!k0h"0,e!`OkOy:“`svvv%yr4jkf0nN)@^,ktmDF+a},WKX`9<Eo;U b|F> %Кy3ulATOCؒp"2sj{45JmSI]KS $r1SQ~j~_η|&%f $h%E)%<ߴ;e8Vgu,EIKM<0UUtm/0̞nL;@ ##&FYh4QF1]m| ],uka6BHp%ӿ?Lx;ueJ“)HddIQɦ-ձ̇ޝ22Ubv5$Y8c[;?1|Ӛ^1 ݣ:)^1k ;(<33(j"V#^n|%[7Wi .+KYTE.ƫ &UkOX8= &h- 7፛0LL@3FaxlD0YVڃ|F\떬TS&aĴ-~(-?.7*#1ȕEZPH/#zhrkalAJ'r:-O *QfpaߦPGV`ɉTf"dQj?F|WgCH@=&I9&L>*\Hm2u{/e~T%cyu˟8P~^ɆDPJe{P.t.l ŭYi1#v8oy 9a•87r~HWLz#E^>VꚦġA#4wRV.df`zptf#{ΒGg.`>ǰ;1ق^AIJT()70\-%֗ʼnH75+zΤRДp/h7ԢT{z . J9n{n$V1Wg,`4n·P:P}< ,(E$Jupuhw* ="=(--(G B 5($F:qz 5)=93VRv2.b@TCAJ19 ~;ނVk&oCc/)95AjLQRZPfuBDoJ%N[];0(GEOڮ5* Viaw$:MӣpmPǗ7"Q)=)VgS| ۂ"{9RLR%)hnuLhu>i`o)qȉtKLXE8s6j+uLCД)M;C$,U 0+U}hxBPu~٤(7e%2J4*ڝ[RCEI#s*0"eMa)!K`Eay+C=ٔ-s6-Ay䪆pIiՕ@^}v%xH/M44Q2mOM{%,0R''`Ӊ'ӤWT>0,F0R h邲Coz2ƴ|x/BLa1ݫ )eo7SZy"GҿۚcOgmxwK.JusL1ʐ:'<>Y%d~i\UB,R2Djeuς9&r׆<,<^QSzp7[R`5]uhΏo@$S6c!6931ՍPZ!G('vɞfn?"2-,n=/(dWP4^ӽB?y<;xZ*QK~GER0y%l)vA"IҦHԄՠ6孊{|GΠI X eE00&c<.U#h.bl<)fLw'1sr]AĠ4|(յ@u!7nMHߤ뜈-(lxEx \gܿb2zAyI7(%Hn@z9w  mCws,5q`wlGG|4j(cN!L6\VNo#S+Q#Nwc8ɸ7մ{|8?_Iݑ(g8JEplHe`r69G䬏vLu桠s^iҹ& (=7CݣWCU2mhSM55CpBYFxG?;1y :G$rF[@=R1 k;0`A> 2gINQ ŕc}1Y .Va]jWdň݀ҏ2n 7]'ٞR;(8(95[AϾSJ0SX,+e3%vƙ\E>ٍ#gUֽbU[P&{4ML9:!;"Z,fuuk'Y[9C#6}k Ǜ *%,e+ %L6M앤42&wi8:}F\,Rk.GaIӡ7'g A pO$Sژܒ-@Tic!>ڠ&U)Jr)q)W b.b )Izpufcb'Io}6nƽ4B֎:Jc)&IUjiRrMhZ$CE?c9:x>'=y%3XV1B &lO۳jeJUryS_fYL[O܃2Gi3=lKFT9zN2*` SBI\v;;:zFJq1\ MtzᰍفF%BfJCY ba w SJ̖ڗ{ٔm6oH sٝ`罜do3dt@ɇ$egslAF@ {1zuwu쎍 #DX!enm@ )6a0CBቊ#/5+!ݰM0@;x<Օ;P}_2ؒͣ ]f&GPnŤwXfhn}٥qݰ#g=.[xͷ (VctY P3VDM]1^Mpb< SQCY$wc4CF܆DgS@/%1>jwҁ7/nS-W߃-h7N^@lFB-1z6"t]5uLSF-wxMyE+աG~`M R.2+,COvz:^q洩pQIl<2W ,c# |K?.7@Ρ{E5@krʌnDU{A/_Gp|2*%+kD/v?-geBKpBr ԺF;BVH _{-eՆ=1URG'SXn= _|x)X=e.D;M~~%믿~wO&?|[?xZ.rpL>iz[`>?_tϻ J#vDU/O`l%~XqwnGmhoBw'''M=\.g(жL8n:^Y% %d>Lq8|S'.é=xg~0 ]{ֺ;ō3ErG!N se2 \!<]s+Q( R$gEu|#<8h6]Tj9(C3[Fr1>/IAs7ࡅ)l9TX-d@v+]>ܲҺ/ooo+qOűZO̳h|0CLWI[tFI+7(agpm6sP(/3dqSlsyyoяܽs]=[ 9T,' qZR?NXY\HRJ;x ʇpw7̜? wUPpbO9ٽl ʽ"OK{uWT[PZ^]ypt$.Ey̼qj&'d6{B,UJC9<o]㣕zsh~yH[kof6,(/h^ ~SKx^M(ÈЫ Ɔ7d|z:b9r*Ea)F&_uɶ]{pW\m^×Ϭ u/țL8(;iF82t}d70wϞw,ށ ߲Pg"MhMՍVdy6]LzNE晁=5v|ڠm 1O*^#HWsVrfIM+|E(Il@Sks}hMIhry ~%۱CAJ&QZi3-|x~+/֜x ޾ kw_>$IDAT"M=Oa\]Ǐ\VR{3P,Dwӌ'uCuC!;cl#qf"{D3b=H7]r3{@@yE6{PRv@_A5c%<2ό?b6^o.)x57_AQ,/Yl}6}JXo bo0H#ҷޮqyR3DnIorс?gIEkR bBu NCφ7Dζ+<2NBa\ni=ӑUO}.e~nȈ*Vj yҏSD-R[ً\RZ2'B4')cR7g9q iM\؁' eFPNʖV0^$Όk^]|%)qب޴öݯ&NeT;n>PS'̚>T utܝ1wo#Y,\ST(!:޸gS!2DIR~Vq2O*ON5bL׮/дvo5oWWQ$bRǓ98KhJ΍5q55T|-Cz@(S|"g\N0IΑʤs99 U*la0AT4L}N"&N"ػ]ܠf\;b`-({/J+#[N&KwKkU%NI8g Z'.ʰԷuR5zKC -#D@9a51CW7zߦ'jbkp.2Xk<#|(e<.ŕ-wziڱ&e])NṊ ʂ3[FPbQu4bz#ZkvK=4αNl ^̣۵(<(+8eBwA^0úʲ!: .`V \!-4C*30-(M5#77%(9;DtЎO .=Е9a6>Ül @+gp9 l:Nl NBĉyPM( 7lȘ!)uo쑢[zʩo,J/zT/ˋXPE.FhR|mmeKBIMёxPKJ;dyAoߔ28a,zxw0HkCWX ?]%ܒ[B ?qWWO;R=x\Ky:~ʚPNGp 4Eh:vvu J3QXcM'.#:z =8q;6 m!N|`/xH1 )ssum9C8A7X4qRQSodG>梂qS˰=(ML'%XgQOdIJ?[-imJj1gPDJ$:UMwہ2o u<>vs@%e2pS0.5 >-?ڀ׏^3N.-SWxewqRl'eT[*ˡ.EW{;DϨx8M419? 79HJIʡojuhB 02TF&˞(w}fA]dv?kvqI .O'-mڙWEHw#hLB\ -AѠ I/ʫARiʼO3I҈~Tz<)Su8ϋ;oG3Dr?/Fa$l]w4#SW Cƒ=x/;C\T7?_-CjeW* Fڌ]W>D6 Xq'^?ȚQ<5@CMk3YGUUGym?zSy<% y8lmmva;lPva;(Qn?{'ҊIENDB`tpxuq~tbclear xfem clear vrsn vrsn.name = 'COMSOL 3.3'; vrsn.ext = 'a'; vrsn.major = 0; vrsn.build = 511; vrsn.rcs = '$Name: $'; vrsn.date = '$Date: 2007/02/02 19:05:58 $'; xfem.version = vrsn; xfem.id = 4; xfem.geomdata = 'geom'; xfem.eqvars = 'on'; xfem.cplbndeq = 'on'; xfem.cplbndsh = 'off'; xfem.drawvalid = 'on'; xfem.geomvalid = 'on'; xfem.solvalid = 'off'; xfem.linshape = 'off'; xfem.linshapetol = 0.1; xfem.meshtime = 't'; clear appl appl.mode.class = 'FlPDEW'; appl.mode.type = 'cartesian'; appl.dim = {'Hrad','Hazi','Haxi','Hrad_t','Hazi_t','Haxi_t'}; appl.sdim = {'x','y','z'}; appl.name = 'Axisymmetric_Dielectric_Resonator_Simulator'; appl.shape = {'shlag(2,''Hrad'')','shlag(2,''Hazi'')','shlag(2,''Haxi'')'}; appl.gporder = 4; appl.cporder = 2; appl.sshape = 2; appl.border = 'off'; appl.assignsuffix = '_Axisymmetric_Dielectric_Resonator_Simulator'; clear prop prop.elemdefault='Lag2'; prop.wave='off'; prop.frame='xy'; clear weakconstr weakconstr.value = 'off'; weakconstr.dim = {'lm1','lm2','lm3','lm4','lm5','lm6'}; prop.weakconstr = weakconstr; appl.prop = prop; clear pnt pnt.weak = {{'0';'0';'0'}}; pnt.dweak = {{'0';'0';'0'}}; pnt.constr = {{'0';'0';'0'}}; pnt.style = {{{'0'},{'0','0','255'}}}; pnt.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, ... 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1]; appl.pnt = pnt; clear bnd bnd.name = {'tangential_H','normal_H','null',''}; bnd.weak = {{'0';'0';'0'},{'0';'0';'0'},{'0';'0';'0'},{'0';'0';'0'}}; bnd.dweak = {{'0';'0';'0'},{'0';'0';'0'},{'0';'0';'0'},{'0';'0';'0'}}; bnd.constr = {{'Hrad*nx+Haxi*ny';'0';'0'},{'Hazi';'-Haxi*nx+Hrad*ny';'0'},{'0'; ... '0';'0'},{'-Hrad';'-Hazi';'-Haxi'}}; bnd.style = {{{'0'},{'0','0','0'},{'solid'}},{{'0'},{'0','0','255'},{'dashed'}}, ... {{'0'},{'0','255','0'},{'dashdot'}},{{'0'},{'255','0','255'},{'dashed'}}}; bnd.ind = [3,2,3,4,3,4,3,4,3,4,3,4,3,4,1,4,2,4,4,4,4,4,4,4,4,4,4,4,4,1, ... 4,2,4,4,4,4,4,4,4,4,4,4,4,4,1,4,2,4,4,4,4,4,4,4,4,4,4,4,4,1,4,2,4,4,4, ... 4,4,4,4,4,4,4,4,4,1,4,2,4,4,4,4,4,4,4,4,4,4,4,4,1,4,2,4,4,4,4,4,4,4,4, ... 4,4,4,4,1,1,1,1,1,1,1,1]; appl.bnd = bnd; clear equ equ.shape = {[1;2;3],[1;2;3],[1;2;3],[1;2;3],[1;2;3],[1;2;3],[1;2;3],[1;2; ... 3],[1;2;3],[1;2;3],[1;2;3]}; equ.gporder = {{1;1;1},{1;1;1},{1;1;1},{1;1;1},{1;1;1},{1;1;1},{1;1;1},{1;1; ... 1},{1;1;1},{1;1;1},{1;1;1}}; equ.cporder = {{1;1;1},{1;1;1},{1;1;1},{1;1;1},{1;1;1},{1;1;1},{1;1;1},{1;1; ... 1},{1;1;1},{1;1;1},{1;1;1}}; equ.init = {{'0';'0';'0';'0';'0';'0'},{'0';'0';'0';'0';'0';'0'},{'0';'0'; ... '0';'0';'0';'0'},{'0';'0';'0';'0';'0';'0'},{'0';'0';'0';'0';'0';'0'},{'0'; ... '0';'0';'0';'0';'0'},{'0';'0';'0';'0';'0';'0'},{'0';'0';'0';'0';'0';'0'}, ... {'0';'0';'0';'0';'0';'0'},{'0';'0';'0';'0';'0';'0'},{'0';'0';'0';'0';'0'; ... '0'}}; equ.usage = {1,1,1,1,1,1,1,1,1,1,1}; equ.name = {'dielectric_0:vacuum','uniax_diel_1','isotrop_diel_1','uniax_diel_2', ... 'uniax_diel_3','uniax_diel_4','uniax_diel_5','isotrop_diel_2','isotrop_diel_3', ... 'isotrop_diel_4','isotrop_diel_5'}; equ.dinit = {{'0';'0';'0';'0';'0';'0'},{'0';'0';'0';'0';'0';'0'},{'0';'0'; ... '0';'0';'0';'0'},{'0';'0';'0';'0';'0';'0'},{'0';'0';'0';'0';'0';'0'},{'0'; ... '0';'0';'0';'0';'0'},{'0';'0';'0';'0';'0';'0'},{'0';'0';'0';'0';'0';'0'}, ... {'0';'0';'0';'0';'0';'0'},{'0';'0';'0';'0';'0';'0'},{'0';'0';'0';'0';'0'; ... '0'}}; equ.weak = {{'((test(Hazi)*Hazi-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*(test(Hrad)*Hrad+test(Haxi)*Haxi))/x+(test(Hazix)*(Hazi-M*Hrad)+Hazix*(test(Hazi)-M*test(Hrad)))-M*(test(Haxi)*Haziy+Haxi*test(Haziy))+x*(test(Hazix)*Hazix+((test(Haxix)-test(Hrady))*(Haxix-Hrady)+Haziy*test(Haziy))))'; ... 'alpha*((test(Hrad)*Hrad-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hazi)*Hazi)/x+(test(Haxiy)+test(Hradx))*(Hrad-M*Hazi)+(test(Hrad)-M*test(Hazi))*(Hradx+Haxiy)+x*(test(Hradx)+test(Haxiy))*(Hradx+Haxiy))'; ... '0'},{'((eperp1*(test(Hazi)*Hazi-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hrad)*Hrad)+epara1*M^2*test(Haxi)*Haxi)/x+eperp1*(test(Hazix)*(Hazi-M*Hrad)+Hazix*(test(Hazi)-M*test(Hrad)))-epara1*M*(test(Haxi)*Haziy+Haxi*test(Haziy))+x*(eperp1*test(Hazix)*Hazix+epara1*((test(Haxix)-test(Hrady))*(Haxix-Hrady)+Haziy*test(Haziy))))/(epara1*eperp1)'; ... 'alpha*((test(Hrad)*Hrad-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hazi)*Hazi)/x+(test(Haxiy)+test(Hradx))*(Hrad-M*Hazi)+(test(Hrad)-M*test(Hazi))*(Hradx+Haxiy)+x*(test(Hradx)+test(Haxiy))*(Hradx+Haxiy))'; ... '0'},{'((test(Hazi)*Hazi-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*(test(Hrad)*Hrad+test(Haxi)*Haxi))/x+(test(Hazix)*(Hazi-M*Hrad)+Hazix*(test(Hazi)-M*test(Hrad)))-M*(test(Haxi)*Haziy+Haxi*test(Haziy))+x*(test(Hazix)*Hazix+((test(Haxix)-test(Hrady))*(Haxix-Hrady)+Haziy*test(Haziy))))/e1'; ... 'alpha*((test(Hrad)*Hrad-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hazi)*Hazi)/x+(test(Haxiy)+test(Hradx))*(Hrad-M*Hazi)+(test(Hrad)-M*test(Hazi))*(Hradx+Haxiy)+x*(test(Hradx)+test(Haxiy))*(Hradx+Haxiy))'; ... '0'},{'((eperp2*(test(Hazi)*Hazi-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hrad)*Hrad)+epara2*M^2*test(Haxi)*Haxi)/x+eperp2*(test(Hazix)*(Hazi-M*Hrad)+Hazix*(test(Hazi)-M*test(Hrad)))-epara2*M*(test(Haxi)*Haziy+Haxi*test(Haziy))+x*(eperp2*test(Hazix)*Hazix+epara2*((test(Haxix)-test(Hrady))*(Haxix-Hrady)+Haziy*test(Haziy))))/(epara2*eperp2)'; ... 'alpha*((test(Hrad)*Hrad-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hazi)*Hazi)/x+(test(Haxiy)+test(Hradx))*(Hrad-M*Hazi)+(test(Hrad)-M*test(Hazi))*(Hradx+Haxiy)+x*(test(Hradx)+test(Haxiy))*(Hradx+Haxiy))'; ... '0'},{'((eperp3*(test(Hazi)*Hazi-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hrad)*Hrad)+epara3*M^2*test(Haxi)*Haxi)/x+eperp3*(test(Hazix)*(Hazi-M*Hrad)+Hazix*(test(Hazi)-M*test(Hrad)))-epara3*M*(test(Haxi)*Haziy+Haxi*test(Haziy))+x*(eperp3*test(Hazix)*Hazix+epara3*((test(Haxix)-test(Hrady))*(Haxix-Hrady)+Haziy*test(Haziy))))/(epara3*eperp3)'; ... 'alpha*((test(Hrad)*Hrad-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hazi)*Hazi)/x+(test(Haxiy)+test(Hradx))*(Hrad-M*Hazi)+(test(Hrad)-M*test(Hazi))*(Hradx+Haxiy)+x*(test(Hradx)+test(Haxiy))*(Hradx+Haxiy))'; ... '0'},{'((eperp4*(test(Hazi)*Hazi-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hrad)*Hrad)+epara4*M^2*test(Haxi)*Haxi)/x+eperp4*(test(Hazix)*(Hazi-M*Hrad)+Hazix*(test(Hazi)-M*test(Hrad)))-epara4*M*(test(Haxi)*Haziy+Haxi*test(Haziy))+x*(eperp4*test(Hazix)*Hazix+epara4*((test(Haxix)-test(Hrady))*(Haxix-Hrady)+Haziy*test(Haziy))))/(epara4*eperp4)'; ... 'alpha*((test(Hrad)*Hrad-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hazi)*Hazi)/x+(test(Haxiy)+test(Hradx))*(Hrad-M*Hazi)+(test(Hrad)-M*test(Hazi))*(Hradx+Haxiy)+x*(test(Hradx)+test(Haxiy))*(Hradx+Haxiy))'; ... '0'},{'((eperp5*(test(Hazi)*Hazi-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hrad)*Hrad)+epara5*M^2*test(Haxi)*Haxi)/x+eperp5*(test(Hazix)*(Hazi-M*Hrad)+Hazix*(test(Hazi)-M*test(Hrad)))-epara5*M*(test(Haxi)*Haziy+Haxi*test(Haziy))+x*(eperp5*test(Hazix)*Hazix+epara5*((test(Haxix)-test(Hrady))*(Haxix-Hrady)+Haziy*test(Haziy))))/(epara5*eperp5)'; ... 'alpha*((test(Hrad)*Hrad-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hazi)*Hazi)/x+(test(Haxiy)+test(Hradx))*(Hrad-M*Hazi)+(test(Hrad)-M*test(Hazi))*(Hradx+Haxiy)+x*(test(Hradx)+test(Haxiy))*(Hradx+Haxiy))'; ... '0'},{'((test(Hazi)*Hazi-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*(test(Hrad)*Hrad+test(Haxi)*Haxi))/x+(test(Hazix)*(Hazi-M*Hrad)+Hazix*(test(Hazi)-M*test(Hrad)))-M*(test(Haxi)*Haziy+Haxi*test(Haziy))+x*(test(Hazix)*Hazix+((test(Haxix)-test(Hrady))*(Haxix-Hrady)+Haziy*test(Haziy))))/e2'; ... 'alpha*((test(Hrad)*Hrad-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hazi)*Hazi)/x+(test(Haxiy)+test(Hradx))*(Hrad-M*Hazi)+(test(Hrad)-M*test(Hazi))*(Hradx+Haxiy)+x*(test(Hradx)+test(Haxiy))*(Hradx+Haxiy))'; ... '0'},{'((test(Hazi)*Hazi-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*(test(Hrad)*Hrad+test(Haxi)*Haxi))/x+(test(Hazix)*(Hazi-M*Hrad)+Hazix*(test(Hazi)-M*test(Hrad)))-M*(test(Haxi)*Haziy+Haxi*test(Haziy))+x*(test(Hazix)*Hazix+((test(Haxix)-test(Hrady))*(Haxix-Hrady)+Haziy*test(Haziy))))/e3'; ... 'alpha*((test(Hrad)*Hrad-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hazi)*Hazi)/x+(test(Haxiy)+test(Hradx))*(Hrad-M*Hazi)+(test(Hrad)-M*test(Hazi))*(Hradx+Haxiy)+x*(test(Hradx)+test(Haxiy))*(Hradx+Haxiy))'; ... '0'},{'((test(Hazi)*Hazi-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*(test(Hrad)*Hrad+test(Haxi)*Haxi))/x+(test(Hazix)*(Hazi-M*Hrad)+Hazix*(test(Hazi)-M*test(Hrad)))-M*(test(Haxi)*Haziy+Haxi*test(Haziy))+x*(test(Hazix)*Hazix+((test(Haxix)-test(Hrady))*(Haxix-Hrady)+Haziy*test(Haziy))))/e4'; ... 'alpha*((test(Hrad)*Hrad-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hazi)*Hazi)/x+(test(Haxiy)+test(Hradx))*(Hrad-M*Hazi)+(test(Hrad)-M*test(Hazi))*(Hradx+Haxiy)+x*(test(Hradx)+test(Haxiy))*(Hradx+Haxiy))'; ... '0'},{'((test(Hazi)*Hazi-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*(test(Hrad)*Hrad+test(Haxi)*Haxi))/x+(test(Hazix)*(Hazi-M*Hrad)+Hazix*(test(Hazi)-M*test(Hrad)))-M*(test(Haxi)*Haziy+Haxi*test(Haziy))+x*(test(Hazix)*Hazix+((test(Haxix)-test(Hrady))*(Haxix-Hrady)+Haziy*test(Haziy))))/e5'; ... 'alpha*((test(Hrad)*Hrad-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hazi)*Hazi)/x+(test(Haxiy)+test(Hradx))*(Hrad-M*Hazi)+(test(Hrad)-M*test(Hazi))*(Hradx+Haxiy)+x*(test(Hradx)+test(Haxiy))*(Hradx+Haxiy))'; ... '0'}}; equ.dweak = {{'fc*x*(Haxitt*test(Haxi)+Hazitt*test(Hazi)+Hradtt*test(Hrad))'; ... '0';'0'},{'fc*x*(Haxitt*test(Haxi)+Hazitt*test(Hazi)+Hradtt*test(Hrad))'; ... '0';'0'},{'fc*x*(Haxitt*test(Haxi)+Hazitt*test(Hazi)+Hradtt*test(Hrad))'; ... '0';'0'},{'fc*x*(Haxitt*test(Haxi)+Hazitt*test(Hazi)+Hradtt*test(Hrad))'; ... '0';'0'},{'fc*x*(Haxitt*test(Haxi)+Hazitt*test(Hazi)+Hradtt*test(Hrad))'; ... '0';'0'},{'fc*x*(Haxitt*test(Haxi)+Hazitt*test(Hazi)+Hradtt*test(Hrad))'; ... '0';'0'},{'fc*x*(Haxitt*test(Haxi)+Hazitt*test(Hazi)+Hradtt*test(Hrad))'; ... '0';'0'},{'fc*x*(Haxitt*test(Haxi)+Hazitt*test(Hazi)+Hradtt*test(Hrad))'; ... '0';'0'},{'fc*x*(Haxitt*test(Haxi)+Hazitt*test(Hazi)+Hradtt*test(Hrad))'; ... '0';'0'},{'fc*x*(Haxitt*test(Haxi)+Hazitt*test(Hazi)+Hradtt*test(Hrad))'; ... '0';'0'},{'fc*x*(Haxitt*test(Haxi)+Hazitt*test(Hazi)+Hradtt*test(Hrad))'; ... '0';'0'}}; equ.constr = {{'0';'0';'0'},{'0';'0';'0'},{'0';'0';'0'},{'0';'0';'0'},{'0'; ... '0';'0'},{'0';'0';'0'},{'0';'0';'0'},{'0';'0';'0'},{'0';'0';'0'},{'0';'0'; ... '0'},{'0';'0';'0'}}; equ.bndweak = {{'0';'0';'0'},{'0';'0';'0'},{'0';'0';'0'},{'0';'0';'0'},{'0'; ... '0';'0'},{'0';'0';'0'},{'0';'0';'0'},{'0';'0';'0'},{'0';'0';'0'},{'0';'0'; ... '0'},{'0';'0';'0'}}; equ.bndgporder = {{1;1;1},{1;1;1},{1;1;1},{1;1;1},{1;1;1},{1;1;1},{1;1;1},{1;1; ... 1},{1;1;1},{1;1;1},{1;1;1}}; equ.style = {{{'0'},{'166','166','166'}},{{'0'},{'0','255','0'}},{{'0'},{'0', ... '0','255'}},{{'0'},{'0','255','255'}},{{'0'},{'255','0','255'}},{{'0'}, ... {'166','166','166'}},{{'0'},{'0','0','255'}},{{'0'},{'0','255','0'}},{{'0'}, ... {'0','255','255'}},{{'0'},{'255','0','255'}},{{'0'},{'166','166','166'}}}; equ.ind = [1,2,1,2,1,2,1,2,2,2,2,2,2,2,1,2,1,2,1,2,1,2,2,2,2,2,2,2,1,2, ... 1,2,1,2,1,2,2,2,2,2,2,2,1,2,1,2,1,2,1]; appl.equ = equ; xfem.appl{1} = appl; xfem.geom = flbinary('fem6','geom','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'); xfem.sdim = {'x','y'}; xfem.frame = {'xy'}; xfem.shape = {'shlag(2,''Hrad'')','shlag(2,''Hazi'')','shlag(2,''Haxi'')'}; xfem.gporder = 4; xfem.cporder = 2; xfem.sshape = 2; xfem.simplify = 'on'; xfem.border = 1; xfem.form = 'coefficient'; clear units; units.basesystem = 'SI'; xfem.units = units; clear equ equ.shape = {[1;2;3],[1;2;3]}; equ.gporder = {{1;1;1},{1;1;1}}; equ.cporder = {{1;1;1},{1;1;1}}; equ.init = {{'0';'0';'0'},{'0';'0';'0'}}; equ.dinit = {{'0';'0';'0'},{'0';'0';'0'}}; equ.weak = {{'(test(Hazi)*Hazi-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*(test(Hrad)*Hrad+test(Haxi)*Haxi))/x+test(Hazix)*(Hazi-M*Hrad)+Hazix*(test(Hazi)-M*test(Hrad))-M*(test(Haxi)*Haziy+Haxi*test(Haziy))+x*(test(Hazix)*Hazix+(test(Haxix)-test(Hrady))*(Haxix-Hrady)+Haziy*test(Haziy))'; ... 'alpha*((test(Hrad)*Hrad-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hazi)*Hazi)/x+(test(Haxiy)+test(Hradx))*(Hrad-M*Hazi)+(test(Hrad)-M*test(Hazi))*(Hradx+Haxiy)+x*(test(Haxiy)+test(Hradx))*(Hradx+Haxiy))'; ... '0'},{'((eperp1*(test(Hazi)*Hazi-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hrad)*Hrad)+epara1*M^2*test(Haxi)*Haxi)/x+eperp1*(test(Hazix)*(Hazi-M*Hrad)+Hazix*(test(Hazi)-M*test(Hrad)))-epara1*M*(test(Haxi)*Haziy+Haxi*test(Haziy))+x*(eperp1*test(Hazix)*Hazix+epara1*((test(Haxix)-test(Hrady))*(Haxix-Hrady)+Haziy*test(Haziy))))/(epara1*eperp1)'; ... 'alpha*((test(Hrad)*Hrad-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hazi)*Hazi)/x+(test(Haxiy)+test(Hradx))*(Hrad-M*Hazi)+(test(Hrad)-M*test(Hazi))*(Hradx+Haxiy)+x*(test(Haxiy)+test(Hradx))*(Hradx+Haxiy))'; ... '0'}}; equ.dweak = {{'fc*x*(Haxitt*test(Haxi)+Hazitt*test(Hazi)+Hradtt*test(Hrad))'; ... '0';'0'},{'fc*x*(Haxitt*test(Haxi)+Hazitt*test(Hazi)+Hradtt*test(Hrad))'; ... '0';'0'}}; equ.constr = {{'0';'0';'0'},{'0';'0';'0'}}; equ.c = {{{'0'},{'0'},{'0'};{'0'},{'0'},{'0'};{'0'},{'0'},{'0'}},{{'0'}, ... {'0'},{'0'};{'0'},{'0'},{'0'};{'0'},{'0'},{'0'}}}; equ.a = {{'0','0','0';'0','0','0';'0','0','0'},{'0','0','0';'0','0','0'; ... '0','0','0'}}; equ.f = {{'0';'0';'0'},{'0';'0';'0'}}; equ.ea = {{'0','0','0';'0','0','0';'0','0','0'},{'0','0','0';'0','0','0'; ... '0','0','0'}}; equ.da = {{'0','0','0';'0','0','0';'0','0','0'},{'0','0','0';'0','0','0'; ... '0','0','0'}}; equ.al = {{{'0';'0'},{'0';'0'},{'0';'0'};{'0';'0'},{'0';'0'},{'0';'0'}; ... {'0';'0'},{'0';'0'},{'0';'0'}},{{'0';'0'},{'0';'0'},{'0';'0'};{'0';'0'}, ... {'0';'0'},{'0';'0'};{'0';'0'},{'0';'0'},{'0';'0'}}}; equ.be = {{{'0';'0'},{'0';'0'},{'0';'0'};{'0';'0'},{'0';'0'},{'0';'0'}; ... {'0';'0'},{'0';'0'},{'0';'0'}},{{'0';'0'},{'0';'0'},{'0';'0'};{'0';'0'}, ... {'0';'0'},{'0';'0'};{'0';'0'},{'0';'0'},{'0';'0'}}}; equ.ga = {{{'0';'0'};{'0';'0'};{'0';'0'}},{{'0';'0'};{'0';'0'};{'0';'0'}}}; equ.sshape = {[1],[1]}; equ.sshapedim = {{1;1;1},{1;1;1}}; equ.ind = [1,2,1,2,1,2,1,2,2,2,2,2,2,2,1,2,1,2,1,2,1,2,2,2,2,2,2,2,1,2, ... 1,2,1,2,1,2,2,2,2,2,2,2,1,2,1,2,1,2,1]; equ.dim = {'Hrad','Hazi','Haxi'}; equ.var = {'absHradx_Axisymmetric_Dielectric_Resonator_Simulator',{'sqrt(Hradx^2+Hrady^2)', ... 'sqrt(Hradx^2+Hrady^2)'}, ... 'abscu1x_Axisymmetric_Dielectric_Resonator_Simulator',{'sqrt(cu1x^2+cu1y^2)', ... 'sqrt(cu1x^2+cu1y^2)'}, ... 'absHazix_Axisymmetric_Dielectric_Resonator_Simulator',{'sqrt(Hazix^2+Haziy^2)', ... 'sqrt(Hazix^2+Haziy^2)'}, ... 'abscu2x_Axisymmetric_Dielectric_Resonator_Simulator',{'sqrt(cu2x^2+cu2y^2)', ... 'sqrt(cu2x^2+cu2y^2)'}, ... 'absHaxix_Axisymmetric_Dielectric_Resonator_Simulator',{'sqrt(Haxix^2+Haxiy^2)', ... 'sqrt(Haxix^2+Haxiy^2)'}, ... 'abscu3x_Axisymmetric_Dielectric_Resonator_Simulator',{'sqrt(cu3x^2+cu3y^2)', ... 'sqrt(cu3x^2+cu3y^2)'}}; equ.expr = {'eperp',{'1','eperp1'}, ... 'epara',{'1','epara1'}}; equ.bnd.weak = {{'0';'0';'0'}}; equ.bnd.gporder = {{1;1;1}}; equ.bnd.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, ... 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1]; equ.bnd.expr = {}; equ.lock = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; equ.mlock = {[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]}; xfem.equ = equ; clear bnd bnd.weak = {{'0';'0';'0'},{'0';'0';'0'},{'0';'0';'0'}}; bnd.dweak = {{'0';'0';'0'},{'0';'0';'0'},{'0';'0';'0'}}; bnd.constr = {{'0';'0';'0'},{'Hazi';'-Haxi*nx+Hrad*ny';'0'},{'Hrad*nx+Haxi*ny'; ... '0';'0'}}; bnd.q = {{'0','0','0';'0','0','0';'0','0','0'},{'0','0','0';'0','0','0'; ... '0','0','0'},{'0','0','0';'0','0','0';'0','0','0'}}; bnd.h = {{'0','0','0';'0','0','0';'0','0','0'},{'0','0','0';'0','0','0'; ... '0','0','0'},{'0','0','0';'0','0','0';'0','0','0'}}; bnd.g = {{'0';'0';'0'},{'0';'0';'0'},{'0';'0';'0'}}; bnd.r = {{'0';'0';'0'},{'0';'0';'0'},{'0';'0';'0'}}; bnd.shape = {[1;2;3],[1;2;3],[1;2;3]}; bnd.sshape = {[1],[1],[1]}; bnd.sshapedim = {{1;1;1},{1;1;1},{1;1;1}}; bnd.gporder = {{1;1;1},{1;1;1},{1;1;1}}; bnd.cporder = {{1;1;1},{1;1;1},{1;1;1}}; bnd.init = {{'';'';''},{'';'';''},{'';'';''}}; bnd.dinit = {{'';'';''},{'';'';''},{'';'';''}}; bnd.ind = [1,2,1,1,1,1,1,1,1,1,1,1,1,1,3,1,2,1,1,1,1,1,1,1,1,1,1,1,1,3, ... 1,2,1,1,1,1,1,1,1,1,1,1,1,1,3,1,2,1,1,1,1,1,1,1,1,1,1,1,1,3,1,2,1,1,1, ... 1,1,1,1,1,1,1,1,1,3,1,2,1,1,1,1,1,1,1,1,1,1,1,1,3,1,2,1,1,1,1,1,1,1,1, ... 1,1,1,1,3,3,3,3,3,3,3,3]; bnd.dim = {'Hrad','Hazi','Haxi'}; bnd.var = {'nx_Axisymmetric_Dielectric_Resonator_Simulator',{'nx','nx','nx'}, ... 'ny_Axisymmetric_Dielectric_Resonator_Simulator',{'ny','ny','ny'}}; bnd.expr = {}; bnd.lock = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; bnd.mlock = {[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]}; xfem.bnd = bnd; clear pnt pnt.weak = {{'0';'0';'0'}}; pnt.dweak = {{'0';'0';'0'}}; pnt.constr = {{'0';'0';'0'}}; pnt.shape = {[1;2;3]}; pnt.sshape = {[1]}; pnt.sshapedim = {{1;1;1}}; pnt.init = {{'';'';''}}; pnt.dinit = {{'';'';''}}; pnt.ind = [1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1, ... 1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1]; pnt.dim = {'Hrad','Hazi','Haxi'}; pnt.var = {}; pnt.expr = {}; pnt.lock = [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; pnt.mlock = {[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]}; xfem.pnt = pnt; xfem.var = {}; xfem.expr = {'DivH','(Hrad-Hazi*M+(Haxiy+Hradx)*x)/x', ... 'MagEnDens','Hrad*Hrad+Hazi*Hazi+Haxi*Haxi', ... 'Drad','(Haxi*M-Haziy*x)/x', ... 'Dazi','-Haxix+Hrady', ... 'Daxi','(Hazi-Hrad*M+Hazix*x)/x', ... 'Erad','Drad/eperp', ... 'Eazi','Dazi/eperp', ... 'Eaxi','Daxi/epara', ... 'ElecMagSqrd','Erad*Erad+Eazi*Eazi+Eaxi*Eaxi', ... 'ElecEnDens','Erad*Drad+Eazi*Dazi+Eaxi*Daxi'}; clear descr descr.expr= {'Daxi','axial component of electric displacement','ElecEnDens','electric energy density','DivH','divergence of magnetic field (should be zero!)','MagEnDens','magnetic energy density','ElecMagSqrd','electric field strength magnitude squared','Erad','radial component of electric field strength','Drad','radial component of electric displacement','Eaxi','axial component of electric field strength','Dazi','azimuthal component of electric displacement','Eazi','azimuthal component of electric field strength'}; xfem.descr = descr; clear draw draw.p.objs = {}; draw.p.name = {}; draw.c.objs = {}; draw.c.name = {}; draw.s.objs = {flbinary('g229','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g227','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g228','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g222','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g221','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g220','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g226','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g225','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g224','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g223','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g202','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g203','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g204','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g216','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g217','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g218','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g219','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g211','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g210','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g213','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g212','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g215','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g250','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g214','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g206','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g205','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g208','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g207','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g209','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g249','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g248','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g247','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g246','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g245','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g244','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g243','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g242','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g241','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g240','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g238','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g239','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g235','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g234','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g237','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g236','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g231','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g230','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g233','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph'),flbinary('g232','draw','bragg_cavity_wall_Q_v7_corrected_reset_c33a.mph')}; draw.s.name = {'s28','s26','s27','s21','s20','s19','s25','s24','s23','s22','s1', ... 's2','s3','s15','s16','s17','s18','s10','s9','s12','s11','s14','s49','s13', ... 's5','s4','s7','s6','s8','s48','s47','s46','s45','s44','s43','s42','s41', ... 's40','s39','s37','s38','s34','s33','s36','s35','s30','s29','s32','s31'}; xfem.draw = draw; xfem.const = {'c','299792458','fc','4*pi^2/c^2','alpha','1','M','0','delta_e','0*1e-3','e1','9.8*(1+delta_e)','e2','1.0','e3','1.0','e4','1.0','e5','1.0','delta_eperp1','0*1e-3','eperp1','9.394*(1+delta_eperp1)','delta_epara1','0*1e-3','epara1','11.593*(1+delta_epara1)','eperp2','1.0','epara2','1.0','eperp3','1.0','epara3','1.0','eperp4','1.0','epara4','1.0','eperp5','1.0','epara5','1.0','e_293K_alumina','9.8','eperp_4K_sapph_UWA','9.2725','epara_4K_sapph_UWA','11.3486','eperp_293K_sapph','9.407','epara_293K_sapph','11.62','eperp_4K_sapph_NPL','9.2848','epara_4K_sapph_NPL','11.3660','eperp_300K_sapph','9.394','epara_300K_sapph','11.593'}; xfem.globalexpr = {}; clear fcns xfem.functions = {}; descr = xfem.descr; descr.const= {'epara_300K_sapph','flory and taber 1993','eperp1','relative permittivity of uniaxial_dielectric_1 perpendicular to cylindrical axis','e_293K_alumina','relative permittivity of alumina at room temperature','c','speed of light (exact!)','delta_epara1','fractional increment (for determining filling factors)','eperp3','etc.','delta_eperp1','fractional increment (for determining filling factors)','eperp2','relative permittivity of uniaxial_dielectric_2 perpendicular to cylindrical axis','M','azimuthal mode order','fc','constant used internally --do not modify','eperp_300K_sapph','flory and taber 1993','eperp_4K_sapph_UWA','UWA values for cryogenic HEMEX sapphire','e3','etc. ...','epara3','etc. ...','delta_e','fractional increment (for determining filling factors)','epara1','relative permittivity of uniaxial_dielectric_1 parallel to cylindrical axis','e1','relative permittivity of isotropic_dielectric_1','e2','ditto for isotropic_dielectric_2','epara2','ditto but parallel to cylindrical axis','eperp_293K_sapph','nominal room temperature values for same','alpha','penalty coefficient on Div H','eperp_4K_sapph_NPL','NPL values: fitted to resonance frequencies of NPL Cs HEMEX resonator'}; xfem.descr = descr; xfem.xmcases = [0]; xfem.mcases = [0]; flbinary clear; xfem.rulingmode = 'Axisymmetric_Dielectric_Resonator_Simulator'; xfem.solform = 'weak'; uq~tgui.solvemodel.toutcomp='off'; gui.solvemodel.currsolver='eig'; gui.solvemodel.solveroption=''; gui.solvemodel.postsolver='lin'; gui.solvemodel.nonlin='auto'; gui.solvemodel.ntol='1.0E-6'; gui.solvemodel.maxiter='25'; gui.solvemodel.manualdamp='off'; gui.solvemodel.hnlin='off'; gui.solvemodel.initstep='1.0'; gui.solvemodel.minstep='1.0E-4'; gui.solvemodel.rstep='10.0'; gui.solvemodel.useaugsolver='off'; gui.solvemodel.autoaugcomp='on'; gui.solvemodel.augcomp=''; gui.solvemodel.augtol='1.0E-6'; gui.solvemodel.augmaxiter='25'; gui.solvemodel.augsolver='lumped'; gui.solvemodel.atol='0.0010'; gui.solvemodel.rtol='0.01'; gui.solvemodel.tlist='0:0.1:1'; gui.solvemodel.tout='tlist'; gui.solvemodel.tsteps='free'; gui.solvemodel.manualreassem='off'; gui.solvemodel.emassconst='on'; gui.solvemodel.massconst='on'; gui.solvemodel.loadconst='on'; gui.solvemodel.constrconst='on'; gui.solvemodel.jacobianconst='on'; gui.solvemodel.constrjacobianconst='on'; gui.solvemodel.manualstep='off'; gui.solvemodel.initialstep='0.0010'; gui.solvemodel.maxorder='5'; gui.solvemodel.maxstep='1.0'; gui.solvemodel.timeusestopcond='off'; gui.solvemodel.paramusestopcond='off'; gui.solvemodel.masssingular='maybe'; gui.solvemodel.consistent='bweuler'; gui.solvemodel.estrat='0'; gui.solvemodel.complex='off'; gui.solvemodel.neigs='2'; gui.solvemodel.shift='10.0e9'; gui.solvemodel.maxeigit='300'; gui.solvemodel.etol='0.0'; gui.solvemodel.krylovdim='0'; gui.solvemodel.eigname='lambda'; gui.solvemodel.eigref='0'; gui.solvemodel.pname=''; gui.solvemodel.plist=''; gui.solvemodel.porder='1'; gui.solvemodel.manualparam='off'; gui.solvemodel.pinitstep='0.0'; gui.solvemodel.pminstep='0.0'; gui.solvemodel.pmaxstep='0.0'; gui.solvemodel.autooldcomp='on'; gui.solvemodel.oldcomp=''; gui.solvemodel.outform='auto'; gui.solvemodel.symmetric='on'; gui.solvemodel.symmhermit='symmetric'; gui.solvemodel.method='eliminate'; gui.solvemodel.nullfun='auto'; gui.solvemodel.blocksize='5000'; gui.solvemodel.uscale='auto'; gui.solvemodel.manscale=''; gui.solvemodel.rowscale='on'; gui.solvemodel.conjugate='off'; gui.solvemodel.complexfun='off'; gui.solvemodel.matherr='on'; gui.solvemodel.solfile='off'; gui.solvemodel.adaptgeom='currgeom'; gui.solvemodel.maxt='10000000'; gui.solvemodel.rmethod='longest'; gui.solvemodel.resmethod='coefficient'; gui.solvemodel.resorder='0'; gui.solvemodel.l2scale='1'; gui.solvemodel.l2staborder='2'; gui.solvemodel.eigselect='1'; gui.solvemodel.tpfun='fltpft'; gui.solvemodel.ngen='2'; gui.solvemodel.tpmult='1.7'; gui.solvemodel.tpworst='0.5'; gui.solvemodel.tpfract='0.5'; gui.solvemodel.autosolver='on'; gui.solvemodel.solcomp='Hazi,Hrad,Haxi'; gui.solvemodel.outcomp='Hazi,Hrad,Haxi'; gui.solvemodel.inittype='init_expr_currsol_radio'; gui.solvemodel.initsolnum='Automatic'; gui.solvemodel.inittime='0'; gui.solvemodel.utype='u_init_radio'; gui.solvemodel.usolnum='Automatic'; gui.solvemodel.utime='0'; gui.solvemodel.scriptcommands=''; gui.solvemodel.usescript='off'; gui.solvemodel.autoscript='off'; gui.solvemodel.sameaxis='off'; gui.solvemodel.linsolvernode.currlinsolver='spooles'; gui.solvemodel.linsolvernode.type='linsolver'; gui.solvemodel.linsolvernode.droptol='0.0'; gui.solvemodel.linsolvernode.thresh='0.1'; gui.solvemodel.linsolvernode.umfalloc='0.7'; gui.solvemodel.linsolvernode.preorder='mmd'; gui.solvemodel.linsolvernode.preroworder='on'; gui.solvemodel.linsolvernode.pardreorder='nd'; gui.solvemodel.linsolvernode.pardrreorder='on'; gui.solvemodel.linsolvernode.pivotperturb='1.0E-8'; gui.solvemodel.linsolvernode.errorchk='on'; gui.solvemodel.linsolvernode.iter='2'; gui.solvemodel.linsolvernode.prefuntype='left'; gui.solvemodel.linsolvernode.iluiter='1'; gui.solvemodel.linsolvernode.itol='1.0E-6'; gui.solvemodel.linsolvernode.rhob='400.0'; gui.solvemodel.linsolvernode.maxlinit='10000'; gui.solvemodel.linsolvernode.itrestart='50'; gui.solvemodel.linsolvernode.seconditer='1'; gui.solvemodel.linsolvernode.relax='1.0'; gui.solvemodel.linsolvernode.amgauto='3'; gui.solvemodel.linsolvernode.mglevels='6'; gui.solvemodel.linsolvernode.mgcycle='v'; gui.solvemodel.linsolvernode.maxcoarsedof='5000'; gui.solvemodel.linsolvernode.oocmemory='128.0'; gui.solvemodel.linsolvernode.oocfilename=''; gui.solvemodel.linsolvernode.modified='off'; gui.solvemodel.linsolvernode.fillratio='2.0'; gui.solvemodel.linsolvernode.respectpattern='on'; gui.solvemodel.linsolvernode.droptype='droptol'; gui.solvemodel.linsolvernode.vankavars=''; gui.solvemodel.linsolvernode.vankasolv='gmres'; gui.solvemodel.linsolvernode.vankatol='0.02'; gui.solvemodel.linsolvernode.vankarestart='100'; gui.solvemodel.linsolvernode.vankarelax='0.8'; gui.solvemodel.linsolvernode.mgauto='shape'; gui.solvemodel.linsolvernode.rmethod='regular'; gui.solvemodel.linsolvernode.coarseassem='on'; gui.solvemodel.linsolvernode.meshscale='2'; gui.solvemodel.linsolvernode.mgautolevels='2'; gui.solvemodel.linsolvernode.mgkeep='off'; gui.solvemodel.linsolvernode.mggeom='Geom1'; gui.solvemodel.linsolvernode.mcase0='on'; gui.solvemodel.linsolvernode.mgassem0='on'; gui.registry.general_currentmodel='Geom1'; gui.registry.general_currmeshcase='0'; gui.registry.general_savedonserver='off'; gui.registry.general_savedchanges='off'; gui.registry.general_rulingmode=''; gui.registry.general_incompletemfilehistory='off'; gui.registry.saved_license='1012177'; gui.registry.saved_version='COMSOL 3.3.0.511'; gui.registry.info_modelname=''; gui.registry.info_author=''; gui.registry.info_company=''; gui.registry.info_department=''; gui.registry.info_reference=''; gui.registry.info_url=''; gui.registry.info_saveddate='1203457666111'; gui.registry.info_creationdate='1141412764044'; gui.registry.info_modelresult=''; gui.reportregistry.report_contents='title=1;toc=1;modelprop=1;const=1;geometry=1;model_Geom1=1;mesh=1;appmode_Axisymmetric_Dielectric_Resonator_Simulator=1;lib_material=1;lib_section=1;coord=1;odesettings=1;intcoup=1;extrcoup=1;projcoup=1;percond=1;idcond=1;functions=1;interp=1;solver=1;post=1;equations=1;variables=1;'; gui.reportregistry.report_outputformat='html'; gui.reportregistry.report_filename='D:\My Documents\Work\Time Programme\Time 2003-2006\E03T 15 A high-flux rubidium\Topics\Finite Element Simulator\AxSim Distributed Bragg Reflector\10th Oct 2006\flory_DBR_10GHz_rep_mo1.html'; gui.reportregistry.report_autoopen='on'; gui.reportregistry.report_paperformat='a4'; gui.reportregistry.report_includedefaults='on'; gui.reportregistry.report_template='none'; gui.reportregistry.report_showemptysections='off'; gui.flmodel{1}.modelname='Geom1'; gui.flmodel{1}.currmode='subdomain'; gui.flmodel{1}.currappl='0'; gui.flmodel{1}.axis.xmin='-0.00822746991971943'; gui.flmodel{1}.axis.xmax='0.06312747165421459'; gui.flmodel{1}.axis.ymin='-0.0033824998885393144'; gui.flmodel{1}.axis.ymax='0.0710324976593256'; gui.flmodel{1}.axis.zmin='-1.0'; gui.flmodel{1}.axis.zmax='1.0'; gui.flmodel{1}.axis.xspacing='0.01'; gui.flmodel{1}.axis.yspacing='0.01'; gui.flmodel{1}.axis.zspacing='0.01'; gui.flmodel{1}.axis.extrax=''; gui.flmodel{1}.axis.extray=''; gui.flmodel{1}.axis.extraz=''; gui.flmodel{1}.camera.xmin='-0.31516103167977205'; gui.flmodel{1}.camera.xmax='0.31516103167977205'; gui.flmodel{1}.camera.ymin='-0.37207498773932457'; gui.flmodel{1}.camera.ymax='0.37207498773932457'; gui.flmodel{1}.camera.camposx='0.27450000867247576'; gui.flmodel{1}.camera.camposy='0.3382499888539314'; gui.flmodel{1}.camera.camposz='3.7207498773932457'; gui.flmodel{1}.camera.camtargetx='0.27450000867247576'; gui.flmodel{1}.camera.camtargety='0.3382499888539314'; gui.flmodel{1}.camera.camtargetz='0.0'; gui.flmodel{1}.camera.camupx='0.0'; gui.flmodel{1}.camera.camupy='1.0'; gui.flmodel{1}.camera.camupz='0.0'; gui.flmodel{1}.lightmodel.headlight.type='point'; gui.flmodel{1}.lightmodel.headlight.name='headlight'; gui.flmodel{1}.lightmodel.headlight.enable='on'; gui.flmodel{1}.lightmodel.headlight.colorr='255'; gui.flmodel{1}.lightmodel.headlight.colorg='255'; gui.flmodel{1}.lightmodel.headlight.colorb='255'; gui.flmodel{1}.lightmodel.headlight.xpos='0.2745'; gui.flmodel{1}.lightmodel.headlight.ypos='0.33824998'; gui.flmodel{1}.lightmodel.headlight.zpos='3.7207499'; gui.flmodel{1}.lightmodel.scenelight{1}.type='directional'; gui.flmodel{1}.lightmodel.scenelight{1}.name='light 1'; gui.flmodel{1}.lightmodel.scenelight{1}.enable='on'; gui.flmodel{1}.lightmodel.scenelight{1}.colorr='255'; gui.flmodel{1}.lightmodel.scenelight{1}.colorg='255'; gui.flmodel{1}.lightmodel.scenelight{1}.colorb='255'; gui.flmodel{1}.lightmodel.scenelight{1}.xdir='1.0'; gui.flmodel{1}.lightmodel.scenelight{1}.ydir='-1.0'; gui.flmodel{1}.lightmodel.scenelight{1}.zdir='1.0'; gui.flmodel{1}.lightmodel.scenelight{2}.type='directional'; gui.flmodel{1}.lightmodel.scenelight{2}.name='light 2'; gui.flmodel{1}.lightmodel.scenelight{2}.enable='on'; gui.flmodel{1}.lightmodel.scenelight{2}.colorr='255'; gui.flmodel{1}.lightmodel.scenelight{2}.colorg='255'; gui.flmodel{1}.lightmodel.scenelight{2}.colorb='255'; gui.flmodel{1}.lightmodel.scenelight{2}.xdir='-1.0'; gui.flmodel{1}.lightmodel.scenelight{2}.ydir='1.0'; gui.flmodel{1}.lightmodel.scenelight{2}.zdir='-1.0'; gui.flmodel{1}.lightmodel.scenelight{3}.type='directional'; gui.flmodel{1}.lightmodel.scenelight{3}.name='light 3'; gui.flmodel{1}.lightmodel.scenelight{3}.enable='on'; gui.flmodel{1}.lightmodel.scenelight{3}.colorr='255'; gui.flmodel{1}.lightmodel.scenelight{3}.colorg='255'; gui.flmodel{1}.lightmodel.scenelight{3}.colorb='255'; gui.flmodel{1}.lightmodel.scenelight{3}.xdir='-1.0'; gui.flmodel{1}.lightmodel.scenelight{3}.ydir='-1.0'; gui.flmodel{1}.lightmodel.scenelight{3}.zdir='1.0'; gui.flmodel{1}.lightmodel.scenelight{4}.type='directional'; gui.flmodel{1}.lightmodel.scenelight{4}.name='light 4'; gui.flmodel{1}.lightmodel.scenelight{4}.enable='on'; gui.flmodel{1}.lightmodel.scenelight{4}.colorr='255'; gui.flmodel{1}.lightmodel.scenelight{4}.colorg='255'; gui.flmodel{1}.lightmodel.scenelight{4}.colorb='255'; gui.flmodel{1}.lightmodel.scenelight{4}.xdir='1.0'; gui.flmodel{1}.lightmodel.scenelight{4}.ydir='-1.0'; gui.flmodel{1}.lightmodel.scenelight{4}.zdir='-1.0'; gui.flmodel{1}.registry.axis_visible='on'; gui.flmodel{1}.registry.axis_auto='on'; gui.flmodel{1}.registry.axis_autoy='on'; gui.flmodel{1}.registry.axis_autoz='on'; gui.flmodel{1}.registry.axis_box='off'; gui.flmodel{1}.registry.axis_equal='on'; gui.flmodel{1}.registry.axis_csys='on'; gui.flmodel{1}.registry.grid_visible='on'; gui.flmodel{1}.registry.grid_auto='on'; gui.flmodel{1}.registry.grid_autoz='on'; gui.flmodel{1}.registry.grid_labels='on'; gui.flmodel{1}.registry.labels_object='on'; gui.flmodel{1}.registry.labels_vertex='off'; gui.flmodel{1}.registry.labels_edge='off'; gui.flmodel{1}.registry.labels_face='off'; gui.flmodel{1}.registry.labels_subdomain='off'; gui.flmodel{1}.registry.symbols_vertexscale='1.0'; gui.flmodel{1}.registry.symbols_edgescale='1.0'; gui.flmodel{1}.registry.symbols_facescale='1.0'; gui.flmodel{1}.registry.select_draw2d='size'; gui.flmodel{1}.registry.select_adj='cycle'; gui.flmodel{1}.registry.light_headlight='off'; gui.flmodel{1}.registry.light_scenelight='off'; gui.flmodel{1}.registry.light_shininess='0.5'; gui.flmodel{1}.registry.camera_mouse='orbit'; gui.flmodel{1}.registry.camera_camconstr='none'; gui.flmodel{1}.registry.camera_mouseconstr='none'; gui.flmodel{1}.registry.camera_perspective='off'; gui.flmodel{1}.registry.camera_moveasbox='off'; gui.flmodel{1}.registry.draw_assembly='off'; gui.flmodel{1}.registry.draw_dialog='off'; gui.flmodel{1}.registry.draw_keepborders='on'; gui.flmodel{1}.registry.draw_keepedges='off'; gui.flmodel{1}.registry.draw_multi='off'; gui.flmodel{1}.registry.draw_snap2grid='on'; gui.flmodel{1}.registry.draw_snap2vtx='on'; gui.flmodel{1}.registry.draw_solid='on'; gui.flmodel{1}.registry.draw_workplane_coordsys='on'; gui.flmodel{1}.registry.draw_workplane_showgeom='on'; gui.flmodel{1}.registry.draw_repair='off'; gui.flmodel{1}.registry.draw_repairtol=''; gui.flmodel{1}.registry.draw_projection='intersection'; gui.flmodel{1}.registry.transparency_value='1.0'; gui.flmodel{1}.registry.mesh_geomdetail='normal'; gui.flmodel{1}.registry.mesh_showquality='off'; gui.flmodel{1}.registry.post_cameraview='2'; gui.flmodel{1}.registry.graphics_scale='10.0'; gui.flmodel{1}.registry.render_mesh='off'; gui.flmodel{1}.registry.render_bndarrow='on'; gui.flmodel{1}.registry.render_vertex='off'; gui.flmodel{1}.registry.render_edge='on'; gui.flmodel{1}.registry.render_face='off'; gui.flmodel{1}.registry.highlight_vertex='off'; gui.flmodel{1}.registry.highlight_edge='on'; gui.flmodel{1}.registry.highlight_face='on'; gui.flmodel{1}.meshparam.hauto='5'; gui.flmodel{1}.meshparam.usehauto='on'; gui.flmodel{1}.meshparam.hmax=''; gui.flmodel{1}.meshparam.hmaxfact='1'; gui.flmodel{1}.meshparam.hcurve='0.3'; gui.flmodel{1}.meshparam.hgrad='1.3'; gui.flmodel{1}.meshparam.hcutoff='0.001'; gui.flmodel{1}.meshparam.hnarrow='1'; gui.flmodel{1}.meshparam.hpnt='10'; gui.flmodel{1}.meshparam.xscale='1.0'; gui.flmodel{1}.meshparam.yscale='1.0'; gui.flmodel{1}.meshparam.jiggle='on'; gui.flmodel{1}.meshparam.mcase='0'; gui.flmodel{1}.meshparam.rmethod='regular'; gui.flmodel{1}.meshparam.hmaxvtx={'','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','',''}; gui.flmodel{1}.meshparam.hgradvtx={'','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','',''}; gui.flmodel{1}.meshparam.hmaxedg={'','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','',''}; gui.flmodel{1}.meshparam.hcutoffedg={'','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','',''}; gui.flmodel{1}.meshparam.hcurveedg={'','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','',''}; gui.flmodel{1}.meshparam.hgradedg={'','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','',''}; gui.flmodel{1}.meshparam.hgradsub={'','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','',''}; gui.flmodel{1}.meshparam.hmaxsub={'','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','','',''}; gui.flmodel{1}.mapmeshparam.edgegroups={{'2','16','4','1'},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{},{}}; gui.flmodel{1}.mapmeshparam.edgelem={1,'10',2,'16',3,'4',5,'10',7,'4',9,'10',11,'4',13,'10',17,'4',32,'8',47,'4',62,'8',77,'4',92,'8'}; gui.flmodel{1}.mapmeshparam.nedges='112'; gui.flmodel{1}.mapmeshparam.nelemtype=[1,1,1,1,1,1,1,1,1,1,1,1,1,1]; gui.flmodel{1}.mapmeshparam.reverse=[0,0,0,0,0,0,0,0,0,0,0,0,0,0]; gui.flmodel{1}.mapmeshparam.ratio={'1','1','1','1','1','1','1','1','1','1','1','1','1','1'}; gui.flmodel{1}.mapmeshparam.hauto='5'; gui.flmodel{1}.postmodel.postplot.triplot='on'; gui.flmodel{1}.postmodel.postplot.tridata={'log10(ElecMagSqrd+3e+1)'}; gui.flmodel{1}.postmodel.postplot.trirangeauto='on'; gui.flmodel{1}.postmodel.postplot.trirangemin='1.4771212547196624'; gui.flmodel{1}.postmodel.postplot.trirangemax='6.627796266763021'; gui.flmodel{1}.postmodel.postplot.tricont='on'; gui.flmodel{1}.postmodel.postplot.triheightdata={'Hrad'}; gui.flmodel{1}.postmodel.postplot.triheightdatacheck='off'; gui.flmodel{1}.postmodel.postplot.trimap='jet'; gui.flmodel{1}.postmodel.postplot.trimapdepth='1024'; gui.flmodel{1}.postmodel.postplot.tribar='on'; gui.flmodel{1}.postmodel.postplot.triusemap='on'; gui.flmodel{1}.postmodel.postplot.tricolor='255,0,0'; gui.flmodel{1}.postmodel.postplot.tricoloring='interp'; gui.flmodel{1}.postmodel.postplot.trifill='fill'; gui.flmodel{1}.postmodel.postplot.contplot='on'; gui.flmodel{1}.postmodel.postplot.contdata={''}; gui.flmodel{1}.postmodel.postplot.contcont='on'; gui.flmodel{1}.postmodel.postplot.contheightdata={'Hrad'}; gui.flmodel{1}.postmodel.postplot.contheightdatacheck='off'; gui.flmodel{1}.postmodel.postplot.contcolordata={'Hrad'}; gui.flmodel{1}.postmodel.postplot.contcolorrangeauto='on'; gui.flmodel{1}.postmodel.postplot.contcolorrangemin=''; gui.flmodel{1}.postmodel.postplot.contcolorrangemax=''; gui.flmodel{1}.postmodel.postplot.contcolordatacheck='off'; gui.flmodel{1}.postmodel.postplot.contmap='cool'; gui.flmodel{1}.postmodel.postplot.contmapdepth='1024'; gui.flmodel{1}.postmodel.postplot.contbar='on'; gui.flmodel{1}.postmodel.postplot.contusemap='off'; gui.flmodel{1}.postmodel.postplot.contcolor='255,255,255'; gui.flmodel{1}.postmodel.postplot.contlevels='8'; gui.flmodel{1}.postmodel.postplot.contvectorlevels=''; gui.flmodel{1}.postmodel.postplot.contisvector='off'; gui.flmodel{1}.postmodel.postplot.contlabel='off'; gui.flmodel{1}.postmodel.postplot.contfill='off'; gui.flmodel{1}.postmodel.postplot.linplot='off'; gui.flmodel{1}.postmodel.postplot.lindata={'Hrad'}; gui.flmodel{1}.postmodel.postplot.linrangeauto='on'; gui.flmodel{1}.postmodel.postplot.linrangemin='-0.7800750712289763'; gui.flmodel{1}.postmodel.postplot.linrangemax='0.7797135542749551'; gui.flmodel{1}.postmodel.postplot.lincont='on'; gui.flmodel{1}.postmodel.postplot.linheightdata={'Hrad'}; gui.flmodel{1}.postmodel.postplot.linheightdatacheck='off'; gui.flmodel{1}.postmodel.postplot.linmap='jet'; gui.flmodel{1}.postmodel.postplot.linmapdepth='1024'; gui.flmodel{1}.postmodel.postplot.linbar='on'; gui.flmodel{1}.postmodel.postplot.linusemap='on'; gui.flmodel{1}.postmodel.postplot.lincolor='255,0,0'; gui.flmodel{1}.postmodel.postplot.lincoloring='interp'; gui.flmodel{1}.postmodel.postplot.arrowplot='off'; gui.flmodel{1}.postmodel.postplot.arrowploton='sub'; gui.flmodel{1}.postmodel.postplot.arrowdata={'Hrad','Haxi'}; gui.flmodel{1}.postmodel.postplot.arrowbnddata={'Hrad','Hazi'}; gui.flmodel{1}.postmodel.postplot.arrowheightdata={'Hrad'}; gui.flmodel{1}.postmodel.postplot.arrowheightdatacheck='off'; gui.flmodel{1}.postmodel.postplot.arrowxspacing='15'; gui.flmodel{1}.postmodel.postplot.arrowxvectorspacing=''; gui.flmodel{1}.postmodel.postplot.arrowxisvector='off'; gui.flmodel{1}.postmodel.postplot.arrowyspacing='15'; gui.flmodel{1}.postmodel.postplot.arrowyvectorspacing=''; gui.flmodel{1}.postmodel.postplot.arrowyisvector='off'; gui.flmodel{1}.postmodel.postplot.arrowtype='arrow'; gui.flmodel{1}.postmodel.postplot.arrowlength='proportional'; gui.flmodel{1}.postmodel.postplot.arrowcolor='255,255,255'; gui.flmodel{1}.postmodel.postplot.arrowautoscale='off'; gui.flmodel{1}.postmodel.postplot.arrowscale='1.4'; gui.flmodel{1}.postmodel.postplot.princplot='off'; gui.flmodel{1}.postmodel.postplot.princdata={'','','','','','','','','','','',''}; gui.flmodel{1}.postmodel.postplot.princheightdata={'Hrad'}; gui.flmodel{1}.postmodel.postplot.princheightdatacheck='off'; gui.flmodel{1}.postmodel.postplot.princxspacing='8'; gui.flmodel{1}.postmodel.postplot.princxvectorspacing=''; gui.flmodel{1}.postmodel.postplot.princxisvector='off'; gui.flmodel{1}.postmodel.postplot.princyspacing='8'; gui.flmodel{1}.postmodel.postplot.princyvectorspacing=''; gui.flmodel{1}.postmodel.postplot.princyisvector='off'; gui.flmodel{1}.postmodel.postplot.princtype='arrow'; gui.flmodel{1}.postmodel.postplot.princlength='proportional'; gui.flmodel{1}.postmodel.postplot.princcolor='0,153,0'; gui.flmodel{1}.postmodel.postplot.princautoscale='on'; gui.flmodel{1}.postmodel.postplot.princscale='1'; gui.flmodel{1}.postmodel.postplot.flowplot='off'; gui.flmodel{1}.postmodel.postplot.flowdata={'Hrad','Hazi'}; gui.flmodel{1}.postmodel.postplot.flowuseexpression='off'; gui.flmodel{1}.postmodel.postplot.flowcolor='255,0,0'; gui.flmodel{1}.postmodel.postplot.flowcolordata={'Hrad'}; gui.flmodel{1}.postmodel.postplot.flowmap='jet'; gui.flmodel{1}.postmodel.postplot.flowmapdepth='1024'; gui.flmodel{1}.postmodel.postplot.flowbar='on'; gui.flmodel{1}.postmodel.postplot.flowheightdata={'Hrad'}; gui.flmodel{1}.postmodel.postplot.flowheightdatacheck='off'; gui.flmodel{1}.postmodel.postplot.flowlines='20'; gui.flmodel{1}.postmodel.postplot.flowstart='sub'; gui.flmodel{1}.postmodel.postplot.flowstartx=''; gui.flmodel{1}.postmodel.postplot.flowstarty=''; gui.flmodel{1}.postmodel.postplot.flowisstartvector='off'; gui.flmodel{1}.postmodel.postplot.flowtol='0.01'; gui.flmodel{1}.postmodel.postplot.flowstattol='0.01'; gui.flmodel{1}.postmodel.postplot.flowlooptol='0.01'; gui.flmodel{1}.postmodel.postplot.flowmaxtime='Inf'; gui.flmodel{1}.postmodel.postplot.flowmaxsteps='5000'; gui.flmodel{1}.postmodel.postplot.flowback='on'; gui.flmodel{1}.postmodel.postplot.flownormal='off'; gui.flmodel{1}.postmodel.postplot.flowdistuniform='0.05'; gui.flmodel{1}.postmodel.postplot.flowlinesvel='20'; gui.flmodel{1}.postmodel.postplot.flowseedmanual='off'; gui.flmodel{1}.postmodel.postplot.flowseed1=''; gui.flmodel{1}.postmodel.postplot.flowseed2=''; gui.flmodel{1}.postmodel.postplot.flowinitref='1'; gui.flmodel{1}.postmodel.postplot.flowignoredist='0.5'; gui.flmodel{1}.postmodel.postplot.flowsat='1.3'; gui.flmodel{1}.postmodel.postplot.flowdistend='0.5'; gui.flmodel{1}.postmodel.postplot.flowdens='none'; gui.flmodel{1}.postmodel.postplot.partplot='off'; gui.flmodel{1}.postmodel.postplot.partmasstype='mass'; gui.flmodel{1}.postmodel.postplot.partplotas='lines'; gui.flmodel{1}.postmodel.postplot.predefforces=''; gui.flmodel{1}.postmodel.postplot.partmass='1'; gui.flmodel{1}.postmodel.postplot.partforce={'','',''}; gui.flmodel{1}.postmodel.postplot.part_massless_flowdata={'Hrad','Hazi'}; gui.flmodel{1}.postmodel.postplot.parttstartauto='on'; gui.flmodel{1}.postmodel.postplot.parttstart=''; gui.flmodel{1}.postmodel.postplot.partvelstart={'0','0','0'}; gui.flmodel{1}.postmodel.postplot.partstartptssel='partstart_explicit'; gui.flmodel{1}.postmodel.postplot.partstartdl=''; gui.flmodel{1}.postmodel.postplot.partstartedim1levels='10'; gui.flmodel{1}.postmodel.postplot.partstartedim1vectorlevels=''; gui.flmodel{1}.postmodel.postplot.partstartedim1isvector='off'; gui.flmodel{1}.postmodel.postplot.explicitcoord={'0','0'}; gui.flmodel{1}.postmodel.postplot.partuseexpression='off'; gui.flmodel{1}.postmodel.postplot.partcolor='255,0,0'; gui.flmodel{1}.postmodel.postplot.partcolordata={'Hrad'}; gui.flmodel{1}.postmodel.postplot.partmap='jet'; gui.flmodel{1}.postmodel.postplot.partmapdepth='1024'; gui.flmodel{1}.postmodel.postplot.partbar='on'; gui.flmodel{1}.postmodel.postplot.partpointcolor='255,0,0'; gui.flmodel{1}.postmodel.postplot.partpointautoscale='on'; gui.flmodel{1}.postmodel.postplot.partpointscale='1'; gui.flmodel{1}.postmodel.postplot.partdroptype='once'; gui.flmodel{1}.postmodel.postplot.partdroptimes=''; gui.flmodel{1}.postmodel.postplot.partdropfreq=''; gui.flmodel{1}.postmodel.postplot.partbnd='stick'; gui.flmodel{1}.postmodel.postplot.partmasslessrtol='0.001'; gui.flmodel{1}.postmodel.postplot.partmasslessatolmanual='off'; gui.flmodel{1}.postmodel.postplot.partmasslessatol={''}; gui.flmodel{1}.postmodel.postplot.partmasslessstepsizemanual='off'; gui.flmodel{1}.postmodel.postplot.partmasslesstendauto='on'; gui.flmodel{1}.postmodel.postplot.partmasslessmaxstepsauto='on'; gui.flmodel{1}.postmodel.postplot.partmasslessedgetol='0.001'; gui.flmodel{1}.postmodel.postplot.partmasslesstvar='partt'; gui.flmodel{1}.postmodel.postplot.partmasslessstatic='off'; gui.flmodel{1}.postmodel.postplot.partmasslessres='5'; gui.flmodel{1}.postmodel.postplot.partrtol='0.001'; gui.flmodel{1}.postmodel.postplot.partatolmanual='off'; gui.flmodel{1}.postmodel.postplot.partatol={'',''}; gui.flmodel{1}.postmodel.postplot.partstepsizemanual='off'; gui.flmodel{1}.postmodel.postplot.parttendauto='on'; gui.flmodel{1}.postmodel.postplot.partmaxstepsauto='on'; gui.flmodel{1}.postmodel.postplot.partedgetol='0.001'; gui.flmodel{1}.postmodel.postplot.partvelvar={'partu','partv','partw'}; gui.flmodel{1}.postmodel.postplot.parttvar='partt'; gui.flmodel{1}.postmodel.postplot.partstatic='off'; gui.flmodel{1}.postmodel.postplot.partres='5'; gui.flmodel{1}.postmodel.postplot.maxminplot='off'; gui.flmodel{1}.postmodel.postplot.maxminsubdata={'Hrad'}; gui.flmodel{1}.postmodel.postplot.maxminsubdatacheck='on'; gui.flmodel{1}.postmodel.postplot.maxminbnddata={'Hrad'}; gui.flmodel{1}.postmodel.postplot.maxminbnddatacheck='off'; gui.flmodel{1}.postmodel.postplot.geom='on'; gui.flmodel{1}.postmodel.postplot.roughplot='off'; gui.flmodel{1}.postmodel.postplot.autorefine='on'; gui.flmodel{1}.postmodel.postplot.refine='3'; gui.flmodel{1}.postmodel.postplot.geomnum={'Geom1'}; gui.flmodel{1}.postmodel.postplot.phase='0'; gui.flmodel{1}.postmodel.postplot.solnum='0'; gui.flmodel{1}.postmodel.postplot.selectvia='stored'; gui.flmodel{1}.postmodel.postplot.autotitle='on'; gui.flmodel{1}.postmodel.postplot.customtitle=''; gui.flmodel{1}.postmodel.postplot.smoothinternal='on'; gui.flmodel{1}.postmodel.postplot.useellogic='off'; gui.flmodel{1}.postmodel.postplot.ellogic=''; gui.flmodel{1}.postmodel.postplot.ellogictype='all'; gui.flmodel{1}.postmodel.postplot.deformplot='off'; gui.flmodel{1}.postmodel.postplot.deformsub='on'; gui.flmodel{1}.postmodel.postplot.deformbnd='on'; gui.flmodel{1}.postmodel.postplot.deformsubdata={'Hrad','Hazi'}; gui.flmodel{1}.postmodel.postplot.deformbnddata={'Hrad','Hazi'}; gui.flmodel{1}.postmodel.postplot.deformautoscale='on'; gui.flmodel{1}.postmodel.postplot.deformscale='0.002997'; gui.flmodel{1}.postmodel.postplot.animate_solnum='0'; gui.flmodel{1}.postmodel.postplot.animate_selectvia='stored'; gui.flmodel{1}.postmodel.postplot.filetype='AVI'; gui.flmodel{1}.postmodel.postplot.width='640'; gui.flmodel{1}.postmodel.postplot.height='480'; gui.flmodel{1}.postmodel.postplot.fps='10'; gui.flmodel{1}.postmodel.postplot.statfunctype='full'; gui.flmodel{1}.postmodel.postplot.statnframes='11'; gui.flmodel{1}.postmodel.postplot.reverse='off'; gui.flmodel{1}.postmodel.postplot.movieinmatlab='off'; gui.flmodel{1}.postmodel.postplot.copyaxis='off'; gui.flmodel{1}.postmodel.intdata{1}.intdata={'Hrad'}; gui.flmodel{1}.postmodel.intdata{1}.phase='0'; gui.flmodel{1}.postmodel.intdata{1}.solnum='0'; gui.flmodel{1}.postmodel.intdata{1}.selectvia='stored'; gui.flmodel{1}.postmodel.intdata{2}.order='4'; gui.flmodel{1}.postmodel.intdata{2}.autoorder='on'; gui.flmodel{1}.postmodel.intdata{2}.multiplyexpr='off'; gui.flmodel{1}.postmodel.intdata{2}.intdata={'2*pi*x*(abs(Hazi)^2+abs((Haxi*nx-Hrad*ny))^2)'}; gui.flmodel{1}.postmodel.intdata{2}.phase='0'; gui.flmodel{1}.postmodel.intdata{2}.solnum='0'; gui.flmodel{1}.postmodel.intdata{2}.selectvia='stored'; gui.flmodel{1}.postmodel.intdata{3}.order='4'; gui.flmodel{1}.postmodel.intdata{3}.autoorder='on'; gui.flmodel{1}.postmodel.intdata{3}.multiplyexpr='off'; gui.flmodel{1}.postmodel.intdata{3}.intdata={'2*pi*x*(Hrad*Hrad+Hazi*Hazi+Haxi*Haxi)'}; gui.flmodel{1}.postmodel.intdata{3}.phase='0'; gui.flmodel{1}.postmodel.intdata{3}.solnum='0'; gui.flmodel{1}.postmodel.intdata{3}.selectvia='stored'; gui.flmodel{1}.postmodel.domainplot.colordata={'Hrad'}; gui.flmodel{1}.postmodel.domainplot.colorrangeauto='on'; gui.flmodel{1}.postmodel.domainplot.colorrangemin=''; gui.flmodel{1}.postmodel.domainplot.colorrangemax=''; gui.flmodel{1}.postmodel.domainplot.colorcont='on'; gui.flmodel{1}.postmodel.domainplot.surfacesurfacemap='jet'; gui.flmodel{1}.postmodel.domainplot.surfacesurfacemapdepth='1024'; gui.flmodel{1}.postmodel.domainplot.surfacesurfacebar='on'; gui.flmodel{1}.postmodel.domainplot.surfacesurfaceusemap='on'; gui.flmodel{1}.postmodel.domainplot.surfacesurfacecolor='255,0,0'; gui.flmodel{1}.postmodel.domainplot.surfacesurfacecoloring='interp'; gui.flmodel{1}.postmodel.domainplot.surfacesurfacefill='fill'; gui.flmodel{1}.postmodel.domainplot.extrusion='off'; gui.flmodel{1}.postmodel.domainplot.lineyaxisdata={'Hrad'}; gui.flmodel{1}.postmodel.domainplot.lineyaxiscont='on'; gui.flmodel{1}.postmodel.domainplot.linexaxisxaxistype='arc'; gui.flmodel{1}.postmodel.domainplot.linexaxisuseexpr='off'; gui.flmodel{1}.postmodel.domainplot.linexaxisdata={'Hrad'}; gui.flmodel{1}.postmodel.domainplot.linelinestyle='solid'; gui.flmodel{1}.postmodel.domainplot.linelinecolor='cyclecolor'; gui.flmodel{1}.postmodel.domainplot.linelinemarker='none'; gui.flmodel{1}.postmodel.domainplot.linelegend='off'; gui.flmodel{1}.postmodel.domainplot.linelinelabels='off'; gui.flmodel{1}.postmodel.domainplot.linecolor='255,0,0'; gui.flmodel{1}.postmodel.domainplot.linesurfacemap='jet'; gui.flmodel{1}.postmodel.domainplot.linesurfacemapdepth='1024'; gui.flmodel{1}.postmodel.domainplot.linesurfacebar='on'; gui.flmodel{1}.postmodel.domainplot.linesurfaceusemap='on'; gui.flmodel{1}.postmodel.domainplot.linesurfacecolor='255,0,0'; gui.flmodel{1}.postmodel.domainplot.linesurfacecoloring='interp'; gui.flmodel{1}.postmodel.domainplot.linesurfacefill='fill'; gui.flmodel{1}.postmodel.domainplot.pointyaxisdata={'Hrad'}; gui.flmodel{1}.postmodel.domainplot.pointxxaxistype=''; gui.flmodel{1}.postmodel.domainplot.pointxuseexpr='off'; gui.flmodel{1}.postmodel.domainplot.pointxdata={'Hrad'}; gui.flmodel{1}.postmodel.domainplot.pointlinestyle='solid'; gui.flmodel{1}.postmodel.domainplot.pointlinecolor='cyclecolor'; gui.flmodel{1}.postmodel.domainplot.pointlinemarker='none'; gui.flmodel{1}.postmodel.domainplot.pointlegend='off'; gui.flmodel{1}.postmodel.domainplot.pointlinelabels='off'; gui.flmodel{1}.postmodel.domainplot.pointcolor='255,0,0'; gui.flmodel{1}.postmodel.domainplot.crossdispcolor='255,0,0'; gui.flmodel{1}.postmodel.domainplot.phase='0'; gui.flmodel{1}.postmodel.domainplot.solnum='0'; gui.flmodel{1}.postmodel.domainplot.selectvia='stored'; gui.flmodel{1}.postmodel.domainplot.autotitle='on'; gui.flmodel{1}.postmodel.domainplot.customtitle=''; gui.flmodel{1}.postmodel.domainplot.autolabelx='on'; gui.flmodel{1}.postmodel.domainplot.customlabelx=''; gui.flmodel{1}.postmodel.domainplot.autolabely='on'; gui.flmodel{1}.postmodel.domainplot.customlabely=''; gui.flmodel{1}.postmodel.domainplot.axistype={'lin','lin'}; gui.flmodel{1}.postmodel.domainplot.smoothinternal='on'; gui.flmodel{1}.postmodel.domainplot.autorefine='on'; gui.flmodel{1}.postmodel.domainplot.refine='3'; gui.flmodel{1}.postmodel.domainplot.plottypeind='0'; gui.flmodel{1}.postmodel.crossplot.extrusion='off'; gui.flmodel{1}.postmodel.crossplot.lineyaxisdata={'Haxi'}; gui.flmodel{1}.postmodel.crossplot.linexaxisxaxistype='x'; gui.flmodel{1}.postmodel.crossplot.linexaxisuseexpr='off'; gui.flmodel{1}.postmodel.crossplot.linexaxisdata={'Hrad'}; gui.flmodel{1}.postmodel.crossplot.linelinestyle='solid'; gui.flmodel{1}.postmodel.crossplot.linelinecolor='cyclecolor'; gui.flmodel{1}.postmodel.crossplot.linelinemarker='none'; gui.flmodel{1}.postmodel.crossplot.linelegend='off'; gui.flmodel{1}.postmodel.crossplot.linelinelabels='off'; gui.flmodel{1}.postmodel.crossplot.linecolor='255,0,0'; gui.flmodel{1}.postmodel.crossplot.linesurfacemap='jet'; gui.flmodel{1}.postmodel.crossplot.linesurfacemapdepth='1024'; gui.flmodel{1}.postmodel.crossplot.linesurfacebar='on'; gui.flmodel{1}.postmodel.crossplot.linesurfaceusemap='on'; gui.flmodel{1}.postmodel.crossplot.linesurfacecolor='255,0,0'; gui.flmodel{1}.postmodel.crossplot.linesurfacecoloring='interp'; gui.flmodel{1}.postmodel.crossplot.linesurfacefill='fill'; gui.flmodel{1}.postmodel.crossplot.lineresolution='200'; gui.flmodel{1}.postmodel.crossplot.linecoord={'0','0.056','0.00','0.00'}; gui.flmodel{1}.postmodel.crossplot.linelevels='5'; gui.flmodel{1}.postmodel.crossplot.linevectorlevels=''; gui.flmodel{1}.postmodel.crossplot.lineisvector='off'; gui.flmodel{1}.postmodel.crossplot.lineactive='off'; gui.flmodel{1}.postmodel.crossplot.pointyaxisdata={'Hrad'}; gui.flmodel{1}.postmodel.crossplot.pointxxaxistype=''; gui.flmodel{1}.postmodel.crossplot.pointxuseexpr='off'; gui.flmodel{1}.postmodel.crossplot.pointxdata={'Hrad'}; gui.flmodel{1}.postmodel.crossplot.pointlinestyle='solid'; gui.flmodel{1}.postmodel.crossplot.pointlinecolor='cyclecolor'; gui.flmodel{1}.postmodel.crossplot.pointlinemarker='none'; gui.flmodel{1}.postmodel.crossplot.pointlegend='off'; gui.flmodel{1}.postmodel.crossplot.pointlinelabels='off'; gui.flmodel{1}.postmodel.crossplot.pointcolor='255,0,0'; gui.flmodel{1}.postmodel.crossplot.pointcoord={'0','0'}; gui.flmodel{1}.postmodel.crossplot.crossdispcolor='255,0,0'; gui.flmodel{1}.postmodel.crossplot.phase='0'; gui.flmodel{1}.postmodel.crossplot.solnum='0'; gui.flmodel{1}.postmodel.crossplot.selectvia='stored'; gui.flmodel{1}.postmodel.crossplot.autotitle='on'; gui.flmodel{1}.postmodel.crossplot.customtitle=''; gui.flmodel{1}.postmodel.crossplot.autolabelx='on'; gui.flmodel{1}.postmodel.crossplot.customlabelx=''; gui.flmodel{1}.postmodel.crossplot.autolabely='on'; gui.flmodel{1}.postmodel.crossplot.customlabely=''; gui.flmodel{1}.postmodel.crossplot.axistype={'lin','lin'}; gui.flmodel{1}.postmodel.crossplot.plottypeind='0'; gui.flmodel{1}.postmodel.dataexport.pntdata={'Hrad'}; gui.flmodel{1}.postmodel.dataexport.pntlocation='element'; gui.flmodel{1}.postmodel.dataexport.pntlagorder='2'; gui.flmodel{1}.postmodel.dataexport.bnddata={'Hrad'}; gui.flmodel{1}.postmodel.dataexport.bndcont='on'; gui.flmodel{1}.postmodel.dataexport.bndlocation='element'; gui.flmodel{1}.postmodel.dataexport.bndlagorder='2'; gui.flmodel{1}.postmodel.dataexport.subdata={'Hrad'}; gui.flmodel{1}.postmodel.dataexport.subcont='on'; gui.flmodel{1}.postmodel.dataexport.subxspacing='10'; gui.flmodel{1}.postmodel.dataexport.subxvectorspacing=''; gui.flmodel{1}.postmodel.dataexport.subxisvector='off'; gui.flmodel{1}.postmodel.dataexport.subyspacing='10'; gui.flmodel{1}.postmodel.dataexport.subyvectorspacing=''; gui.flmodel{1}.postmodel.dataexport.subyisvector='off'; gui.flmodel{1}.postmodel.dataexport.sublocation='element'; gui.flmodel{1}.postmodel.dataexport.sublagorder='2'; gui.flmodel{1}.postmodel.dataexport.phase='0'; gui.flmodel{1}.postmodel.dataexport.solnum='0'; gui.flmodel{1}.postmodel.dataexport.selectvia='stored'; gui.flmodel{1}.postmodel.dataexport.exportformat='ptd'; gui.flmodel{1}.postmodel.dataexport.exportedim='2'; gui.flmodel{1}.postmodel.dataexport.struct='off'; gui.flmodel{1}.postmodel.datadisplay.fullprecision='off'; gui.flmodel{1}.postmodel.datadisplay.phase='0'; gui.flmodel{1}.postmodel.datadisplay.solnum='0'; gui.flmodel{1}.postmodel.datadisplay.selectvia='stored'; gui.flmodel{1}.postmodel.datadisplay.interpdata={'Hrad'}; gui.flmodel{1}.postmodel.datadisplay.coord={'0','0'}; gui.flmodel{1}.postmodel.globalplot.globalyaxisexprs={}; gui.flmodel{1}.postmodel.globalplot.globalyaxisexprsdisp={}; gui.flmodel{1}.postmodel.globalplot.linexaxisxaxistype=''; gui.flmodel{1}.postmodel.globalplot.linexaxisuseexpr='off'; gui.flmodel{1}.postmodel.globalplot.globallinestyle='solid'; gui.flmodel{1}.postmodel.globalplot.globallinecolor='cyclecolor'; gui.flmodel{1}.postmodel.globalplot.globallinemarker='none'; gui.flmodel{1}.postmodel.globalplot.globallegend='off'; gui.flmodel{1}.postmodel.globalplot.globallinelabels='off'; gui.flmodel{1}.postmodel.globalplot.globalcolor='255,0,0'; gui.flmodel{1}.postmodel.globalplot.autotitle='on'; gui.flmodel{1}.postmodel.globalplot.customtitle=''; gui.flmodel{1}.postmodel.globalplot.autolabelx='on'; gui.flmodel{1}.postmodel.globalplot.customlabelx=''; gui.flmodel{1}.postmodel.globalplot.autolabely='on'; gui.flmodel{1}.postmodel.globalplot.customlabely=''; gui.flmodel{1}.postmodel.globalplot.axistype={'lin','lin'}; gui.flmodel{1}.postmodel.globalplot.phase='0'; gui.flmodel{1}.postmodel.globalplot.solnum='0'; gui.flmodel{1}.postmodel.globalplot.selectvia='stored'; gui.flmodel{1}.postmodel.globalplot.plotimported='off'; gui.flmodel{1}.geommodel.suppressed{1}=[]; gui.flmodel{1}.geommodel.suppressed{2}=[]; gui.flmodel{1}.geommodel.suppressed{3}=[]; gui.flmodel{1}.workplane.type='0'; gui.flmodel{1}.workplane.wrkpln='0,1,0,0,0,1,0,0,0'; gui.flmodel{1}.workplane.localsys='0,1,0,0,0,0,1,0,0,0,0,1'; gui.flmodel{1}.workplane.model2d='Geom1'; gui.flmodel{1}.workplane.quicktype='10'; gui.flmodel{1}.workplane.parameter='0'; gui.flmodel{1}.workplane.zdir='up'; gui.flmodel{1}.meshmodel.meshplot.subplot='on'; gui.flmodel{1}.meshmodel.meshplot.meshsubuseexpression='off'; gui.flmodel{1}.meshmodel.meshplot.meshsubcolor='128,128,128'; gui.flmodel{1}.meshmodel.meshplot.meshsubbordercheck='off'; gui.flmodel{1}.meshmodel.meshplot.meshwiresubplot='on'; gui.flmodel{1}.meshmodel.meshplot.meshwiresubcolor='192,192,192'; gui.flmodel{1}.meshmodel.meshplot.bndplot='on'; gui.flmodel{1}.meshmodel.meshplot.meshbndcolor='0,0,0'; gui.flmodel{1}.meshmodel.meshplot.useellogic='off'; gui.flmodel{1}.meshmodel.meshplot.ellogic=''; gui.flmodel{1}.meshmodel.meshplot.ellogictype='all'; gui.flmodel{1}.meshmodel.meshplot.meshkeepfraction='1'; gui.flmodel{1}.meshmodel.meshplot.meshkeeptype='random'; gui.flmodel{1}.meshmodel.meshplot.meshmap='jet'; gui.flmodel{1}.meshmodel.meshplot.meshmapdepth='1024'; gui.flmodel{1}.meshmodel.meshplot.meshbar='on'; srcom.femlab.geom.Geom2Vc/Oxrcom.femlab.geom.Geom͹6{L geomAssoctLcom/femlab/geom/GeomAssoc;Lversionq~xpwq~wduq~NdGeom2|=-C6?1du礧?p= ף?1du礧?s?Q?p= ף?Q?s????? BezierCurve1du礧?p= ף??1du礧?s?? BezierCurve1du礧?p= ף??Q?p= ף?? BezierCurve1du礧?s??Q?s?? BezierCurveQ?p= ף??Q?s?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?5?p= ף?5?s?V_?p= ף?V_?s????? BezierCurve5?p= ף??5?s?? BezierCurve5?p= ף??V_?p= ף?? BezierCurve5?s??V_?s?? BezierCurveV_?p= ף??V_?s?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?V_?p= ף?V_?s?1du礧?p= ף?1du礧?s????? BezierCurveV_?p= ף??V_?s?? BezierCurveV_?p= ף??1du礧?p= ף?? BezierCurveV_?s??1du礧?s?? BezierCurve1du礧?p= ף??1du礧?s?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?1du礧? -?1du礧?p= ף?Q? -?Q?p= ף????? BezierCurve1du礧? -??1du礧?p= ף?? BezierCurve1du礧? -??Q? -?? BezierCurve1du礧?p= ף??Q?p= ף?? BezierCurveQ? -??Q?p= ף?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?V_? -?V_?p= ף?1du礧? -?1du礧?p= ף????? BezierCurveV_? -??V_?p= ף?? BezierCurveV_? -??1du礧? -?? BezierCurveV_?p= ף??1du礧?p= ף?? BezierCurve1du礧? -??1du礧?p= ף?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?5? -?5?p= ף?V_? -?V_?p= ף????? BezierCurve5? -??5?p= ף?? BezierCurve5? -??V_? -?? BezierCurve5?p= ף??V_?p= ף?? BezierCurveV_? -??V_?p= ף?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?p= ף?p= ף?p= ף?s?5?p= ף?5?s????? BezierCurvep= ף?p= ף??p= ף?s?? BezierCurvep= ף?p= ף??5?p= ף?? BezierCurvep= ף?s??5?s?? BezierCurve5?p= ף??5?s?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?[[%?p= ף?[[%?s?p= ף?p= ף?p= ף?s????? BezierCurve[[%?p= ף??[[%?s?? BezierCurve[[%?p= ף??p= ף?p= ף?? BezierCurve[[%?s??p= ף?s?? BezierCurvep= ף?p= ף??p= ף?s?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?^K=?p= ף?^K=?s?[[%?p= ף?[[%?s????? BezierCurve]K=?p= ף??]K=?s?? BezierCurve]K=?p= ף??[[%?p= ף?? BezierCurve]K=?s??[[%?s?? BezierCurve[[%?p= ף??[[%?s?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?p= ף?s?^K=?p= ף?^K=?s????? BezierCurvep= ף??s?? BezierCurvep= ף??]K=?p= ף?? BezierCurves??]K=?s?? BezierCurve]K=?p= ף??]K=?s?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?Q?^K=?^K=?Q????? BezierCurve?Q?? BezierCurve?^K=?? BezierCurveQ??^K=?Q?? BezierCurve^K=??^K=?Q?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?^K=?^K=?Q?[[%?[[%?Q????? BezierCurve]K=??]K=?Q?? BezierCurve]K=??[[%?? BezierCurve]K=?Q??[[%?Q?? BezierCurve[[%??[[%?Q?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?[[%?[[%?Q?p= ף?p= ף?Q????? BezierCurve[[%??[[%?Q?? BezierCurve[[%??p= ף?? BezierCurve[[%?Q??p= ף?Q?? BezierCurvep= ף??p= ף?Q?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6? -?p= ף?^K=? -?^K=?p= ף????? BezierCurve -??p= ף?? BezierCurve -??\K=? -?? BezierCurvep= ף??\K=?p= ף?? BezierCurve\K=? -??\K=?p= ף?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?^K=? -?^K=?p= ף?[[%? -?[[%?p= ף????? BezierCurve\K=? -??\K=?p= ף?? BezierCurve\K=? -??[[%? -?? BezierCurve\K=?p= ף??[[%?p= ף?? BezierCurve[[%? -??[[%?p= ף?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?[[%? -?[[%?p= ף?p= ף? -?p= ף?p= ף????? BezierCurve[[%? -??[[%?p= ף?? BezierCurve[[%? -??p= ף? -?? BezierCurve[[%?p= ף??p= ף?p= ף?? BezierCurvep= ף? -??p= ף?p= ף?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?p= ף? -?p= ף?p= ף?5? -?5?p= ף????? BezierCurvep= ף? -??p= ף?p= ף?? BezierCurvep= ף? -??5? -?? BezierCurvep= ף?p= ף??5?p= ף?? BezierCurve5? -??5?p= ף?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?[[%?Q?[[%? -?p= ף?Q?p= ף? -????? BezierCurve[[%?Q??[[%? -?? BezierCurve[[%?Q??p= ף?Q?? BezierCurve[[%? -??p= ף? -?? BezierCurvep= ף?Q??p= ף? -?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?^K=?Q?^K=? -?[[%?Q?[[%? -????? BezierCurve]K=?Q??]K=? -?? BezierCurve]K=?Q??[[%?Q?? BezierCurve]K=? -??[[%? -?? BezierCurve[[%?Q??[[%? -?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?5?Q?5? -?V_?Q?V_? -????? BezierCurve5?Q??5? -?? BezierCurve5?Q??V_?Q?? BezierCurve5? -??V_? -?? BezierCurveV_?Q??V_? -?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?p= ף?Q?p= ף? -?5?Q?5? -????? BezierCurvep= ף?Q??p= ף? -?? BezierCurvep= ף?Q??5?Q?? BezierCurvep= ף? -??5? -?? BezierCurve5?Q??5? -?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?1du礧?Q?1du礧? -?Q?Q?Q? -????? BezierCurve1du礧?Q??1du礧? -?? BezierCurve1du礧?Q??Q?Q?? BezierCurve1du礧? -??Q? -?? BezierCurveQ?Q??Q? -?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?1du礧?D?1du礧? Q?Q?D?Q? Q????? BezierCurve1du礧?D??1du礧? Q?? BezierCurve1du礧?D??Q?D?? BezierCurve1du礧? Q??Q? Q?? BezierCurveQ?D??Q? Q?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?V_?Q?V_? -?1du礧?Q?1du礧? -????? BezierCurveV_?Q??V_? -?? BezierCurveV_?Q??1du礧?Q?? BezierCurveV_? -??1du礧? -?? BezierCurve1du礧?Q??1du礧? -?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?5?5?Q?V_?V_?Q????? BezierCurve5??5?Q?? BezierCurve5??V_?? BezierCurve5?Q??V_?Q?? BezierCurveV_??V_?Q?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?p= ף?p= ף?Q?5?5?Q????? BezierCurvep= ף??p= ף?Q?? BezierCurvep= ף??5?? BezierCurvep= ף?Q??5?Q?? BezierCurve5??5?Q?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?1du礧?1du礧?Q?Q?Q?Q????? BezierCurve1du礧??1du礧?Q?? BezierCurve1du礧??Q?? BezierCurve1du礧?Q??Q?Q?? BezierCurveQ??Q?Q?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?V_?V_?Q?1du礧?1du礧?Q????? BezierCurveV_??V_?Q?? BezierCurveV_??1du礧?? BezierCurveV_?Q??1du礧?Q?? BezierCurve1du礧??1du礧?Q?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?Q? -?^K=?Q?^K=? -????? BezierCurveQ?? -?? BezierCurveQ??^K=?Q?? BezierCurve -??^K=? -?? BezierCurve^K=?Q??^K=? -?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?V_?D?V_? Q?1du礧?D?1du礧? Q????? BezierCurveV_?D??V_? Q?? BezierCurveV_?D??1du礧?D?? BezierCurveV_? Q??1du礧? Q?? BezierCurve1du礧?D??1du礧? Q?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?5?D?5? Q?V_?D?V_? Q????? BezierCurve5?D??5? Q?? BezierCurve5?D??V_?D?? BezierCurve5? Q??V_? Q?? BezierCurveV_?D??V_? Q?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?p= ף?D?p= ף? Q?5?D?5? Q????? BezierCurvep= ף?D??p= ף? Q?? BezierCurvep= ף?D??5?D?? BezierCurvep= ף? Q??5? Q?? BezierCurve5?D??5? Q?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?[[%?D?[[%? Q?p= ף?D?p= ף? Q????? BezierCurve[[%?D??[[%? Q?? BezierCurve[[%?D??p= ף?D?? BezierCurve[[%? Q??p= ף? Q?? BezierCurvep= ף?D??p= ף? Q?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?^K=?D?^K=? Q?[[%?D?[[%? Q????? BezierCurve^K=?D??^K=? Q?? BezierCurve^K=?D??[[%?D?? BezierCurve^K=? Q??[[%? Q?? BezierCurve[[%?D??[[%? Q?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?D? Q?^K=?D?^K=? Q????? BezierCurveD?? Q?? BezierCurveD??^K=?D?? BezierCurve Q??^K=? Q?? BezierCurve^K=?D??^K=? Q?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?1du礧?' ?1du礧?D?Q?' ?Q?D????? BezierCurve1du礧?' ??1du礧?D?? BezierCurve1du礧?' ??Q?' ?? BezierCurve1du礧?D??Q?D?? BezierCurveQ?' ??Q?D?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?V_?' ?V_?D?1du礧?' ?1du礧?D????? BezierCurveV_?' ??V_?D?? BezierCurveV_?' ??1du礧?' ?? BezierCurveV_?D??1du礧?D?? BezierCurve1du礧?' ??1du礧?D?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?5?' ?5?D?V_?' ?V_?D????? BezierCurve5?' ??5?D?? BezierCurve5?' ??V_?' ?? BezierCurve5?D??V_?D?? BezierCurveV_?' ??V_?D?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?p= ף?' ?p= ף?D?5?' ?5?D????? BezierCurvep= ף?' ??p= ף?D?? BezierCurvep= ף?' ??5?' ?? BezierCurvep= ף?D??5?D?? BezierCurve5?' ??5?D?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?^K=?' ?^K=?D?[[%?' ?[[%?D????? BezierCurve^K=?' ??^K=?D?? BezierCurve^K=?' ??[[%?' ?? BezierCurve^K=?D??[[%?D?? BezierCurve[[%?' ??[[%?D?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?[[%?' ?[[%?D?p= ף?' ?p= ף?D????? BezierCurve[[%?' ??[[%?D?? BezierCurve[[%?' ??p= ף?' ?? BezierCurve[[%?D??p= ף?D?? BezierCurvep= ף?' ??p= ף?D?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?V_?s?V_?' ?1du礧?s?1du礧?' ????? BezierCurveV_?s??V_?' ?? BezierCurveV_?s??1du礧?s?? BezierCurveV_?' ??1du礧?' ?? BezierCurve1du礧?s??1du礧?' ?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?5?s?5?' ?V_?s?V_?' ????? BezierCurve5?s??5?' ?? BezierCurve5?s??V_?s?? BezierCurve5?' ??V_?' ?? BezierCurveV_?s??V_?' ?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?' ?D?^K=?' ?^K=?D????? BezierCurve' ??D?? BezierCurve' ??^K=?' ?? BezierCurveD??^K=?D?? BezierCurve^K=?' ??^K=?D?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?1du礧?s?1du礧?' ?Q?s?Q?' ????? BezierCurve1du礧?s??1du礧?' ?? BezierCurve1du礧?s??Q?s?? BezierCurve1du礧?' ??Q?' ?? BezierCurveQ?s??Q?' ?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?^K=?s?^K=?' ?[[%?s?[[%?' ????? BezierCurve^K=?s??^K=?' ?? BezierCurve^K=?s??[[%?s?? BezierCurve^K=?' ??[[%?' ?? BezierCurve[[%?s??[[%?' ?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?s?' ?^K=?s?^K=?' ????? BezierCurves??' ?? BezierCurves??^K=?s?? BezierCurve' ??^K=?' ?? BezierCurve^K=?s??^K=?' ?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?p= ף?s?p= ף?' ?5?s?5?' ????? BezierCurvep= ף?s??p= ף?' ?? BezierCurvep= ף?s??5?s?? BezierCurvep= ף?' ??5?' ?? BezierCurve5?s??5?' ?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wduq~NdGeom2|=-C6?[[%?s?[[%?' ?p= ף?s?p= ף?' ????? BezierCurve[[%?s??[[%?' ?? BezierCurve[[%?s??p= ף?s?? BezierCurve[[%?' ??p= ף?' ?? BezierCurvep= ף?s??p= ף?' ?? AssocAttrib VectorInt&@p22xwxsq~Uwq~wRuq~NRGeom2|=-C6?@Q? -?p= ף?s?' ?D? Q?^K=?^K=?Q?^K=? -?^K=?p= ף?^K=?s?^K=?' ?^K=?D?^K=? Q?[[%?[[%?Q?[[%? -?[[%?p= ף?[[%?s?[[%?' ?[[%?D?[[%? Q?p= ף?p= ף?Q?p= ף? -?p= ף?p= ף?p= ף?s?p= ף?' ?p= ף?D?p= ף? Q?5?5?Q?5? -?5?p= ף?5?s?5?' ?5?D?5? Q?V_?V_?Q?V_? -?V_?p= ף?V_?s?V_?' ?V_?D?V_? Q?1du礧?1du礧?Q?1du礧? -?1du礧?p= ף?1du礧?s?1du礧?' ?1du礧?D?1du礧? Q?Q?Q?Q?Q? -?Q?p= ף?Q?s?Q?' ?Q?D?Q? Q?p? ?? ?? ?? ??  ? ? ? ? ?? ? ? ?  ?  ?  ?  ?  ?  ?  ? ? ? ?? ??? ? !?"? #?$? %?&? '?(? )?*?+?, ?-?.!?/?0"?1?2#?3?4$?5?6%?7?8&?9 ?:'?; (?<!"?=!)?>"#??"*?@#$?A#+?B$%? C$,? D%&?!E%-?! F&'?"G&.?"!H'(?#I'/?#"J(0?#K)*?$L)1?$M*+?%N*2?%$O+,?&P+3?&%Q,-? 'R,4?'&S-.?!(T-5?('U./?")V.6?)(W/0?#*X/7?*)Y08?*Z12?$+[19?+\23?%,]2:?,+^34?&-_3;?-,`45?'.a4<?.-b56?(/c5=?/.d67?)0e6>?0/f78?*1g7??10h8@?1i9:?+j:;?,k;<?-l<=?.m=>?/n>??0o?@?1pp BezierCurve?Q?? BezierCurve?^K=?? BezierCurveQ?? -?? BezierCurveQ??^K=?Q?? BezierCurve -??p= ף?? BezierCurve -??^K=? -?? BezierCurvep= ף??s?? BezierCurvep= ף??\K=?p= ף?? BezierCurves??' ?? BezierCurves??]K=?s?? BezierCurve' ??D?? BezierCurve' ??^K=?' ?? BezierCurveD?? Q?? BezierCurveD??^K=?D?? BezierCurve Q??^K=? Q?? BezierCurve^K=??^K=?Q?? BezierCurve]K=??[[%?? BezierCurve^K=?Q??^K=? -?? BezierCurve]K=?Q??[[%?Q?? BezierCurve\K=? -??\K=?p= ף?? BezierCurve]K=? -??[[%? -?? BezierCurve]K=?p= ף??]K=?s?? BezierCurve\K=?p= ף??[[%?p= ף?? BezierCurve^K=?s??^K=?' ?? BezierCurve]K=?s??[[%?s?? BezierCurve^K=?' ??^K=?D?? BezierCurve^K=?' ??[[%?' ?? BezierCurve^K=?D??^K=? Q?? BezierCurve^K=?D??[[%?D?? BezierCurve^K=? Q??[[%? Q?? BezierCurve[[%??[[%?Q?? BezierCurve[[%??p= ף?? BezierCurve[[%?Q??[[%? -?? BezierCurve[[%?Q??p= ף?Q?? BezierCurve[[%? -??[[%?p= ף?? BezierCurve[[%? -??p= ף? -?? BezierCurve[[%?p= ף??[[%?s?? BezierCurve[[%?p= ף??p= ף?p= ף?? BezierCurve[[%?s??[[%?' ?? BezierCurve[[%?s??p= ף?s?? BezierCurve[[%?' ??[[%?D?? BezierCurve[[%?' ??p= ף?' ?? BezierCurve[[%?D??[[%? Q?? BezierCurve[[%?D??p= ף?D?? BezierCurve[[%? Q??p= ף? Q?? BezierCurvep= ף??p= ף?Q?? BezierCurvep= ף??5?? BezierCurvep= ף?Q??p= ף? -?? BezierCurvep= ף?Q??5?Q?? BezierCurvep= ף? -??p= ף?p= ף?? BezierCurvep= ף? -??5? -?? BezierCurvep= ף?p= ף??p= ף?s?? BezierCurvep= ף?p= ף??5?p= ף?? BezierCurvep= ף?s??p= ף?' ?? BezierCurvep= ף?s??5?s?? BezierCurvep= ף?' ??p= ף?D?? BezierCurvep= ף?' ??5?' ?? BezierCurvep= ף?D??p= ף? Q?? BezierCurvep= ף?D??5?D?? BezierCurvep= ף? Q??5? Q?? BezierCurve5??5?Q?? BezierCurve5??V_?? BezierCurve5?Q??5? -?? BezierCurve5?Q??V_?Q?? BezierCurve5? -??5?p= ף?? BezierCurve5? -??V_? -?? BezierCurve5?p= ף??5?s?? BezierCurve5?p= ף??V_?p= ף?? BezierCurve5?s??5?' ?? BezierCurve5?s??V_?s?? BezierCurve5?' ??5?D?? BezierCurve5?' ??V_?' ?? BezierCurve5?D??5? Q?? BezierCurve5?D??V_?D?? BezierCurve5? Q??V_? Q?? BezierCurveV_??V_?Q?? BezierCurveV_??1du礧?? BezierCurveV_?Q??V_? -?? BezierCurveV_?Q??1du礧?Q?? BezierCurveV_? -??V_?p= ף?? BezierCurveV_? -??1du礧? -?? BezierCurveV_?p= ף??V_?s?? BezierCurveV_?p= ף??1du礧?p= ף?? BezierCurveV_?s??V_?' ?? BezierCurveV_?s??1du礧?s?? BezierCurveV_?' ??V_?D?? BezierCurveV_?' ??1du礧?' ?? BezierCurveV_?D??V_? Q?? BezierCurveV_?D??1du礧?D?? BezierCurveV_? Q??1du礧? Q?? BezierCurve1du礧??1du礧?Q?? BezierCurve1du礧??Q?? BezierCurve1du礧?Q??1du礧? -?? BezierCurve1du礧?Q??Q?Q?? BezierCurve1du礧? -??1du礧?p= ף?? BezierCurve1du礧? -??Q? -?? BezierCurve1du礧?p= ף??1du礧?s?? BezierCurve1du礧?p= ף??Q?p= ף?? BezierCurve1du礧?s??1du礧?' ?? BezierCurve1du礧?s??Q?s?? BezierCurve1du礧?' ??1du礧?D?? BezierCurve1du礧?' ??Q?' ?? BezierCurve1du礧?D??1du礧? Q?? BezierCurve1du礧?D??Q?D?? BezierCurve1du礧? Q??Q? Q?? BezierCurveQ??Q?Q?? BezierCurveQ?Q??Q? -?? BezierCurveQ? -??Q?p= ף?? BezierCurveQ?p= ף??Q?s?? BezierCurveQ?s??Q?' ?? BezierCurveQ?' ??Q?D?? BezierCurveQ?D??Q? Q??xwxq~Xq~Zq~\q~^q~`q~bq~dq~fq~hq~jq~lq~nq~pq~rq~tq~vq~xq~zq~|q~~q~q~q~q~q~q~q~q~q~q~q~q~q~q~q~q~q~q~q~q~q~q~q~q~q~q~q~q~q~srcom.femlab.api.client.MFileInfo3$$LfemNameq~[historyt[Ljava/lang/String;[mfileTagsAndTypest[[Ljava/lang/String;[ resetHistoryq~[ storedNamesq~Lversionq~xpwq~uq~t7% COMSOL Multiphysics Model M-file % Generated by COMSOL 3.3a (COMSOL 3.3.0.511, $Date: 2007/02/02 19:05:58 $) % Some geometry objects are stored in a separate file. % The name of this file is given by the variable 'flbinaryfile'. flclear fem % COMSOL version clear vrsn vrsn.name = 'COMSOL 3.3'; vrsn.ext = 'a'; vrsn.major = 0; vrsn.build = 511; vrsn.rcs = '$Name: $'; vrsn.date = '$Date: 2007/02/02 19:05:58 $'; fem.version = vrsn; flbinaryfile=''; % Constants fem.const = {'c','299792458', ... 'fc','4*pi^2/c^2', ... 'alpha','1', ... 'M','0', ... 'delta_e','0*1e-3', ... 'e1','9.8*(1+delta_e)', ... 'e2','1.0', ... 'e3','1.0', ... 'e4','1.0', ... 'e5','1.0', ... 'delta_eperp1','0*1e-3', ... 'eperp1','9.394*(1+delta_eperp1)', ... 'delta_epara1','0*1e-3', ... 'epara1','11.593*(1+delta_epara1)', ... 'eperp2','1.0', ... 'epara2','1.0', ... 'eperp3','1.0', ... 'epara3','1.0', ... 'eperp4','1.0', ... 'epara4','1.0', ... 'eperp5','1.0', ... 'epara5','1.0', ... 'e_293K_alumina','9.8', ... 'eperp_4K_sapph_UWA','9.2725', ... 'epara_4K_sapph_UWA','11.3486', ... 'eperp_293K_sapph','9.407', ... 'epara_293K_sapph','11.62', ... 'eperp_4K_sapph_NPL','9.2848', ... 'epara_4K_sapph_NPL','11.3660', ... 'eperp_300K_sapph','9.394', ... 'epara_300K_sapph','11.593'}; % Geometry clear draw g229=flbinary('g229','draw',flbinaryfile); g227=flbinary('g227','draw',flbinaryfile); g228=flbinary('g228','draw',flbinaryfile); g222=flbinary('g222','draw',flbinaryfile); g221=flbinary('g221','draw',flbinaryfile); g220=flbinary('g220','draw',flbinaryfile); g226=flbinary('g226','draw',flbinaryfile); g225=flbinary('g225','draw',flbinaryfile); g224=flbinary('g224','draw',flbinaryfile); g223=flbinary('g223','draw',flbinaryfile); g202=flbinary('g202','draw',flbinaryfile); g203=flbinary('g203','draw',flbinaryfile); g204=flbinary('g204','draw',flbinaryfile); g216=flbinary('g216','draw',flbinaryfile); g217=flbinary('g217','draw',flbinaryfile); g218=flbinary('g218','draw',flbinaryfile); g219=flbinary('g219','draw',flbinaryfile); g211=flbinary('g211','draw',flbinaryfile); g210=flbinary('g210','draw',flbinaryfile); g213=flbinary('g213','draw',flbinaryfile); g212=flbinary('g212','draw',flbinaryfile); g215=flbinary('g215','draw',flbinaryfile); g250=flbinary('g250','draw',flbinaryfile); g214=flbinary('g214','draw',flbinaryfile); g206=flbinary('g206','draw',flbinaryfile); g205=flbinary('g205','draw',flbinaryfile); g208=flbinary('g208','draw',flbinaryfile); g207=flbinary('g207','draw',flbinaryfile); g209=flbinary('g209','draw',flbinaryfile); g249=flbinary('g249','draw',flbinaryfile); g248=flbinary('g248','draw',flbinaryfile); g247=flbinary('g247','draw',flbinaryfile); g246=flbinary('g246','draw',flbinaryfile); g245=flbinary('g245','draw',flbinaryfile); g244=flbinary('g244','draw',flbinaryfile); g243=flbinary('g243','draw',flbinaryfile); g242=flbinary('g242','draw',flbinaryfile); g241=flbinary('g241','draw',flbinaryfile); g240=flbinary('g240','draw',flbinaryfile); g238=flbinary('g238','draw',flbinaryfile); g239=flbinary('g239','draw',flbinaryfile); g235=flbinary('g235','draw',flbinaryfile); g234=flbinary('g234','draw',flbinaryfile); g237=flbinary('g237','draw',flbinaryfile); g236=flbinary('g236','draw',flbinaryfile); g231=flbinary('g231','draw',flbinaryfile); g230=flbinary('g230','draw',flbinaryfile); g233=flbinary('g233','draw',flbinaryfile); g232=flbinary('g232','draw',flbinaryfile); draw.s.objs = {g229,g227,g228,g222,g221,g220,g226,g225,g224,g223,g202, ... g203,g204,g216,g217,g218,g219,g211,g210,g213,g212,g215,g250, ... g214,g206,g205,g208,g207,g209,g249,g248,g247,g246,g245,g244, ... g243,g242,g241,g240,g238,g239,g235,g234,g237,g236,g231,g230, ... g233,g232}; draw.s.name = {'s28','s26','s27','s21','s20','s19','s25','s24','s23','s22','s1', ... 's2','s3','s15','s16','s17','s18','s10','s9','s12','s11','s14','s49','s13', ... 's5','s4','s7','s6','s8','s48','s47','s46','s45','s44','s43','s42','s41', ... 's40','s39','s37','s38','s34','s33','s36','s35','s30','s29','s32','s31'}; draw.s.tags = {'g229','g227','g228','g222','g221','g220','g226','g225','g224', ... 'g223','g202','g203','g204','g216','g217','g218','g219','g211','g210','g213', ... 'g212','g215','g250','g214','g206','g205','g208','g207','g209','g249','g248', ... 'g247','g246','g245','g244','g243','g242','g241','g240','g238','g239','g235', ... 'g234','g237','g236','g231','g230','g233','g232'}; fem.draw = draw; fem.geom = geomcsg(fem); % (Default values are not included) % Application mode 1 clear appl appl.mode.class = 'FlPDEW'; appl.dim = {'Hrad','Hazi','Haxi','Hrad_t','Hazi_t','Haxi_t'}; appl.name = 'Axisymmetric_Dielectric_Resonator_Simulator'; appl.gporder = 4; appl.cporder = 2; appl.assignsuffix = '_Axisymmetric_Dielectric_Resonator_Simulator'; clear bnd bnd.name = {'tangential_H','normal_H','null',''}; bnd.constr = {'Hrad*nx+Haxi*ny',{'Hazi';'-Haxi*nx+Hrad*ny'},0,{'-Hrad';'-Hazi'; ... '-Haxi'}}; bnd.ind = [3,2,3,4,3,4,3,4,3,4,3,4,3,4,1,4,2,4,4,4,4,4,4,4,4,4,4,4,4,1, ... 4,2,4,4,4,4,4,4,4,4,4,4,4,4,1,4,2,4,4,4,4,4,4,4,4,4,4,4,4,1,4,2,4,4,4, ... 4,4,4,4,4,4,4,4,4,1,4,2,4,4,4,4,4,4,4,4,4,4,4,4,1,4,2,4,4,4,4,4,4,4,4, ... 4,4,4,4,1,1,1,1,1,1,1,1]; appl.bnd = bnd; clear equ equ.dweak = 'fc*x*(Haxitt*test(Haxi)+Hazitt*test(Hazi)+Hradtt*test(Hrad))'; equ.name = {'dielectric_0:vacuum','uniax_diel_1','isotrop_diel_1','uniax_diel_2', ... 'uniax_diel_3','uniax_diel_4','uniax_diel_5','isotrop_diel_2','isotrop_diel_3', ... 'isotrop_diel_4','isotrop_diel_5'}; equ.weak = {{'((test(Hazi)*Hazi-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*(test(Hrad)*Hrad+test(Haxi)*Haxi))/x+(test(Hazix)*(Hazi-M*Hrad)+Hazix*(test(Hazi)-M*test(Hrad)))-M*(test(Haxi)*Haziy+Haxi*test(Haziy))+x*(test(Hazix)*Hazix+((test(Haxix)-test(Hrady))*(Haxix-Hrady)+Haziy*test(Haziy))))'; ... 'alpha*((test(Hrad)*Hrad-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hazi)*Hazi)/x+(test(Haxiy)+test(Hradx))*(Hrad-M*Hazi)+(test(Hrad)-M*test(Hazi))*(Hradx+Haxiy)+x*(test(Hradx)+test(Haxiy))*(Hradx+Haxiy))'}, ... {'((eperp1*(test(Hazi)*Hazi-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hrad)*Hrad)+epara1*M^2*test(Haxi)*Haxi)/x+eperp1*(test(Hazix)*(Hazi-M*Hrad)+Hazix*(test(Hazi)-M*test(Hrad)))-epara1*M*(test(Haxi)*Haziy+Haxi*test(Haziy))+x*(eperp1*test(Hazix)*Hazix+epara1*((test(Haxix)-test(Hrady))*(Haxix-Hrady)+Haziy*test(Haziy))))/(epara1*eperp1)'; ... 'alpha*((test(Hrad)*Hrad-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hazi)*Hazi)/x+(test(Haxiy)+test(Hradx))*(Hrad-M*Hazi)+(test(Hrad)-M*test(Hazi))*(Hradx+Haxiy)+x*(test(Hradx)+test(Haxiy))*(Hradx+Haxiy))'}, ... {'((test(Hazi)*Hazi-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*(test(Hrad)*Hrad+test(Haxi)*Haxi))/x+(test(Hazix)*(Hazi-M*Hrad)+Hazix*(test(Hazi)-M*test(Hrad)))-M*(test(Haxi)*Haziy+Haxi*test(Haziy))+x*(test(Hazix)*Hazix+((test(Haxix)-test(Hrady))*(Haxix-Hrady)+Haziy*test(Haziy))))/e1'; ... 'alpha*((test(Hrad)*Hrad-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hazi)*Hazi)/x+(test(Haxiy)+test(Hradx))*(Hrad-M*Hazi)+(test(Hrad)-M*test(Hazi))*(Hradx+Haxiy)+x*(test(Hradx)+test(Haxiy))*(Hradx+Haxiy))'}, ... {'((eperp2*(test(Hazi)*Hazi-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hrad)*Hrad)+epara2*M^2*test(Haxi)*Haxi)/x+eperp2*(test(Hazix)*(Hazi-M*Hrad)+Hazix*(test(Hazi)-M*test(Hrad)))-epara2*M*(test(Haxi)*Haziy+Haxi*test(Haziy))+x*(eperp2*test(Hazix)*Hazix+epara2*((test(Haxix)-test(Hrady))*(Haxix-Hrady)+Haziy*test(Haziy))))/(epara2*eperp2)'; ... 'alpha*((test(Hrad)*Hrad-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hazi)*Hazi)/x+(test(Haxiy)+test(Hradx))*(Hrad-M*Hazi)+(test(Hrad)-M*test(Hazi))*(Hradx+Haxiy)+x*(test(Hradx)+test(Haxiy))*(Hradx+Haxiy))'}, ... {'((eperp3*(test(Hazi)*Hazi-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hrad)*Hrad)+epara3*M^2*test(Haxi)*Haxi)/x+eperp3*(test(Hazix)*(Hazi-M*Hrad)+Hazix*(test(Hazi)-M*test(Hrad)))-epara3*M*(test(Haxi)*Haziy+Haxi*test(Haziy))+x*(eperp3*test(Hazix)*Hazix+epara3*((test(Haxix)-test(Hrady))*(Haxix-Hrady)+Haziy*test(Haziy))))/(epara3*eperp3)'; ... 'alpha*((test(Hrad)*Hrad-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hazi)*Hazi)/x+(test(Haxiy)+test(Hradx))*(Hrad-M*Hazi)+(test(Hrad)-M*test(Hazi))*(Hradx+Haxiy)+x*(test(Hradx)+test(Haxiy))*(Hradx+Haxiy))'}, ... {'((eperp4*(test(Hazi)*Hazi-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hrad)*Hrad)+epara4*M^2*test(Haxi)*Haxi)/x+eperp4*(test(Hazix)*(Hazi-M*Hrad)+Hazix*(test(Hazi)-M*test(Hrad)))-epara4*M*(test(Haxi)*Haziy+Haxi*test(Haziy))+x*(eperp4*test(Hazix)*Hazix+epara4*((test(Haxix)-test(Hrady))*(Haxix-Hrady)+Haziy*test(Haziy))))/(epara4*eperp4)'; ... 'alpha*((test(Hrad)*Hrad-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hazi)*Hazi)/x+(test(Haxiy)+test(Hradx))*(Hrad-M*Hazi)+(test(Hrad)-M*test(Hazi))*(Hradx+Haxiy)+x*(test(Hradx)+test(Haxiy))*(Hradx+Haxiy))'}, ... {'((eperp5*(test(Hazi)*Hazi-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hrad)*Hrad)+epara5*M^2*test(Haxi)*Haxi)/x+eperp5*(test(Hazix)*(Hazi-M*Hrad)+Hazix*(test(Hazi)-M*test(Hrad)))-epara5*M*(test(Haxi)*Haziy+Haxi*test(Haziy))+x*(eperp5*test(Hazix)*Hazix+epara5*((test(Haxix)-test(Hrady))*(Haxix-Hrady)+Haziy*test(Haziy))))/(epara5*eperp5)'; ... 'alpha*((test(Hrad)*Hrad-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hazi)*Hazi)/x+(test(Haxiy)+test(Hradx))*(Hrad-M*Hazi)+(test(Hrad)-M*test(Hazi))*(Hradx+Haxiy)+x*(test(Hradx)+test(Haxiy))*(Hradx+Haxiy))'}, ... {'((test(Hazi)*Hazi-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*(test(Hrad)*Hrad+test(Haxi)*Haxi))/x+(test(Hazix)*(Hazi-M*Hrad)+Hazix*(test(Hazi)-M*test(Hrad)))-M*(test(Haxi)*Haziy+Haxi*test(Haziy))+x*(test(Hazix)*Hazix+((test(Haxix)-test(Hrady))*(Haxix-Hrady)+Haziy*test(Haziy))))/e2'; ... 'alpha*((test(Hrad)*Hrad-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hazi)*Hazi)/x+(test(Haxiy)+test(Hradx))*(Hrad-M*Hazi)+(test(Hrad)-M*test(Hazi))*(Hradx+Haxiy)+x*(test(Hradx)+test(Haxiy))*(Hradx+Haxiy))'}, ... {'((test(Hazi)*Hazi-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*(test(Hrad)*Hrad+test(Haxi)*Haxi))/x+(test(Hazix)*(Hazi-M*Hrad)+Hazix*(test(Hazi)-M*test(Hrad)))-M*(test(Haxi)*Haziy+Haxi*test(Haziy))+x*(test(Hazix)*Hazix+((test(Haxix)-test(Hrady))*(Haxix-Hrady)+Haziy*test(Haziy))))/e3'; ... 'alpha*((test(Hrad)*Hrad-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hazi)*Hazi)/x+(test(Haxiy)+test(Hradx))*(Hrad-M*Hazi)+(test(Hrad)-M*test(Hazi))*(Hradx+Haxiy)+x*(test(Hradx)+test(Haxiy))*(Hradx+Haxiy))'}, ... {'((test(Hazi)*Hazi-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*(test(Hrad)*Hrad+test(Haxi)*Haxi))/x+(test(Hazix)*(Hazi-M*Hrad)+Hazix*(test(Hazi)-M*test(Hrad)))-M*(test(Haxi)*Haziy+Haxi*test(Haziy))+x*(test(Hazix)*Hazix+((test(Haxix)-test(Hrady))*(Haxix-Hrady)+Haziy*test(Haziy))))/e4'; ... 'alpha*((test(Hrad)*Hrad-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hazi)*Hazi)/x+(test(Haxiy)+test(Hradx))*(Hrad-M*Hazi)+(test(Hrad)-M*test(Hazi))*(Hradx+Haxiy)+x*(test(Hradx)+test(Haxiy))*(Hradx+Haxiy))'}, ... {'((test(Hazi)*Hazi-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*(test(Hrad)*Hrad+test(Haxi)*Haxi))/x+(test(Hazix)*(Hazi-M*Hrad)+Hazix*(test(Hazi)-M*test(Hrad)))-M*(test(Haxi)*Haziy+Haxi*test(Haziy))+x*(test(Hazix)*Hazix+((test(Haxix)-test(Hrady))*(Haxix-Hrady)+Haziy*test(Haziy))))/e5'; ... 'alpha*((test(Hrad)*Hrad-M*(test(Hazi)*Hrad+Hazi*test(Hrad))+M^2*test(Hazi)*Hazi)/x+(test(Haxiy)+test(Hradx))*(Hrad-M*Hazi)+(test(Hrad)-M*test(Hazi))*(Hradx+Haxiy)+x*(test(Hradx)+test(Haxiy))*(Hradx+Haxiy))'}}; equ.ind = [1,2,1,2,1,2,1,2,2,2,2,2,2,2,1,2,1,2,1,2,1,2,2,2,2,2,2,2,1,2, ... 1,2,1,2,1,2,2,2,2,2,2,2,1,2,1,2,1,2,1]; appl.equ = equ; fem.appl{1} = appl; fem.border = 1; clear units; units.basesystem = 'SI'; fem.units = units; % Subdomain settings clear equ equ.ind = [1,2,1,2,1,2,1,2,2,2,2,2,2,2,1,2,1,2,1,2,1,2,2,2,2,2,2,2,1,2, ... 1,2,1,2,1,2,2,2,2,2,2,2,1,2,1,2,1,2,1]; equ.dim = {'Hrad','Hazi','Haxi'}; % Subdomain expressions equ.expr = {'eperp',{1,'eperp1'}, ... 'epara',{1,'epara1'}}; fem.equ = equ; % Scalar expressions fem.expr = {'DivH','(Hrad-Hazi*M+(Haxiy+Hradx)*x)/x', ... 'MagEnDens','Hrad*Hrad+Hazi*Hazi+Haxi*Haxi', ... 'Drad','(Haxi*M-Haziy*x)/x', ... 'Dazi','-Haxix+Hrady', ... 'Daxi','(Hazi-Hrad*M+Hazix*x)/x', ... 'Erad','Drad/eperp', ... 'Eazi','Dazi/eperp', ... 'Eaxi','Daxi/epara', ... 'ElecMagSqrd','Erad*Erad+Eazi*Eazi+Eaxi*Eaxi', ... 'ElecEnDens','Erad*Drad+Eazi*Dazi+Eaxi*Daxi'}; % Descriptions clear descr descr.expr= {'Daxi','axial component of electric displacement','ElecEnDens','electric energy density','DivH','divergence of magnetic field (should be zero!)','MagEnDens','magnetic energy density','ElecMagSqrd','electric field strength magnitude squared','Erad','radial component of electric field strength','Drad','radial component of electric displacement','Eaxi','axial component of electric field strength','Dazi','azimuthal component of electric displacement','Eazi','azimuthal component of electric field strength'}; fem.descr = descr; % Descriptions descr = fem.descr; descr.const= {'epara_300K_sapph','flory and taber 1993','eperp1','relative permittivity of uniaxial_dielectric_1 perpendicular to cylindrical axis','e_293K_alumina','relative permittivity of alumina at room temperature','c','speed of light (exact!)','delta_epara1','fractional increment (for determining filling factors)','eperp3','etc.','delta_eperp1','fractional increment (for determining filling factors)','eperp2','relative permittivity of uniaxial_dielectric_2 perpendicular to cylindrical axis','M','azimuthal mode order','fc','constant used internally --do not modify','eperp_300K_sapph','flory and taber 1993','eperp_4K_sapph_UWA','UWA values for cryogenic HEMEX sapphire','e3','etc. ...','epara3','etc. ...','delta_e','fractional increment (for determining filling factors)','epara1','relative permittivity of uniaxial_dielectric_1 parallel to cylindrical axis','e1','relative permittivity of isotropic_dielectric_1','e2','ditto for isotropic_dielectric_2','epara2','ditto but parallel to cylindrical axis','eperp_293K_sapph','nominal room temperature values for same','alpha','penalty coefficient on Div H','eperp_4K_sapph_NPL','NPL values: fitted to resonance frequencies of NPL Cs HEMEX resonator'}; fem.descr = descr; % Multiphysics fem=multiphysics(fem); uq~ptclear mfile clear vrsn vrsn.name = 'COMSOL 3.3'; vrsn.ext = 'a'; vrsn.major = 0; vrsn.build = 511; vrsn.rcs = '$Name: $'; vrsn.date = '$Date: 2007/02/02 19:05:58 $'; mfile.version=vrsn; mfile.fem='fem'; mfile.stored={'fem0','fem1'}; mfile.tags={'g229','g227','g228','g222','g221','g220','g226','g225','g224','g223','g202','g203','g204','g216','g217','g218','g219','g211','g210','g213','g212','g250','g215','g214','g206','g205','g208','g207','g209','g249','g248','g247','g246','g245','g244','g243','g242','g241','g240','g238','g239','g235','g234','g237','g236','g231','g230','g233','g232'}; mfile.types={'draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw','draw'}; x