for altura in 1000 1500 2000; do
    grdfft MAG.grd -V -M -L -C"$altura" -GMAG-ContPraCima${altura}.grd
    grd2cpt MAG-ContPraCima${altura}.grd -Crainbow -Z > MAG-ContPraCima${altura}.cpt
    grdimage MAG-ContPraCima${altura}.grd -JM8 -CMAG-ContPraCima${altura}.cpt -K -V \
    > MAG-ContPraCima${altura}.ps
    pscoast -B.25 -N2 -R -JM8 -O -K >> MAG-ContPraCima${altura}.ps
    psscale -E -D11/8/16/1.0 -CMAG-ContPraCima${altura}.cpt                            \
            -B100a50:"Continuacao pra cima, ${altura}m (nT)": -O                    \
            >> MAG-ContPraCima${altura}.ps

    # Mapa residual
    grdmath MAG.grd MAG-ContPraCima${altura}.grd SUB = MAG-ContPraCima${altura}-res.grd
    grd2cpt MAG-ContPraCima${altura}-res.grd -Crainbow -Z \
    > MAG-ContPraCima${altura}-res.cpt
    grdimage MAG-ContPraCima${altura}-res.grd -JM8 -CMAG-ContPraCima${altura}-res.cpt \
    -K -V > MAG-ContPraCima${altura}-residual.ps
    pscoast  -B.25 -N2 -R -JM8 -O -K >> MAG-ContPraCima-${altura}-residual.ps
    psscale -D11/8/16/1.0 -CMAG-ContPraCima${altura}-res.cpt                         \
            -B100a50:"Continuacao pra cima, ${altura}m -- Residual (nT)": -O         \
            >> MAG-ContPraCima${altura}-residual.ps
done