function plot(args) time1=subwrd(args, 1) *time1=12Z9dec2000 say time1 'reinit' 'open /export-4/sgi109/pacs/mcidas/data/global_ir/global_ir_geos_merged_half-deg/global_0.5deg_0.5hourly_tmp-lnx.ctl' 'set gxout shaded' 'set grid off' *'set mpdset mres' *'run bw.gs' *'run /export/lnx196/yzhou/CWlink/IR/rgbset.gs' 'set grads off' 'set map 0 1 1' 'set lat -50 50' 'set lon -180 180' 'set time 'time1 'q dims' line5=sublin(result,5) tmp2=subwrd(line5,9) say 'tmp2='tmp2 time1=tmp2-720+1 time2=tmp2 timem=(time1+time2)/2.0 'set t 'timem 'enable print /disk2/irstat/web_plots/meta/am_ir_seasonal.meta' 'draw string 5.15 1.6 Degrees K' 'run /disk2/irstat/web_plots/web.specs.seasonal' say 'd ave(tmpm,t='time1',t='time2',1,-b)/10' 'd ave(tmpm,t='time1',t='time2',1,-b)/10' *'d ave(tmpm,t-1,t+0)/10' 'run /disk2/irstat/web_plots/cbar1.gs' *'run /export/lnx196/yzhou/CWlink/IR/date.gs' 'draw string 5.5 6.3 'mydate() say mydate() 'print' 'disable print' 'quit' function mydate 'query time' sres = subwrd(result,3) i = 1 while (substr(sres,i,1)!='Z') i = i + 1 endwhile hour = substr(sres,1,i) isav = i i = i + 1 while (substr(sres,i,1)>='0' & substr(sres,i,1)<='9') i = i + 1 endwhile day = substr(sres,isav+1,i-isav-1) month = substr(sres,i,3) year = substr(sres,i+3,4) return (day' 'month' 'year)