#!/bin/ksh 
export LANG=en_US.UTF-8
export LC_ALL=en_US.UTF-8

datai=2022120100
dataf=2023030100

out=/share/ioper/tempo/MERGE/GOES/OBS/USP

while [ $datai -le $dataf ] ; do
 echo $datai
 mkdir -p ${out}/${datai:0:6}
 cd ${out}/${datai:0:6}
 wget -q -N http://compmodel.iag.usp.br/dados/superficie/${datai:0:6}/is${datai:0:4}-${datai:4:2}-${datai:6:2}-${datai:8:2}00
 sed '1,10d' is${datai:0:4}-${datai:4:2}-${datai:6:2}-${datai:8:2}00 |grep "${datai:0:4} ${datai:4:2} ${datai:6:2} ${datai:8:2}00 A"| awk -v nstid=INMET '{ if($26 != -999.0 && $26 != "" ) print $5"  "nstid"  "$7"  "$6"  "$26 }' | sort -u > PREC_USP_$datai.txt
 rm -f is${datai:0:4}-${datai:4:2}-${datai:6:2}-${datai:8:2}00
 datai=$(date -u -d "${datai:0:8} ${datai:8:2}:00 1 hours  " "+%Y%m%d%H")
done
exit
