startday=today(3)
starttime=0.0
ndays=7
forward=-1
thetav=480
source='NMC'
; generate grid
latgrid=findgen(59)+30.
longrid=findgen(355)
latin=[0]
lonin=[0]
for j=0,n_elements(latgrid)-1 do begin
   lonin=[lonin,longrid]
   latin=[latin,replicate(latgrid(j),n_elements(longrid))]
   endfor
latin=latin(1:*)
lonin=lonin(1:*)

trajdriver,startday=startday,starttime=starttime,ndays=ndays,latin=latin,$
 lonin=lonin,theta=thetav,forward=forward,/noxdr,$
 latsave=latsave,lonsave=lonsave,psave=psave,temsave=temsave,tsave=tsave,$
 qsave=qsave,source='ASM',save_memory=50,grid='GG1X1',forecast='F-FLK-01'

; now resort the results into the right grid points
s=size(qsave)
pv_rdf=fltarr(n_elements(longrid),n_elements(latgrid))
n=where(pv_rdf eq 0.)
q1=qsave(*,s(2)-1)*1.0e5
sn=n_elements(latgrid)
sb=n_elements(longrid)
for j=0,sn-1 do pv_rdf(0,j)=q1(j*sb:(j+1)*sb-1)
dir='/misc/mrs13/'
file='rdf_'+str(thetav)+'_'+startday+'.ps
tvzinit
lat0=90.
lon0=0.
proj=1
tgcoltab,'spectrum2a'
viewgraph
font
scl=[0,max(pv_rdf)]
putmap,pv_rdf,latgrid,longrid,proj,lat0,lon0,/wrap,scl
gettheta,startday,'T',thetav,temp_data,lat=lats,lon=lons,source=source,hemis=0
putmap,temp_data,lats,lons,proj,lat0,lon0,/wrap,/contour,c_levels=180.+findgen(20)*10,c_col=0
putmap,temp_data,lats,lons,proj,lat0,lon0,/wrap,/contour,c_levels=[188,195],c_thick=5,c_col=0
colbar,'PV ',scl,[450,80,20,300]
xyouts,/norm,.2,.9,'7 day PV RDF and T, '+str(thetav)+'K, '+nicedate(startday),size=1.2
spawn,'cp /app/netpub/www/httpd/htdocs/916/Personnel/people/Schoeberl,_Mark_R./rdf.gif /app/netpub/www/httpd/htdocs/916/Personnel/people/Schoeberl,_Mark_R./rdf2.gif
make_gif,'/app/netpub/www/httpd/htdocs/916/Personnel/people/Schoeberl,_Mark_R./rdf.gif'
; store the gif files
dayno=julday(startday)
if dayno lt 10 then stdayno='00'+str(dayno)
if dayno ge 10 and dayno lt 100 then stdayno='0'+str(dayno)
if dayno ge  100 then stdayno=str(dayno)
print,stdayno
spawn,'cp /app/netpub/www/httpd/htdocs/916/Personnel/people/Schoeberl,_Mark_R./rdf.gif /misc/mrs03/rdf_nh'+strmid(startday,0,2)+stdayno+'.gif'
end