;;---------------------------------------
;  TRACE
;
;  **** DO NOT MODIFY ****
;
;  set parameters in plot.def then type:
;
;  @trace.idl
;
;;---------------------------------------

field1='' & field2='' & winds=''
backup_data='no' & already_data='no'  	
datafolder='./' & plotfolder='./'
topo=0
extract='no'

SPAWN, '\rm user.idl'
SPAWN, 'grep -v grep plot.def > user.idl'

;; to avoid dumb and useless warnings from IDL
space=0. & altitude=0. & lon=0. & lat=0.  
what_I_plot=0. & title_user='' & missing_value=0.

@user.idl
if (nam5 eq " interp_method = 0 ") then coord='model_level'
if ((nam5 eq " interp_method = 1 ") and (tabnam(0) lt tabnam(1))) then coord='height'

.compile call_arwpost
if (extract eq 'yes') then call_arwpost, nam1, nam2, nam3, nam4, nam5, tabnam

denom=plot+'_'+coord+'_'+field1
if ((field2 ne '') and (topo eq 0)) then denom=denom+'_'+field2
if (topo eq 1) then denom=denom+'_HGT'
if (winds(0) ne '') then denom=denom+'_'+winds(0)+winds(1)

if (already_data eq 'yes') then SPAWN, 'ln -sf  '+datafolder+'/'+coord+'.input.ctl input.ctl'
if (already_data eq 'yes') then SPAWN, 'ln -sf  '+datafolder+'/'+coord+'.input.dat input.dat'

SPAWN, '\rm param_plot.idl'
SPAWN, 'grep grepALL plot.def > param_plot.idl'
if ((plot eq 'meridional') or (plot eq 'zonal')) then SPAWN, 'grep grepSEC plot.def >> param_plot.idl'
if (plot eq 'map') then SPAWN, 'grep grepMAP plot.def >> param_plot.idl'
.compile map_latlon
.compile section

!p.charthick = 2.0
!p.thick = 3.0
!x.thick = 2.0
!y.thick = 2.0

out_wrf, $
        plot=plot, $
        field1=field1, $
        field2=field2, $
        when=1, $
        winds=winds, $
        level=level,$
        save_data=yeah,$        ;; data is loaded only once
        topo=topo,$
;        range=range,$
        colors=colors, $
        nlon=nlon, $
        nlat=nlat, $
        save_ps=denom+string(1001,'(I0)')
;
; display the first one as a sample
;
SPAWN, 'display ./'+denom+'1001.ps &'

if (num ge 2) then for i=2,num do $
	out_wrf, $
	plot=plot, $
	field1=field1, $
	field2=field2, $
	when=i, $
	winds=winds, $
	level=level,$
        save_data=yeah,$	;; data is loaded only once
	topo=topo,$
	range=range,$
	colors=colors, $
	nlon=nlon, $
	nlat=nlat, $
	save_ps=denom+string(1000+i,'(I0)')

SPAWN, 'mv -f ./'+denom+'????.ps '+plotfolder+'/'
if (backup_data eq 'yes') then SPAWN, 'cp -f namelist.ARWpost '+datafolder+'/'+coord+'.namelist.ARWpost'
if (backup_data eq 'yes') then SPAWN, 'cp -f input.ctl '+datafolder+'/'+coord+'.input.ctl'
if (backup_data eq 'yes') then SPAWN, 'cp -f input.dat '+datafolder+'/'+coord+'.input.dat'

SPAWN, 'cp -f plot.def '+plotfolder+'/'+denom+'.def'