;+ ; NAME: ; Routine_name ; ; PURPOSE: ; ; CATEGORY: ; ; CALLING SEQUENCE: ; ; INPUTS: ; ; OPTIONAL INPUTS: ; ; KEYWORD PARAMETERS: ; ; OUTPUTS: ; ; OPTIONAL OUTPUTS: ; ; COMMON BLOCKS: ; ; SIDE EFFECTS: ; ; RESTRICTIONS: ; ; PROCEDURE: ; ; EXAMPLE: ; ; MODIFICATION HISTORY: ; Written by: Wayne Feltz, CIMSS/SSEC ; ; $Date:$ ; $Id:$ ; $Log:$ ; ; ;- dummy='' date='' sondfile=pickfile(path='/home3/waynef/public_html/camexiii/archive/sondes/walviz/',filter='*.da*') result=str_sep(sondfile,'/') lastrec=n_elements(result)-1 file=result(lastrec) filepart=str_sep(file,'.') sdate='98'+strmid(filepart(0),0,4) stime=strmid(filepart(0),4,9) data=ddread(sondfile,offset=4,/countallrows) nl=n_elements(data(1,*)) sonde_pr=fltarr(nl) sonde_t=fltarr(nl) sonde_td=fltarr(nl) sonde_rh=fltarr(nl) sonde_mr=fltarr(nl) sonde_alt=fltarr(nl) n_columns=n_elements(data(*,0)) if n_columns eq 6 then begin sonde_pr(*)=data(1,*) sonde_alt(*)=data(2,*) sonde_t(*)=data(3,*) sonde_td(*)=data(3,*)-data(4,*) sonde_rh(*)=data(5,*) sonde_mr=float(nl) sonde_mr=wvmr(sonde_pr,sonde_t+273.13,sonde_td+273.13) endif if n_columns eq 5 then begin sonde_pr(*)=data(1,*) sonde_t(*)=data(2,*) sonde_td(*)=data(3,*) sonde_rh(*)=data(4,*) sonde_mr=float(nl) sonde_mr=wvmr(sonde_pr,sonde_t+273.13,sonde_td+273.13) endif if n_columns eq 8 then begin sonde_pr(*)=data(1,*) sonde_t(*)=data(2,*) sonde_rh(*)=data(3,*) result=where(sonde_pr gt 0.0 ,nl) sonde_pr=sonde_pr(result) sonde_t=sonde_t(result) sonde_rh=sonde_rh(result) sonde_mr=fltarr(nl) sonde_td=fltarr(nl) sonde_mr(*)=mrrh(sonde_pr(*),sonde_t(*),sonde_rh(*)) sonde_td(*)=dewpt(sonde_pr(*),sonde_t(*),sonde_mr(*),nl) endif if n_columns eq 6 then begin stability_indices2_sondes,nl, reverse(sonde_pr), reverse(sonde_alt), reverse(sonde_t), reverse(sonde_td), reverse(sonde_mr), li, totalt, kindex,$ pres_lcl,sw,pres_lfc, pres_eq, cape, cin endif title='Wallops VIZ Sounding '+sdate+' '+stime+' UTC' skewt, [-30, 40], everyT=10, everyDA=10, everySA=2, everyW=2,title=title plot_skewt, sonde_t(*),sonde_td(*),sonde_pr(*),0 if n_columns eq 6 then begin xyouts,.02, .96, 'CAPE = '+strmid(strtrim(string(cape),1),0,4), color=0,/normal,charsize=.7 xyouts,.02, .93, 'CIN = '+strmid(strtrim(string(cin),1),0,4), color=0,/normal,charsize=.7 xyouts,.02, .90, 'K Index = '+strmid(strtrim(string(kindex),1),0,2), color=0,/normal,charsize=.7 xyouts,.02, .87, 'LI = '+strmid(strtrim(string(li),1),0,5), color=0,/normal,charsize=.7 xyouts,.02, .84, 'Showalter = '+strmid(strtrim(string(sw),1),0,4), color=0,/normal,charsize=.7 xyouts,.02, .81, 'LCL = '+strmid(strtrim(string(pres_lcl),1),0,3)+' mb', color=0,/normal,charsize=.7 xyouts,.02, .78, 'LFC = '+strmid(strtrim(string(pres_lfc),1),0,3)+' mb', color=0,/normal,charsize=.7 xyouts,.02, .75, 'EL = '+strmid(strtrim(string(pres_eq),1),0,3)+' mb', color=0,/normal,charsize=.7 endif tvlct,r,g,b,/Get giffile=tvrd() write_gif, '~/public_html/camexiii/archive/quicklooks/'+sdate+'/'+sdate+'_'+stime+'_walviz_sonde.gif', giffile, r, g, b end