Jose, tu script es útil pero tiene dos detalles que sería bueno modificar. El primero es que la línea que define el valor NA es opcional en el formato y por ejemplo QGIS no la pone si es que no tienes NA cuando exportas al formato ascii. Entonces un if (~line) fprintf('NODATA_value -9999\n'). Lo otro es que la DEM en el archivo con formato GRASS no incluye la zona UTM, ya que tampoco viene en el formato ascii, entonces una cajita para agregar la zona UTM al output en formato GRASS no te cuesta nada!

Como has estado? Un saludo desde Escocia!

Caco