;SET_PLOT,'PS' ;DEVICE, FILENAME='../plots/thor/tuv19z18utnew.ps',/landscape DEVICE, SET_FONT='Helvetica Bold', /TT_FONT SET_CHARACTER_SIZE=[70,90] ;device,decomposed=0 ;device,/color loadct,13 ; ****************************************************** nmf2=fltarr(73,36,216) nmf215=fltarr(73,36,216) nmf2q=fltarr(73,36) nmf2x=fltarr(73,36,432) lonx=fltarr(73,36) latx=fltarr(73,36) un=fltarr(73,36) vn=fltarr(73,36) scal=fltarr(5,36) xscal=fltarr(5,36) yscal=fltarr(5,36) labels=strarr(5) outbound=fltarr(24) inbound=fltarr(24) alata=fltarr(121) alona=fltarr(121) xlt=fltarr(121) alatn=fltarr(40,24) alonn=fltarr(40,24) time=fltarr(40,24) alatnx=fltarr(40*5,24*5) alonnx=fltarr(40*5,24*5) timex=fltarr(40*5,24*5) spv=-999999.9 dum=fltarr(73,36) ;openr,8,'../out/solmax1/nmf2.dat' ;readf,8,'$(10E11.4)',nmf2 ;close,8 ;openr,8,'../out/solmax2/nmf2.dat' ;readf,8,'$(10E11.4)',nmf215 ;close,8 ;openr,8,'fluxout.dat' openr,8,'alatn.dat' readf,8,'$(10E11.4)',alatn close,8 openr,8,'outbound.dat' readf,8,'$(10E11.4)',outbound close,8 openr,8,'inbound.dat' readf,8,'$(10E11.4)',inbound close,8 openr,8,'alonn.dat' readf,8,'$(10E11.4)',alonn close,8 openr,8,'hintm.dat' readf,8,'$(10E11.4)',time close,8 openr,8,'traj1xlt.dat' readf,8,'$(10E11.4)',xlt close,8 ; ****************************************************** ; fill lat and long arrays ; ****************************************************** for j = 0,35,1 do begin for i = 0,72,1 do begin lonx(i,j) = -180. + i*5. latx(i,j) = -87.5 + j*5. endfor endfor ; ****************************************************** ; set a zero longitude to be plotted ; ****************************************************** zlon=0. day = 14 ; ****************************************************** ; ****************************************************** ;loop top ;for nnn = 0,431,1 do begin nnn=108 ; ****************************************************** ; ****************************************************** ; set a time to be plotted (0 -> 141) from 6 0 10 to 6 23 40 ; ****************************************************** hour = nnn/3 day = hour/24 + 1 if hour gt 23 then begin hour = hour - 24 endif if hour gt 23 then begin hour = hour - 24 endif if hour gt 23 then begin hour = hour - 24 endif if hour gt 23 then begin hour = hour - 24 endif if hour gt 23 then begin hour = hour - 24 endif print,'day hour ',day,hour minute = (nnn - nnn/3*3)*20. itime = hour*6 + minute/20 jtime = hour nnnq = hour*3 + minute/20 rottime= float(hour) + float(minute)/60. ; ****************************************************** ; ****************************************************** ; set a rotation to be plotted ; ****************************************************** rot=-rottime*15. ;rot=0. print,rot,hour,minute ;stop ; if this does not work try -rot plat=1. tlon= rot + zlon if plat ge 0 then tlon = tlon+180. if tlon gt 180. then tlon = tlon-360. if tlon lt -180. then tlon = tlon+360. xlon = rot + zlon if xlon gt 180. then xlon = xlon-360. if xlon lt -180. then xlon = xlon+360. ; ; ; ****************************************************** ; find or set a maximum and minimum ; ****************************************************** ; ; max=1. ; min=1. xmax=-99999999999.9 xmin=99999999999.9 for j = 0,23,1 do begin for i = 0,39,1 do begin if time(i,j) gt xmax then xmax=time(i,j) if time(i,j) lt xmin then xmin=time(i,j) ; if alatn(i,j) eq 0. then time(i,j)=0. if time(i,j) gt 0. then begin print,time(i,j),alatn(i,j),alonn(i,j),i,j endif endfor endfor max=xmax min=xmin ; max=80. max=1.e8 min=0. print,'xmax xmin',xmax,xmin for j = 0,23,1 do begin for i = 0,39,1 do begin if time(i,j) gt max then time(i,j)=max if time(i,j) lt min then time(i,j)=min ; time(i,j)=i endfor endfor n=5 xn=5. for j = 0,23,1 do begin for i = 0,39,1 do begin alatnx(i*n,j*n) = alatn(i,j) alonnx(i*n,j*n) = alonn(i,j) timex(i*n,j*n) = time(i,j) endfor endfor for j = 0,23,1 do begin for i = 0,38,1 do begin for ii = 1,n-1,1 do begin alatnx(i*n+ii,j*n) = alatnx(i*n,j*n)*(xn-ii)/xn + alatnx((i+1)*n,j*n)*ii/xn alonnx(i*n+ii,j*n) = alonnx(i*n,j*n)*(xn-ii)/xn + alonnx((i+1)*n,j*n)*ii/xn timex(i*n+ii,j*n) = timex(i*n,j*n)*(xn-ii)/xn + timex((i+1)*n,j*n)*ii/xn endfor endfor endfor for j = 0,22,1 do begin for ii = 0,39*n,1 do begin for jj = 1,n-1,1 do begin alatnx(ii,j*n+jj) = alatnx(ii,j*n)*(xn-jj)/xn + alatnx(ii,(j+1)*n)*jj/xn alonnx(ii,j*n+jj) = alonnx(ii,j*n)*(xn-jj)/xn + alonnx(ii,(j+1)*n)*jj/xn timex(ii,j*n+jj) = timex(ii,j*n)*(xn-jj)/xn + timex(ii,(j+1)*n)*jj/xn endfor endfor endfor !p.position=[0.1,0.1,0.80,0.80] ;map_set,0.,0.,0.,/mollweide,/isotropic,/noborder,/horizon, $ ;color=255*8/10,mlinestyle=0 map_set,90.,zlon,rot,/orthographic,/isotropic, $ limit=[20,-180,90,180],/noborder,color=255.*8./10.,mlinestyle=0 contour,timex,alonnx,alatnx,/overplot,/cell_fill $ ,max_val=max,min_val=min+0.5,levels=min+findgen(30)*(max-min)/30. ;map_grid,/label,color=255 map_continents,color=255 for i=1,360,1 do begin plots,[-180.+i,-180.+i+ 1.],[20.,20.],thick=2.,color=255.*8./10. endfor for i=0,22,1 do begin plots,[-180.+i*15.,-180.+i*15.+ 15.],[outbound(i),outbound(i+1)], $ thick=2.,color=255.*8./10. endfor plots,[-180.+23*15.,-180.+0.*15.],[outbound(23),outbound(0)], $ thick=2.,color=255.*8./10. for i=0,22,1 do begin plots,[-180.+i*15.,-180.+i*15.+ 15.],[inbound(i),inbound(i+1)], $ thick=2.,color=255.*8./10. endfor plots,[-180.+23*15.,-180.+0.*15.],[inbound(23),inbound(0)], $ thick=2.,color=255.*8./10. ; ****************************************************** ; label graph ; ****************************************************** xval1=-180.-rottime*15. if xval1 lt -180. then xval1=xval1+360. xval2=0.-rottime*15. if xval2 lt -180. then xval2=xval2+360. xval3=90.-rottime*15. if xval3 lt -180. then xval3=xval3+360. xval4=-90.-rottime*15. if xval4 lt -180. then xval4=xval4+360. ;print,xval1,xval2 xyouts, xval1,20.,'12.',charsize=1.75,alignment=.5,color=255.*8./10.,font=1 xyouts, xval2,20.,'0.',charsize=1.75,alignment=.5,color=255.*8./10.,font=1 xyouts, xval4,20.,'18.',charsize=1.75,alignment=.5,color=255.*8./10.,font=1 xyouts, xval3,20.,'6.',charsize=1.75,alignment=.5,color=255.*8./10.,font=1 toplabel='Height int Fric in aur /mass UT ' + string(hour,format='(i2)') $ + ':' + string(minute,format='(i2)') + $ ' May 1997 ' + ' Day ' + string(day,format='(i2)') toplabel2=strtrim(toplabel,2) ;xyouts,xval1, 8., toplabel2,alignment=.5,charsize=1.4,color=0 xyouts,0.45,0.85, toplabel2,alignment=.5,charsize=1.9,color=255*8/10, $ /normal,font=1 ; ****************************************************** ; Draw label bar ; ****************************************************** ; bw is bar width in normalized coordinates (nc), bh is bar height (nc) ; xo, yo are the origin points in nc,labels is an array of character strings ; labcolor 0-254 ; charsize leave as default ; lableft put0 ;labtop 0 ;laboffset=0.02 ;labcolor0-254, try 0,1,254,255 bw=0.05 bh=0.8 xo=0.825 yo=0.1 botcolor=1 topcolor=254 lableft=0 labtop=0 laboffset=0.0 labcolor=255*8/10 charsize=1.25 labels(0)=string(min,format='(E10.2)') labels(1)=string(min+1./4.*(max-min),format='(E10.2)') labels(2)=string(min + 1./2.*(max-min),format='(E10.2)') labels(3)=string(min + 3./4.*(max-min),format='(E10.2)') labels(4)=string(max,format='(E10.2)') !p.color=0 clrbarg,bw,bh,xo,yo,labels,BOTCOLOR=botcolor,TOPCOLOR=topcolor,$ LABCOLOR=labcolor,CHARSIZE=charsize,LABLEFT=lableft,LABTOP=labtop,$ LABOFFSET=laboffset ; ****************************************************** filename='plots/htjoum' + $ strcompress(string(nnn,format="(i3)"),/remove_all) + '.jpg' dump=tvrd(true=3) write_jpeg,filename,dump,true=3,quality=100 ; ****************************************************** ; ****************************************************** ;end do loop ;endfor ; ****************************************************** ; ****************************************************** ;DEVICE,/CLOSE stop end ; ******************************************************