Revisión de 11:06 4 ene 2023

Creación de un nuevo domino de simulación

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

AMBA ideal

Creando el directorio de intentar correr el caso ideal con las emisiones creados por nosotres a partir de un laburo de Sol

Directorio de trabajo:

/home/solange.luque//MODELOS/SMOKE/DATA/AMBAideal

Desplegando SMOKE en el directorio copiando del contendio del SMOKE-TestCase

cp /home/solange.luque//MODELOS/SMOKE/DATA/SMOKE-TestCase/* /home/solange.luque//MODELOS/SMOKE/DATA/AMBAideal

Modificamos la configuración:

cd AMBAideal/scripts
diff ../directory_definitions.csh /home/solange.luque/MODELOS/SMOKE/DATA/SMOKE-TestCase/
scripts/directory_definitions.csh
19c19
< setenv MET_ROOT "/home/solange.luque/MODELOS/SMOKE/DATA/AMBAideal/inputs/biogenics"
---
> setenv MET_ROOT "/home/solange.luque/MODELOS/SMOKE/DATA/SMOKE-TestCase/inputs/biogenics"
24c24
< setenv MET_ROOT_3D "/home/solange.luque/MODELOS/SMOKE/DATA/AMBAideal/inputs/biogenics"
---
> setenv MET_ROOT_3D "/home/solange.luque/MODELOS/SMOKE/DATA/SMOKE-TestCase/inputs/biogenics"
30c30
< setenv CASE "AMBAideal"
---
> setenv CASE "SMOKE-TestCase"
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
46,47c46
< #setenv SMOKE_LOCATION "/home/solange.luque/MODELOS/SMOKE/gcc/Linux2_x86_64"
< setenv SMOKE_LOCATION "/home/solange.luque/MODELOS/SMOKE_dbg/gcc/Linux2_x86_64"
---
> setenv SMOKE_LOCATION "/home/solange.luque/MODELOS/SMOKE/gcc/Linux2_x86_64"

cd nonpoint
diff Annual_np_oilgas_12US1_2017gb_17j_TestCase.csh /home/solange.luque/MODELOS/SMOKE/DATA/SMOKE-TestCase/
scripts/nonpoint/Annual_np_oilgas_12US1_2017gb_17j_TestCase.csh 
7d6
< setenv GRIDMASK "/home/solange.luque/estudios/CMAQ-ideal/DATA/GRIDMAKS/geocode_AMBA.ncf"
33c32
< setenv RUN_MONTHS "11"
---
> setenv RUN_MONTHS "7"
39,41c38,40
< setenv BASE_YEAR "2012"
< setenv EPI_STDATE_TIME "${BASE_YEAR}-11-05 00:00:00.0"
< setenv EPI_ENDATE_TIME "${BASE_YEAR}-11-05 23:59:00.0"
---
> setenv BASE_YEAR "2017"
> setenv EPI_STDATE_TIME "${BASE_YEAR}-01-01 00:00:00.0"
> setenv EPI_ENDATE_TIME "${BASE_YEAR}-12-31 23:59:00.0"
56c55
< setenv COSTCY "${GE_DAT}/costcy_AMBA.txt"
---
> setenv COSTCY "${GE_DAT}/costcy_for_2017platform_24apr2020_nf_v1.txt"
58c57
< setenv GRIDDESC "${GE_DAT}/gridding/griddesc_AMBA.txt"
---
> setenv GRIDDESC "${GE_DAT}/gridding/griddesc_lambertonly_18jan2019_v7.txt"
88c87
< setenv EMISINV_A "$CASEINPUTS/np_oilgas/np_oilgas_AMBA.csv"
---
> setenv EMISINV_A "$CASEINPUTS/np_oilgas/np_oilgas_2017NEI_NONPOINT_20200501_04may2020_nf_v1.csv"

Al archivo /home/solange.luque//MODELOS/SMOKE/DATA/smoke4.7/scripts/emf/smk_ar_annual_emf.cs

setenv IMPORT_GRDNETCDF_YN Y  # esto es en teoría para que tome el inventario pre-gridded
setenv HOUR_SPECIFIC_YN    N  # esto en teoría para que sepa el inventario es horario
setenv DAY_SPECIFIC_YN     N  #esto si es por dia
#setenv USE_EXP_GEO_CODES   Y  #esto para que use los files GEOCODE_LEVELi en vez de costcy

setenv NETCDF_POL_UNIT "kg m-2 s-1"
setenv NETCDF_INV_YEAR "2012"

Al correr tenemos un segmentation fault:

./Annual_np_oilgas_12US1_2017gb_17j_TestCase.csh >& run_nonpoint.log
vim run_nonpoint.log
(...)
Backtrace for this error:
#0  0x1519b474abd0 in ???
#1  0x1519b4749e25 in ???
#2  0x1519b4428d5f in ???
#3  0x1519b4494e99 in ???
#4  0x55c093c2f155 in ???
#5  0x55c093a6fb5b in rdgrdncf_
        at /home/solange.luque/MODELOS/SMOKE_dbg/gcc/src/smkinven/rdgrdncf.f:389
#6  0x55c093aa29b5 in smkinven
        at /home/solange.luque/MODELOS/SMOKE_dbg/gcc/src/smkinven/smkinven.f:368
#7  0x55c093a268c8 in main
        at /home/solange.luque/MODELOS/SMOKE_dbg/gcc/src/smkinven/smkinven.f:61
Segmentation fault
0.044u 0.016s 0:00.06 83.3%     0+0k 0+24io 0pf+0w

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

Cambiando 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 sección Chapter 5. Surrogate Tool

Datos AMBA

Fuente de datos de capital en: data.buenosaires.gob.ar

Callejero CABA: callejero.zip

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