Ficheros de costas con provincias para Metview

Hola

Hasta ahora he estado usando los ficheros de costas con provincias para Metview de reciente creación mediante una colaboración entre Alberto Fernández y Alfons Callado, utilizando la opción siguiente, dentro de la macro correspondiente:

MAP_COASTLINE_RESOLUTION : 'HIGH', (o 'MEDIUM')
MAP_COASTLINE_PATH : '/home/tmp/stm/metview/bin'

En la carpeta a que hace referencia el path guardaba los archivos NOAAH_r8.bin y NOAAM_r8.bin, correspondientes a los mapas con límites provinciales o de C.C.A.A.

Ahora estoy empezando a trabajar en sureste y he comprobado que esto no me funciona, sea cual sea el path que use para guardar los ficheros de provincias. Por defecto los toma de usr/sprod/metview/etc_magics/coast, si no me equivoco, pero los que hay allí, aunque tienen resolución diferente (alta, media y baja) no tienen límites provinciales o de comunidades.

¿Me podéis decir alguno como resolver esto?

Muchas gracias de antemano

Alex

Ya está resuelto este problema. No hay ningún problema con el COASTLINE_PATH, pero los ficheros NOAAH_r8.bin y NOAAM_r8.bin que estaba usando no servían, porque estaban compilados para otra arquitectura diferente de la del cluster. Ya me ha pasado Carlos Santos los ficheros buenos y todo funciona correctamente.

He comprobado que algo no va bien del todo con ellos.
Al menos en la versión 3.11.5 de Metview instalada en nuestros PC/Linux/Ubuntu.
Cuando usamos estos ficheros, la opción LAND_SHADE deja de funcionar.
No se si en el Cluster Linux ocurre lo mismo.

Un saludo,
Ernesto


P.D.: Alex, quizá no deberías haber cerrado tu hilo...

Pues lo vuelvo a abrir y continuamos el debate allí. Voy a intentar fusionar los temas, a ver si lo consigo.

P.S. Ya está hecho.

Los siguientes documentos podrían serviros de guia:

El último es más técnico y trata de la generación de los bin con un 'make' que se adapta al sistema operativo.
Lo ha escrito Alfons Callado y posiblemente él os pueda guiar en los detalles.

manual_pdf
make_pdf

Gracias Alberto y Alfonso por los enlaces facilitados.
He usado vuestra aplicación para reconstruir los ficheros de líneas de costa en PC/Linux/Ubuntu/PGI, pero el problema persiste.
Os cuento mi experiencia.

Estoy interesado en areas geográficas relativamente pequeñas en las que la resolución de los contornos costeros resulta relevante.
Este es el caso de la salidas de los modelos de predicción de oleaje. GMT utiliza la misma fuente de datos que Magics/Metview (GSHHS) para las lineas de costa, sin embargo el primero permite plotear contornos con muy alta resolución sin grandes problemas.

El fichero de costas de alta resolución ('NOAAH_r8.bin') que viene por defecto en Metview 3.11.5 tiene un factor de reducción de polígonos de 0.04. El resultado con los contornos insulares resulta muy pobre:


He usado vuestro programa para aumentar la resolución mediante la reducción del filtrado. Los resultados han sido algo desconcertantes.
El límite parece estar alrededor de 0.017. Para este valor se obtiene una mejora en los contornos insulares pero sigue siendo deficiente:


Por debajo de 0.017 los contornos costeros dejan de dibujarse y Metview se queja amargamente en la terminal.

Hay que citar además que con reducción de polígonos por debajo de 0.02, el LAND_SHADE a escala global también da problemas:



Curiosamente este fenómeno parece depender del área representada porque desaparece cuando realizamos un zoom.

Resumiendo: parece que existe un límite efectivo a la resolución de los contornos costeros en Magics/Metview y que éste no está relacionado con la resolución real de los datos de origen.

Cualquier sugerencia será bienvenida.

Un saludo,
Ernesto

Hola Ernesto,

Vamos por partes.

Efectivamente, MetView tiene una limitación en
el número máximo de polígonos a dibujar, alrededor
de los 23.000, que se corresponden al famoso
0.017 ... se supone que algún día el ECMWF solucionará
este problema.

Te adjunto "provimundo.coast", una base en formato
Magics/Asccii que creamos conjuntamente entre Alberto
y quién te escribe. Tiene las islas Canarias, Baleares,
y la península con provincias a una resolución bastante
más buena a la que muestras en el foro, y el resto del
mundo al famosillo 0.017. Te adjunto las imágenes
"canarias_*.1.jpeg" para que lo veas.

Si esa resolución no es suficiente, y sólo deseas las
islas Canarias, tienes otra solución:

(1) Coge la base original ascii a muy buena resolución
(ej. de la NOAA: GSHHS).

(2) Elimina TODOS los polígonos menos los de las islas
Canarias, y adecua la cabecera (primera línia del
fichero) a los cambios:

NOAA_Canarias.ascii

(3) Transfórmalo con "t_convnoaa_AEMET" al formato
binario que necesita MetView sin reducción
alguna (-e 0.0):

t_convnoaa_AEMET -fNE -e 0.0 NOAA_Canarias.ascii
NOAA_Canarias.bin

(4) Y ya puedes plotear con MetView las Canarias
a la misma alta resolución original NOAA con
"NOAA_Canarias.bin". La salvedad es
que el resto del mundo será océano.

Si lo creéis oportuno, podéis subir este correo
al foro.

Saludos,

Alfons Callado Pallarès

Muchísimas gracias por tan detallada y eficaz respuesta.
No conocía el límite de los 23.000 polígonos en Metview y me pregunto si será posible modificar el dimensionado en el código.

a) He probado la solución de eliminar polígonos dejando solo el archipiélago canario sin filtrado y funciona perfectamente tal y como sugieres. Eso resuelve el problema que planteaba ya que a la hora de plotear salidas de modelos, la extensión geográfica está prefijada
y resultaría sencillo construir ficheros específicos que contienen solo los polígonos que se muestran, a resolución plena.

b) Sin embargo, no he tenido buenos resultados con el fichero 'provimundo.coast' que facilitas. ¿Quizá no estoy usándolo correctamente?. He generado una versión binaria para Metview con él y los problemas son exactamente los mismos que he comentado en el otro post (LAND_SHADE erróneo dependiendo de la escala geográfica, desaparición de líneas de costas y mensaje de 'Magics Engine Crashed').


Un saludo,
Ernesto