/NuevoDominio

De Wikicima
(Diferencias entre revisiones)
Saltar a: navegación, buscar
(→‎Instalación de MSA)
(→‎AMBA ideal)
 
(No se muestran 9 ediciones intermedias realizadas por un usuario)
Línea 2: Línea 2:
   
 
Siguiendo las instrucciones del manual de [https://www.cmascenter.org/smoke/documentation/4.5/manual_smokev45.pdf SMOKE v4.5] página 146
  
Siguiendo las instrucciones del manual de [https://www.cmascenter.org/smoke/documentation/4.5/manual_smokev45.pdf 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:
  +
<pre>
  +
/home/solange.luque//MODELOS/SMOKE/DATA/AMBAideal
  +
</pre>
  +
  +
Desplegando SMOKE en el directorio copiando del contendio del SMOKE-TestCase
  +
<pre style="shell">
  +
cp /home/solange.luque//MODELOS/SMOKE/DATA/SMOKE-TestCase/* /home/solange.luque//MODELOS/SMOKE/DATA/AMBAideal
  +
</pre>
  +
  +
Modificamos la configuración:
  +
<pre style="shell">
  +
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"
  +
</pre>
  +
  +
<pre style="shell">
  +
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"
  +
</pre>
  +
  +
Al archivo <code>/home/solange.luque//MODELOS/SMOKE/DATA/smoke4.7/scripts/emf/smk_ar_annual_emf.csv</code>
  +
  +
<pre>
  +
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"
  +
</pre>
  +
Cargamos al SMOKE las emisiones regrilladas:
  +
<pre style="shell">
  +
ln -s /home/solange.luque/estudios/CMAQ-ideal/DATA/emisiones/*nc ./
  +
</pre>
  +
  +
Archivo con la lista de las emisiones:
  +
<pre>
  +
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
  +
</pre>
  +
  +
  +
Al correr tenemos un segmentation fault:
  +
<pre>
  +
./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
  +
</pre>
   
 
= Agregar proyección =
  
= Agregar proyección =
Línea 50: Línea 174:
 
</pre>
  
</pre>
   
−
= Cambiand archivos de entrada espaciales =
 +
= Cambiando archivos de entrada espaciales =
   
 
== SRGDESC ==
  
== SRGDESC ==
Línea 77: Línea 201:
 
==== Datos AMBA ====
  
==== Datos AMBA ====
 
Fuente de datos de capital en: [https://data.buenosaires.gob.ar/dataset/calles/resource/juqdkmgo-302-resource data.buenosaires.gob.ar]
  
Fuente de datos de capital en: [https://data.buenosaires.gob.ar/dataset/calles/resource/juqdkmgo-302-resource data.buenosaires.gob.ar]
  +
  +
Callejero CABA: [https://cdn.buenosaires.gob.ar/datosabiertos/datasets/jefatura-de-gabinete-de-ministros/calles/callejero.zip callejero.zip]
   
 
Para e MSA necesitamos shape files de las emsisiones
  
Para e MSA necesitamos shape files de las emsisiones

Ú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

Contenido

[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
Obtenido de «https://wiki.cima.fcen.uba.ar/mediawiki/index.php?title=/NuevoDominio&oldid=2254»
Herramientas personales