Error in daily ptfire

SMOKE Emissions Inventory
1

Hello everyone,

I’m running SMOKE4.5 with NEI-2015 platform. After successful running of one-time ptfire sector, I face an error in daily ptfire section (in smkinven log file) about GENPDOUT subroutine. This is lat 10 lines of the log file:

 Value for HOURLY_TO_PROFILE not defined;returning default:   FALSE
 Reading day  -specific data...
 Successful open for emissions file: 
 /scratch/se75432/workDir/emis_2015/2015fd_cb6_15j/inputs/ptfire/ptday_ptfire_2015_large_split_ff10_10oct2017_v0
 Successful open for emissions file: 
 /scratch/se75432/workDir/emis_2015/2015fd_cb6_15j/inputs/ptfire/ptday_ptfire_2015_large_split_ff10_prevdec_10oct2017_v0
 
 *** ERROR ABORT in subroutine GENPDOUT
 No period-specific sources found in input file

Does anyone have any idea how should I deal with it?
Thank you in advance for your help.

Best,
Soroush

2

Just a reminder!
I’d appreciate it if someone can give me a hint that what might be wrong resulted in the above error.

3

@Soroush. Sorry for the late response. This error normally occurs when there is no matched sources between annual (PTINV) and daily (PTDAY) fire inventories based on FIPS, SCC, FireID and other source characteristics. I would suggest to check those characteristics between the files and make sure they are shared correctly.

PTINV: https://www.cmascenter.org/smoke/documentation/4.7/html/ch08s02s08.html
PTDAY:https://www.cmascenter.org/smoke/documentation/4.7/html/ch08s02s06.html

4

Thank you for the reply @bbaek.
I checked my annual and daily inventories. As expected the headers are different, but SCC and other characteristics are the same in both inventories.
Below is the first few lines of annual and daily inventories:

PTINV:

#DESC Apportioned large fires to 12km grid cells based on GeoMAC polygon with circle fallback
#EXPORT_DATE=Wed Feb 14 14:35:18 EST 2018
#EXPORT_VERSION_NAME=Initial Version
#EXPORT_VERSION_NUMBER=0
country_cd,region_cd,tribal_code,facility_id,unit_id,rel_point_id,process_id,agy_facility_id,agy_unit_id,agy_rel_point_id,agy_process_id,scc,poll,ann_value,ann_pct_red,facility_name,erptype,stkhgt,stkdiam,stktemp,stkflow,stkvel,naics,longitude,latitude,ll_datum,horiz_coll_mthd,design_capacity,design_capacity_units,reg_codes,fac_source_type,unit_type_code,control_ids,control_measures,current_cost,cumulative_cost,projection_factor,submitter_id,calc_method,data_set_id,facil_category_code,oris_facility_code,oris_boiler_id,ipm_yn,calc_year,date_updated,fug_height,fug_width_ydim,fug_length_xdim,fug_angle,zipcode,annual_avg_hours_per_year,jan_value,feb_value,mar_value,apr_value,may_value,jun_value,jul_value,aug_value,sep_value,oct_value,nov_value,dec_value,jan_pctred,feb_pctred,mar_pctred,apr_pctred,may_pctred,jun_pctred,jul_pctred,aug_pctred,sep_pctred,oct_pctred,nov_pctred,dec_pctred,comment
“US”,“01001”,“C12915437539758390”,“180”,“b”,“c”,“2811015001”,“ACRESBURNED”,121.999999700000004,“SF11E1373567”,1,1,1,1,1,-86.8620000000000232,32.6180000000000021,“2015”,“Hms006400 Fire”
“US”,“01001”,“C12915437539758390”,“180”,“b”,“c”,“2811015001”,“CO”,16.1761624000000026,“SF11E1373567”,1,1,1,1,1,-86.8620000000000232,32.6180000000000021,“2015”,“Hms006400 Fire”
“US”,“01001”,“C12915437539758390”,“180”,“b”,“c”,“2811015001”,“CO2”,107.246757790000004,“SF11E1373567”,1,1,1,1,1,-86.8620000000000232,32.6180000000000021,“2015”,“Hms006400 Fire”
“US”,“01001”,“C12915437539758390”,“180”,“b”,“c”,“2811015001”,“HFLUX”,1000,“SF11E1373567”,1,1,1,1,1,-86.8620000000000232,32.6180000000000021,“2015”,“Hms006400 Fire”
“US”,“01001”,“C12915437539758390”,“180”,“b”,“c”,“2811015001”,“NH3”,0.262563029999999975,“SF11E1373567”,1,1,1,1,1,-86.8620000000000232,32.6180000000000021,“2015”,“Hms006400 Fire”

PTDAY:

#FORMAT=FF10_DAILY_POINT
#COUNTRY=US
#YEAR=2015
#DESC Wildland fires for 2015 provided by STI
#DESC Recalculated smoldering SCCs to include only residual emissions. Moved balance of smoldering emissions to flaming.
#DESC Apportioned large fires to 12km grid cells based on GeoMAC polygon with circle fallback
#EXPORT_DATE=Wed Feb 14 14:35:06 EST 2018
#EXPORT_VERSION_NAME=Initial Version
#EXPORT_VERSION_NUMBER=0
country_cd,region_cd,tribal_code,facility_id,unit_id,rel_point_id,process_id,scc,poll,op_type_cd,calc_method,date_updated,monthnum,monthtot,dayval1,dayval2,dayval3,dayval4,dayval5,dayval6,dayval7,dayval8,dayval9,dayval10,dayval11,dayval12,dayval13,dayval14,dayval15,dayval16,dayval17,dayval18,dayval19,dayval20,dayval21,dayval22,dayval23,dayval24,dayval25,dayval26,dayval27,dayval28,dayval29,dayval30,dayval31,comment
“US”,“01001”,“C12915437539758390”,“180”,“b”,“c”,“2811015001”,“ACRESBURNED”,1,121.999999700000004,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,121.999999700000004,0,0,0,0,
“US”,“01001”,“C12915437539758390”,“180”,“b”,“c”,“2811015001”,“CO”,1,16.1761624000000026,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,16.1761624000000026,0,0,0,0,
“US”,“01001”,“C12915437539758390”,“180”,“b”,“c”,“2811015001”,“CO2”,1,107.246757790000004,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,107.246757790000004,0,0,0,0,
“US”,“01001”,“C12915437539758390”,“180”,“b”,“c”,“2811015001”,“HFLUX”,1,1000,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1000,0,0,0,0,
“US”,“01001”,“C12915437539758390”,“180”,“b”,“c”,“2811015001”,“NH3”,1,0.262563029999999975,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.262563029999999975,0,0,0,0,

Do you know what else might be wrong in this case?

Thanks,
Soroush

5

Hi @bbaek.
I’m still confused about the error I’m receiving, and couldn’t find a solution yet.
I also tried turning off smk_reports for ptfire, but it didn’t help.
Do you have any idea what else might be the reason for this error?

Thanks,
Soroush

6

To debug this, we would need the following information:

The run script and the entire smkinven log file that contains the ERROR message.

Also the ptinv file that was posted is the FORMAT, COUNTRY, and YEAR part of the header.

You should verify that the FORMAT is FF10_POINT, COUNTRY is US, and YEAR is 2015.

image001.jpg

7

You need to make sure that your daily inventory data covers your modeling period. If not, Smkinven will also give you the same kind of error message you are getting

8

Thanks for the reply.
Please find run script and smkinven log file attached (the format of log file is changed to text file to be uploaded).
Also, I checked the ptinv headers. FORMAT, COUNTRY, and YEAR is set to values you mentioned (it just was not copied in the earlier post).

Annual_ptfire_daily_12US2_2015fd_cb6_15j.csh (6.0 KB)
smkinven_ptfire_jul_2015fd_cb6_15j.log.txt (804.4 KB)

9

Thanks @bbaek.
I checked daily inventory data. It does cover my modeling period. So, I don’t think data coverage be the problem.

10

I just checked the format of your inventory data (PTINV and PTDAY). There is a missing data in your data for “tribal_code” from both files. If there is no data to enter, then make sure you put a comma between the columns instead of skipping it since these input files are in CSV format. Please put them in excel to make sure all data columns are in right place.

11

I checked my inventories again. The original inventory files are correct (values belong to correct columns in excel), but I don’t know why when I copy them here with blockquote, the format changes.
Thank you though.

12

Perhaps the file needs to be emailed.

13

Just to update this topic. The problem is resolved after switching from self-compiled to pre-compiled SMOKE executables.