Senhores, boa tarde!

Outra opção:

intervalos = seq( trunc(min(ColumnAmountO3)), trunc(max(ColumnAmountO3)), 5 )
map( xlim=c(-100,-10), ylim=c(-60,10)) -> tmp
filled.contour(sort(lon), lat, media_ColumnAmountO3[order(lon),], 
               levels=intervalos, lwd=2, col=terrain.colors(14), 
               xlim=c(-100,-10), ylim=c(-60,10),
               plot.axes={axis(1); axis(2); lines(tmp)})
title( main="Campo médio de ozonio janeiros" )   # título do gráfico


​

​
================================================
Éder Comunello
Agronomist (UEM), MSc in Environ. Sciences (UEM)
DSc in Agricultural Systems Engineering (USP/Esalq)
Brazilian Agricultural Research Corporation (Embrapa)
Dourados, MS, Brazil |<O>|
================================================
GEO, -22.2752, -54.8182, 408m
UTC-04:00 / DST: UTC-03:00