'reinit'
res4=write('obs.txt', 'Stn' ' ' 'P.Clm' ' ' 'Day1' ' ' 'Day2' ' ' 'Day3' ' ' 'Day4' ' ' 'Day5' ' ' 'Day6' ' ' 'Day7')
m = 1
while(m<=24)
res= read('senegal.txt'); res1=sublin(res,2); nam=subwrd(res1,1);ltt=subwrd(res1,2);lnn=subwrd(res1,3)
'open /cpc/africawrf/ebekele/data/obs//tmax/output/obs_daily_tmax_percentile.ctl'
'set lat 12 17'
'set lon -18 -11'
'set time 26Apr'
'set lev 80'
'define clm = tmax'
'set lat 'ltt
'set lon 'lnn
'define clim=clm'
'd clim'; res2=sublin(result,1);clmm=subwrd(res2,4)
'close 1'
'open /cpc/data/cwlinks/temp/GLOBAL/hi_res/y2023/CPC_GLOBAL_T_V0.x_10min.lnx.2023.ctl'
'set lat 12 17'
'set lon -18 -11'
'set time 26Apr2023'; 'q dims'; res1=sublin(result,5); tin=subwrd(res1,9);ten=tin+6

'set lat 'ltt
'set lon 'lnn
'set t 'tin+0
'define tmm = tmax'
'd tmm'; res2=sublin(result,1);ress1=subwrd(res2,4)
'set t 'tin+1
'define tmm = tmax'
'd tmm'; res2=sublin(result,1);ress2=subwrd(res2,4)
'set t 'tin+2
'define tmm = tmax'
'd tmm'; res2=sublin(result,1);ress3=subwrd(res2,4)
'set t 'tin+3
'define tmm = tmax'
'd tmm'; res2=sublin(result,1);ress4=subwrd(res2,4)
'set t 'tin+4
'define tmm = tmax'
'd tmm'; res2=sublin(result,1);ress5=subwrd(res2,4)
'set t 'tin+5
'define tmm = tmax'
'd tmm'; res2=sublin(result,1);ress6=subwrd(res2,4)
'set t 'tin+6
'define tmm = tmax'
'd tmm'; res7=sublin(result,1);ress7=subwrd(res2,4)

tt0=math_format('%-6.2f',clmm)
tt1=math_format('%-6.2f',ress1)
tt2=math_format('%-6.2f',ress2)
tt3=math_format('%-6.2f',ress3)
tt4=math_format('%-6.2f',ress4)
tt5=math_format('%-6.2f',ress5)
tt6=math_format('%-6.2f',ress6)
tt7=math_format('%-6.2f',ress7)

res4=write('obs.txt', nam ' ' tt0 ' ' tt1 ' ' tt2 ' ' tt3 ' ' tt4 ' ' tt5 ' ' tt6 ' ' tt7)  
'close 1'
m = m + 1
endwhile

'quit'