Última revisión de 09:36 11 ene 2023

Creación de un nuevo domino de simulación

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

[editar] 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"
88,89c87
< setenv EMISINV_A "$CASEINPUTS/np_oilgas/NOxAMBAtot_anual.nc"
< setenv EMISINV_B "$CASEINPUTS/np_oilgas/VOCAMBAtot_anual.nc"
---
> 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.csv

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"

Cargamos al SMOKE las emisiones regrilladas:

ln -s /home/solange.luque/estudios/CMAQ-ideal/DATA/emisiones/*nc ./

Archivo con la lista de las emisiones:

cat arinv_np/home/solange.luque/MODELOS/SMOKE/DATA/AMBAideal/inputs/np_oilgas/arinv_np_oilgas_AMBAideal.lst
#LIST GRID
/home/solange.luque/MODELOS/SMOKE/DATA/AMBAideal/inputs/np_oilgas/NOxTOTideal1x1.nc
/home/solange.luque/MODELOS/SMOKE/DATA/AMBAideal/inputs/np_oilgas/VOCTOTideal1x1.nc

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

[editar] 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

[editar] Cambiando archivos de entrada espaciales

[editar] 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

[editar] 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

[editar] 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

[editar] 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