PRO nv

@compile_opt.pro	;

dirs = findallsubdirs(/recursive, count=count)

; Don't process the output directories

i = where(strpos(dirs,'ipsd_vuinsitu_output') NE -1,complement=j)
dirs = dirs[j]

i = where(strpos(dirs,'ipsdt_vuinsitu_output') NE -1,complement=j)
dirs = dirs[j]

body = jpl_body(/mars)

count = n_elements(dirs)

cd, current=current
FOR i=0,count-1 DO BEGIN

	cd, dirs[i]
	files = file_search('nv3d*.*',count=nn)

	FOR j=0,nn-1 DO BEGIN

		tmp = GetFileSpec(files[j],from='name',upto='type')

		hdr = vu_select(files[j],/read, ff=ff)
		vu_insitu, hdr, ff, body=body, /noplot, ascii='../ipsd_vuinsitu_output/'+tmp+'.txt_new'

	ENDFOR

ENDFOR


FOR i=0,count-1 DO BEGIN

	cd, dirs[i]
	files = file_search('nv3f*.*',count=nn)

	IF nn NE 0 THEN BEGIN

		tmp = dirs[i]
		tmp = strmid(tmp,0,strlen(tmp)-1)		; Remove trialing \
		tmp = GetFileSpec(tmp, part='name')		; Subdirectory name

		hdr = vu_select(files,/read, ff=ff)
		vu_insitu, hdr, ff, body=body, /noplot, ascii='../ipsdt_vuinsitu_output/'+tmp+'.txt_new'

	ENDIF
	
ENDFOR

cd, current

RETURN  &  END