/NuevoDominio

Creación de un nuevo domino de simulación

Siguiendo las instrucciones del manual de SMOKE v4.5 página 146

Agregar proyección

Añadir en el archivo GRIDDESC (del .csh a utilizar) se tiene que añadir la proyección

Partimos de un archivo existente y creamos uno para AMBA (en nuestro caso ${GE_DAT}=/home/solange.luque/MODELOS/SMOKE/DATA/):

cp ${GE_DAT}/gridding/griddesc_lambertonly_18jan2019_v7.txt ${GE_DAT}/gridding/griddesc_AMBA.txt

Creando nueva proyección con la configuración de 3-dominios llamada `ChemGBsAs_cmaq` (del namelist.wps de WRF):

&geogrid
 parent_id         =   1,   1,  2,
 parent_grid_ratio =   1,   5,  3,
 i_parent_start    =   1,  43, 40,
 j_parent_start    =   1,  42, 25,
 e_we              = 111, 121, 142,
 e_sn              = 101, 101, 127,
 geog_data_res     = '2m', 'modis_30s+2m', 'modis_30s+30s',
 dx                = 15000,
 dy                = 15000,
 map_proj          = 'lambert',
 ref_lat           = -34.6,
 ref_lon           = -58.4,
 truelat1          = -35.6,
 truelat2          = -33.6,
 stand_lon         = -58.4,
 geog_data_path    = '/media/lluis/ExtDiskA_ext4/DATA/WRF/geog'
 opt_geogrid_tbl_path = './'
/

Copiando archivo para saber cambios introducidos (partiendo de la versión gcc!)

cp ../directory_definitions.csh ../directory_definitions.csh.old
diff ../directory_definitions.csh ../directory_definitions.csh.old
34,36c34,36
< setenv REGION "AMBA 1km"
< setenv REGION_ABBREV "AMBA1" # affects filename labeling
< setenv REGION_IOAPI_GRIDNAME "AMBA_1k" # should match GRIDDESC
---
> setenv REGION "Continental US 12km large"
> setenv REGION_ABBREV "12US1" # affects filename labeling
> setenv REGION_IOAPI_GRIDNAME "12US1_459X299" # should match GRIDDESC

Cambiand archivos de entrada espaciales

SRGDESC

Archivo de descripción de las sustitutas espaciales.

Cómo no existen archivos para AMBA, tendremos que generarlos desde cero. Siguiendo el manual, utilizando el `MIMS Spatial Allocator` ( MSA) provisto desde CMAS

Instalación de MSA

En hydra en /home/solange.luque/libraries, siguiendo instalación v4.4

git clone https://github.com/CMASCenter/Spatial-Allocator.git
cd Spatial-Allocator/
cp bin/sa_setup.csh  bin/sa_setup.csh.old
diff bin/sa_setup.csh bin/sa_setup.csh.old 
2c2
< setenv SA_HOME /home/solange.luque/libraries/Spatial-Allocator
---
> setenv SA_HOME /proj/ie/apps/longleaf/sallocator/Spatial-Allocator

Siguiendo los pasos descritos en el Users' Guide

Para e MSA necesitamos shape files de las emsisiones

Emisiones NOx y VOC AMBA

Ejemplo de transformación de archivos de emisiones en netCDF a shape-file:

gdalwarp -of GTiff NETCDF:NOxAMBAtot.nc NOxAMBAtot.shp
gdalwarp -of GTiff NETCDF:VOCAMBAtot.nc VOCAMBAtot.shp