GModule electrohydro

5 4 GfsElectroHydro GfsBox GfsGEdge { x = 0.5 } {

    Global {
      #define Volt 1.0
    }
    VariableTracer Cpos { scheme = none }
    VariableTracer Cneg { scheme = none } 
    AdaptGradient { istep = 5 } { cmax = 0.02 minlevel = 3 maxlevel = 5 } Cneg 

    Init {} {
       Phi = Volt*(1.-x/5.)
       Cpos =  1. 
       Cneg =  1. 
    }
    Time { end = 3.5 dtmax = 0.01 }
   
    SourceDiffusion Cpos 1.0
    SourceDiffusionExplicit Cpos Cpos Phi
    
    SourceDiffusion Cneg 1.0
    SourceDiffusionExplicit  Cneg -Cneg Phi 
      
    EventStop { istep = 10 } Phi 1e-4 DPhi

#    OutputSimulation { istep = 5 } stdout 
    OutputLocation { start = end } { 
	awk '{ if ($1 != "#") print $2, $9, $12, $13; }' > profile 
    } points
} {
    # Electric parameters
    perm = 1.0
    charge = (Cpos - Cneg)
}
GfsBox {
    top = Boundary 
    bottom = Boundary
    left = Boundary { 
	BcDirichlet Phi Volt 
	BcDirichlet Cpos exp(-Volt) 
	BcDirichlet Cneg exp(Volt) 
    }
}
GfsBox {
    top = Boundary 
    bottom = Boundary
}
GfsBox {
    top = Boundary 
    bottom = Boundary
}
GfsBox {
    top = Boundary 
    bottom = Boundary
}
GfsBox {
    top = Boundary 
    bottom = Boundary
    right = Boundary { 
	BcDirichlet Phi 0.
        BcDirichlet Cpos 1.
        BcDirichlet Cneg 1. 
    }
}
1 2 right
2 3 right
3 4 right
4 5 right