<div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif">Olá, <br><br></div><div class="gmail_default" style="font-family:verdana,sans-serif">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.<br></div><div class="gmail_default" style="font-family:verdana,sans-serif">Já consegui gerar um arquivo .txt, porém o programa que utilizo não lê arquivos modo texto.<br>abaixo o arquivo de entrada também .nc<br>Tentei utilizar a função create da biblioteca net<br><a href="https://www.dropbox.com/s/ntvmtx7wonwo839/teste.nc?dl=0">https://www.dropbox.com/s/ntvmtx7wonwo839/teste.nc?dl=0</a><br></div><div class="gmail_default" style="font-family:verdana,sans-serif"><br>library(maps)<br>library(ncdf4)<br><br><br> dados <- nc_open('<a href="http://teste.nc">teste.nc</a>')<br> # lendo coordenadas espaço-temporal<br> lat <- ncvar_get( dados, 'lat' )<br> lon <- ncvar_get( dados, 'lon' )<br> time <- ncvar_get( dados, 'time' )<br><br><br> # lendo dados coluna total de Ozônio<br> ColumnAmountO3 <- ncvar_get( dados, 'ColumnAmountO3' )<br><br> # dimensoes da variavel ColumnAmountO3<br> dims_ColumnAmountO3 <- dim(ColumnAmountO3)<br><br> # tornando o arranjo 3D (ColumnAmountO3) em um 2D, organizado em ptos de grade X tempo<br><br> dim(ColumnAmountO3) <- c(dims_ColumnAmountO3[1]*dims_ColumnAmountO3[2], dims_ColumnAmountO3[3] )<br><br> # calculando a média e retornado-a em 2D<br><br> media_ColumnAmountO3 <- rowMeans( ColumnAmountO3)<br> dim(media_ColumnAmountO3) <- c( dims_ColumnAmountO3[1],dims_ColumnAmountO3[2] )<br><br> <br> # longitude varia de 0 a 360, convertendo para -180 a 180, essaconversão é feita para plotagem sobre o mapa<br> for (i in 1:dim(lon)) { if (lon[i] > 180) { lon[i] <- lon[i]-360 } }<br><br><br></div><div class="gmail_default" style="font-family:verdana,sans-serif">desde já agradeço<br clear="all"></div><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><font face="tahoma, sans-serif">____________________________________________________________________________</font></div><font face="tahoma, sans-serif"><div dir="ltr"><font face="tahoma, sans-serif"><br></font></div>MATEUS DIAS NUNES</font><div><div><font face="tahoma, sans-serif">MESTRANDO DO PROGRAMA DE PÓS-GRADUAÇÃO EM METEOROLOGIA - PPGMET</font></div><div><span style="font-family:tahoma,sans-serif">UNIVERSIDADE FEDERAL DE PELOTAS - UFPEL</span><font face="tahoma, sans-serif"><br></font></div><div><font face="tahoma, sans-serif"><span style="background-color:rgb(255,255,255)"><span style="color:rgb(0,0,0)">TELEFONE: +55 (53) 81125154           </span>  </span><br>____________________________________________<span style="background-color:rgb(255,255,255)"><span></span></span>________________________________</font></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
</div>