pro field_map,pv_arr,t_arr,m_arr,lons,lats,theta,trop_hgt,lat00,lon00,xrdf=xrdf,no_trop_hgt=no_trop_hgt
     
; map  mpv,T, trop height and rdf fields
; rdf is normally not done on the call
; /xrdf means that the rdf is done.

pvlevels=findgen(20)*.15
if lat00 gt 45. then pvlevels=findgen(25)*.1+1.0
if keyword_set(xrdf) then goto,xxx
    contour,pv_arr,/cell_fill,lons,lats,/over,levels=pvlevels
 ; temperature     
     tnat=nat_theta(10.e-9,4.e-6,theta,type=1)
     tice=frost_theta(4.e-6,theta)
     contour,t_arr,/follow,lons,lats,/over,levels=[185,190,192.5,195,196,197,198,200,205,210,220,230],c_labels=[1,1,1,1,1,1,1,1,1],c_col=0,thick=2
     contour,t_arr,/follow,lons,lats,/over,levels=[tice,tnat],c_col=0,c_line=1,c_thick=3,c_annotat=['ICE','NAT']
; streamlines
     m_arr=m_arr*1.e-5
     mlevels=findgen(15)*(max(m_arr)-min(m_arr))/15.+min(m_arr)
     contour,m_arr,/follow,lons,lats,/over,levels=mlevels,c_col=0,c_line=2,c_labels=0*mlevels
     range_ring,lat00,lon00,3000,360,bearing,latp,lonp
     oplot,lonp,latp,col=75
; trop_ height
if keyword_set(no_trop_hgt) then return
     if max(trop_hgt) ne min(trop_hgt) then begin
    
     contour,trop_hgt*1.0e-3,/follow,lons,lats,/over,$
        levels=[8.,10.,12.,14.],c_labels=[1,1,1,1],c_col=180,thick=2
     
     endif
xxx:
     
return
end