; Dec 14, 1994 Alice Lecinski with help from Joan Burkepile. ; ; This idl procedure can be used to read in and contour ; ANY MK____.[EW]ib file. ; psit='y' dirpath=' ' roti=' ' mlmb=' ' myhei=' ' height='1.74' print,'Enter in the directory path to the MK????.??ib file' print,'eg. mk3_synoptic/heightL/ ' read,dirpath print,'Enter the rotation, eg. 1845 ' read,roti print,'Enter the limb, eg. E ' read,mlmb print,'Enter the height, eg. L or M ' read,myhei & if(myhei eq 'L')then height='1.36' if (psit eq 'y') then begin ;{ !P.MULTI = [ 0, 1 , 1 ] psfilenm=roti+myhei+mlmb+'.ps' set_plot,'ps' device,/landscape,/inches,xsize=10.0,ysize=3.75,$ xoffset=0.5,yoffset=10.5,filename=psfilenm endif else begin;}{ window,0,xsize=360*3+50,ysize=181*3+50 endelse;} ahdr=fltarr(8) ckey=bytarr(30) cpad=bytarr(2) histmap=fltarr(30,37) bigmap=fltarr(390,181) ;;; /puu/ftpmk3/map/history/mk3_synoptic/heightL/MK1841.LWib filenm=dirpath+'MK'+roti+'.'+myhei+mlmb+'ib' close,1 openr,1,filenm & readu,1,ahdr & readu,1,ckey readu,1,cpad readu,1,histmap close,1 rot = fix(ahdr(0)) ahdr(1) = ahdr(1) + 1900. doy = 0 for i=0,389,13 do begin ;{ One day is ~13 degrees... lat = 0 spa=12 bigmap(i:i+spa,0:2) = histmap(doy,lat) lat = lat + 1 for j=3,181-4,5 do begin ;{ bigmap(i:i+spa,j:j+4) = histmap(doy,lat) lat = lat + 1 endfor ;} bigmap(i:i+spa,178:180) = histmap(doy,lat) if (ckey(doy) eq 48 and histmap(doy , 19) le -1000.) then ckey(doy)=0 doy = doy +1 endfor ;} doy=doy-1 bigmap(351:359,178:180) = histmap(doy,lat) if(ahdr(4) eq 2.) then bigmap=rotate(bigmap,7) ; Flip the west limb map. bigmap=smooth(bigmap,5) lvl2= [ 1000., 1200., 1400., 1600.] & lvl2=lvl2/100. clr2= [120,140,180,220] lvl= [200., 400., 600., 800., 1000., 1200., 1400., 1600.] lvl=lvl/100. mln= [1, 0, 2, 0, 2, 0, 0, 0] mln= [1, 0, 0, 0, 0, 0, 0, 0] thk= [1, 1, 2, 3, 4, 2, 2, 3] myyticks=6 & myytickv=findgen(7)*30 fedoy=ahdr(3) & edoy = fix(fedoy) fbdoy=ahdr(2) & bdoy = fix(fbdoy) rot='Rotation '+string(fix(ahdr(0))) limb=' East Limb' & if (ahdr(4) eq 2.) then limb=' West Limb' bdate='Bdoy '+string(fbdoy) & bdate=strcompress(bdate) edate='Edoy '+string(fedoy) & edate=strcompress(edate) year='Year '+string(fix(ahdr(1))) & year=strcompress(year) byr=fix(ahdr(1)) leap= long((byr mod 4 eq 0) and (byr mod 100 ne 0) or $ (byr mod 400 eq 0)) fedoy = fedoy + (365+leap)*(edoy lt bdoy)*1.0 edoy = edoy + (365+leap)*(edoy lt bdoy) numday = edoy-bdoy myxticks=numday+1 myxtickv=fltarr(numday+2) deltax = 13. offset = 0. myxname=replicate(' X ',numday+2) for ii=0,myxticks do begin ;{ if(ckey(ii+0) eq 48) then begin ;{ myxname(ii) = string(format='(i8)',edoy-ii) endif;} myxtickv(ii) = offset + ii*deltax endfor ;} myyname=replicate(' ',19) myyname(0)='-90' myyname(6)=' 90' myyname(3)=' 0' contour,bigmap,/follow,levels=lvl,$ c_linestyle=mln,$ c_thick=thk,$ xstyle=1,ystyle=1,zstyle=1,charsize=1.0,$ xcharsize=0.6,$ xtitle='DOY', ytitle='LATITUDE',$ title=rot+limb+' '+year+' '+bdate+' '+edate ,$ subtitle=height, $ xtickname=myxname,$ ytickname=myyname,$ xticks=myxticks,xtickv=myxtickv, $ yticks=myyticks,ytickv=myytickv, $ ticklen= -0.015 , xmargin=[10,10] , ymargin=[4,4] contour,bigmap,/follow,levels=lvl2,c_color=clr2,overplot=1,fill=1 espot=(edoy*1.0-fedoy)*13. +13. oplot,[espot,espot],[0,180],thick=3 for ii=0,5 do begin ;{ espot = espot-0.2 oplot,[espot,espot],[0,180],thick=3 endfor;} bspot=(edoy*1.0-fbdoy)*13. +13. oplot,[bspot,bspot],[0,180],thick=3 for ii=0,5 do begin ;{ bspot=bspot+0.2 oplot,[bspot,bspot],[0,180],thick=3 endfor;} if (psit eq 'y') then begin ;{ device,/close set_plot,'x' endif ;} end