took usual twolayers code and changed upper voltage BC, only depends on x for now
added Vb potential in same form as Radu's paper but no time dependence yet. Can't see any difference in results
made the applied voltage Vb wider and now results are different although q is super weird
made ii->yf in Vb, added T dependence in Vb now it's roughly doing the right thing but q is very unhappy
added Vb in IC for q but still wiggly q
added Vb in IC for q but still wiggly q
got rid of Vb in IC because wasn't helping with squiggly q
added fluQ to be able to change value of fluxw
added fluQ to the steady part in epsbig code too./a.out
got rid of charge (set q=0) so as to match Radu's stuff
changed del=1d-1 to del=1 to have surface tension at higher order but not much seems to change, back to 1d-1 now
set length to be 4pi and fluQ=0 to get initial condtion for my Vbxt stuff
tried modifying case B to add Vbxt inside Vbxtepsbig but for some reason doesn't work at all
cp twolayersvbxt into vbxtepsbig and now it seems to work
now reading the hini.dat file which I manually got to look like I wanted
applying voltage simimlar to perturbation using IC from previous calculation with no voltage, but this doesn't drag any waves
saving Vbxt parameters to Vbxt.par
added xsi1 and xsi2 in Vb to make periodic but turns out there's a problem with how I define it