INPUT:
rmin= minimum R
rmax= maximum R
nrs= grid in R
zmin= minimum z
zmax= maximum z
nzs= grid in z
ncontours= number of contours
justcontours= (False) if True, just plot contours
savefilename= save to or restore from this savefile (pickle)
log= if True, plot the log density