[gradsusr] extract to ascii nesting while

[RafaelCB]

Espero que ajude!

cit=read('o_estac.dat')
lin=sublin(cit,2)
_sig=subwrd(lin,1)
#say "cit="cit
#say "lin="lin
#say "sig="_sig
la=subwrd(lin,2)
lo=subwrd(lin,3)

'd gr2stn(t2.1,min.2)-273.15'
t2m=sublin(result,_n2)
#say "t2m="t2m
anl=sublin(result,_n3)


'd gr2stn(rh.1,min.2)'
rh=sublin(result,_n2)
#say "rh="rh
'd gr2stn(u10.1,min.2)'
u1=sublin(result,_n2)
#say "u1="u1
'd gr2stn(v10.1,min.2)'
v1=sublin(result,_n2)
#say "v1="v1
'd gr2stn(slvl.1,min.2)'
slp1=sublin(result,_n2)
#say "slp1="slp1

#'define chva=(RAINC.1(t='n')+RAINNC.1(t='n'))'
#'define chvb=(RAINC.1(t='n-4')+RAINNC.1(t='n-4'))'
#'d gr2stn((rainc.1(t='n')+RAINNC.1(t='n')),min.2)'
'd
gr2stn((rainc.1(t='n')+RAINNC.1(t='n'))-(RAINC.1(t='n-4')+RAINNC.1(t='n-4')),min.2)'
#'d gr2stn((chva-chvb),min.2)'
chv1=sublin(result,_n2)
chv1=subwrd(chv1,1)
say "chv1="chv1

if(chv1<0.1)
 chv1=0
endif
say "chv1if="chv1


#'!printf "%4.4s  %2.2s  %2.2s  %2.2s     %5.5s %6.6s %6.6s %6.6s %6.6s
%6.6s %6.6s %6.6s\n" 'ano' 'mes' 'dia' 'hora' '_sig' 'lo ' 'la' 't2m' 'rh '
'u1 ' 'v1 ' 'slp1' >> verific.dat'
'!echo
'ano','mes','dia','hora','_sig','lo','la','t2m','rh','u1','v1','slp1','chv1'
| sed s@\ @@g >> verific.dat'
*'!printf "\n"                        >> verific.dat'

-- 
Rafael Censi Borges
Téc.Meteorologia - Epagri/Ciram
+55 (48)3239-8064/8066
http://ciram.epagri.sc.gov.br

2011/1/5 claudiomet <claudiomet at gmail.com>

>
> Greetings GrADS users !!!
>
> I need to extract values from reanalysis grib file into ascii file
> the format of the ascii file must be
>
> date,lat,lon,lev,var1,var2,var3
>
> varying date and lev (nesting while) for each lat and lon
>
> ¿How to do this?
> --
> Claudio Cortes
> +56 (2) 2994121
>
> Meteorologo
> Laboratorio de Investigación, Innovación e Informática Aplicada (L3IA)
> Unidad de Modelacion y Gestion de la Calidad del Aire (UMGCA)
> Centro Nacional del Medio Ambiente (CENMA)
>
> --
> Claudio Cortes
> +56 (2) 2994121
>
> Meteorologist
> Laboratory Research, Innovation and Applied Informatics
> Modeling and Air Quality Management Unit
> National Enviroment Center, Chile (CENMA)
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110105/7f07b4fc/attachment-0003.html