g\set /def g\set exp .8 def char dt*12[2] logi*8 def real ut[2] def inte ndim m k o scan1 scan2 scan3 ! sic user let logi 'sys_info' ! find get cx_num[1] let dt[1] 'date_observed' let ut[1] utobs|pi*12 get cx_num[found] let dt[2] 'date_observed' let ut[2] utobs|pi*12 ! if logi.eq."OBSERVER" then sic output meteo.1 say 'dt[1]' 'dt[2]' /format a12 a12 sic output ! sic output meteo.2 say "typeset -l date1 date2" say "rm -f meteo.3 meteo.4" say "read -r date1 date2 < meteo.1" say "dat1=${date1%19*}; dat2=${date1##*19}" say "cat /users1/oper/operator/METEO/data/$dat1$dat2.met > meteo.3" say "cat /users1/oper/operator/METEO/data/$dat1$dat2.met > meteo.5" say "if [[ $date1 != $date2 ]] then" say " dat1=${date2%19*}; dat2=${date2##*19}" say " cat /users1/oper/operator/METEO/data/$dat1$dat2.met > meteo.4" say " cat /users1/oper/operator/METEO/data/$dat1$dat2.met >> meteo.5" say "fi" sic output ! sys ". meteo.2; rm -f meteo.1 meteo.2" ! col x 1 /file meteo.5 define real time[nxy] col x 1 /file meteo.3 for i 1 to nxy let ndim nxy let time[i] x[i]-2 next sic search meteo.4 if sic$exist then let ut[2] ut[2]+24 col x 1 /file meteo.4 let x x+24 for i 1 to nxy let m i+ndim let time[m] x[i]-2 next endif ! col y 2 /file meteo.5 let x time g\set box 3.5 12 13 18 lim ut[1] ut[2] * * cl pl box n conn ! col y 3 g\set box 3.5 12 7.8 12.8 lim ut[1] ut[2] * * box n conn ! col y 4 g\set box 15.5 24 13 18 lim ut[1] ut[2] 0 100 box n conn ! col y 6 g\set box 15.5 24 7.8 12.8 lim ut[1] ut[2] * * box label "Wind speed (m/s)"/y label "Time (UT)"/x conn endif ! set /def set rece 2 find/proc cali if found.eq.0 then set rece 1 find/proc cali endif for i 1 to found-1 get cx_num[i] let scan1 scan let m i+1 get cx_num[m] let scan2 scan if i+2.lt.found then ! cal cold let o i+2 get cx_num[o] let scan3 scan endif if scan1.eq.scan2.and.scan2.ne.scan3 then get cx_num[m] endif if scan1.eq.scan2 then def real ws[nant] mws rws for j 1 to nant let ws[j] water[j] next compute mws mean ws compute rws rms ws let k k+1 def real nmwt[k] nrwt[k] ntim[k] npr[k] nhu[k] nte[k] if k.ne.1 then for l 1 to k-1 let nmwt[l] omwt[l] let nrwt[l] orwt[l] let ntim[l] otim[l] let npr[l] opr[l] let nhu[l] ohu[l] let nte[l] ote[l] next del /var omwt orwt otim opr ohu ote endif let nmwt[k] mws let nrwt[k] rws let ntim[k] utobs let npr[k] pressure let nhu[k] humidity let nte[k] ambiant_t-273.15 def real omwt[k] orwt[k] otim[k] opr[k] ohu[k] ote[k] let omwt nmwt let orwt nrwt let otim ntim let opr npr let ohu nhu let ote nte del /var nmwt nrwt ntim ws mws rws npr nhu nte endif next ! for i 1 to k-1 let m i+1 if (otim[m].lt.otim[i]) let otim[m] otim[m]+2*pi endif next let otim otim*12|pi ! g\set box 3.5 12 2.6 7.6 lim /var otim omwt lim ut[1] ut[2] = = box label "Water (mm) / Rec."'receiver'/y label "Time (UT)"/x g\set mar 20 0 0.1 point otim omwt g\set mar 0 0 0 err y otim omwt orwt ! g\set mar 20 0 0.1 g\set box 3.5 12 13 18 if (logi.eq."OBSERVER") then col y 2 /file meteo.5 lim ut[1] ut[2] * * else lim /var otim opr box n endif point otim opr label "Pressure (mbar)"/y ! g\set box 3.5 12 7.8 12.8 if (logi.eq."OBSERVER") then col y 3 /file meteo.5 lim ut[1] ut[2] * * else lim /var otim ote box n endif point otim ote label "Temperature (C\Uo)"/y ! g\set box 15.5 24 13 18 if (logi.eq."OBSERVER") then col y 4 /file meteo.5 lim ut[1] ut[2] 0 100 else lim /var otim ohu box n endif point otim ohu label "Humidity (%)"/y ! g\set mar 0 0 0.3 g\set exp 1. sys "rm -f meteo.ps meteo.3 meteo.4 meteo.5" hard meteo.ps /dev ps fast