/NuevoDominio
Creación de un nuevo domino de simulación
Siguiendo las instrucciones del manual de SMOKE v4.5 página 146
Contenido |
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.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"
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
Cargamos al SMOKE las emisiones regrilladas:
ln -s /home/solange.luque/estudios/CMAQ-ideal/DATA/emisiones/*nc ./
Archivo con la lista de las emsisiones
cat arinv_np/home/solange.luque/MODELOS/SMOKE/DATA/AMBAideal/inputs/np_oilgas/_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
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
