i have runed below system of PDES as:

c := 3*10^8;
                           
hbar := 0.105457148e-33;
                             
kB := 0.13806503e-22;
                                
epsilon0 := 0.885418782e-11;
                                                        
timeunit := 1;
                             
g := 1*timeunit;
                            
t0 := 0.10e-4/g;
                           
k := 5/(0.1e10*g);
                                                          
td := 0.10e-4/g;
B := 1;
                            
L := 0.1e-2;
                           
OD := 0.1e7;
              
eta := g*OD/(2*L);
                                 
                       
Omegap0 := .1*g;
      

Omegac(z):=(1e7*B*g)/(sqrt(1+((z-L)/(0.2*L))^(2)));

                 
PDE1 := diff(rho31(t, z), t) = -(1/2)*g*rho31(t, z)+(.5*I)*Omegac(z)*rho21(t, z)+(.5*I)*Omegap(t, z);
                        
PDE2 := diff(rho21(t, z), t) = (.5*I)*Omegac(z)*rho31(t, z);
                                                          
         
PDE3 := diff(Omegap(t, z), t) = -c*(diff(Omegap(t, z), z))+I*c*eta*rho31(t, z);
       

IBC := {Omegap(0, z) = 0, Omegap(t, 0) = Omegap0*exp(-((t-t0)/k)^2), rho21(0, z) = 0, rho31(0, z) = 0};
 
                              
pds := pdsolve({PDE1, PDE2, PDE3}, IBC, numeric, time = t, range = 0 .. 0.1e-2);
                      

every things is ok but I don't know how can i plot rho31, rho21 and Omegap in 2D or 3D plot.

I have tried

but i see the bellow error:

Error, (in pdsolve/numeric/plot) unable to compute solution for t>INFO["failtime"]:
unable to store -37500000000000000*I when datatype=float[8]

Please help me to solve the error?

Thanks