; rdwdjro.pro 1029-?? min to run now in tiegcm_crit!! ; JPG, JPEG, PNG all fail set_plot, 'PS' if15sep=1 if (if15sep eq 1) then begin device, filename="jro_viz_15sep05.ps",yoff=1.0, ysize=5.0, /inches, /color endif else begin device, filename="jro_viz_13-22sep05.ps",yoff=1.0, ysize=5.0, /inches, /color endelse ifhd = 0 ifh2 = 0 year1 = 2005 day1 = 256 day2 = 266 ntmpda = 24*3 y1=-75 y2=85 lenda = day2 - day1 y90= y1 + 0.90*(y2-y1) y110 = y1 + 1.05*(y2-y1) ; Read IMF file used for TIEGCM run of 15-min aves every 5 min (t=0 from t=-5 to -20m) fileimf = '/hao/aim1/tgcm/data/imf_OMNI_2005001-2005365.nc' nc = ncdf_open(fileimf) ncdf_varget, nc, 'date', yrdayfracimf ncdf_varget, nc, 'bx', yrbximf ncdf_varget, nc, 'by', yrbyimf ncdf_varget, nc, 'bz', yrbzimf ncdf_varget, nc, 'swvel', yrvswimf ntimf=size(yrdayfracimf,/dimensions) ncdf_close, nc ; ntimf=105122 but expected 105121 since 365*24*12=105120 times every 5 min nt1 = (day1-1)*288 nt2 = day2*288 + 1 print,ntimf,day1,nt1,day2,nt2,yrdayfracimf(nt1),yrdayfracimf(nt2) ; Why is nt1=7904 and nt2=11073 when they should be 73440 and 76609?? nt1=73440 nt2=76609 ntot = nt2-nt1 print,ntot,yrdayfracimf(nt1),yrdayfracimf(nt2) dafracimf = fltarr(ntot) bzimf=fltarr(ntot) byz=fltarr(ntot) byimf=fltarr(ntot) bximf=fltarr(ntot) btmag=fltarr(ntot) vsw=fltarr(ntot) for n=nt1,nt2-1 do begin ; Fails because is a long?? ; dafracimf(n-nt1) = yrdayfracimf(n)-year1*1000 dafracimf(n-nt1) = day1 + (n-nt1)*5./288. bzimf(n-nt1) = yrbzimf(n) byimf(n-nt1) = yrbyimf(n) bximf(n-nt1) = yrbximf(n) byz(n-nt1) = abs(yrbyimf(n)/yrbzimf(n)) btmag(n-nt1) = sqrt(yrbximf(n)^2+yrbyimf(n)^2+yrbzimf(n)^2) vsw(n-nt1) = yrvswimf(n) endfor n=0 print,dafracimf(n),bzimf(n),byimf(n),byz(n),btmag(n),vsw(n) close,9 char8=string(8) jro_viz=fltarr(30000) yday=fltarr(30000) jro_viz_20min=fltarr(800) for n=0,799 do begin jro_viz_20min(n) = -999.99 endfor i=0 ;openr,9,'/hao/polar1/emery/tarkin_d/hermes_i_bozo/jro/2005drift/cedar_avg_drifts_2005/jro_av_viz_2005' ;nskip = 2250 openr,9,'jro_av_viz_2005' nskip=2271 ; Read over nskip lines for n=1,nskip do begin readf,9,'$(a8)',char8 ;if (n gt 2250) then print,char8 endfor iday0 = 256 minendp = 55 while ~ eof(9) do begin readf,9,'$(19x,i2,8x,i2,1x,2f7.2,f10.4)',minbeg,minend,jro_Fav_viz,err_viz,dayfrac yday(i)= dayfrac jro_viz(i) = jro_Fav_viz ; JRO Obs are mostly every 5 min, but have gaps nummin = fix( (dayfrac - iday0)*24.*60. - 2.49) ; This works if minbeg=minendp or if minbeg function of 20 min ; - fails for minendp function of 20 min followed by a gap m20 = fix(minbeg/20) if (m20*20 eq minbeg) then begin im1 = i m20e = fix(minendp/20) if (m20e*20 eq minendp) then im1 = i-1 i20 = fix((nummin+1)/20) jro_viz_20min(i20) = 0.5*(jro_viz(i)+jro_viz(im1)) ;print,'dayfrac i im1 i20 viz20m ',dayfrac,i,im1,i20,jro_viz_20min(i20) endif minendp = minend i=i+1 endwhile i=i-1 print,'i yday, viz',i,yday(i),jro_viz(i) for i2=i+1,30000-1 do begin yday(i2) = yday(i) jro_viz(i2) = jro_viz(i) endfor close,9 ; Jicamarca coordinates (files from Anderson and Chau) jj=15 ii=21 ; F reg aves are for 24 ranges from 218-577 km in the JRO Viz data (02Apr ilev=13-24) ; JULIA 150-km drifts are ~140-170 km (ilev=11-13? 200~15?) with kindat=1950 (<2006 and avg >=2006) ht1 = 218.e+5 ht2 = 577.e+5 ntim = lenda*ntmpda ; no values on first time so remove ;ntim = ntim-1 jroav_wconcrit=fltarr(ntim) jroav_wdcrit15=fltarr(ntim) timewcon=fltarr(ntim) timewd15=fltarr(ntim) itwcon = 0 itwd15 = 0 ;for num=1,37 do begin for num=1,36 do begin filewcon=string(num,format='("/hao/aim2/emery/tiegcm_crit/EMERY.tiegcm_concrit.s_sep05_",i3.3,".nc")') nc = ncdf_open(filewcon) ncdf_varget, nc, 'mtime', timed timew=(timed(0,*)-day1)*24+timed(1,*)+timed(2,*)/60. ncdf_varget, nc, 'ZG', zwd ncdf_varget, nc, 'WI_VEL', wi_wcon n4dd=size(wi_wcon,/dimensions) ngd=n4dd(0) & nLd=n4dd(1) & nzd=n4dd(2) & ntd= n4dd(3) ncdf_close, nc ; print,'first,last time =',num,itwcon,itwd15,year1,day1,timed(0,0),timewcon(0),timewcon(ntd-1)/24. ; Skip first model time if not from a spin-up (it1=1 instead of it1=0) it1=1 it1=0 for it=it1,ntd-1 do begin tem=0. nk = 0 for kk=1,nzd-1 do begin if (zwd(ii,jj,kk,it) gt ht1 and zwd(ii,jj,kk,it) lt ht2) then begin tem = tem + wi_wcon(ii,jj,kk,it) nk = nk + 1 endif endfor ; Convert from cm/s to m/s jroav_wconcrit(itwcon+it-it1)=tem*0.01/nk timewcon(itwcon+it-it1) = timew(it) endfor ; print,'wcon nht last ave =',nk,timewcon(itwcon+ntd-1-it1),jroav_wconcrit(itwcon+ntd-1-it1) itwcon = itwcon + ntd-it1 filewd15=string(num,format='("/hao/aim2/emery/tiegcm_crit/EMERY.tiegcm_dcrit15.s_sep05_",i3.3,".nc")') nc = ncdf_open(filewd15) ncdf_varget, nc, 'mtime', timeded timewe=(timeded(0,*)-day1)*24+timeded(1,*)+timeded(2,*)/60. ncdf_varget, nc, 'ZG', zwded ncdf_varget, nc, 'WI_VEL', wi_wd15 n4dd=size(wi_wd15,/dimensions) ngd=n4dd(0) & nLd=n4dd(1) & nzd=n4dd(2) & ntded= n4dd(3) ncdf_close, nc ; print,'wd15 first day =',timeded(0,0) ; Skip first model time if not from a spin-up (it1=1 instead of it1=0) it1=1 it1=0 for it=it1,ntded-1 do begin tem=0. nk = 0 for kk=1,nzd-1 do begin if (zwded(ii,jj,kk,it) gt ht1 and zwded(ii,jj,kk,it) lt ht2) then begin tem = tem + wi_wd15(ii,jj,kk,it) nk = nk + 1 endif endfor ; Convert from cm/s to m/s jroav_wdcrit15(itwd15+it-it1)=tem*0.01/nk timewd15(itwd15+it-it1) = timewe(it) endfor ; print,'wd15 nht last ave =',nk,timewd15(itwd15+ntded-1-it1),jroav_wdcrit15(itwd15+ntded-1-it1) itwd15 = itwd15 + ntded-it1 endfor ; end of for num=1,37 jroav_w2d10=fltarr(ntim*10) timew2d10=fltarr(ntim*10) itw2d10 = 0 jroav_wdcrit10=fltarr(ntim) timewd10=fltarr(ntim) itwd10 = 0 ;for num=1,13 do begin for num=1,12 do begin filewd10=string(num,format='("/hao/aim2/emery/tiegcm_crit/EMERY.tiegcm_dcrit10.swi2_sep05_",i3.3,".nc")') nc = ncdf_open(filewd10) ncdf_varget, nc, 'mtime', timed10 timew10=(timed10(0,*)-day1)*24+timed10(1,*)+timed10(2,*)/60. ncdf_varget, nc, 'ZG', zwd10 ncdf_varget, nc, 'WI_VEL', wi_wd10 n4dd=size(wi_wd10,/dimensions) ngd=n4dd(0) & nLd=n4dd(1) & nzd=n4dd(2) & ntd= n4dd(3) ncdf_close, nc print,'wd10 last time =',num,ntd,timew10(ntd-1)/24. ; Skip first model time if not from a spin-up (it1=1 instead of it1=0) it1=1 it1=0 for it=it1,ntd-1 do begin tem=0. nk = 0 for kk=1,nzd-1 do begin if (zwd10(ii,jj,kk,it) gt ht1 and zwd10(ii,jj,kk,it) lt ht2) then begin tem = tem + wi_wd10(ii,jj,kk,it) nk = nk + 1 endif endfor ; Convert from cm/s to m/s jroav_w2d10(itw2d10+it-it1)=tem*0.01/nk timew2d10(itw2d10+it-it1) = timew10(it) endfor ; save now every 20 min instead of every 2 min for it=it1,ntd-1,10 do begin jroav_wdcrit10(itwd10)=jroav_w2d10(itw2d10+it-it1) timewd10(itwd10) = timew2d10(itw2d10+it-it1) itwd10 = itwd10 + 1 endfor itw2d10 = itw2d10 + ntd-it1 print,itwd10,itw2d10,timewd10(0),timewd10(itwd10-1) endfor ; end of for num=1,12 for 2min saves of crit2-1=10deg jroav_wdcrit05=fltarr(ntim) timewd05=fltarr(ntim) itwd05 = 0 jroav_w10d05=fltarr(ntim*2) timew10d05=fltarr(ntim*2) itw10d05 = 0 for num=1,14 do begin filewd05=string(num,format='("/hao/aim2/emery/tiegcm_crit/EMERY.tiegcm_dcrit05.swi10_sep05_",i3.3,".nc")') nc = ncdf_open(filewd05) ncdf_varget, nc, 'mtime', timed05 timew05=(timed05(0,*)-day1)*24+timed05(1,*)+timed05(2,*)/60. ncdf_varget, nc, 'ZG', zwd05 ncdf_varget, nc, 'WI_VEL', wi_wd05 n4dd=size(wi_wd05,/dimensions) ngd=n4dd(0) & nLd=n4dd(1) & nzd=n4dd(2) & ntd= n4dd(3) ncdf_close, nc print,'wd05 last time =',num,ntd,timew05(ntd-1)/24. ; Skip first model time if not from a spin-up (it1=1 instead of it1=0) it1=1 it1=0 for it=it1,ntd-1 do begin tem=0. nk = 0 for kk=1,nzd-1 do begin if (zwd05(ii,jj,kk,it) gt ht1 and zwd05(ii,jj,kk,it) lt ht2) then begin tem = tem + wi_wd05(ii,jj,kk,it) nk = nk + 1 endif endfor ; Convert from cm/s to m/s jroav_w10d05(itw10d05+it-it1)=tem*0.01/nk timew10d05(itw10d05+it-it1) = timew05(it) endfor ; save now every 20 min instead of every 10 min for it=it1,ntd-1,2 do begin jroav_wdcrit05(itwd05)=jroav_w10d05(itw10d05+it-it1) timewd05(itwd05) = timew10d05(itw10d05+it-it1) itwd05 = itwd05 + 1 endfor itw10d05 = itw10d05 + ntd-it1 print,itwd05,itw10d05,timewd05(0),timewd05(itwd05-1) endfor ; end of for num=1,15 for 10 min saves of crit2-1=05deg jroav_wd2p5=fltarr(ntim) timewd2p5=fltarr(ntim) itwd2p5 = 0 for num=1,6 do begin filewd2p5=string(num,format='("/hao/aim2/emery/tiegcm_crit/EMERY.tiegcm_dc1_2p5_mn10_diff15.s_sep05_",i3.3,".nc")') nc = ncdf_open(filewd2p5) ncdf_varget, nc, 'mtime', timed2p5 timew2p5=(timed2p5(0,*)-day1)*24+timed2p5(1,*)+timed2p5(2,*)/60. ncdf_varget, nc, 'ZG', zwd2p5 ncdf_varget, nc, 'WI_VEL', wi_wd2p5 n4dd=size(wi_wd2p5,/dimensions) ngd=n4dd(0) & nLd=n4dd(1) & nzd=n4dd(2) & ntd= n4dd(3) ncdf_close, nc print,'wd2p5 last time =',num,ntd,timew2p5(ntd-1)/24. ; Skip first model time if not from a spin-up (it1=1 instead of it1=0) it1=1 it1=0 for it=it1,ntd-1 do begin tem=0. nk = 0 for kk=1,nzd-1 do begin if (zwd2p5(ii,jj,kk,it) gt ht1 and zwd2p5(ii,jj,kk,it) lt ht2) then begin tem = tem + wi_wd2p5(ii,jj,kk,it) nk = nk + 1 endif endfor ; Convert from cm/s to m/s jroav_wd2p5(itwd2p5+it-it1)=tem*0.01/nk timewd2p5(itwd2p5+it-it1) = timew2p5(it) endfor itwd2p5 = itwd2p5 + ntd-it1 endfor ; end of for num=1,7 for 20 min saves of crit2-1=15deg, c1=crb+2.5 jroav_wd7p5=fltarr(ntim) timewd7p5=fltarr(ntim) itwd7p5 = 0 for num=1,6 do begin filewd7p5=string(num,format='("/hao/aim2/emery/tiegcm_crit/EMERY.tiegcm_dc1_7p5_mn10_diff15.s_sep05_",i3.3,".nc")') nc = ncdf_open(filewd7p5) ncdf_varget, nc, 'mtime', timed7p5 timew7p5=(timed7p5(0,*)-day1)*24+timed7p5(1,*)+timed7p5(2,*)/60. ncdf_varget, nc, 'ZG', zwd7p5 ncdf_varget, nc, 'WI_VEL', wi_wd7p5 n4dd=size(wi_wd7p5,/dimensions) ngd=n4dd(0) & nLd=n4dd(1) & nzd=n4dd(2) & ntd= n4dd(3) ncdf_close, nc print,'wd7p5 last time =',num,ntd,timew7p5(ntd-1)/24. ; Skip first model time if not from a spin-up (it1=1 instead of it1=0) it1=1 it1=0 for it=it1,ntd-1 do begin tem=0. nk = 0 for kk=1,nzd-1 do begin if (zwd7p5(ii,jj,kk,it) gt ht1 and zwd7p5(ii,jj,kk,it) lt ht2) then begin tem = tem + wi_wd7p5(ii,jj,kk,it) nk = nk + 1 endif endfor ; Convert from cm/s to m/s jroav_wd7p5(itwd7p5+it-it1)=tem*0.01/nk timewd7p5(itwd7p5+it-it1) = timew7p5(it) endfor itwd7p5 = itwd7p5 + ntd-it1 endfor ; end of for num=1,7 for 20 min saves of crit2-1=15deg, c1=crb+7.5 jroav_wd17p5=fltarr(ntim) timewd17p5=fltarr(ntim) itwd17p5 = 0 for num=1,6 do begin filewd17p5=string(num,format='("/hao/aim2/emery/tiegcm_crit/EMERY.tiegcm_dyncrit.ssep05_crb_17p5_05_",i3.3,".nc")') nc = ncdf_open(filewd17p5) ncdf_varget, nc, 'mtime', timed17p5 timew17p5=(timed17p5(0,*)-day1)*24+timed17p5(1,*)+timed17p5(2,*)/60. ncdf_varget, nc, 'ZG', zwd17p5 ncdf_varget, nc, 'WI_VEL', wi_wd17p5 n4dd=size(wi_wd17p5,/dimensions) ngd=n4dd(0) & nLd=n4dd(1) & nzd=n4dd(2) & ntd= n4dd(3) ncdf_close, nc print,'wd17p5 last time =',num,ntd,timew17p5(ntd-1)/24. ; Skip first model time if not from a spin-up (it1=1 instead of it1=0) it1=1 it1=0 for it=it1,ntd-1 do begin tem=0. nk = 0 for kk=1,nzd-1 do begin if (zwd17p5(ii,jj,kk,it) gt ht1 and zwd17p5(ii,jj,kk,it) lt ht2) then begin tem = tem + wi_wd17p5(ii,jj,kk,it) nk = nk + 1 endif endfor ; Convert from cm/s to m/s jroav_wd17p5(itwd17p5+it-it1)=tem*0.01/nk timewd17p5(itwd17p5+it-it1) = timew17p5(it) endfor itwd17p5 = itwd17p5 + ntd-it1 endfor ; end of for num=1,7 for 20 min saves of crit2-1=5deg, c1=crb+17.5 ifwr20=1 if (ifwr20 eq 1) then begin ; Print out values of Wi from TIEGCM to a file ncwrfil = 'wi_7b_tie_imf_sep05' openw,wrunit,ncwrfil,/get_lun printf,wrunit,format='("JRO Viz (VperpN) 20 min aves and Weimer05/TIE-GCM Wi values in m/s every 20 min for different crit1,2")' printf,wrunit,format='("CON=15,+15; CRB17P5_5=crb+17.5,+5; CRB7P5_15=crb+7.5,+15 etc listing from most equatorward dynamic crit2 to least equatorward")' printf,wrunit,format='("year daynum uthr uthincr Bx By Bz Bt Vsw VizJRO CON CRB17P5_5 7P5_15 CRB5_15 2P5_15 CRB5_10 CRB5_5")' ida0 = fix(timewd15(0)/24.) ; i=0 at 0 UT 256 has -NaN WI_VEL at start, so start with i=1 at 0020 UT ;for i=1,itwd15-1 do begin for i=0,itwd15-1 do begin iday = fix(timewd15(i)/24.+day1) utincr = timewd15(i) - ida0*24. ndays = fix(utincr/24.) utda = utincr - ndays*24 i4=i*4 ; JRO is i+1 because i=0 is 0020 for models instead of 0000 with it1=1 ; printf,wrunit,iday,utincr,utda,jroav_wconcrit(i),jroav_wdcrit15(i),jro_viz_20min(i+1),format='("2005 ",i03," ",5f8.2)' printf,wrunit,iday,utincr,utda,bximf(i4),byimf(i4),bzimf(i4),btmag(i4),vsw(i4),jro_viz_20min(i),jroav_wconcrit(i),jroav_wd17p5(i),jroav_wd7p5(i),jroav_wdcrit15(i),jroav_wd2p5(i),jroav_wdcrit10(i),jroav_wdcrit05(i),format='("2005 ",i03," ",15f8.2)' endfor endif ; Read values at every 20 min and calculate root mean square openr,9,'tie_jro_wi_sep05' timejro20=fltarr(720) jro20=fltarr(720) sqwcon=fltarr(720) sqwd15=fltarr(720) sqwd10=fltarr(720) sqwd05=fltarr(720) sqwd2p5=fltarr(720) sqwd7p5=fltarr(720) sqwd17p5=fltarr(720) iend = 720 - 1 j20=intarr(10) dmswcon=fltarr(10) dmswd15=fltarr(10) dmswd10=fltarr(10) dmswd05=fltarr(10) dmswd2p5=fltarr(10) dmswd7p5=fltarr(10) dmswd17p5=fltarr(10) i20 = 0 readf,9,'$(a8)',char8 rmswcon = 0. rmswd15 = 0. rmswd10 = 0. rmswd05 = 0. rmswd2p5 = 0. rmswd7p5 = 0. rmswd17p5 = 0. for j=0,9 do begin j20(j) = 0 dmswcon(j) = 0. dmswd15(j) = 0. dmswd10(j) = 0. dmswd05(j) = 0. dmswd2p5(j) = 0. dmswd7p5(j) = 0. dmswd17p5(j) = 0. endfor for i=0,iend-1 do begin readf,9,'$(8x,f9.2,24x,f8.2)',timej,av20 ;print,i,timej,av20 if (av20 gt -990.) then begin timejro20(i20)=timej jro20(i20)=av20 ; calculate rms sqwcon(i20) = (av20-jroav_wconcrit(i))^2 sqwd15(i20) = (av20-jroav_wdcrit15(i))^2 sqwd10(i20) = (av20-jroav_wdcrit10(i))^2 sqwd05(i20) = (av20-jroav_wdcrit05(i))^2 sqwd2p5(i20) = (av20-jroav_wd2p5(i))^2 sqwd7p5(i20) = (av20-jroav_wd7p5(i))^2 sqwd17p5(i20) = (av20-jroav_wd17p5(i))^2 rmswcon = rmswcon + sqwcon(i20) rmswd15 = rmswd15 + sqwd15(i20) rmswd10 = rmswd10 + sqwd10(i20) rmswd05 = rmswd05 + sqwd05(i20) rmswd2p5 = rmswd2p5 + sqwd2p5(i20) rmswd7p5 = rmswd7p5 + sqwd7p5(i20) rmswd17p5 = rmswd17p5 + sqwd17p5(i20) ; calculate for each day ida = fix(timej/24.) dmswcon(ida) = dmswcon(ida) + sqwcon(i20) dmswd15(ida) = dmswd15(ida) + sqwd15(i20) dmswd10(ida) = dmswd10(ida) + sqwd10(i20) dmswd05(ida) = dmswd05(ida) + sqwd05(i20) dmswd2p5(ida) = dmswd2p5(ida) + sqwd2p5(i20) dmswd7p5(ida) = dmswd7p5(ida) + sqwd7p5(i20) dmswd17p5(ida) = dmswd17p5(ida) + sqwd17p5(i20) j20(ida) = j20(ida) + 1 i20 = i20 + 1 endif endfor ; Convert to m/s rms: rmswcon = sqrt(rmswcon/i20) rmswd15 = sqrt(rmswd15/i20) rmswd10 = sqrt(rmswd10/i20) rmswd05 = sqrt(rmswd05/i20) rmswd2p5 = sqrt(rmswd2p5/i20) rmswd7p5 = sqrt(rmswd7p5/i20) rmswd17p5 = sqrt(rmswd17p5/i20) print,rmswcon,rmswd15,rmswd10,rmswd05,rmswd2p5,rmswd7p5,rmswd17p5 for j=0,9 do begin dmswcon(j) = sqrt(dmswcon(j)/j20(j)) dmswd15(j) = sqrt(dmswd15(j)/j20(j)) dmswd10(j) = sqrt(dmswd10(j)/j20(j)) dmswd05(j) = sqrt(dmswd05(j)/j20(j)) dmswd2p5(j) = sqrt(dmswd2p5(j)/j20(j)) dmswd7p5(j) = sqrt(dmswd7p5(j)/j20(j)) dmswd17p5(j) = sqrt(dmswd17p5(j)/j20(j)) print,j,dmswcon(j),dmswd15(j),dmswd10(j),dmswd05(j),dmswd2p5(j),dmswd7p5(j),dmswd17p5(j) endfor for i2=i20-1,iend-1 do begin timejro20(i2) = timejro20(i20-2) jro_viz(i2) = jro20(i20-2) endfor close,9 loadct,39 ; 40=dark blue, 75=med blue, 100=light blue, 125=blue-green, 150=green, 175=yellow-green, 200=gold, 225=orange-red, 250=red, 275=blank (>255!) imfn=-10 imfx=10 yi110 = 1.1*(imfx-imfn) + imfn yi10 = 0.05*(imfx-imfn) + imfn if (if15sep eq 1) then begin sday1 = 258 sday2 = 259 plot,yday,jro_viz,xrange=[sday1,sday2],yrange=[y1,y2],xstyle=1,ystyle=1,position=[0.05,0.05,0.99,0.85],ytitle='ISR F-region Viz (m/s)',xtitle='daynumber' oplot,timewcon/24.+day1,jroav_wconcrit,color=150 oplot,timewd15/24.+day1,jroav_wdcrit15,color=250 oplot,timewd10/24.+day1,jroav_wdcrit10,color=100 oplot,timewd05/24.+day1,jroav_wdcrit05,color=190 oplot,timewd2p5/24.+day1,jroav_wd2p5,color=40 oplot,timewd7p5/24.+day1,jroav_wd7p5,color=175 oplot,timewd17p5/24.+day1,jroav_wd17p5,color=220 oplot,timew2d10/24.+day1,jroav_w2d10,color=100,linestyle=1 oplot,timew10d05/24.+day1,jroav_w10d05,color=190,linestyle=1 oplot,[sday1,sday2],[0,0] xyouts,sday1,y110,'Jicamarca (12S,76.8W)',charsize=1.5 xyouts,0.01*(sday2-sday1)+sday1,y90,'WCONCRIT',color=150 xyouts,0.15*(sday2-sday1)+sday1,y90,'WDCRIT15',color=250 xyouts,0.29*(sday2-sday1)+sday1,y90,'WDCRIT10',color=100 xyouts,0.43*(sday2-sday1)+sday1,y90,'WDCRIT05',color=190 xyouts,0.57*(sday2-sday1)+sday1,y90,'CRB_2.5_15',color=40 xyouts,0.70*(sday2-sday1)+sday1,y90,'CRB_7.5_15',color=175 xyouts,0.83*(sday2-sday1)+sday1,y90,'CRB_17.5_5',color=220 xyouts,0.50*(sday2-sday1)+sday1,y110,'TIEGCM W05 15 Sep 2005',size=1.5 plot,dafracimf,bzimf,xrange=[sday1,sday2],yrange=[imfn,imfx],xstyle=1,ystyle=1,position=[0.05,0.05,0.99,0.85],ytitle='IMF (nT)',xtitle='daynumber' oplot,dafracimf,-vsw/100.,color=100 oplot,dafracimf,byimf,color=150 oplot,dafracimf,btmag oplot,dafracimf,byz,color=250 oplot,[sday1,sday2],[0,0] xyouts,0.04*(sday2-sday1)+sday1,yi10,'-Vsw/100',color=100 xyouts,0.23*(sday2-sday1)+sday1,yi10,'Bz,Bt' xyouts,0.42*(sday2-sday1)+sday1,yi10,'By',color=150 xyouts,0.61*(sday2-sday1)+sday1,yi10,'|By/Bz|',color=250 xyouts,0.61*(sday2-sday1)+sday1,yi110,'15 Sep 2005',size=1.5 endif else begin plot,yday,jro_viz,xrange=[day1,day2],yrange=[y1,y2],xstyle=1,ystyle=1,position=[0.05,0.05,0.99,0.85],ytitle='ISR F-region Viz (m/s)',xtitle='daynumber' ;plot,yday,jro_viz,xrange=[day1,day2],yrange=[y1,y2],xstyle=1,ystyle=1,position=[0.01,0.01,0.99,0.47],ytitle='Viz (m/s)',xtitle='daynumber' oplot,timewcon/24.+day1,jroav_wconcrit,color=150 ;oplot,timejro20/24.+day1,jro20 oplot,timewd15/24.+day1,jroav_wdcrit15,color=250 ;oplot,timewd10/24.+day1,jroav_wdcrit10,color=125 ;oplot,timewd05/24.+day1,jroav_wdcrit05,color=225 oplot,timewd2p5/24.+day1,jroav_wd2p5,color=40 oplot,timewd7p5/24.+day1,jroav_wd7p5,color=175 oplot,timewd17p5/24.+day1,jroav_wd17p5,color=220 oplot,timew2d10/24.+day1,jroav_w2d10,color=100,linestyle=1 oplot,timew10d05/24.+day1,jroav_w10d05,color=190,linestyle=1 oplot,[day1,day2],[0,0] xyouts,day1,y110,'Jicamarca (12S,76.8W)',charsize=1.5 xyouts,0.01*(day2-day1)+day1,y90,'WCONCRIT',color=150 xyouts,0.15*(day2-day1)+day1,y90,'WDCRIT15',color=250 xyouts,0.29*(day2-day1)+day1,y90,'WDCRIT10',color=100 xyouts,0.43*(day2-day1)+day1,y90,'WDCRIT05',color=190 xyouts,0.57*(day2-day1)+day1,y90,'CRB_2.5_15',color=40 xyouts,0.70*(day2-day1)+day1,y90,'CRB_7.5_15',color=175 xyouts,0.83*(day2-day1)+day1,y90,'CRB_17.5_5',color=220 xyouts,0.50*(day2-day1)+day1,y110,'TIEGCM W05 13-22 Sep 2005',size=1.5 plot,dafracimf,bzimf,xrange=[day1,day2],yrange=[imfn,imfx],xstyle=1,ystyle=1,position=[0.05,0.05,0.99,0.85],ytitle='IMF (nT)',xtitle='daynumber' oplot,dafracimf,-vsw/100.,color=100 oplot,dafracimf,byimf,color=150 oplot,dafracimf,btmag oplot,dafracimf,byz,color=250 oplot,[day1,day2],[0,0] xyouts,0.04*(day2-day1)+day1,yi10,'-Vsw/100',color=100 xyouts,0.23*(day2-day1)+day1,yi10,'Bz,Bt' xyouts,0.42*(day2-day1)+day1,yi10,'By',color=150 xyouts,0.61*(day2-day1)+day1,yi10,'|By/Bz|',color=250 xyouts,0.61*(day2-day1)+day1,yi110,'13-22 Sep 2005',size=1.5 endelse device, /close end