Using PDE to model 3D Maxwell equations for wavelength generation

Hi everyone

I want to use PDE to model 3D Maxwell equations for wavelength generation by simulating the ultrafast laser beam interacting with gasous material in a gas cell. The equations in cylindrical coordinate are:
∇2 E(ρ, z, t)− 1/c2* ∂2E(ρ, z, t)/∂ t2 =ω2(ρ, z, t)/c2* E(ρ, z, t)+μ* ∂ J(ρ, z, t)/∂ t
∇2 E(ρ, z, t)− 1/c2* ∂2E(ρ, z, t)/∂ t2 = ω2(ρ, z, t)/c2* E(ρ, z, t)+μ* ∂2P(ρ, z, t)/∂ t2 .

Initial E(ρ, z, t) can be a simple Gauss equation,
ω(ρ, z, t) is plasma frequency and well known
J(ρ, z, t) is absorption term and known.
P(ρ, z, t) is nonlinear polarization and known.
C and μ are constants.

I can use 3.5 and 4.0 Comsol.
May i ask whether it is possible to use coefficient Form PDE?
How to set up COMSol PDE mode to model above equations?
the intial laser pulse will input from one side of gas cell and output new wavelength from the other side, how to integrate the geometry with PDE?
Is that possible to model propagation process?

thank you very much for any of your help in advance.

Johansen