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 ____________________________________________________________________________