Senhores, bom dia!
Mateus, seria bom você disponibilizar o arquivo de dados ('
novembros20052014.nc') pra podermos testar. Mas de antemão posso te dizer que uma opção é transformar a matriz
media_ColumnAmountO3 em um objeto {raster} e salvar num dos formatos possíveis (ascii, geotiff, ...).
Sugiro adicionar o código abaixo no teu script e testar. Não garanto que vá funcionar, pois sem os dados é um tiro no escuro. Retorne o resultado na lista...
### <code r>
...
require(raster)
tmp <- raster(media_ColumnAmountO3)
r <- t(flip(tmp, 1)) ### para colocar na posição correta!
plot(r)
extent(r) <- c(range(lon), range(lat))
plot(r)
writeFormats() # formatos possíveis
writeRaster(r, "nomearq.tif", "GTiff")
test <- raster("nomearq.tif")
plot(test)
### </code>