[R-br] Gerar arquivo de saída ascii

[Mateus Dias Nunes]

Olá preciso gerar um arquivo de saída ascii que contenha as informações de
lat e lon para este campo meteorológico.
Já consegui gerar o arquivo ascii, porém ele fica sem as informações de lat
(-60S,10N) e lon(-120W,10E)
abaixo o script.

library(maps)

library(ncdf4)


 dados <- nc_open('teste.nc')


lat <- ncvar_get( dados, 'lat' )

lon <- ncvar_get( dados, 'lon' )

time <- ncvar_get( dados, 'time' )


# lendo dados coluna total de Ozônio

ColumnAmountO3 <- ncvar_get( dados, 'ColumnAmountO3' )



# dimensoes da variavel ColumnAmountO3

dims_ColumnAmountO3 <- dim(ColumnAmountO3)



# tornando o arranjo 3D (ColumnAmountO3) em um 2D, organizado em ptos de
grade X tempo

dim(ColumnAmountO3) <- c( dims_ColumnAmountO3[1]*dims_ColumnAmountO3[2],
dims_ColumnAmountO3[3] )

# calculando a média e retornado-a em 2D

media_ColumnAmountO3 <- rowMeans( ColumnAmountO3)

dim(media_ColumnAmountO3) <- c( dims_ColumnAmountO3[1],
dims_ColumnAmountO3[2] )



#==========================================================================================================

# longitude varia de 0 a 360, convertendo para -180 a 180, essa conversão é
feita para plotagem sobre o mapa

for (i in 1:dim(lon)) { if (lon[i]>180) { lon[i] <- lon[i]-360 } }

plotando mapa da America do Sul

map(xlim=c(-90,-30), ylim=c(-60,10))

map.axes()              # plotando eixos

title( main="Campo médio de Ozônio " )   # título do gráfico


# definindo intervalo de 5 Dobson Units (DU)

intervalos = seq( trunc(min(ColumnAmountO3)), trunc(max(ColumnAmountO3)), 5
)


#PARA PLOTAR GRAFICO CONTORNOS EM PRETO E BRANCO

contour( sort(lon), lat, media_ColumnAmountO3[ order(lon), ], add=T,
levels=intervalos, lwd=2, labcex=1.3  )

teste=data.frame(media_ColumnAmountO3)

write.table(teste,file="teste.txt",row.names=TRUE,sep="",quote=FALSE)
____________________________________________________________________________

MATEUS DIAS NUNES
MESTRANDO DO PROGRAMA DE PÓS-GRADUAÇÃO EM METEOROLOGIA - PPGMET
UNIVERSIDADE FEDERAL DE PELOTAS - UFPEL
TELEFONE: +55 (53) 81125154
____________________________________________________________________________
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160221/3dedadde/attachment.html>