[R-br] Gerar arquivo de saída netCDF no R

[Thiago V. dos Santos]

A resposta já foi postada antes - veja a sua própria pergunta intitulada "Criar arquvo de saída a partir de dados NetCDF". 

Nenhuma das respostas funcionou?

Em que ponto do seu código está a função "create" para nós tentarmos te ajudar com o erro?
Greetings,
 -- Thiago V. dos Santos

PhD student
Land and Atmospheric Science
University of Minnesota


On Thursday, February 25, 2016 1:00 PM, Mateus Dias Nunes <nunes.mateusdias em gmail.com> wrote:



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             
____________________________________________________________________________
_______________________________________________
R-br mailing list
R-br em listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne� c�igo m�imo reproduz�el.