Skip to content
Snippets Groups Projects
Commit 5425d68f authored by Guillaume Demesy's avatar Guillaume Demesy
Browse files

Revert "normalized units : 2*pi*c=a so we find exactly"

This reverts commit 44ad8846.
parent e4b435c0
No related branches found
No related tags found
No related merge requests found
......@@ -3,6 +3,10 @@
//// _data.geo ////
///////////////////////////////////
nm = 1.;
epsilon0 = 8.854187817e-3*nm;
mu0 = 400.*Pi*nm;
cel = 1.0/(Sqrt[epsilon0 * mu0]);
deg2rad = Pi/180;
pp0 = "1Geometry/0";
......@@ -23,6 +27,7 @@ colorpp5 = "Ivory";
DefineConstant[ a_lat = {50 , Name StrCat[pp0 , "1grating period d [nm]"] , Highlight Str[colorpp0] , Closed close_menu} ];
// normalization factor
norm = a_lat/(2.*Pi*cel);
DefineConstant[
d_sq = {0.806 , Name StrCat[pp0 , "2sq [d]"] , Highlight Str[colorpp0] , Closed close_menu} ,
......@@ -53,15 +58,6 @@ DefineConstant[
flag_outEigvec = {1 , Name StrCat[pp4, "output eigenvector?"], Choices{0,1}}
];
cel = a_lat/(2*Pi);
epsf = 8.854187817e-3;
muf = 400.*Pi;
nm = 2*Pi/(a_lat*Sqrt[epsf*muf]);
epsilon0 = epsf*nm;
mu0 = muf*nm;
norm = a_lat/(2.*Pi*cel);
// Normalization
d_sq = d_sq * a_lat;
space2pml = space2pml * a_lat;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment