how to create contour plots?

I have an expression

Bz1=((1j*h)/r)*((R2*exp((alpha1/r)*z))+(S2*exp((-alpha1/r)*z)))*(exp(1j*h*(((2*pi*N/60)*t)-theta)))

I have used the cylindrical coordinate system for my analysis where: Bz1: Flux density in the z direction. h:harmonic order and it will take values from a for loop (inner loop). r:radius and this will also take its values from a for loop(outer loop). R2,S2,alpha1 are calculated during the iteration. t-time period and it is given a constant value. I need to create a contour of Bz1 in the (theta,z)plane with theta=-(pi/p):0.1:(pi/p),p-pole pair and z=0:0.005:(g+e),g-air gap and e thickness of the disc. Since am struck up here , am not able to proceed my analysis please help.