function main(args)

'reinit'
mes=subwrd(args,1)
df=subwrd(args,2)
var=subwrd(args,3)
ctl=subwrd(args,4)
out=subwrd(args,5)
ano=subwrd(args,6)

'open /share/ioper/tempo/MERGE/GPM/DAILY/MERGE_CPTEC_2022.ctl'
'sdfopen  'ctl

*while (ano<=ano)
'reset'
 if (mes=feb) 
   if ( (math_mod(ano,4)=0) & (math_mod(ano,400)!=0) | (math_mod(ano,400)=0) )  
     df=29
   else
     df=28
   endif
 endif

 'set dfile 2'
 'set lon -81.95  -34.65'
 'set lat -56.15   12.85'
 'set t 1' 
 'define tmt=ttmed(t=1,z=1)'
 'set dfile 1'
 'set lon -81.95  -34.65'
 'set lat -56.15   12.85' 
 'set time 12Z01'mes%ano
 'tm=lterp(tmt,prec.1(t=1))'
 'set gxout stat'
 'd tm'
 n_undef=sublin(result,7)
 n_valid=subwrd(n_undef,4)
 say n_valid
 if (n_valid >=171427)
  say ano"  "mes"  "df
 'set time 12Z01'mes%ano
 'set sdfwrite -nc4 -zip  -3dt 'out
 'sdfwrite tm'
 endif
 ano=ano+1
*endwhile
'quit'
return