I have one question about the EPIC Site Information Generation Tool. For the BELD4 file for the domain (beld4_cmaq12km_2006.nc), is it packaged in the database of Spatial Allocator or I need to generate that using landuseTool_WRFCMAQ_BELD4.csh? I try to find this file in data folder of Spatial Allocator and also FEST-C data but I didn’t find that. Could you offer more information?
Also, can I use the NLCD 2011 (2016 edition, the new version) to generate the files using Spatial Allocator for FESTC since the NLCD website only offers the new edition data, not the 2011 edition (I couldn’t find a place to download the old version of NLCD 2011)?
Any feedback you can give would be greatly appreciated.
Thanks.
I don’t think that beld4_cmaq12km_2006.nc is packed in Spatial Allocator package.
In festc, this file can be created under a scenario you created ( based on your research domain ) in festc interface.
If you already downloaded 2011 NLCD and MODIS data, you can use them to create the *nc file, doesn’t matter 2016 version or older version. The utility in festc interface is the first one "BELD4 Data Generation. " The program used is: $SA_HOME/bin/64bits/computeGridLandUse_beld4.exe
Thank you for your answer. I have already downloaded 2011 NLCD and MODISv6 data for creating the beld4 NC file for the CONUS domain and used them to run computeGridLandUse_beld4.exe. However, I got some problem (I also see some people posted a similar question before but without any solution after that):
Projecting Shapefile…
Projecting 12US1.shp to +proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=23 +lon_0=-96 +x_0=0 +y_0=0 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs and stored it in a temp file: tmp62Input file size is 161190, 104424
ERROR 3: Failed to write scanline 7743 to file.
ERROR 1: An error occured while writing a dirty block
ERROR 4: Unable to open tmp1093988704_img.bil to obtain file list.
0…6063271_prj.shp
/home/ll70/models/Spatial-Allocator/src/libs/gdal-2.0.2/local/bin/ogr2ogr -t_srs "+proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=23 +lon_0=-96 +x_0=0 +y_0=0 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs " tmp626063271_prj.shp 12US1.shp
Completed in projecting the shapefile.
Rasterizing projected Shapefile…
Rasterized grid Shapefile is stored in file: tmp1093988704_img.bil
/home/ll70/models/Spatial-Allocator/src/libs/gdal-2.0.2/local/bin/gdal_translate -ot UInt32 -a_ullr -2651655.000000 3745995.000000 2916795.000000 136185.000000 -scale 0 255 0 0 -outsize 185615 120327 -of EHdr /projects/dsc1/ll70/model/Spatial-Allocator/data/sat/NLCD2011/NLCD_2011_Land_Cover_L48_20190424.img tmp1093988704_img.bil
Error in system call: /home/ll70/models/Spatial-Allocator/src/libs/gdal-2.0.2/local/bin/gdal_translate -ot UInt32 -a_ullr -2651655.000000 3745995.000000 2916795.000000 136185.000000 -scale 0 255 0 0 -outsize 185615 120327 -of EHdr /projects/dsc1/ll70/model/Spatial-Allocator/data/sat/NLCD2011/NLCD_2011_Land_Cover_L48_20190424.img tmp1093988704_img.bil
I have no idea about the reason for this problem. Could you give some suggestion to solve this problem?
Any feedback you can give me would be greatly appreciated.
Thanks,
Lina