;.rnew read_rs_downs.pro ; Sebastian Hoch - April 2020 ;;read_netCDF,'/Users/shoch/CFOG/DATA/RS_DOWNS/netCDF/CFOG_Ferryland_20180906_085834.nc',netdat,att,st ;filename='/Users/shoch/CFOG/DATA/RS_DOWNS/netCDF/CFOG_Ferryland_20180906_085834.nc' ;;------------------------------------------------ ;; Open netcdf file and get an ID for file (ncid) ;;------------------------------------------------ ;ncid=ncdf_open(filename,/nowrite) ;;------------------------------------------------ ;; Find the variable ID ;;------------------------------------------------ ;varid=ncdf_varid(ncid,var_name) ;;------------------------------------------------ ;; Get variable ;;------------------------------------------------ ;ncdf_varget,ncid,varid,var ;;------------------------------------------------ ;; close netcdf file and end routine ;;------------------------------------------------ ;ncdf_close,ncid output='/Users/shoch/CFOG/DATA/RS_DOWNS/PLOTS/cfog_soundings_downs.ps' outfolder='/Users/shoch/CFOG/DATA/RS_DOWNS/DATA-OUT/' ;/Users/shoch/bin/GruanToolRs92/CFOG_Ferryland_20181006_051758_SynchronizedSoundingData.csv prefix='CFOG_Ferryland_' endfix='_SynchronizedSoundingData.csv' folder='/Users/shoch/bin/GruanToolRs92/' cd, folder sstr=FILE_SEARCH( prefix+'2018*'+endfix) & files=sstr(*) & nfiles=n_elements(files) ;sstr=FILE_SEARCH( prefix+'20180928*'+endfix) & files=sstr(*) & nfiles=n_elements(files) ;nfiles=10 reads,STRMID(sstr(*),strlen(prefix)+00,4),yy,format='(i4.4)' & yystr=STRMID(sstr(*),strlen(prefix)+00,4) reads,STRMID(sstr(*),strlen(prefix)+04,2),mo,format='(i4.4)' & mostr=STRMID(sstr(*),strlen(prefix)+04,2) reads,STRMID(sstr(*),strlen(prefix)+06,2),da,format='(i4.4)' & dastr=STRMID(sstr(*),strlen(prefix)+06,2) reads,STRMID(sstr(*),strlen(prefix)+09,2),hr,format='(i4.4)' & hrstr=STRMID(sstr(*),strlen(prefix)+09,2) reads,STRMID(sstr(*),strlen(prefix)+11,2),mi,format='(i4.4)' & mistr=STRMID(sstr(*),strlen(prefix)+11,2) reads,STRMID(sstr(*),strlen(prefix)+12,2),se,format='(i4.4)' & sestr=STRMID(sstr(*),strlen(prefix)+12,2) jultime_UTC=julday(mo(*),da(*),yy(*),hr(*),mi(*),se(*)) datetimestr=strarr(nfiles) for i=0,nfiles-1 do datetimestr(i)=yystr(i)+'/'+mostr(i)+'/'+dastr(i)+' '+hrstr(i)+''+mistr(i)+' UTC' rs_id=datetimestr(*) rs_id_UTC=STRMID(rs_id,11,8) rs_id_DAY_UTC=STRMID(rs_id,8,2)+'/'+STRMID(rs_id,11,8) nlines=3000 sdatarr=fltarr(nfiles,20,nlines) & sdatarr(*)=!VALUES.F_NAN ; SONDDATA(0,*) : SEC NACH START ; SONDDATA(1,*) : HOEHE KILO-!!!-METER ; SONDDATA(2,*) : PRESSURE ; SONDDATA(3,*) : TEMPERATURE ; SONDDATA(4,*) : RELATIVE HUMIDITY ; SONDDATA(5,*) : DEWPOINT TEMPERATURE ; SONDDATA(6,*) : WINDDIRECTION -/? s.o./ COS ; SONDDATA(7,*) : WINDSPEED ; SONDDATA(8,*) : q ; SONDDATA(9,*) : MISCHUNGSVERHAELTNIS for snr=0,nfiles-1 do begin d=read_csv(sstr(snr),HEADER=RSHeader) print,RSHeader ; 01 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 ;SoundingIdPk RadioRxTimePk [s] DataSrvTime Pressure [hPa] Temperature Humidity [%] WindDir [deg] WindSpeed [m/s] WindNorth [m/s] WindEast [m/s] Height [gpm] PtuStatus WindInterpolated Latitude [deg] Longitude [deg] North [m] East [m] Up [m] Altitude [m] Dropping ; FIELD01 STRING Array[976] ; FIELD02 DOUBLE Array[976] ; FIELD03 STRING Array[976] ; FIELD04 DOUBLE Array[976] ; FIELD05 DOUBLE Array[976] ; FIELD06 DOUBLE Array[976] ; FIELD07 DOUBLE Array[976] ; FIELD08 DOUBLE Array[976] ; FIELD09 DOUBLE Array[976] ; FIELD10 DOUBLE Array[976] ; FIELD11 DOUBLE Array[976] ; FIELD12 LONG Array[976] ; FIELD13 STRING Array[976] ; FIELD14 DOUBLE Array[976] ; FIELD15 DOUBLE Array[976] ; FIELD16 DOUBLE Array[976] ; FIELD17 DOUBLE Array[976] ; FIELD18 DOUBLE Array[976] ; FIELD19 DOUBLE Array[976] ; FIELD20 LONG Array[976] ;FUNCTION HUMIDITY_CONVERSIONS, inputstring, INVAR, P, TK humarr=HUMIDITY_CONVERSIONS('RH', d.field06, d.field04, d.field05) qv1=humarr(3,*) index=where(d.field20 eq 0,cnt) linesfile=n_elements(d.field01(index)) if linesfile ge nlines then begin sdatarr(snr,0,*)=d.field02(0:nlines-1)-d.field02(0) sdatarr(snr,1,*)=d.field11(0:nlines-1)/1000. sdatarr(snr,2,*)=d.field04(0:nlines-1) sdatarr(snr,3,*)=d.field05(0:nlines-1)-273.15 sdatarr(snr,4,*)=d.field06(0:nlines-1) sdatarr(snr,5,*)=humarr(1, 0:nlines-1)-273.15 sdatarr(snr,6,*)=d.field07(0:nlines-1) sdatarr(snr,7,*)=d.field08(0:nlines-1) ;sdatarr(snr,8,*)=humarr(3, 0:nlines-1) sdatarr(snr,9,*)=humarr(3, 0:nlines-1) endif else begin lin=linesfile sdatarr(snr,0,0:lin-1)=d.field02(index)-d.field02(0) sdatarr(snr,1,0:lin-1)=d.field11(index)/1000. sdatarr(snr,2,0:lin-1)=d.field04(index) sdatarr(snr,3,0:lin-1)=d.field05(index)-273.15 sdatarr(snr,4,0:lin-1)=d.field06(index) sdatarr(snr,5,0:lin-1)=humarr(1, index) -273.15 sdatarr(snr,6,0:lin-1)=d.field07(index) sdatarr(snr,7,0:lin-1)=d.field08(index) ;sdatarr(snr,8,0:lin-1)=humarr(1, *) sdatarr(snr,9,0:lin-1)=humarr(3, index) endelse endfor PRINT, 'READ ALL DATA' ; PLOTTING ;------------------------------------------------------ for snr=0,nfiles-1 do begin close,1 openw,1,outfolder+'RS_CFOG_DOWNS_'+STRMID(sstr(snr),15,15)+'.dat' ;CFOG_Ferryland_20181006_051758_SynchronizedSoundingData.csv ; SONDDATA(0,*) : SEC NACH START ; SONDDATA(1,*) : HOEHE KILO-!!!-METER ; SONDDATA(2,*) : PRESSURE ; SONDDATA(3,*) : TEMPERATURE ; SONDDATA(4,*) : RELATIVE HUMIDITY ; SONDDATA(5,*) : DEWPOINT TEMPERATURE ; SONDDATA(6,*) : WINDDIRECTION -/? s.o./ COS ; SONDDATA(7,*) : WINDSPEED ; SONDDATA(8,*) : q ; SONDDATA(9,*) : MISCHUNGSVERHAELTNIS outvars=[0,1,2,3,4,5,6,7,9] printf,1,'SECOND[s] ALT[km] P[hPa] T['+STRING("260B)+'C] RH[%] DPT['+STRING("260B)+'C] WDIR['+STRING("260B)+'] WSP[m/s] MR[kg/kg]' printf,1,sdatarr(snr,outvars,*),format='(f9.1,f9.4,f9.2,f9.2,f7.1,f9.2,f7.1,f9.2,f10.3)' close,1 endfor ; PLOTTING ;------------------------------------------------------ set_plot, 'ps' device, /color, portrait=1 device, filename= output device, xsize=21-2, ysize=27 device, xoffset=1, yoffset=1 !p.font=1 ;!p.font=0 ;device, set_font='Helvetica', /bold,/tt_font , font_size=18 ;------------------------------------------------------ ;device, set_font='Helvetica', /tt_font , font_size=28 device, font_size=18 A = FINDGEN(17) * (!PI*2/16.) & USERSYM, COS(A), SIN(A), /FILL !p.thick=3 !x.thick=3 !y.thick=3 xm1=10 & xm2=10 ym1= 4 & ym2= 6 ym1= 2 & ym2= 2 xm1= 8 & xm2=13 ;ym1= 4 & ym2= 6 ;------------------------------------------------------------------------------------------------ ; Red Green Blue DEF TVLCT, 0, 0, 0, 1 ; black, schwarz ; TVLCT, 127, 0, 0, 2 ; brown, braun ; TVLCT, 255, 0, 0, 3 ; red, rot ; TVLCT, 255, 127, 0, 4 ; lightorange, helles orange ; TVLCT, 255, 255, 0, 5 ; yellow, gelb ; TVLCT, 0, 255, 0, 6 ; green, gruen ; TVLCT, 0, 255, 255, 7 ; lightblue, hellblau ; TVLCT, 0, 0, 255, 8 ; blue, blau ; TVLCT, 255, 255, 255, 10 ; white, weiss ; ;------------------------------------------------------------------------------------------------ ;TVLCT, 255, 0, 255, 9 ; pink TVLCT, 128, 0, 128, 9 ; purple ;------------------------------------------------------ tvlct, 0, 0, 128, 100 ; darkblue tvlct, 0, 0, 255, 101 ; blue tvlct, 0, 128, 255, 102 ; lightblue tvlct, 0, 192, 255, 103 ; tvlct, 0, 255, 255, 104 ; cyan tvlct, 255, 255, 0, 105 ; yellow tvlct, 255, 192, 0, 106 ; tvlct, 255, 128, 0, 107 ; orange tvlct, 255, 92, 0, 108 ; orange tvlct, 255, 0, 0, 109 ; red tvlct, 192, 0, 0, 110 ; darkred tvlct, 128, 0, 0, 111 ; darkred tvlct, 96, 0, 0, 112 ; darkred tvlct, 64, 0, 0, 113 ; darkred tvlct, 32, 0, 0, 114 ; darkred tvlct, 0, 0, 0, 115 ; darkred ;------------------------------------------------------ TVLCT, 41, 41, 41, 20 ; grey ; TVLCT, 51, 51, 51, 21 ; light grey ; TVLCT, 255-30, 255-30, 255-30, 21 ; light grey ; TVLCT, 255-15, 255-15, 255-15, 20 ; light grey ; ;------------------------------------------------------ ;------------------------------------------------------ nullstrarr=strarr(20) nullstrarr(*)=' ' !P.MULTI=[0,4,5] for snr=0,nfiles-1 do begin ;------------------------------------------------------ ;------------------------------------------------------ print,sstr(snr) comparr=indgen(1) comparr(0)=snr ym1= 2 & ym2=4 xm1= 4 & xm2=4 altitude= 0.032 max_alt = 3.2;2.0 hlkmas =max_alt col=[0,1,2,3,4,5,6,7,8] rawlines=100 index=where(sdatarr(comparr(0),1,*) gt (altitude + hlkmas),count) if count eq 0 then begin ;index= [0,rawlines -2] index= [ rawlines -2] endif ; Plotting of temperature and dewpoint ;------------------------------------- plot, sdatarr(comparr(0),3,*),sdatarr(comparr(0),1,*),/nodata, $ ;ytitle='Height above sea level (km)', $ ytitle='Height (km a.s.l.)', $ yrange=[altitude, max_alt], $ ystyle=1+8-8, $ xtitle='T and T!Dd!N ('+STRING("260B)+'C)' , $ ;xtitle='T ('+STRING("260B)+'C)' , $ xmargin=[xm1,xm2], $ ymargin=[ym1,ym2], $ title= exp_name,charsize=1.3,$ ;xrange=[6,19],xstyle=1,xtickinterval=2,xminor=2;,$ ;xrange=[-80.,max(sdatarr(comparr(0),3,*))] xrange=[min(sdatarr(comparr(*),3,0:index(0))), $ max(sdatarr(comparr(*),3,0:index(0)))] ;axis,yaxis=1,yrange=(!Y.CRANGE-altitude), $ ; ystyle=1, ycharsize=1.3, $ ; ytitle='Height above ground (km)' for i=0,n_elements(comparr)-1 do begin oplot, sdatarr(comparr(i),3,*),sdatarr(comparr(i),1,*),color=col(i),psym=-3 oplot, sdatarr(comparr(i),5,*),sdatarr(comparr(i),1,*),color=col(i),psym=-3,line=1,NSUM=10 endfor xshift= 1 & yshift= -1.2 for i=0,n_elements(comparr)-1 do begin xyouts,!x.crange(0) + (xshift ) *(!x.crange(1)-!x.crange(0))/15.,$ !y.crange(1) - ((i+yshift)+1./4) *(!y.crange(1)-!y.crange(0))/20.,$ strcompress(string(comparr(i),format='(i2)')) +' '+rs_id(comparr(i)), color= col(i),charsize=0.65 endfor ; Plotting of relative humidity ;-------------------------------- plot, sdatarr(comparr(0),4,*),sdatarr(comparr(0),1,*),/nodata, $ ;ytitle='Height above sea level (km)', $ ytitle='Height (km a.s.l.)', $ yrange=[altitude, max_alt], $ ystyle=1+8-8, $ xmargin=[xm1,xm2], $ ymargin=[ym1,ym2], $ xtitle='RH (%)', $ title= sitestring,charsize=1.3,$ xrange=[20,100] ;axis,yaxis=1,yrange=(!Y.CRANGE-altitude), $ ; ystyle=1, ycharsize=1.3, $ ; ytitle='Height above ground (km)' for i=0,n_elements(comparr)-1 do begin oplot, sdatarr(comparr(i),4,*),sdatarr(comparr(i),1,*),color=col(i),psym=-3 endfor ; Plotting of windspeed ;-------------------------------- plot, sdatarr(comparr(0),7,*),sdatarr(comparr(0),1,*),/nodata, $ ;ytitle='Height above sea level (km)', $ ytitle='Height (km a.s.l)', $ yrange=[altitude, max_alt], $ ystyle=1+8-8, $ xmargin=[xm1,xm2], $ ymargin=[ym1,ym2], $ xtitle='WS (m s!U-1!N)', $ title= location,charsize=1.3 ;axis,yaxis=1,yrange=(!Y.CRANGE-altitude), $ ; ystyle=1, ycharsize=1.3, $ ; ytitle='Height above ground (km)' for i=0,n_elements(comparr)-1 do begin oplot, sdatarr(comparr(i),7,*),sdatarr(comparr(i),1,*),color=col(i),psym=-3 endfor ; Plotting of wind direction ;-------------------------------- plot, sdatarr(comparr(0),6,*),sdatarr(comparr(0),1,*),psym=3,/nodata, $ ;ytitle='Height above sea level (km)', $ ytitle='Height (km a.s.l.)', $ yrange=[altitude, max_alt], $ ystyle=1+8-8, $ xmargin=[xm1,xm2], $ ymargin=[ym1,ym2], $ xtitle='WD ('+STRING("260B)+')', $ xrange=[0,360],xstyle = 1, xticks=4,xminor=3, $ title= filename,charsize=1.3 ;axis,yaxis=1,yrange=(!Y.CRANGE-altitude), $ ; ystyle=1, ycharsize=1.3, $ ; ytitle='Height above ground (km)' for i=0,n_elements(comparr)-1 do begin oplot, sdatarr(comparr(i),6,*),sdatarr(comparr(i),1,*),color=col(i),psym=3 endfor ;------------------------------------------------------ ;------------------------------------------------------ endfor ;------------------------------------------------------ ;------------------------------------------------------ ;------------------------------------------------------ ;EVENT 2 / 12 Sept ;------------------------------------------------------ !P.MULTI=[0,4,1] ;------------------------------------------------------ ;PLOTTING COMPS ;------------------------------------------------------ comparr=[47,49,51,54,58,61,64,66,68] col= [1,2,3,4,5,6,7,8,9] comparr=[18,19] col= [3,8] ym1= 2 & ym2=4 xm1= 4 & xm2=4 altitude= 0.032 max_alt = 3.2;2.0 max_alt = 5.2;2.0 hlkmas =max_alt rawlines=100 index=where(sdatarr(comparr(0),1,*) gt (altitude + hlkmas),count) if count eq 0 then begin ;index= [0,rawlines -2] index= [ rawlines -2] endif ; Plotting of temperature and dewpoint ;------------------------------------- plot, sdatarr(comparr(0),3,*),sdatarr(comparr(0),1,*),/nodata, $ position=[0.1,0.6,0.3,0.9],$ ;ytitle='Height above sea level (km)', $ ytitle='Height (km a.s.l.)', $ yrange=[altitude, max_alt], $ ystyle=1+8-8, $ xtitle='T and T!Dd!N ('+STRING("260B)+'C)' , $ ;xtitle='T ('+STRING("260B)+'C)' , $ xmargin=[xm1,xm2], $ ymargin=[ym1,ym2], $ title= exp_name,charsize=1.3,$ ;xrange=[-80.,max(sdatarr(comparr(0),3,*))] xrange=[5-15-5,19],xstyle=1,xtickinterval=5,xminor=5;,$ ;xrange=[min(sdatarr(comparr(*),3,0:index(0))), $ ; max(sdatarr(comparr(*),3,0:index(0)))-0.1] ;axis,yaxis=1,yrange=(!Y.CRANGE-altitude), $ ; ystyle=1, ycharsize=1.3, $ ; ytitle='Height above ground (km)' for i=0,n_elements(comparr)-1 do begin oplot, sdatarr(comparr(i),3,*),sdatarr(comparr(i),1,*),color=col(i),psym=-3 oplot, sdatarr(comparr(i),5,*),sdatarr(comparr(i),1,*),color=col(i),psym=-3,line=1,NSUM=10 endfor xshift= 1 & yshift= -1.2 for i=0,n_elements(comparr)-1 do begin ;xyouts,!x.crange(0) + (xshift ) *(!x.crange(1)-!x.crange(0))/15.,$ ; !y.crange(1) - ((i+yshift)+1./4) *(!y.crange(1)-!y.crange(0))/20.,$ ; ;strcompress(string(comparr(i),format='(i2)')) +' '+rs_id(comparr(i)), color= col(i),charsize=0.65 ; ;rs_id(comparr(i)), color= col(i),charsize=0.65 ; rs_id_UTC(comparr(i)), color= col(i),charsize=0.65 endfor ; Plotting of relative humidity ;-------------------------------- plot, sdatarr(comparr(0),4,*),sdatarr(comparr(0),1,*),/nodata, $ position=[0.3,0.6,0.5,0.9],ytickname=nullstrarr,$ ;ytitle='Height above sea level (km)', $ ;ytitle='Height (km ASL)', $ yrange=[altitude, max_alt], $ ystyle=1+8-8, $ xmargin=[xm1,xm2], $ ymargin=[ym1,ym2], $ xtitle='RH (%)', $ title= sitestring,charsize=1.3,$ xrange=[20,100] ;axis,yaxis=1,yrange=(!Y.CRANGE-altitude), $ ; ystyle=1, ycharsize=1.3, $ ; ytitle='Height above ground (km)' for i=0,n_elements(comparr)-1 do begin oplot, sdatarr(comparr(i),4,*),sdatarr(comparr(i),1,*),color=col(i),psym=-3 endfor ; Plotting of windspeed ;-------------------------------- plot, sdatarr(comparr(0),7,*),sdatarr(comparr(0),1,*),/nodata, $ position=[0.5,0.6,0.7,0.9],ytickname=nullstrarr,$ xrange=[0,14.99+5],xstyle=1,xtickinterval=5,xminor=5,$ ;ytitle='Height above sea level (km)', $ ;ytitle='Height (km ASL)', $ yrange=[altitude, max_alt], $ ystyle=1+8-8, $ xmargin=[xm1,xm2], $ ymargin=[ym1,ym2], $ xtitle='WS (m s!U-1!N)', $ title= location,charsize=1.3 ;axis,yaxis=1,yrange=(!Y.CRANGE-altitude), $ ; ystyle=1, ycharsize=1.3, $ ; ytitle='Height above ground (km)' for i=0,n_elements(comparr)-1 do begin oplot, sdatarr(comparr(i),7,*),sdatarr(comparr(i),1,*),color=col(i),psym=-3 endfor ; Plotting of wind direction ;-------------------------------- plot, sdatarr(comparr(0),6,*),sdatarr(comparr(0),1,*),psym=3,/nodata, $ position=[0.7,0.6,0.9,0.9],ytickname=nullstrarr,$ ;ytitle='Height above sea level (km)', $ ;ytitle='Height (km ASL)', $ yrange=[altitude, max_alt], $ ystyle=1+8-8, $ xmargin=[xm1,xm2], $ ymargin=[ym1,ym2], $ xtitle='WD ('+STRING("260B)+')', $ xrange=[0,360],xstyle = 1, xticks=4,xminor=3, $ title= filename,charsize=1.3 ;axis,yaxis=1,yrange=(!Y.CRANGE-altitude), $ ; ystyle=1, ycharsize=1.3, $ ; ytitle='Height above ground (km)' for i=0,n_elements(comparr)-1 do begin oplot, sdatarr(comparr(i),6,*),sdatarr(comparr(i),1,*),color=col(i),psym=3 endfor xshift= 15.25 & yshift= 1 for i=0,n_elements(comparr)-1 do begin xyouts,!x.crange(0) + (xshift ) *(!x.crange(1)-!x.crange(0))/15.,$ !y.crange(1) - ((i+yshift)+1./4) *(!y.crange(1)-!y.crange(0))/15.,$ ;strcompress(string(comparr(i),format='(i2)')) +' '+rs_id(comparr(i)), color= col(i),charsize=0.65 ;rs_id(comparr(i)), color= col(i),charsize=0.65 rs_id_DAY_UTC(comparr(i)), color= col(i),charsize=0.65 endfor ;------------------------------------------------------ ;EVENT 11 / SUPER-IOP ;------------------------------------------------------ !P.MULTI=[0,4,1] ;------------------------------------------------------ ;PLOTTING COMPS ;------------------------------------------------------ comparr=[47,49,51,54,58,61,64,66,68] ;comparr=[2,3] col= [1,2,3,4,5,6,7,8,9] ym1= 2 & ym2=4 xm1= 4 & xm2=4 altitude= 0.032 max_alt = 3.2;2.0 hlkmas =max_alt rawlines=100 index=where(sdatarr(comparr(0),1,*) gt (altitude + hlkmas),count) if count eq 0 then begin ;index= [0,rawlines -2] index= [ rawlines -2] endif ; Plotting of temperature and dewpoint ;------------------------------------- plot, sdatarr(comparr(0),3,*),sdatarr(comparr(0),1,*),/nodata, $ position=[0.1,0.6,0.3,0.9],$ ;ytitle='Height above sea level (km)', $ ytitle='Height (km ASL)', $ yrange=[altitude, max_alt], $ ystyle=1+8-8, $ ;xtitle='T and T!Dd!N ('+STRING("260B)+'C)' , $ xtitle='T ('+STRING("260B)+'C)' , $ xmargin=[xm1,xm2], $ ymargin=[ym1,ym2], $ title= exp_name,charsize=1.3,$ ;xrange=[-80.,max(sdatarr(comparr(0),3,*))] xrange=[5,19],xstyle=1,xtickinterval=5,xminor=5;,$ ;xrange=[min(sdatarr(comparr(*),3,0:index(0))), $ ; max(sdatarr(comparr(*),3,0:index(0)))-0.1] ;axis,yaxis=1,yrange=(!Y.CRANGE-altitude), $ ; ystyle=1, ycharsize=1.3, $ ; ytitle='Height above ground (km)' for i=0,n_elements(comparr)-1 do begin oplot, sdatarr(comparr(i),3,*),sdatarr(comparr(i),1,*),color=col(i),psym=-3 ; oplot, sdatarr(comparr(i),5,*),sdatarr(comparr(i),1,*),color=col(i),psym=-3,line=1,NSUM=10 endfor xshift= 1 & yshift= -1.2 for i=0,n_elements(comparr)-1 do begin ;xyouts,!x.crange(0) + (xshift ) *(!x.crange(1)-!x.crange(0))/15.,$ ; !y.crange(1) - ((i+yshift)+1./4) *(!y.crange(1)-!y.crange(0))/20.,$ ; ;strcompress(string(comparr(i),format='(i2)')) +' '+rs_id(comparr(i)), color= col(i),charsize=0.65 ; ;rs_id(comparr(i)), color= col(i),charsize=0.65 ; rs_id_UTC(comparr(i)), color= col(i),charsize=0.65 endfor ; Plotting of relative humidity ;-------------------------------- plot, sdatarr(comparr(0),4,*),sdatarr(comparr(0),1,*),/nodata, $ position=[0.3,0.6,0.5,0.9],ytickname=nullstrarr,$ ;ytitle='Height above sea level (km)', $ ;ytitle='Height (km ASL)', $ yrange=[altitude, max_alt], $ ystyle=1+8-8, $ xmargin=[xm1,xm2], $ ymargin=[ym1,ym2], $ xtitle='RH (%)', $ title= sitestring,charsize=1.3,$ xrange=[20,100] ;axis,yaxis=1,yrange=(!Y.CRANGE-altitude), $ ; ystyle=1, ycharsize=1.3, $ ; ytitle='Height above ground (km)' for i=0,n_elements(comparr)-1 do begin oplot, sdatarr(comparr(i),4,*),sdatarr(comparr(i),1,*),color=col(i),psym=-3 endfor ; Plotting of windspeed ;-------------------------------- plot, sdatarr(comparr(0),7,*),sdatarr(comparr(0),1,*),/nodata, $ position=[0.5,0.6,0.7,0.9],ytickname=nullstrarr,$ xrange=[0,14.99+5],xstyle=1,xtickinterval=5,xminor=5,$ ;ytitle='Height above sea level (km)', $ ;ytitle='Height (km ASL)', $ yrange=[altitude, max_alt], $ ystyle=1+8-8, $ xmargin=[xm1,xm2], $ ymargin=[ym1,ym2], $ xtitle='WS (m/s)', $ title= location,charsize=1.3 ;axis,yaxis=1,yrange=(!Y.CRANGE-altitude), $ ; ystyle=1, ycharsize=1.3, $ ; ytitle='Height above ground (km)' for i=0,n_elements(comparr)-1 do begin oplot, sdatarr(comparr(i),7,*),sdatarr(comparr(i),1,*),color=col(i),psym=-3 endfor ; Plotting of wind direction ;-------------------------------- plot, sdatarr(comparr(0),6,*),sdatarr(comparr(0),1,*),psym=3,/nodata, $ position=[0.7,0.6,0.9,0.9],ytickname=nullstrarr,$ ;ytitle='Height above sea level (km)', $ ;ytitle='Height (km ASL)', $ yrange=[altitude, max_alt], $ ystyle=1+8-8, $ xmargin=[xm1,xm2], $ ymargin=[ym1,ym2], $ xtitle='WD ('+STRING("260B)+')', $ xrange=[0,360],xstyle = 1, xticks=4,xminor=3, $ title= filename,charsize=1.3 ;axis,yaxis=1,yrange=(!Y.CRANGE-altitude), $ ; ystyle=1, ycharsize=1.3, $ ; ytitle='Height above ground (km)' for i=0,n_elements(comparr)-1 do begin oplot, sdatarr(comparr(i),6,*),sdatarr(comparr(i),1,*),color=col(i),psym=3 endfor xshift= 15.25 & yshift= 1 for i=0,n_elements(comparr)-1 do begin xyouts,!x.crange(0) + (xshift ) *(!x.crange(1)-!x.crange(0))/15.,$ !y.crange(1) - ((i+yshift)+1./4) *(!y.crange(1)-!y.crange(0))/15.,$ ;strcompress(string(comparr(i),format='(i2)')) +' '+rs_id(comparr(i)), color= col(i),charsize=0.65 ;rs_id(comparr(i)), color= col(i),charsize=0.65 rs_id_DAY_UTC(comparr(i)), color= col(i),charsize=0.65 endfor ;------------------------------------------------------ ; PAPER CLIVE ;------------------------------------------------------ !P.MULTI=[0,4,1] ;------------------------------------------------------ ;PLOTTING COMPS ;------------------------------------------------------ comparr=[53,55,58,59,61,63,65,66,67] ;comparr=[ 55,58,59,61,63,65,66,67] ;comparr=[2,3] col= [1,2,3,4,5,6,7,8,9] lin= [1,0,0,0,0,0,0,0,0] nnsum= [15,1,1,1,1,1,1,1,1] ym1= 2 & ym2=4 xm1= 4 & xm2=4 altitude= 0.032 max_alt = 3.2;2.0 hlkmas =max_alt rawlines=100 index=where(sdatarr(comparr(0),1,*) gt (altitude + hlkmas),count) if count eq 0 then begin ;index= [0,rawlines -2] index= [ rawlines -2] endif ; Plotting of temperature and dewpoint ;------------------------------------- plot, sdatarr(comparr(0),3,*),sdatarr(comparr(0),1,*),/nodata, $ position=[0.1,0.6,0.3,0.9],$ ;ytitle='Height above sea level (km)', $ ytitle='Height (km a.s.l.)', $ yrange=[altitude, max_alt], $ ystyle=1+8-8, $ ;xtitle='T and T!Dd!N ('+STRING("260B)+'C)' , $ xtitle='T ('+STRING("260B)+'C)' , $ xmargin=[xm1,xm2], $ ymargin=[ym1,ym2], $ title= exp_name,charsize=1.3,$ ;xrange=[-80.,max(sdatarr(comparr(0),3,*))] xrange=[5,19],xstyle=1,xtickinterval=5,xminor=5;,$ ;xrange=[min(sdatarr(comparr(*),3,0:index(0))), $ ; max(sdatarr(comparr(*),3,0:index(0)))-0.1] ;axis,yaxis=1,yrange=(!Y.CRANGE-altitude), $ ; ystyle=1, ycharsize=1.3, $ ; ytitle='Height above ground (km)' for i=0,n_elements(comparr)-1 do begin oplot, sdatarr(comparr(i),3,*),sdatarr(comparr(i),1,*),color=col(i),psym=-3,line=lin(i),NSUM=nnsum(i) ; oplot, sdatarr(comparr(i),5,*),sdatarr(comparr(i),1,*),color=col(i),psym=-3,line=1,NSUM=10 endfor xshift= 1 & yshift= -1.2 for i=0,n_elements(comparr)-1 do begin ;xyouts,!x.crange(0) + (xshift ) *(!x.crange(1)-!x.crange(0))/15.,$ ; !y.crange(1) - ((i+yshift)+1./4) *(!y.crange(1)-!y.crange(0))/20.,$ ; ;strcompress(string(comparr(i),format='(i2)')) +' '+rs_id(comparr(i)), color= col(i),charsize=0.65 ; ;rs_id(comparr(i)), color= col(i),charsize=0.65 ; rs_id_UTC(comparr(i)), color= col(i),charsize=0.65 endfor ; Plotting of relative humidity ;-------------------------------- plot, sdatarr(comparr(0),4,*),sdatarr(comparr(0),1,*),/nodata, $ position=[0.3,0.6,0.5,0.9],ytickname=nullstrarr,$ ;ytitle='Height above sea level (km)', $ ;ytitle='Height (km ASL)', $ yrange=[altitude, max_alt], $ ystyle=1+8-8, $ xmargin=[xm1,xm2], $ ymargin=[ym1,ym2], $ xtitle='RH (%)', $ title= sitestring,charsize=1.3,$ xrange=[20,100] ;axis,yaxis=1,yrange=(!Y.CRANGE-altitude), $ ; ystyle=1, ycharsize=1.3, $ ; ytitle='Height above ground (km)' for i=0,n_elements(comparr)-1 do begin oplot, sdatarr(comparr(i),4,*),sdatarr(comparr(i),1,*),color=col(i),psym=-3,line=lin(i),NSUM=nnsum(i) endfor ; Plotting of windspeed ;-------------------------------- plot, sdatarr(comparr(0),7,*),sdatarr(comparr(0),1,*),/nodata, $ position=[0.5,0.6,0.7,0.9],ytickname=nullstrarr,$ xrange=[0,14.99+5],xstyle=1,xtickinterval=5,xminor=5,$ ;ytitle='Height above sea level (km)', $ ;ytitle='Height (km ASL)', $ yrange=[altitude, max_alt], $ ystyle=1+8-8, $ xmargin=[xm1,xm2], $ ymargin=[ym1,ym2], $ xtitle='WS (m s!U-1!N)', $ title= location,charsize=1.3 ;axis,yaxis=1,yrange=(!Y.CRANGE-altitude), $ ; ystyle=1, ycharsize=1.3, $ ; ytitle='Height above ground (km)' for i=0,n_elements(comparr)-1 do begin oplot, sdatarr(comparr(i),7,*),sdatarr(comparr(i),1,*),color=col(i),psym=-3,line=lin(i),NSUM=nnsum(i) endfor ; Plotting of wind direction ;-------------------------------- plot, sdatarr(comparr(0),6,*),sdatarr(comparr(0),1,*),psym=3,/nodata, $ position=[0.7,0.6,0.9,0.9],ytickname=nullstrarr,$ ;ytitle='Height above sea level (km)', $ ;ytitle='Height (km ASL)', $ yrange=[altitude, max_alt], $ ystyle=1+8-8, $ xmargin=[xm1,xm2], $ ymargin=[ym1,ym2], $ xtitle='WD ('+STRING("260B)+')', $ xrange=[0,360],xstyle = 1, xticks=4,xminor=3, $ title= filename,charsize=1.3 ;axis,yaxis=1,yrange=(!Y.CRANGE-altitude), $ ; ystyle=1, ycharsize=1.3, $ ; ytitle='Height above ground (km)' for i=0,n_elements(comparr)-1 do begin oplot, sdatarr(comparr(i),6,*),sdatarr(comparr(i),1,*),color=col(i),psym=-3,line=lin(i),NSUM=nnsum(i) endfor xshift= 15.25 & yshift= 1 for i=0,n_elements(comparr)-1 do begin xyouts,!x.crange(0) + (xshift ) *(!x.crange(1)-!x.crange(0))/15.,$ !y.crange(1) - ((i+yshift)+1./4) *(!y.crange(1)-!y.crange(0))/15.,$ ;strcompress(string(comparr(i),format='(i2)')) +' '+rs_id(comparr(i)), color= col(i),charsize=0.65 ;rs_id(comparr(i)), color= col(i),charsize=0.65 rs_id_DAY_UTC(comparr(i)), color= col(i),charsize=0.65 endfor DEVICE,/CLOSE END