Olá, eu preciso gerar um arquivo de saída .nc (netcdf) dos resultados da média da coluna de ozônio que carregue as informações de lat, lon e valor em cada ponto de grade da região analisada. Já consegui gerar um arquivo .txt, porém o programa que utilizo não lê arquivos modo texto. abaixo o arquivo de entrada também .nc Tentei utilizar a função create da biblioteca net https://www.dropbox.com/s/ntvmtx7wonwo839/teste.nc?dl=0 library(maps) library(ncdf4) dados <- nc_open('teste.nc') # lendo coordenadas espaço-temporal 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, essaconversão é feita para plotagem sobre o mapa for (i in 1:dim(lon)) { if (lon[i] > 180) { lon[i] <- lon[i]-360 } } desde já agradeço ____________________________________________________________________________ MATEUS DIAS NUNES MESTRANDO DO PROGRAMA DE PÓS-GRADUAÇÃO EM METEOROLOGIA - PPGMET UNIVERSIDADE FEDERAL DE PELOTAS - UFPEL TELEFONE: +55 (53) 81125154 ____________________________________________________________________________