******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-29 23:05:06.975527 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-29 18:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 1 it: 0 ndt: 1 idate: 2019-10-29 18:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-29 23:15:07.200494 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-29 18:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 1 it: 0 ndt: 1 idate: 2019-10-29 18:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-29 23:25:11.217855 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-29 18:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 1 it: 0 ndt: 1 idate: 2019-10-29 18:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-29 23:35:06.042653 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-29 18:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 1 it: 0 ndt: 1 idate: 2019-10-29 18:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-29 23:45:07.564234 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-29 18:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 1 it: 0 ndt: 1 idate: 2019-10-29 18:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-29 23:55:11.417950 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-29 18:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 1 it: 0 ndt: 1 idate: 2019-10-29 18:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-30 00:05:06.506912 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-29 18:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 1 it: 0 ndt: 1 idate: 2019-10-29 18:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-30 00:15:06.545977 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-29 18:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 1 it: 0 ndt: 1 idate: 2019-10-29 18:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== Creating output file that is 413P x 346L. Processing input file /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif. Using internal nodata values (e.g. -3.40282e+38) for image /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 413P x 346L. Processing input file /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_72h.tif. Using internal nodata values (e.g. -3.40282e+38) for image /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_72h.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-30 00:25:11.713051 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-29 18:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 1 it: 0 ndt: 1 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist GDACS/1000624/1_HWRF 1000624 ndt: 1 it: 0 ndt: 1 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_72h.tif input file exist home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret 0 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-10-29 18:00:00 currdate 2019-10-29 18:00:00 ndt: 0 delta: 6 nt 72 alldate: DatetimeIndex(['2019-10-29 18:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files lon0,lon1,lat0,lat1 99.34 120.0 5.0 22.3 processing curr bull... 2019-10-29 18:00:00 99.34 120.0 99.34 120.0 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/ /mnt/input/grib_HWRF/2019/1000624/ input file /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind_stept0.tif -r bilinear -te 99.34 5.0 120.0 22.3 -dstnodata 0 Resample tif file /usr/bin/gdalwarp -r bilinear -te 99.34 5.0 120.0 22.3 -dstnodata 0 "/mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind_stept0.tif" outFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind_stept0.tif dd: 0.0500242130751 >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019102918/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019102918/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019102918/all_inpData.txt 99.3650121065 119.974987893 5.01663438257 22.2749878935 lon 413 lat 346 read TC track /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019102918/all_inpData.txt currDate 2019-10-29 18:00:00 f 1 nt-1 5 nt 6 currDate1 2019-10-29 18:00:00 2019-10-30 00:00:00 sdate 2019-10-30 00:00:00 2019-10-30 00:00:00 1 CAT ['TS', 'TS'] savemap /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind_stept0.jpg 99.34 120.0 99.34 120.0 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/ /mnt/input/grib_HWRF/2019/1000624/ input file /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_72h.tif /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind.tif -r bilinear -te 99.34 5.0 120.0 22.3 -dstnodata 0 Resample tif file /usr/bin/gdalwarp -r bilinear -te 99.34 5.0 120.0 22.3 -dstnodata 0 "/mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_72h.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind.tif" outFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind.tif dd: 0.0500242130751 >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019102918/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019102918/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019102918/all_inpData.txt 99.3650121065 119.974987893 5.01663438257 22.2749878935 lon 413 lat 346 read TC track /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019102918/all_inpData.txt currDate 2019-10-29 18:00:00 f 1 nt-1 71 nt 72 currDate1 2019-10-29 18:00:00 2019-11-01 18:00:00 sdate 2019-11-01 18:00:00 2019-11-01 18:00:00 12 CAT ['TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD'] savemap /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind.jpg processing all past bull only if Past=True... False no past >> 2.2. Process final data files compfile /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/20191029.18_Final_completed_wind.txt FINAL alldate: DatetimeIndex(['2019-10-29 18:00:00'], dtype='datetime64[ns]', freq='6H') alldate.shape 1 **FIRST cp /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind.tif /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif max file created /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif stormname 1000624/23w stormname 1000624/23w TCbullsource 23w >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019102918/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019102918/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019102918/all_inpData.txt 99.3900242131 120.0 4.99162227603 22.2499757869 lon 413 lat 346 read TC track /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019102918/all_inpData.txt currDate 2019-10-29 18:00:00 f 1 nt-1 71 nt 72 currDate1 2019-10-29 18:00:00 2019-11-01 18:00:00 sdate 2019-11-01 18:00:00 2019-11-01 18:00:00 12 CAT ['TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD'] savemap /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_FINAL.jpg ret: 0 ============================================ 3. Classify meteo + GDACS index score ============================================ >> 3.1. Classify curr + past forecast Creating output file that is 2479P x 2077L. Processing input file /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind_stept0.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind_stept0.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind_stept0.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18//wind10m_res_stept0.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/lspop20141.tif. Using internal nodata values (e.g. -2.14748e+09) for image /mnt/output/GDACS/DATA/lspop20141.tif. Copying nodata values from source /mnt/output/GDACS/DATA/lspop20141.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18//wind10m_popfile_stept0_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/countries.tif. Using internal nodata values (e.g. -32768) for image /mnt/output/GDACS/DATA/countries.tif. Copying nodata values from source /mnt/output/GDACS/DATA/countries.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18//wind10m_countryfile_stept0_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind_stept0.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18// PopFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18//wind10m_popfile_stept0_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18//wind10m_countryfile_stept0_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18/wind_popDensValues_stept0.xml outDir created popfile: LandScan popCellSize= 0.00833333333333 projection= GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]] >> 1. resample the tif file to the resolution and proj of pop density 0.00833333333333 deg >> 2. read the charactristics of the input file >> 3a. extract a piece of pop. file corresponding to the required bounding box split,rastershift False -1 99.3650121065 4.96665456014 120.02334544 22.2749878935 HWRF /usr/bin/gdalwarp -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18//wind10m_popfile_stept0_clipped.tif" >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density gdalwarp -r near -tr 0.00833333333333 0.00833333333333 -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18//wind10m_countryfile_stept0_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.0083333333334 cellsize 0.00833333333333 >> 5. count the popolation in each cell and assign to the class and write to output >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18/wind_popDensValues_stept0.xml >> 7. remove files done t0 completed Creating output file that is 2479P x 2077L. Processing input file /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18//wind10m_res_all.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/lspop20141.tif. Using internal nodata values (e.g. -2.14748e+09) for image /mnt/output/GDACS/DATA/lspop20141.tif. Copying nodata values from source /mnt/output/GDACS/DATA/lspop20141.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18//wind10m_popfile_all_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/countries.tif. Using internal nodata values (e.g. -32768) for image /mnt/output/GDACS/DATA/countries.tif. Copying nodata values from source /mnt/output/GDACS/DATA/countries.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18//wind10m_countryfile_all_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18// PopFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18//wind10m_popfile_all_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18//wind10m_countryfile_all_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18/wind_popDensValues_all.xml popfile: LandScan popCellSize= 0.00833333333333 projection= GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]] >> 1. resample the tif file to the resolution and proj of pop density 0.00833333333333 deg >> 2. read the charactristics of the input file >> 3a. extract a piece of pop. file corresponding to the required bounding box split,rastershift False -1 99.3650121065 4.96665456014 120.02334544 22.2749878935 HWRF /usr/bin/gdalwarp -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18//wind10m_popfile_all_clipped.tif" >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density gdalwarp -r near -tr 0.00833333333333 0.00833333333333 -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18//wind10m_countryfile_all_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.0083333333334 cellsize 0.00833333333333 >> 5. count the popolation in each cell and assign to the class and write to output >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18/wind_popDensValues_all.xml >> 7. remove files done Creating output file that is 2479P x 2077L. Processing input file /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_res_final.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/lspop20141.tif. Using internal nodata values (e.g. -2.14748e+09) for image /mnt/output/GDACS/DATA/lspop20141.tif. Copying nodata values from source /mnt/output/GDACS/DATA/lspop20141.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_popfile_final_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/countries.tif. Using internal nodata values (e.g. -32768) for image /mnt/output/GDACS/DATA/countries.tif. Copying nodata values from source /mnt/output/GDACS/DATA/countries.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_countryfile_final_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. t0 completed ============================================ 4. Copy files ============================================ copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues_t0.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind.jpg >> !!!! Check Copy !!!! , listfiles_copy , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml , listfiles_all , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml all files copied copycheck: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/20191029.18_copy_completed_wind.txt >> 3.2. Classify final folder /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final/20191029.18_final_completed_wind.txt input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final// PopFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_popfile_final_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_countryfile_final_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final/wind_popDensValues_final.xml outDir created popfile: LandScan popCellSize= 0.00833333333333 projection= GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]] >> 1. resample the tif file to the resolution and proj of pop density 0.00833333333333 deg >> 2. read the charactristics of the input file >> 3a. extract a piece of pop. file corresponding to the required bounding box split,rastershift False -1 99.3650121065 4.96665456014 120.02334544 22.2749878935 HWRF /usr/bin/gdalwarp -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_popfile_final_clipped.tif" >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density gdalwarp -r near -tr 0.00833333333333 0.00833333333333 -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_countryfile_final_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.0083333333334 cellsize 0.00833333333333 >> 5. count the popolation in each cell and assign to the class and write to output >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final/wind_popDensValues_final.xml >> 7. remove files done copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final/wind_popDensValues_final.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/final/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_FINAL.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/final/wind.jpg ******************************************************* * Case Completed * ******************************************************* Case starterd at: 2019-10-30 00:25:11.713051 UTC Case completed at: 2019-10-30 00:26:32.754195 UTC **********E N D O F J O B*********************** Creating output file that is 413P x 346L. Processing input file /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif. Using internal nodata values (e.g. -3.40282e+38) for image /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 413P x 346L. Processing input file /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_72h.tif. Using internal nodata values (e.g. -3.40282e+38) for image /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_72h.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. 0 .. 10 .. 20 .. 30 .. 40 .. 50 .. 60 .. 70 .. 80 .. 90 .. 100 - Done ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-30 08:05:06.094168 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-30 00:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 2 it: 0 ndt: 2 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 2 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist GDACS/1000624/1_HWRF 1000624 ndt: 2 it: 0 ndt: 2 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_72h.tif input file exist it: 1 ndt: 2 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_72h.tif input file exist home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret 0 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-10-29 18:00:00 currdate 2019-10-30 00:00:00 ndt: 6 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-10-29 18:00:00', '2019-10-30 00:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files lon0,lon1,lat0,lat1 99.34 120.0 5.0 22.3 processing curr bull... 2019-10-30 00:00:00 99.34 120.0 99.34 120.0 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/ /mnt/input/grib_HWRF/2019/1000624/ input file /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/20191030.00_wind_stept0.tif -r bilinear -te 99.34 5.0 120.0 22.3 -dstnodata 0 Resample tif file /usr/bin/gdalwarp -r bilinear -te 99.34 5.0 120.0 22.3 -dstnodata 0 "/mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/20191030.00_wind_stept0.tif" outFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/20191030.00_wind_stept0.tif dd: 0.0500242130751 >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103000/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103000/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103000/all_inpData.txt 99.3650121065 119.974987893 5.01663438257 22.2749878935 lon 413 lat 346 read TC track /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103000/all_inpData.txt currDate 2019-10-30 00:00:00 f 2 nt-1 5 nt 6 currDate1 2019-10-30 00:00:00 2019-10-30 06:00:00 sdate 2019-10-30 06:00:00 2019-10-30 06:00:00 2 CAT ['TS', 'TS', 'TS'] savemap /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/20191030.00_wind_stept0.jpg 99.34 120.0 99.34 120.0 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/ /mnt/input/grib_HWRF/2019/1000624/ input file /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_72h.tif /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/20191030.00_wind.tif -r bilinear -te 99.34 5.0 120.0 22.3 -dstnodata 0 Resample tif file /usr/bin/gdalwarp -r bilinear -te 99.34 5.0 120.0 22.3 -dstnodata 0 "/mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_72h.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/20191030.00_wind.tif" outFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/20191030.00_wind.tif dd: 0.0500242130751 >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103000/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103000/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103000/all_inpData.txt 99.3650121065 119.974987893 5.01663438257 22.2749878935 lon 413 lat 346 read TC track /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103000/all_inpData.txt currDate 2019-10-30 00:00:00 f 2 nt-1 71 nt 72 currDate1 2019-10-30 00:00:00 2019-11-02 00:00:00 sdate 2019-11-02 00:00:00 2019-11-02 00:00:00 13 CAT ['TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD'] savemap /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/20191030.00_wind.jpg processing all past bull only if Past=True... True itdate, istime 2019-10-29 18:00:00 20191029.18 meteo-processing past forecast already completed >> 2.2. Process final data files compfile /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/20191030.00_Final_completed_wind.txt FINAL remove maxtif alldate: DatetimeIndex(['2019-10-29 18:00:00', '2019-10-30 00:00:00'], dtype='datetime64[ns]', freq='6H') alldate.shape 2 date: 2019-10-30 00:00:00 python /usr/bin/gdal_calc.py -A /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind_stept0.tif -B /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/20191030.00_wind.tif --outfile=/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif --calc="maximum(A,B)" max file created /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif stormname 1000624/23w stormname 1000624/23w TCbullsource 23w >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103000/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103000/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103000/all_inpData.txt 99.3900242131 120.0 4.99162227603 22.2499757869 lon 413 lat 346 read TC track /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103000/all_inpData.txt currDate 2019-10-30 00:00:00 f 2 nt-1 71 nt 72 currDate1 2019-10-30 00:00:00 2019-11-02 00:00:00 sdate 2019-11-02 00:00:00 2019-11-02 00:00:00 13 CAT ['TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD'] savemap /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_FINAL.jpg ret: 0 ============================================ 3. Classify meteo + GDACS index score ============================================ >> 3.1. Classify curr + past forecast Creating output file that is 2479P x 2077L. Processing input file /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/20191030.00_wind_stept0.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/20191030.00_wind_stept0.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/20191030.00_wind_stept0.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00//wind10m_res_stept0.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/lspop20141.tif. Using internal nodata values (e.g. -2.14748e+09) for image /mnt/output/GDACS/DATA/lspop20141.tif. Copying nodata values from source /mnt/output/GDACS/DATA/lspop20141.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00//wind10m_popfile_stept0_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/countries.tif. Using internal nodata values (e.g. -32768) for image /mnt/output/GDACS/DATA/countries.tif. Copying nodata values from source /mnt/output/GDACS/DATA/countries.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00//wind10m_countryfile_stept0_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/20191030.00_wind_stept0.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00// PopFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00//wind10m_popfile_stept0_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00//wind10m_countryfile_stept0_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00/wind_popDensValues_stept0.xml outDir created popfile: LandScan popCellSize= 0.00833333333333 projection= GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]] >> 1. resample the tif file to the resolution and proj of pop density 0.00833333333333 deg >> 2. read the charactristics of the input file >> 3a. extract a piece of pop. file corresponding to the required bounding box split,rastershift False -1 99.3650121065 4.96665456014 120.02334544 22.2749878935 HWRF /usr/bin/gdalwarp -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00//wind10m_popfile_stept0_clipped.tif" >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density gdalwarp -r near -tr 0.00833333333333 0.00833333333333 -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00//wind10m_countryfile_stept0_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.0083333333334 cellsize 0.00833333333333 >> 5. count the popolation in each cell and assign to the class and write to output >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00/wind_popDensValues_stept0.xml >> 7. remove files done t0 completed Creating output file that is 2479P x 2077L. Processing input file /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/20191030.00_wind.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/20191030.00_wind.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/20191030.00_wind.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00//wind10m_res_all.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/lspop20141.tif. Using internal nodata values (e.g. -2.14748e+09) for image /mnt/output/GDACS/DATA/lspop20141.tif. Copying nodata values from source /mnt/output/GDACS/DATA/lspop20141.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00//wind10m_popfile_all_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/countries.tif. Using internal nodata values (e.g. -32768) for image /mnt/output/GDACS/DATA/countries.tif. Copying nodata values from source /mnt/output/GDACS/DATA/countries.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00//wind10m_countryfile_all_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/20191030.00_wind.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00// PopFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00//wind10m_popfile_all_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00//wind10m_countryfile_all_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00/wind_popDensValues_all.xml popfile: LandScan popCellSize= 0.00833333333333 projection= GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]] >> 1. resample the tif file to the resolution and proj of pop density 0.00833333333333 deg >> 2. read the charactristics of the input file >> 3a. extract a piece of pop. file corresponding to the required bounding box split,rastershift False -1 99.3650121065 4.96665456014 120.02334544 22.2749878935 HWRF /usr/bin/gdalwarp -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00//wind10m_popfile_all_clipped.tif" >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density gdalwarp -r near -tr 0.00833333333333 0.00833333333333 -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00//wind10m_countryfile_all_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.0083333333334 cellsize 0.00833333333333 >> 5. count the popolation in each cell and assign to the class and write to output >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00/wind_popDensValues_all.xml >> 7. remove files done Creating output file that is 2479P x 2077L. Processing input file /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif. Using internal nodata values (e.g. 3.40282e+38) for image /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_res_final.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/lspop20141.tif. Using internal nodata values (e.g. -2.14748e+09) for image /mnt/output/GDACS/DATA/lspop20141.tif. Copying nodata values from source /mnt/output/GDACS/DATA/lspop20141.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_popfile_final_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/countries.tif. Using internal nodata values (e.g. -32768) for image /mnt/output/GDACS/DATA/countries.tif. Copying nodata values from source /mnt/output/GDACS/DATA/countries.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_countryfile_final_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. t0 completed ============================================ 4. Copy files ============================================ copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues_t0.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind.jpg >> !!!! Check Copy !!!! , listfiles_copy , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml , listfiles_all , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml all files copied copycheck: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/20191029.18_copy_completed_wind.txt copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues_t0.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/20191030.00_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind.jpg >> !!!! Check Copy !!!! , listfiles_copy , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues.xml , listfiles_all , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues.xml all files copied copycheck: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/20191030.00_copy_completed_wind.txt >> 3.2. Classify final folder /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final/20191030.00_final_completed_wind.txt input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final// PopFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_popfile_final_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_countryfile_final_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final/wind_popDensValues_final.xml xml file exists...REMOVE popfile: LandScan popCellSize= 0.00833333333333 projection= GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]] >> 1. resample the tif file to the resolution and proj of pop density 0.00833333333333 deg >> 2. read the charactristics of the input file >> 3a. extract a piece of pop. file corresponding to the required bounding box split,rastershift False -1 99.3650121065 4.96665456014 120.02334544 22.2749878935 HWRF /usr/bin/gdalwarp -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_popfile_final_clipped.tif" >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density gdalwarp -r near -tr 0.00833333333333 0.00833333333333 -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_countryfile_final_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.0083333333334 cellsize 0.00833333333333 >> 5. count the popolation in each cell and assign to the class and write to output >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final/wind_popDensValues_final.xml >> 7. remove files done copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final/wind_popDensValues_final.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/final/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_FINAL.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/final/wind.jpg ******************************************************* * Case Completed * ******************************************************* Case starterd at: 2019-10-30 08:05:06.094168 UTC Case completed at: 2019-10-30 08:06:35.615436 UTC **********E N D O F J O B*********************** ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-30 08:15:07.130805 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-30 00:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 2 it: 0 ndt: 2 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 2 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist GDACS/1000624/1_HWRF 1000624 ndt: 2 it: 0 ndt: 2 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_72h.tif input file exist it: 1 ndt: 2 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_72h.tif input file exist home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret 0 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-10-29 18:00:00 currdate 2019-10-30 00:00:00 ndt: 6 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-10-29 18:00:00', '2019-10-30 00:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files lon0,lon1,lat0,lat1 99.34 120.0 5.0 22.3 processing curr bull... 2019-10-30 00:00:00 ********** no new data check copy copy ok copy ok no new data ******************************************************* * Case Completed * ******************************************************* Case starterd at: 2019-10-30 08:15:07.130805 UTC Case completed at: 2019-10-30 08:15:07.231226 UTC **********E N D O F J O B*********************** ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-30 08:25:11.047211 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-30 00:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 2 it: 0 ndt: 2 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 2 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist GDACS/1000624/1_HWRF 1000624 ndt: 2 it: 0 ndt: 2 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_72h.tif input file exist it: 1 ndt: 2 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_72h.tif input file exist home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret 0 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-10-29 18:00:00 currdate 2019-10-30 00:00:00 ndt: 6 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-10-29 18:00:00', '2019-10-30 00:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files lon0,lon1,lat0,lat1 99.34 120.0 5.0 22.3 processing curr bull... 2019-10-30 00:00:00 ********** no new data check copy copy ok copy ok no new data ******************************************************* * Case Completed * ******************************************************* Case starterd at: 2019-10-30 08:25:11.047211 UTC Case completed at: 2019-10-30 08:25:11.164265 UTC **********E N D O F J O B*********************** ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-30 12:35:08.253778 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-30 06:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 3 it: 0 ndt: 3 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 3 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 3 idate: 2019-10-30 06:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== Creating output file that is 413P x 346L. Processing input file /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif. Using internal nodata values (e.g. -3.40282e+38) for image /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 413P x 346L. Processing input file /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_72h.tif. Using internal nodata values (e.g. -3.40282e+38) for image /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_72h.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. 0 .. 10 .. 20 .. 30 .. 40 .. 50 .. 60 .. 70 .. 80 .. 90 .. 100 - Done 0 .. 10 .. 20 .. 30 .. 40 .. 50 .. 60 .. 70 .. 80 .. 90 .. 100 - Done ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-30 12:45:10.473925 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-30 06:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 3 it: 0 ndt: 3 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 3 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 3 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif input file exist GDACS/1000624/1_HWRF 1000624 ndt: 3 it: 0 ndt: 3 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_72h.tif input file exist it: 1 ndt: 3 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_72h.tif input file exist it: 2 ndt: 3 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_72h.tif input file exist home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret 0 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-10-29 18:00:00 currdate 2019-10-30 06:00:00 ndt: 12 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-10-29 18:00:00', '2019-10-30 00:00:00', '2019-10-30 06:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files lon0,lon1,lat0,lat1 99.34 120.0 5.0 22.3 processing curr bull... 2019-10-30 06:00:00 99.34 120.0 99.34 120.0 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.06/ /mnt/input/grib_HWRF/2019/1000624/ input file /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.06/20191030.06_wind_stept0.tif -r bilinear -te 99.34 5.0 120.0 22.3 -dstnodata 0 Resample tif file /usr/bin/gdalwarp -r bilinear -te 99.34 5.0 120.0 22.3 -dstnodata 0 "/mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.06/20191030.06_wind_stept0.tif" outFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.06/20191030.06_wind_stept0.tif dd: 0.0500242130751 >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103006/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103006/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103006/all_inpData.txt 99.3650121065 119.974987893 5.01663438257 22.2749878935 lon 413 lat 346 read TC track /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103006/all_inpData.txt currDate 2019-10-30 06:00:00 f 3 nt-1 5 nt 6 currDate1 2019-10-30 06:00:00 2019-10-30 12:00:00 sdate 2019-10-30 12:00:00 2019-10-30 12:00:00 3 CAT ['TS', 'TS', 'TS', 'TS'] savemap /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.06/20191030.06_wind_stept0.jpg 99.34 120.0 99.34 120.0 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.06/ /mnt/input/grib_HWRF/2019/1000624/ input file /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_72h.tif /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.06/20191030.06_wind.tif -r bilinear -te 99.34 5.0 120.0 22.3 -dstnodata 0 Resample tif file /usr/bin/gdalwarp -r bilinear -te 99.34 5.0 120.0 22.3 -dstnodata 0 "/mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_72h.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.06/20191030.06_wind.tif" outFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.06/20191030.06_wind.tif dd: 0.0500242130751 >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103006/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103006/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103006/all_inpData.txt 99.3650121065 119.974987893 5.01663438257 22.2749878935 lon 413 lat 346 read TC track /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103006/all_inpData.txt currDate 2019-10-30 06:00:00 f 3 nt-1 71 nt 72 currDate1 2019-10-30 06:00:00 2019-11-02 06:00:00 sdate 2019-11-02 06:00:00 2019-11-02 06:00:00 14 CAT ['TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD'] savemap /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.06/20191030.06_wind.jpg processing all past bull only if Past=True... True itdate, istime 2019-10-29 18:00:00 20191029.18 meteo-processing past forecast already completed itdate, istime 2019-10-30 00:00:00 20191030.00 meteo-processing past forecast already completed >> 2.2. Process final data files compfile /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/20191030.06_Final_completed_wind.txt FINAL remove maxtif alldate: DatetimeIndex(['2019-10-29 18:00:00', '2019-10-30 00:00:00', '2019-10-30 06:00:00'], dtype='datetime64[ns]', freq='6H') alldate.shape 3 date: 2019-10-30 00:00:00 python /usr/bin/gdal_calc.py -A /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind_stept0.tif -B /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/20191030.00_wind_stept0.tif --outfile=/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif --calc="maximum(A,B)" date: 2019-10-30 06:00:00 python /usr/bin/gdal_calc.py -A /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif -B /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.06/20191030.06_wind.tif --outfile=/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif --calc="maximum(A,B)" max file created /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif stormname 1000624/23w stormname 1000624/23w TCbullsource 23w >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103006/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103006/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103006/all_inpData.txt 99.3900242131 120.0 4.99162227603 22.2499757869 lon 413 lat 346 read TC track /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103006/all_inpData.txt currDate 2019-10-30 06:00:00 f 3 nt-1 71 nt 72 currDate1 2019-10-30 06:00:00 2019-11-02 06:00:00 sdate 2019-11-02 06:00:00 2019-11-02 06:00:00 14 CAT ['TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD'] savemap /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_FINAL.jpg ret: 0 ============================================ 3. Classify meteo + GDACS index score ============================================ >> 3.1. Classify curr + past forecast Creating output file that is 2479P x 2077L. Processing input file /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.06/20191030.06_wind_stept0.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.06/20191030.06_wind_stept0.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.06/20191030.06_wind_stept0.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06//wind10m_res_stept0.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/lspop20141.tif. Using internal nodata values (e.g. -2.14748e+09) for image /mnt/output/GDACS/DATA/lspop20141.tif. Copying nodata values from source /mnt/output/GDACS/DATA/lspop20141.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06//wind10m_popfile_stept0_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/countries.tif. Using internal nodata values (e.g. -32768) for image /mnt/output/GDACS/DATA/countries.tif. Copying nodata values from source /mnt/output/GDACS/DATA/countries.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06//wind10m_countryfile_stept0_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.06/20191030.06_wind_stept0.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06// PopFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06//wind10m_popfile_stept0_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06//wind10m_countryfile_stept0_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06/wind_popDensValues_stept0.xml outDir created popfile: LandScan popCellSize= 0.00833333333333 projection= GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]] >> 1. resample the tif file to the resolution and proj of pop density 0.00833333333333 deg >> 2. read the charactristics of the input file >> 3a. extract a piece of pop. file corresponding to the required bounding box split,rastershift False -1 99.3650121065 4.96665456014 120.02334544 22.2749878935 HWRF /usr/bin/gdalwarp -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06//wind10m_popfile_stept0_clipped.tif" >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density gdalwarp -r near -tr 0.00833333333333 0.00833333333333 -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06//wind10m_countryfile_stept0_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.0083333333334 cellsize 0.00833333333333 >> 5. count the popolation in each cell and assign to the class and write to output >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06/wind_popDensValues_stept0.xml >> 7. remove files done t0 completed Creating output file that is 2479P x 2077L. Processing input file /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.06/20191030.06_wind.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.06/20191030.06_wind.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.06/20191030.06_wind.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06//wind10m_res_all.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/lspop20141.tif. Using internal nodata values (e.g. -2.14748e+09) for image /mnt/output/GDACS/DATA/lspop20141.tif. Copying nodata values from source /mnt/output/GDACS/DATA/lspop20141.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06//wind10m_popfile_all_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/countries.tif. Using internal nodata values (e.g. -32768) for image /mnt/output/GDACS/DATA/countries.tif. Copying nodata values from source /mnt/output/GDACS/DATA/countries.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06//wind10m_countryfile_all_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.06/20191030.06_wind.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06// PopFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06//wind10m_popfile_all_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06//wind10m_countryfile_all_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06/wind_popDensValues_all.xml popfile: LandScan popCellSize= 0.00833333333333 projection= GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]] >> 1. resample the tif file to the resolution and proj of pop density 0.00833333333333 deg >> 2. read the charactristics of the input file >> 3a. extract a piece of pop. file corresponding to the required bounding box split,rastershift False -1 99.3650121065 4.96665456014 120.02334544 22.2749878935 HWRF /usr/bin/gdalwarp -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06//wind10m_popfile_all_clipped.tif" >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density gdalwarp -r near -tr 0.00833333333333 0.00833333333333 -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06//wind10m_countryfile_all_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.0083333333334 cellsize 0.00833333333333 >> 5. count the popolation in each cell and assign to the class and write to output >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06/wind_popDensValues_all.xml >> 7. remove files done Creating output file that is 2479P x 2077L. Processing input file /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif. Using internal nodata values (e.g. 3.40282e+38) for image /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_res_final.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/lspop20141.tif. Using internal nodata values (e.g. -2.14748e+09) for image /mnt/output/GDACS/DATA/lspop20141.tif. Copying nodata values from source /mnt/output/GDACS/DATA/lspop20141.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_popfile_final_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/countries.tif. Using internal nodata values (e.g. -32768) for image /mnt/output/GDACS/DATA/countries.tif. Copying nodata values from source /mnt/output/GDACS/DATA/countries.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_countryfile_final_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. t0 completed ============================================ 4. Copy files ============================================ copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues_t0.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind.jpg >> !!!! Check Copy !!!! , listfiles_copy , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml , listfiles_all , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml all files copied copycheck: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/20191029.18_copy_completed_wind.txt copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues_t0.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/20191030.00_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind.jpg >> !!!! Check Copy !!!! , listfiles_copy , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues.xml , listfiles_all , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues.xml all files copied copycheck: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/20191030.00_copy_completed_wind.txt copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues_t0.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.06/20191030.06_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind.jpg >> !!!! Check Copy !!!! , listfiles_copy , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues.xml , listfiles_all , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues.xml all files copied copycheck: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/20191030.06_copy_completed_wind.txt >> 3.2. Classify final folder /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final/20191030.06_final_completed_wind.txt input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final// PopFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_popfile_final_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_countryfile_final_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final/wind_popDensValues_final.xml xml file exists...REMOVE popfile: LandScan popCellSize= 0.00833333333333 projection= GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]] >> 1. resample the tif file to the resolution and proj of pop density 0.00833333333333 deg >> 2. read the charactristics of the input file >> 3a. extract a piece of pop. file corresponding to the required bounding box split,rastershift False -1 99.3650121065 4.96665456014 120.02334544 22.2749878935 HWRF /usr/bin/gdalwarp -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_popfile_final_clipped.tif" >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density gdalwarp -r near -tr 0.00833333333333 0.00833333333333 -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_countryfile_final_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.0083333333334 cellsize 0.00833333333333 >> 5. count the popolation in each cell and assign to the class and write to output >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final/wind_popDensValues_final.xml >> 7. remove files done copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final/wind_popDensValues_final.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/final/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_FINAL.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/final/wind.jpg ******************************************************* * Case Completed * ******************************************************* Case starterd at: 2019-10-30 12:45:10.473925 UTC Case completed at: 2019-10-30 12:46:34.698904 UTC **********E N D O F J O B*********************** Creating output file that is 413P x 346L. Processing input file /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif. Using internal nodata values (e.g. -3.40282e+38) for image /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 413P x 346L. Processing input file /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_72h.tif. Using internal nodata values (e.g. -3.40282e+38) for image /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_72h.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. 0 .. 10 .. 20 .. 30 .. 40 .. 50 .. 60 .. 70 .. 80 .. 90 .. 100 - Done 0 .. 10 .. 20 .. 30 .. 40 .. 50 .. 60 .. 70 .. 80 .. 90 .. 100 - Done 0 .. 10 .. 20 .. 30 .. 40 .. 50 .. 60 .. 70 .. 80 .. 90 .. 100 - Done ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-30 21:55:14.649854 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-30 12:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 4 it: 0 ndt: 4 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 4 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 4 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif input file exist it: 3 ndt: 4 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif input file exist GDACS/1000624/1_HWRF 1000624 ndt: 4 it: 0 ndt: 4 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_72h.tif input file exist it: 1 ndt: 4 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_72h.tif input file exist it: 2 ndt: 4 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_72h.tif input file exist it: 3 ndt: 4 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_72h.tif input file exist home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret 0 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-10-29 18:00:00 currdate 2019-10-30 12:00:00 ndt: 18 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-10-29 18:00:00', '2019-10-30 00:00:00', '2019-10-30 06:00:00', '2019-10-30 12:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files lon0,lon1,lat0,lat1 99.34 120.0 5.0 22.3 processing curr bull... 2019-10-30 12:00:00 99.34 120.0 99.34 120.0 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.12/ /mnt/input/grib_HWRF/2019/1000624/ input file /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.12/20191030.12_wind_stept0.tif -r bilinear -te 99.34 5.0 120.0 22.3 -dstnodata 0 Resample tif file /usr/bin/gdalwarp -r bilinear -te 99.34 5.0 120.0 22.3 -dstnodata 0 "/mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.12/20191030.12_wind_stept0.tif" outFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.12/20191030.12_wind_stept0.tif dd: 0.0500242130751 >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103012/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103012/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103012/all_inpData.txt 99.3650121065 119.974987893 5.01663438257 22.2749878935 lon 413 lat 346 read TC track /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103012/all_inpData.txt currDate 2019-10-30 12:00:00 f 4 nt-1 5 nt 6 currDate1 2019-10-30 12:00:00 2019-10-30 18:00:00 sdate 2019-10-30 18:00:00 2019-10-30 18:00:00 4 CAT ['TS', 'TS', 'TS', 'TS', 'TS'] savemap /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.12/20191030.12_wind_stept0.jpg 99.34 120.0 99.34 120.0 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.12/ /mnt/input/grib_HWRF/2019/1000624/ input file /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_72h.tif /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.12/20191030.12_wind.tif -r bilinear -te 99.34 5.0 120.0 22.3 -dstnodata 0 Resample tif file /usr/bin/gdalwarp -r bilinear -te 99.34 5.0 120.0 22.3 -dstnodata 0 "/mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_72h.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.12/20191030.12_wind.tif" outFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.12/20191030.12_wind.tif dd: 0.0500242130751 >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103012/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103012/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103012/all_inpData.txt 99.3650121065 119.974987893 5.01663438257 22.2749878935 lon 413 lat 346 read TC track /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103012/all_inpData.txt currDate 2019-10-30 12:00:00 f 4 nt-1 71 nt 72 currDate1 2019-10-30 12:00:00 2019-11-02 12:00:00 sdate 2019-11-02 12:00:00 2019-11-02 12:00:00 15 CAT ['TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD'] savemap /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.12/20191030.12_wind.jpg processing all past bull only if Past=True... True itdate, istime 2019-10-29 18:00:00 20191029.18 meteo-processing past forecast already completed itdate, istime 2019-10-30 00:00:00 20191030.00 meteo-processing past forecast already completed itdate, istime 2019-10-30 06:00:00 20191030.06 meteo-processing past forecast already completed >> 2.2. Process final data files compfile /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/20191030.12_Final_completed_wind.txt FINAL remove maxtif alldate: DatetimeIndex(['2019-10-29 18:00:00', '2019-10-30 00:00:00', '2019-10-30 06:00:00', '2019-10-30 12:00:00'], dtype='datetime64[ns]', freq='6H') alldate.shape 4 date: 2019-10-30 00:00:00 python /usr/bin/gdal_calc.py -A /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind_stept0.tif -B /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/20191030.00_wind_stept0.tif --outfile=/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif --calc="maximum(A,B)" date: 2019-10-30 06:00:00 python /usr/bin/gdal_calc.py -A /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif -B /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.06/20191030.06_wind_stept0.tif --outfile=/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif --calc="maximum(A,B)" date: 2019-10-30 12:00:00 python /usr/bin/gdal_calc.py -A /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif -B /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.12/20191030.12_wind.tif --outfile=/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif --calc="maximum(A,B)" max file created /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif stormname 1000624/23w stormname 1000624/23w TCbullsource 23w >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103012/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103012/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103012/all_inpData.txt 99.3900242131 120.0 4.99162227603 22.2499757869 lon 413 lat 346 read TC track /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103012/all_inpData.txt currDate 2019-10-30 12:00:00 f 4 nt-1 71 nt 72 currDate1 2019-10-30 12:00:00 2019-11-02 12:00:00 sdate 2019-11-02 12:00:00 2019-11-02 12:00:00 15 CAT ['TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD'] savemap /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_FINAL.jpg ret: 0 ============================================ 3. Classify meteo + GDACS index score ============================================ >> 3.1. Classify curr + past forecast Creating output file that is 2479P x 2077L. Processing input file /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.12/20191030.12_wind_stept0.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.12/20191030.12_wind_stept0.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.12/20191030.12_wind_stept0.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.12//wind10m_res_stept0.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/lspop20141.tif. Using internal nodata values (e.g. -2.14748e+09) for image /mnt/output/GDACS/DATA/lspop20141.tif. Copying nodata values from source /mnt/output/GDACS/DATA/lspop20141.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.12//wind10m_popfile_stept0_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/countries.tif. Using internal nodata values (e.g. -32768) for image /mnt/output/GDACS/DATA/countries.tif. Copying nodata values from source /mnt/output/GDACS/DATA/countries.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.12//wind10m_countryfile_stept0_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.12/20191030.12_wind_stept0.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.12// PopFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.12//wind10m_popfile_stept0_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.12//wind10m_countryfile_stept0_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.12/wind_popDensValues_stept0.xml popfile: LandScan popCellSize= 0.00833333333333 projection= GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]] >> 1. resample the tif file to the resolution and proj of pop density 0.00833333333333 deg >> 2. read the charactristics of the input file >> 3a. extract a piece of pop. file corresponding to the required bounding box split,rastershift False -1 99.3650121065 4.96665456014 120.02334544 22.2749878935 HWRF /usr/bin/gdalwarp -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.12//wind10m_popfile_stept0_clipped.tif" >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density gdalwarp -r near -tr 0.00833333333333 0.00833333333333 -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.12//wind10m_countryfile_stept0_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.0083333333334 cellsize 0.00833333333333 >> 5. count the popolation in each cell and assign to the class and write to output >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.12/wind_popDensValues_stept0.xml >> 7. remove files done t0 completed Creating output file that is 2479P x 2077L. Processing input file /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.12/20191030.12_wind.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.12/20191030.12_wind.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.12/20191030.12_wind.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.12//wind10m_res_all.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/lspop20141.tif. Using internal nodata values (e.g. -2.14748e+09) for image /mnt/output/GDACS/DATA/lspop20141.tif. Copying nodata values from source /mnt/output/GDACS/DATA/lspop20141.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.12//wind10m_popfile_all_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/countries.tif. Using internal nodata values (e.g. -32768) for image /mnt/output/GDACS/DATA/countries.tif. Copying nodata values from source /mnt/output/GDACS/DATA/countries.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.12//wind10m_countryfile_all_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.12/20191030.12_wind.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.12// PopFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.12//wind10m_popfile_all_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.12//wind10m_countryfile_all_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.12/wind_popDensValues_all.xml popfile: LandScan popCellSize= 0.00833333333333 projection= GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]] >> 1. resample the tif file to the resolution and proj of pop density 0.00833333333333 deg >> 2. read the charactristics of the input file >> 3a. extract a piece of pop. file corresponding to the required bounding box split,rastershift False -1 99.3650121065 4.96665456014 120.02334544 22.2749878935 HWRF /usr/bin/gdalwarp -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.12//wind10m_popfile_all_clipped.tif" >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density gdalwarp -r near -tr 0.00833333333333 0.00833333333333 -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.12//wind10m_countryfile_all_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.0083333333334 cellsize 0.00833333333333 >> 5. count the popolation in each cell and assign to the class and write to output >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.12/wind_popDensValues_all.xml >> 7. remove files done Creating output file that is 2479P x 2077L. Processing input file /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif. Using internal nodata values (e.g. 3.40282e+38) for image /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_res_final.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/lspop20141.tif. Using internal nodata values (e.g. -2.14748e+09) for image /mnt/output/GDACS/DATA/lspop20141.tif. Copying nodata values from source /mnt/output/GDACS/DATA/lspop20141.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_popfile_final_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/countries.tif. Using internal nodata values (e.g. -32768) for image /mnt/output/GDACS/DATA/countries.tif. Copying nodata values from source /mnt/output/GDACS/DATA/countries.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_countryfile_final_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. t0 completed ============================================ 4. Copy files ============================================ check t0 xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues_t0.xml check all xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind.jpg >> !!!! Check Copy !!!! , listfiles_copy , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml , listfiles_all , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml all files copied copycheck: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/20191029.18_copy_completed_wind.txt check t0 xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues_t0.xml check all xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/20191030.00_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind.jpg >> !!!! Check Copy !!!! , listfiles_copy , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues.xml , listfiles_all , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues.xml all files copied copycheck: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/20191030.00_copy_completed_wind.txt check t0 xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues_t0.xml check all xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.06/20191030.06_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind.jpg >> !!!! Check Copy !!!! , listfiles_copy , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues.xml , listfiles_all , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues.xml all files copied copycheck: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/20191030.06_copy_completed_wind.txt check t0 xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.12/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind_popDensValues_t0.xml check all xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.12/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.12/20191030.12_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind.jpg >> !!!! Check Copy !!!! , listfiles_copy , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind_popDensValues.xml , listfiles_all , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind_popDensValues.xml all files copied copycheck: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/20191030.12_copy_completed_wind.txt >> 3.2. Classify final folder /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final/20191030.12_final_completed_wind.txt input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final// PopFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_popfile_final_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_countryfile_final_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final/wind_popDensValues_final.xml xml file exists...REMOVE popfile: LandScan popCellSize= 0.00833333333333 projection= GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]] >> 1. resample the tif file to the resolution and proj of pop density 0.00833333333333 deg >> 2. read the charactristics of the input file >> 3a. extract a piece of pop. file corresponding to the required bounding box split,rastershift False -1 99.3650121065 4.96665456014 120.02334544 22.2749878935 HWRF /usr/bin/gdalwarp -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_popfile_final_clipped.tif" >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density gdalwarp -r near -tr 0.00833333333333 0.00833333333333 -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_countryfile_final_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.0083333333334 cellsize 0.00833333333333 >> 5. count the popolation in each cell and assign to the class and write to output >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final/wind_popDensValues_final.xml >> 7. remove files done copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final/wind_popDensValues_final.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/final/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_FINAL.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/final/wind.jpg ******************************************************* * Case Completed * ******************************************************* Case starterd at: 2019-10-30 21:55:14.649854 UTC Case completed at: 2019-10-30 21:56:45.283726 UTC **********E N D O F J O B*********************** ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-30 22:05:06.594759 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-30 12:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 4 it: 0 ndt: 4 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 4 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 4 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif input file exist it: 3 ndt: 4 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif input file exist GDACS/1000624/1_HWRF 1000624 ndt: 4 it: 0 ndt: 4 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_72h.tif input file exist it: 1 ndt: 4 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_72h.tif input file exist it: 2 ndt: 4 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_72h.tif input file exist it: 3 ndt: 4 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_72h.tif input file exist home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret 0 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-10-29 18:00:00 currdate 2019-10-30 12:00:00 ndt: 18 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-10-29 18:00:00', '2019-10-30 00:00:00', '2019-10-30 06:00:00', '2019-10-30 12:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files lon0,lon1,lat0,lat1 99.34 120.0 5.0 22.3 processing curr bull... 2019-10-30 12:00:00 ********** no new data check copy copy ok copy ok copy ok copy ok no new data ******************************************************* * Case Completed * ******************************************************* Case starterd at: 2019-10-30 22:05:06.594759 UTC Case completed at: 2019-10-30 22:05:06.957813 UTC **********E N D O F J O B*********************** ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-30 22:15:11.342880 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-30 12:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 4 it: 0 ndt: 4 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 4 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 4 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif input file exist it: 3 ndt: 4 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif input file exist GDACS/1000624/1_HWRF 1000624 ndt: 4 it: 0 ndt: 4 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_72h.tif input file exist it: 1 ndt: 4 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_72h.tif input file exist it: 2 ndt: 4 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_72h.tif input file exist it: 3 ndt: 4 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_72h.tif input file exist home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret 0 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-10-29 18:00:00 currdate 2019-10-30 12:00:00 ndt: 18 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-10-29 18:00:00', '2019-10-30 00:00:00', '2019-10-30 06:00:00', '2019-10-30 12:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files lon0,lon1,lat0,lat1 99.34 120.0 5.0 22.3 processing curr bull... 2019-10-30 12:00:00 ********** no new data check copy copy ok copy ok copy ok copy ok no new data ******************************************************* * Case Completed * ******************************************************* Case starterd at: 2019-10-30 22:15:11.342880 UTC Case completed at: 2019-10-30 22:15:11.700300 UTC **********E N D O F J O B*********************** ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-31 00:35:06.366082 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-30 18:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 5 it: 0 ndt: 5 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 5 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 5 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif input file exist it: 3 ndt: 5 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif input file exist it: 4 ndt: 5 idate: 2019-10-30 18:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== Creating output file that is 413P x 346L. Processing input file /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_6h.tif. Using internal nodata values (e.g. -3.40282e+38) for image /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_6h.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 413P x 346L. Processing input file /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_72h.tif. Using internal nodata values (e.g. -3.40282e+38) for image /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_72h.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. 0 .. 10 .. 20 .. 30 .. 40 .. 50 .. 60 .. 70 .. 80 .. 90 .. 100 - Done 0 .. 10 .. 20 .. 30 .. 40 .. 50 .. 60 .. 70 .. 80 .. 90 .. 100 - Done 0 .. 10 .. 20 .. 30 .. 40 .. 50 .. 60 .. 70 .. 80 .. 90 .. 100 - Done 0 .. 10 .. 20 .. 30 .. 40 .. 50 .. 60 .. 70 .. 80 .. 90 .. 100 - Done ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-31 00:45:07.285041 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-30 18:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 5 it: 0 ndt: 5 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 5 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 5 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif input file exist it: 3 ndt: 5 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif input file exist it: 4 ndt: 5 idate: 2019-10-30 18:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_6h.tif input file exist GDACS/1000624/1_HWRF 1000624 ndt: 5 it: 0 ndt: 5 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_72h.tif input file exist it: 1 ndt: 5 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_72h.tif input file exist it: 2 ndt: 5 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_72h.tif input file exist it: 3 ndt: 5 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_72h.tif input file exist it: 4 ndt: 5 idate: 2019-10-30 18:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_72h.tif input file exist home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret 0 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-10-29 18:00:00 currdate 2019-10-30 18:00:00 ndt: 24 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-10-29 18:00:00', '2019-10-30 00:00:00', '2019-10-30 06:00:00', '2019-10-30 12:00:00', '2019-10-30 18:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files lon0,lon1,lat0,lat1 99.34 120.0 5.0 22.3 processing curr bull... 2019-10-30 18:00:00 99.34 120.0 99.34 120.0 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.18/ /mnt/input/grib_HWRF/2019/1000624/ input file /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_6h.tif /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.18/20191030.18_wind_stept0.tif -r bilinear -te 99.34 5.0 120.0 22.3 -dstnodata 0 Resample tif file /usr/bin/gdalwarp -r bilinear -te 99.34 5.0 120.0 22.3 -dstnodata 0 "/mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_6h.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.18/20191030.18_wind_stept0.tif" outFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.18/20191030.18_wind_stept0.tif dd: 0.0500242130751 >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103018/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103018/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103018/all_inpData.txt 99.3650121065 119.974987893 5.01663438257 22.2749878935 lon 413 lat 346 read TC track /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103018/all_inpData.txt currDate 2019-10-30 18:00:00 f 5 nt-1 5 nt 6 currDate1 2019-10-30 18:00:00 2019-10-31 00:00:00 sdate 2019-10-31 00:00:00 2019-10-31 00:00:00 5 CAT ['TS', 'TS', 'TS', 'TS', 'TS', 'TS'] savemap /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.18/20191030.18_wind_stept0.jpg 99.34 120.0 99.34 120.0 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.18/ /mnt/input/grib_HWRF/2019/1000624/ input file /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_72h.tif /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.18/20191030.18_wind.tif -r bilinear -te 99.34 5.0 120.0 22.3 -dstnodata 0 Resample tif file /usr/bin/gdalwarp -r bilinear -te 99.34 5.0 120.0 22.3 -dstnodata 0 "/mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_72h.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.18/20191030.18_wind.tif" outFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.18/20191030.18_wind.tif dd: 0.0500242130751 >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103018/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103018/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103018/all_inpData.txt 99.3650121065 119.974987893 5.01663438257 22.2749878935 lon 413 lat 346 read TC track /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103018/all_inpData.txt currDate 2019-10-30 18:00:00 f 5 nt-1 71 nt 72 currDate1 2019-10-30 18:00:00 2019-11-02 18:00:00 sdate 2019-11-02 18:00:00 2019-11-02 18:00:00 16 CAT ['TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD'] savemap /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.18/20191030.18_wind.jpg processing all past bull only if Past=True... True itdate, istime 2019-10-29 18:00:00 20191029.18 meteo-processing past forecast already completed itdate, istime 2019-10-30 00:00:00 20191030.00 meteo-processing past forecast already completed itdate, istime 2019-10-30 06:00:00 20191030.06 meteo-processing past forecast already completed itdate, istime 2019-10-30 12:00:00 20191030.12 meteo-processing past forecast already completed >> 2.2. Process final data files compfile /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/20191030.18_Final_completed_wind.txt FINAL remove maxtif alldate: DatetimeIndex(['2019-10-29 18:00:00', '2019-10-30 00:00:00', '2019-10-30 06:00:00', '2019-10-30 12:00:00', '2019-10-30 18:00:00'], dtype='datetime64[ns]', freq='6H') alldate.shape 5 date: 2019-10-30 00:00:00 python /usr/bin/gdal_calc.py -A /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind_stept0.tif -B /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/20191030.00_wind_stept0.tif --outfile=/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif --calc="maximum(A,B)" date: 2019-10-30 06:00:00 python /usr/bin/gdal_calc.py -A /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif -B /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.06/20191030.06_wind_stept0.tif --outfile=/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif --calc="maximum(A,B)" date: 2019-10-30 12:00:00 python /usr/bin/gdal_calc.py -A /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif -B /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.12/20191030.12_wind_stept0.tif --outfile=/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif --calc="maximum(A,B)" date: 2019-10-30 18:00:00 python /usr/bin/gdal_calc.py -A /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif -B /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.18/20191030.18_wind.tif --outfile=/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif --calc="maximum(A,B)" max file created /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif stormname 1000624/23w stormname 1000624/23w TCbullsource 23w >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103018/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103018/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103018/all_inpData.txt 99.3900242131 120.0 4.99162227603 22.2499757869 lon 413 lat 346 read TC track /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103018/all_inpData.txt currDate 2019-10-30 18:00:00 f 5 nt-1 71 nt 72 currDate1 2019-10-30 18:00:00 2019-11-02 18:00:00 sdate 2019-11-02 18:00:00 2019-11-02 18:00:00 16 CAT ['TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD'] savemap /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_FINAL.jpg ret: 0 ============================================ 3. Classify meteo + GDACS index score ============================================ >> 3.1. Classify curr + past forecast Creating output file that is 2479P x 2077L. Processing input file /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.18/20191030.18_wind_stept0.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.18/20191030.18_wind_stept0.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.18/20191030.18_wind_stept0.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.18//wind10m_res_stept0.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/lspop20141.tif. Using internal nodata values (e.g. -2.14748e+09) for image /mnt/output/GDACS/DATA/lspop20141.tif. Copying nodata values from source /mnt/output/GDACS/DATA/lspop20141.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.18//wind10m_popfile_stept0_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/countries.tif. Using internal nodata values (e.g. -32768) for image /mnt/output/GDACS/DATA/countries.tif. Copying nodata values from source /mnt/output/GDACS/DATA/countries.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.18//wind10m_countryfile_stept0_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.18/20191030.18_wind_stept0.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.18// PopFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.18//wind10m_popfile_stept0_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.18//wind10m_countryfile_stept0_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.18/wind_popDensValues_stept0.xml popfile: LandScan popCellSize= 0.00833333333333 projection= GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]] >> 1. resample the tif file to the resolution and proj of pop density 0.00833333333333 deg >> 2. read the charactristics of the input file >> 3a. extract a piece of pop. file corresponding to the required bounding box split,rastershift False -1 99.3650121065 4.96665456014 120.02334544 22.2749878935 HWRF /usr/bin/gdalwarp -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.18//wind10m_popfile_stept0_clipped.tif" >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density gdalwarp -r near -tr 0.00833333333333 0.00833333333333 -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.18//wind10m_countryfile_stept0_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.0083333333334 cellsize 0.00833333333333 >> 5. count the popolation in each cell and assign to the class and write to output >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.18/wind_popDensValues_stept0.xml >> 7. remove files done t0 completed Creating output file that is 2479P x 2077L. Processing input file /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.18/20191030.18_wind.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.18/20191030.18_wind.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.18/20191030.18_wind.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.18//wind10m_res_all.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/lspop20141.tif. Using internal nodata values (e.g. -2.14748e+09) for image /mnt/output/GDACS/DATA/lspop20141.tif. Copying nodata values from source /mnt/output/GDACS/DATA/lspop20141.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.18//wind10m_popfile_all_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/countries.tif. Using internal nodata values (e.g. -32768) for image /mnt/output/GDACS/DATA/countries.tif. Copying nodata values from source /mnt/output/GDACS/DATA/countries.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.18//wind10m_countryfile_all_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.18/20191030.18_wind.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.18// PopFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.18//wind10m_popfile_all_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.18//wind10m_countryfile_all_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.18/wind_popDensValues_all.xml popfile: LandScan popCellSize= 0.00833333333333 projection= GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]] >> 1. resample the tif file to the resolution and proj of pop density 0.00833333333333 deg >> 2. read the charactristics of the input file >> 3a. extract a piece of pop. file corresponding to the required bounding box split,rastershift False -1 99.3650121065 4.96665456014 120.02334544 22.2749878935 HWRF /usr/bin/gdalwarp -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.18//wind10m_popfile_all_clipped.tif" >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density gdalwarp -r near -tr 0.00833333333333 0.00833333333333 -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.18//wind10m_countryfile_all_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.0083333333334 cellsize 0.00833333333333 >> 5. count the popolation in each cell and assign to the class and write to output >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.18/wind_popDensValues_all.xml >> 7. remove files done Creating output file that is 2479P x 2077L. Processing input file /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif. Using internal nodata values (e.g. 3.40282e+38) for image /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_res_final.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/lspop20141.tif. Using internal nodata values (e.g. -2.14748e+09) for image /mnt/output/GDACS/DATA/lspop20141.tif. Copying nodata values from source /mnt/output/GDACS/DATA/lspop20141.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_popfile_final_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/countries.tif. Using internal nodata values (e.g. -32768) for image /mnt/output/GDACS/DATA/countries.tif. Copying nodata values from source /mnt/output/GDACS/DATA/countries.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_countryfile_final_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. t0 completed ============================================ 4. Copy files ============================================ check t0 xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues_t0.xml check all xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind.jpg >> !!!! Check Copy !!!! , listfiles_copy , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml , listfiles_all , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml all files copied copycheck: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/20191029.18_copy_completed_wind.txt check t0 xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues_t0.xml check all xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/20191030.00_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind.jpg >> !!!! Check Copy !!!! , listfiles_copy , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues.xml , listfiles_all , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues.xml all files copied copycheck: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/20191030.00_copy_completed_wind.txt check t0 xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues_t0.xml check all xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.06/20191030.06_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind.jpg >> !!!! Check Copy !!!! , listfiles_copy , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues.xml , listfiles_all , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues.xml all files copied copycheck: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/20191030.06_copy_completed_wind.txt check t0 xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.12/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind_popDensValues_t0.xml check all xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.12/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.12/20191030.12_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind.jpg >> !!!! Check Copy !!!! , listfiles_copy , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind_popDensValues.xml , listfiles_all , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind_popDensValues.xml all files copied copycheck: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/20191030.12_copy_completed_wind.txt check t0 xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.18/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.18/wind_popDensValues_t0.xml check all xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.18/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.18/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.18/20191030.18_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.18/wind.jpg >> !!!! Check Copy !!!! , listfiles_copy , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.18/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.18/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.18/wind_popDensValues.xml , listfiles_all , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.18/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.18/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.18/wind_popDensValues.xml all files copied copycheck: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.18/20191030.18_copy_completed_wind.txt >> 3.2. Classify final folder /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final/20191030.18_final_completed_wind.txt input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final// PopFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_popfile_final_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_countryfile_final_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final/wind_popDensValues_final.xml xml file exists...REMOVE popfile: LandScan popCellSize= 0.00833333333333 projection= GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]] >> 1. resample the tif file to the resolution and proj of pop density 0.00833333333333 deg >> 2. read the charactristics of the input file >> 3a. extract a piece of pop. file corresponding to the required bounding box split,rastershift False -1 99.3650121065 4.96665456014 120.02334544 22.2749878935 HWRF /usr/bin/gdalwarp -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_popfile_final_clipped.tif" >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density gdalwarp -r near -tr 0.00833333333333 0.00833333333333 -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_countryfile_final_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.0083333333334 cellsize 0.00833333333333 >> 5. count the popolation in each cell and assign to the class and write to output >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final/wind_popDensValues_final.xml >> 7. remove files done copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final/wind_popDensValues_final.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/final/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_FINAL.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/final/wind.jpg ******************************************************* * Case Completed * ******************************************************* Case starterd at: 2019-10-31 00:45:07.285041 UTC Case completed at: 2019-10-31 00:46:37.913835 UTC **********E N D O F J O B*********************** Creating output file that is 413P x 346L. Processing input file /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_6h.tif. Using internal nodata values (e.g. -3.40282e+38) for image /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_6h.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 413P x 346L. Processing input file /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_72h.tif. Using internal nodata values (e.g. -3.40282e+38) for image /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_72h.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. 0 .. 10 .. 20 .. 30 .. 40 .. 50 .. 60 .. 70 .. 80 .. 90 .. 100 - Done 0 .. 10 .. 20 .. 30 .. 40 .. 50 .. 60 .. 70 .. 80 .. 90 .. 100 - Done 0 .. 10 .. 20 .. 30 .. 40 .. 50 .. 60 .. 70 .. 80 .. 90 .. 100 - Done 0 .. 10 .. 20 .. 30 .. 40 .. 50 .. 60 .. 70 .. 80 .. 90 .. 100 - Done 0 .. 10 .. 20 .. 30 .. 40 .. 50 .. 60 .. 70 .. 80 .. 90 .. 100 - Done ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-31 07:55:10.236888 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-31 00:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 6 it: 0 ndt: 6 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 6 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 6 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif input file exist it: 3 ndt: 6 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif input file exist it: 4 ndt: 6 idate: 2019-10-30 18:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_6h.tif input file exist it: 5 ndt: 6 idate: 2019-10-31 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_6h.tif input file exist GDACS/1000624/1_HWRF 1000624 ndt: 6 it: 0 ndt: 6 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_72h.tif input file exist it: 1 ndt: 6 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_72h.tif input file exist it: 2 ndt: 6 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_72h.tif input file exist it: 3 ndt: 6 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_72h.tif input file exist it: 4 ndt: 6 idate: 2019-10-30 18:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_72h.tif input file exist it: 5 ndt: 6 idate: 2019-10-31 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_72h.tif input file exist home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret 0 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-10-29 18:00:00 currdate 2019-10-31 00:00:00 ndt: 30 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-10-29 18:00:00', '2019-10-30 00:00:00', '2019-10-30 06:00:00', '2019-10-30 12:00:00', '2019-10-30 18:00:00', '2019-10-31 00:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files lon0,lon1,lat0,lat1 99.34 120.0 5.0 22.3 processing curr bull... 2019-10-31 00:00:00 99.34 120.0 99.34 120.0 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191031.00/ /mnt/input/grib_HWRF/2019/1000624/ input file /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_6h.tif /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191031.00/20191031.00_wind_stept0.tif -r bilinear -te 99.34 5.0 120.0 22.3 -dstnodata 0 Resample tif file /usr/bin/gdalwarp -r bilinear -te 99.34 5.0 120.0 22.3 -dstnodata 0 "/mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_6h.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191031.00/20191031.00_wind_stept0.tif" outFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191031.00/20191031.00_wind_stept0.tif dd: 0.0500242130751 >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103100/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103100/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103100/all_inpData.txt 99.3650121065 119.974987893 5.01663438257 22.2749878935 lon 413 lat 346 read TC track /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103100/all_inpData.txt currDate 2019-10-31 00:00:00 f 6 nt-1 5 nt 6 currDate1 2019-10-31 00:00:00 2019-10-31 06:00:00 sdate 2019-10-31 06:00:00 2019-10-31 06:00:00 6 CAT ['TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TD'] savemap /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191031.00/20191031.00_wind_stept0.jpg 99.34 120.0 99.34 120.0 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191031.00/ /mnt/input/grib_HWRF/2019/1000624/ input file /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_72h.tif /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191031.00/20191031.00_wind.tif -r bilinear -te 99.34 5.0 120.0 22.3 -dstnodata 0 Resample tif file /usr/bin/gdalwarp -r bilinear -te 99.34 5.0 120.0 22.3 -dstnodata 0 "/mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_72h.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191031.00/20191031.00_wind.tif" outFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191031.00/20191031.00_wind.tif dd: 0.0500242130751 >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103100/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103100/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103100/all_inpData.txt 99.3650121065 119.974987893 5.01663438257 22.2749878935 lon 413 lat 346 read TC track /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103100/all_inpData.txt currDate 2019-10-31 00:00:00 f 6 nt-1 71 nt 72 currDate1 2019-10-31 00:00:00 2019-11-03 00:00:00 sdate 2019-11-03 00:00:00 2019-11-03 00:00:00 17 CAT ['TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD'] savemap /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191031.00/20191031.00_wind.jpg processing all past bull only if Past=True... True itdate, istime 2019-10-29 18:00:00 20191029.18 meteo-processing past forecast already completed itdate, istime 2019-10-30 00:00:00 20191030.00 meteo-processing past forecast already completed itdate, istime 2019-10-30 06:00:00 20191030.06 meteo-processing past forecast already completed itdate, istime 2019-10-30 12:00:00 20191030.12 meteo-processing past forecast already completed itdate, istime 2019-10-30 18:00:00 20191030.18 meteo-processing past forecast already completed >> 2.2. Process final data files compfile /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/20191031.00_Final_completed_wind.txt FINAL remove maxtif alldate: DatetimeIndex(['2019-10-29 18:00:00', '2019-10-30 00:00:00', '2019-10-30 06:00:00', '2019-10-30 12:00:00', '2019-10-30 18:00:00', '2019-10-31 00:00:00'], dtype='datetime64[ns]', freq='6H') alldate.shape 6 date: 2019-10-30 00:00:00 python /usr/bin/gdal_calc.py -A /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind_stept0.tif -B /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/20191030.00_wind_stept0.tif --outfile=/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif --calc="maximum(A,B)" date: 2019-10-30 06:00:00 python /usr/bin/gdal_calc.py -A /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif -B /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.06/20191030.06_wind_stept0.tif --outfile=/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif --calc="maximum(A,B)" date: 2019-10-30 12:00:00 python /usr/bin/gdal_calc.py -A /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif -B /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.12/20191030.12_wind_stept0.tif --outfile=/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif --calc="maximum(A,B)" date: 2019-10-30 18:00:00 python /usr/bin/gdal_calc.py -A /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif -B /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.18/20191030.18_wind_stept0.tif --outfile=/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif --calc="maximum(A,B)" date: 2019-10-31 00:00:00 python /usr/bin/gdal_calc.py -A /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif -B /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191031.00/20191031.00_wind.tif --outfile=/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif --calc="maximum(A,B)" max file created /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif stormname 1000624/23w stormname 1000624/23w TCbullsource 23w >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103100/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103100/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103100/all_inpData.txt 99.3900242131 120.0 4.99162227603 22.2499757869 lon 413 lat 346 read TC track /mnt/output/GDACS/TC/TCTrack/1000624/HWRF/2019103100/all_inpData.txt currDate 2019-10-31 00:00:00 f 6 nt-1 71 nt 72 currDate1 2019-10-31 00:00:00 2019-11-03 00:00:00 sdate 2019-11-03 00:00:00 2019-11-03 00:00:00 17 CAT ['TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TD'] savemap /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_FINAL.jpg ret: 0 ============================================ 3. Classify meteo + GDACS index score ============================================ >> 3.1. Classify curr + past forecast Creating output file that is 2479P x 2077L. Processing input file /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191031.00/20191031.00_wind_stept0.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191031.00/20191031.00_wind_stept0.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191031.00/20191031.00_wind_stept0.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191031.00//wind10m_res_stept0.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/lspop20141.tif. Using internal nodata values (e.g. -2.14748e+09) for image /mnt/output/GDACS/DATA/lspop20141.tif. Copying nodata values from source /mnt/output/GDACS/DATA/lspop20141.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191031.00//wind10m_popfile_stept0_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/countries.tif. Using internal nodata values (e.g. -32768) for image /mnt/output/GDACS/DATA/countries.tif. Copying nodata values from source /mnt/output/GDACS/DATA/countries.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191031.00//wind10m_countryfile_stept0_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191031.00/20191031.00_wind_stept0.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191031.00// PopFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191031.00//wind10m_popfile_stept0_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191031.00//wind10m_countryfile_stept0_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191031.00/wind_popDensValues_stept0.xml outDir created popfile: LandScan popCellSize= 0.00833333333333 projection= GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]] >> 1. resample the tif file to the resolution and proj of pop density 0.00833333333333 deg >> 2. read the charactristics of the input file >> 3a. extract a piece of pop. file corresponding to the required bounding box split,rastershift False -1 99.3650121065 4.96665456014 120.02334544 22.2749878935 HWRF /usr/bin/gdalwarp -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191031.00//wind10m_popfile_stept0_clipped.tif" >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density gdalwarp -r near -tr 0.00833333333333 0.00833333333333 -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191031.00//wind10m_countryfile_stept0_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.0083333333334 cellsize 0.00833333333333 >> 5. count the popolation in each cell and assign to the class and write to output >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191031.00/wind_popDensValues_stept0.xml >> 7. remove files done t0 completed Creating output file that is 2479P x 2077L. Processing input file /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191031.00/20191031.00_wind.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191031.00/20191031.00_wind.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191031.00/20191031.00_wind.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191031.00//wind10m_res_all.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/lspop20141.tif. Using internal nodata values (e.g. -2.14748e+09) for image /mnt/output/GDACS/DATA/lspop20141.tif. Copying nodata values from source /mnt/output/GDACS/DATA/lspop20141.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191031.00//wind10m_popfile_all_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/countries.tif. Using internal nodata values (e.g. -32768) for image /mnt/output/GDACS/DATA/countries.tif. Copying nodata values from source /mnt/output/GDACS/DATA/countries.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191031.00//wind10m_countryfile_all_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191031.00/20191031.00_wind.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191031.00// PopFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191031.00//wind10m_popfile_all_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191031.00//wind10m_countryfile_all_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191031.00/wind_popDensValues_all.xml popfile: LandScan popCellSize= 0.00833333333333 projection= GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]] >> 1. resample the tif file to the resolution and proj of pop density 0.00833333333333 deg >> 2. read the charactristics of the input file >> 3a. extract a piece of pop. file corresponding to the required bounding box split,rastershift False -1 99.3650121065 4.96665456014 120.02334544 22.2749878935 HWRF /usr/bin/gdalwarp -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191031.00//wind10m_popfile_all_clipped.tif" >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density gdalwarp -r near -tr 0.00833333333333 0.00833333333333 -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191031.00//wind10m_countryfile_all_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.0083333333334 cellsize 0.00833333333333 >> 5. count the popolation in each cell and assign to the class and write to output >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191031.00/wind_popDensValues_all.xml >> 7. remove files done Creating output file that is 2479P x 2077L. Processing input file /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif. Using internal nodata values (e.g. 3.40282e+38) for image /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_res_final.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/lspop20141.tif. Using internal nodata values (e.g. -2.14748e+09) for image /mnt/output/GDACS/DATA/lspop20141.tif. Copying nodata values from source /mnt/output/GDACS/DATA/lspop20141.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_popfile_final_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2479P x 2077L. Processing input file /mnt/output/GDACS/DATA/countries.tif. Using internal nodata values (e.g. -32768) for image /mnt/output/GDACS/DATA/countries.tif. Copying nodata values from source /mnt/output/GDACS/DATA/countries.tif to destination /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_countryfile_final_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. t0 completed ============================================ 4. Copy files ============================================ check t0 xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues_t0.xml check all xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191029.18/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191029.18/20191029.18_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind.jpg >> !!!! Check Copy !!!! , listfiles_copy , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml , listfiles_all , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/wind_popDensValues.xml all files copied copycheck: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191029.18/20191029.18_copy_completed_wind.txt check t0 xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues_t0.xml check all xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.00/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.00/20191030.00_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind.jpg >> !!!! Check Copy !!!! , listfiles_copy , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues.xml , listfiles_all , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/wind_popDensValues.xml all files copied copycheck: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.00/20191030.00_copy_completed_wind.txt check t0 xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues_t0.xml check all xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.06/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.06/20191030.06_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind.jpg >> !!!! Check Copy !!!! , listfiles_copy , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues.xml , listfiles_all , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/wind_popDensValues.xml all files copied copycheck: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.06/20191030.06_copy_completed_wind.txt check t0 xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.12/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind_popDensValues_t0.xml check all xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.12/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.12/20191030.12_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind.jpg >> !!!! Check Copy !!!! , listfiles_copy , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind_popDensValues.xml , listfiles_all , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/wind_popDensValues.xml all files copied copycheck: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.12/20191030.12_copy_completed_wind.txt check t0 xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.18/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.18/wind_popDensValues_t0.xml check all xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191030.18/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.18/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191030.18/20191030.18_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.18/wind.jpg >> !!!! Check Copy !!!! , listfiles_copy , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.18/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.18/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.18/wind_popDensValues.xml , listfiles_all , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.18/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.18/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.18/wind_popDensValues.xml all files copied copycheck: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191030.18/20191030.18_copy_completed_wind.txt check t0 xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191031.00/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191031.00/wind_popDensValues_t0.xml check all xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/20191031.00/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191031.00/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/20191031.00/20191031.00_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191031.00/wind.jpg >> !!!! Check Copy !!!! , listfiles_copy , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191031.00/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191031.00/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191031.00/wind_popDensValues.xml , listfiles_all , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191031.00/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191031.00/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191031.00/wind_popDensValues.xml all files copied copycheck: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/20191031.00/20191031.00_copy_completed_wind.txt >> 3.2. Classify final folder /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final/20191031.00_final_completed_wind.txt input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_final.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final// PopFile: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_popfile_final_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_countryfile_final_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final/wind_popDensValues_final.xml xml file exists...REMOVE popfile: LandScan popCellSize= 0.00833333333333 projection= GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]] >> 1. resample the tif file to the resolution and proj of pop density 0.00833333333333 deg >> 2. read the charactristics of the input file >> 3a. extract a piece of pop. file corresponding to the required bounding box split,rastershift False -1 99.3650121065 4.96665456014 120.02334544 22.2749878935 HWRF /usr/bin/gdalwarp -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_popfile_final_clipped.tif" >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density gdalwarp -r near -tr 0.00833333333333 0.00833333333333 -te 99.3650121065 4.96665456014 120.02334544 22.2749878935 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final//wind10m_countryfile_final_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.0083333333334 cellsize 0.00833333333333 >> 5. count the popolation in each cell and assign to the class and write to output >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final/wind_popDensValues_final.xml >> 7. remove files done copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/class/final/wind_popDensValues_final.xml in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/final/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/tif/final/wind_FINAL.jpg in: /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/delft3d/final/wind.jpg ******************************************************* * Case Completed * ******************************************************* Case starterd at: 2019-10-31 07:55:10.236888 UTC Case completed at: 2019-10-31 07:56:40.901756 UTC **********E N D O F J O B*********************** ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-31 08:05:07.915999 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-31 00:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 6 it: 0 ndt: 6 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 6 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 6 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif input file exist it: 3 ndt: 6 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif input file exist it: 4 ndt: 6 idate: 2019-10-30 18:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_6h.tif input file exist it: 5 ndt: 6 idate: 2019-10-31 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_6h.tif input file exist GDACS/1000624/1_HWRF 1000624 ndt: 6 it: 0 ndt: 6 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_72h.tif input file exist it: 1 ndt: 6 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_72h.tif input file exist it: 2 ndt: 6 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_72h.tif input file exist it: 3 ndt: 6 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_72h.tif input file exist it: 4 ndt: 6 idate: 2019-10-30 18:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_72h.tif input file exist it: 5 ndt: 6 idate: 2019-10-31 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_72h.tif input file exist home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret 0 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-10-29 18:00:00 currdate 2019-10-31 00:00:00 ndt: 30 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-10-29 18:00:00', '2019-10-30 00:00:00', '2019-10-30 06:00:00', '2019-10-30 12:00:00', '2019-10-30 18:00:00', '2019-10-31 00:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files lon0,lon1,lat0,lat1 99.34 120.0 5.0 22.3 processing curr bull... 2019-10-31 00:00:00 ********** no new data check copy copy ok copy ok copy ok copy ok copy ok copy ok no new data ******************************************************* * Case Completed * ******************************************************* Case starterd at: 2019-10-31 08:05:07.915999 UTC Case completed at: 2019-10-31 08:05:08.280406 UTC **********E N D O F J O B*********************** ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-31 12:35:07.637284 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-31 06:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 7 it: 0 ndt: 7 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 7 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 7 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif input file exist it: 3 ndt: 7 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif input file exist it: 4 ndt: 7 idate: 2019-10-30 18:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_6h.tif input file exist it: 5 ndt: 7 idate: 2019-10-31 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_6h.tif input file exist it: 6 ndt: 7 idate: 2019-10-31 06:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-31 12:45:06.700377 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-31 06:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 7 it: 0 ndt: 7 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 7 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 7 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif input file exist it: 3 ndt: 7 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif input file exist it: 4 ndt: 7 idate: 2019-10-30 18:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_6h.tif input file exist it: 5 ndt: 7 idate: 2019-10-31 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_6h.tif input file exist it: 6 ndt: 7 idate: 2019-10-31 06:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-31 12:55:11.252375 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-31 06:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 7 it: 0 ndt: 7 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 7 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 7 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif input file exist it: 3 ndt: 7 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif input file exist it: 4 ndt: 7 idate: 2019-10-30 18:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_6h.tif input file exist it: 5 ndt: 7 idate: 2019-10-31 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_6h.tif input file exist it: 6 ndt: 7 idate: 2019-10-31 06:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-31 13:05:06.868404 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-31 06:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 7 it: 0 ndt: 7 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 7 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 7 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif input file exist it: 3 ndt: 7 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif input file exist it: 4 ndt: 7 idate: 2019-10-30 18:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_6h.tif input file exist it: 5 ndt: 7 idate: 2019-10-31 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_6h.tif input file exist it: 6 ndt: 7 idate: 2019-10-31 06:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-31 13:15:07.494090 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-31 06:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 7 it: 0 ndt: 7 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 7 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 7 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif input file exist it: 3 ndt: 7 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif input file exist it: 4 ndt: 7 idate: 2019-10-30 18:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_6h.tif input file exist it: 5 ndt: 7 idate: 2019-10-31 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_6h.tif input file exist it: 6 ndt: 7 idate: 2019-10-31 06:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-31 13:25:11.695448 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-31 06:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 7 it: 0 ndt: 7 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 7 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 7 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif input file exist it: 3 ndt: 7 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif input file exist it: 4 ndt: 7 idate: 2019-10-30 18:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_6h.tif input file exist it: 5 ndt: 7 idate: 2019-10-31 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_6h.tif input file exist it: 6 ndt: 7 idate: 2019-10-31 06:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-31 15:05:07.120597 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-31 06:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 7 it: 0 ndt: 7 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 7 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 7 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif input file exist it: 3 ndt: 7 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif input file exist it: 4 ndt: 7 idate: 2019-10-30 18:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_6h.tif input file exist it: 5 ndt: 7 idate: 2019-10-31 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_6h.tif input file exist it: 6 ndt: 7 idate: 2019-10-31 06:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-31 15:15:09.124729 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-31 06:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 7 it: 0 ndt: 7 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 7 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 7 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif input file exist it: 3 ndt: 7 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif input file exist it: 4 ndt: 7 idate: 2019-10-30 18:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_6h.tif input file exist it: 5 ndt: 7 idate: 2019-10-31 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_6h.tif input file exist it: 6 ndt: 7 idate: 2019-10-31 06:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-31 15:25:13.505528 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-31 06:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 7 it: 0 ndt: 7 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 7 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 7 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif input file exist it: 3 ndt: 7 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif input file exist it: 4 ndt: 7 idate: 2019-10-30 18:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_6h.tif input file exist it: 5 ndt: 7 idate: 2019-10-31 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_6h.tif input file exist it: 6 ndt: 7 idate: 2019-10-31 06:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-31 15:35:07.559955 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-31 06:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 7 it: 0 ndt: 7 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 7 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 7 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif input file exist it: 3 ndt: 7 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif input file exist it: 4 ndt: 7 idate: 2019-10-30 18:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_6h.tif input file exist it: 5 ndt: 7 idate: 2019-10-31 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_6h.tif input file exist it: 6 ndt: 7 idate: 2019-10-31 06:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-31 15:45:06.996857 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-31 06:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 7 it: 0 ndt: 7 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 7 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 7 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif input file exist it: 3 ndt: 7 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif input file exist it: 4 ndt: 7 idate: 2019-10-30 18:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_6h.tif input file exist it: 5 ndt: 7 idate: 2019-10-31 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_6h.tif input file exist it: 6 ndt: 7 idate: 2019-10-31 06:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-31 15:55:09.927637 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-31 06:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 7 it: 0 ndt: 7 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 7 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 7 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif input file exist it: 3 ndt: 7 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif input file exist it: 4 ndt: 7 idate: 2019-10-30 18:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_6h.tif input file exist it: 5 ndt: 7 idate: 2019-10-31 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_6h.tif input file exist it: 6 ndt: 7 idate: 2019-10-31 06:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-31 16:05:06.240254 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-31 06:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 7 it: 0 ndt: 7 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 7 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 7 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif input file exist it: 3 ndt: 7 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif input file exist it: 4 ndt: 7 idate: 2019-10-30 18:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_6h.tif input file exist it: 5 ndt: 7 idate: 2019-10-31 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_6h.tif input file exist it: 6 ndt: 7 idate: 2019-10-31 06:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-31 16:15:07.580256 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-31 06:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 7 it: 0 ndt: 7 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 7 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 7 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif input file exist it: 3 ndt: 7 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif input file exist it: 4 ndt: 7 idate: 2019-10-30 18:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_6h.tif input file exist it: 5 ndt: 7 idate: 2019-10-31 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_6h.tif input file exist it: 6 ndt: 7 idate: 2019-10-31 06:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-31 16:25:11.502027 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-31 06:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 7 it: 0 ndt: 7 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 7 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 7 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif input file exist it: 3 ndt: 7 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif input file exist it: 4 ndt: 7 idate: 2019-10-30 18:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_6h.tif input file exist it: 5 ndt: 7 idate: 2019-10-31 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_6h.tif input file exist it: 6 ndt: 7 idate: 2019-10-31 06:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-31 16:35:07.103027 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-31 06:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 7 it: 0 ndt: 7 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 7 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 7 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif input file exist it: 3 ndt: 7 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif input file exist it: 4 ndt: 7 idate: 2019-10-30 18:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_6h.tif input file exist it: 5 ndt: 7 idate: 2019-10-31 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_6h.tif input file exist it: 6 ndt: 7 idate: 2019-10-31 06:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-31 16:45:14.414799 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-31 06:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 7 it: 0 ndt: 7 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 7 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 7 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif input file exist it: 3 ndt: 7 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif input file exist it: 4 ndt: 7 idate: 2019-10-30 18:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_6h.tif input file exist it: 5 ndt: 7 idate: 2019-10-31 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_6h.tif input file exist it: 6 ndt: 7 idate: 2019-10-31 06:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-10-31 16:55:10.247308 UTC inp1= 20191029.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000624/23w submitting calc 2019-10-29 18:00:00 2019-10-31 06:00:00 delft3d 99.34 120.0 5.0 22.3 4.0 HWRF 72 15 True GDACS/1000624/1_HWRF 6 1 False False 5 aa.txt 20191029.18 1000624/23w 6 wind False *************---------------------****************** GDACS/1000624/1_HWRF 1000624 ndt: 7 it: 0 ndt: 7 idate: 2019-10-29 18:00:00 /mnt/input/grib_HWRF/2019/1000624/twenty-th23w.2019102918.swath.wind_6h.tif input file exist it: 1 ndt: 7 idate: 2019-10-30 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103000.swath.wind_6h.tif input file exist it: 2 ndt: 7 idate: 2019-10-30 06:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103006.swath.wind_6h.tif input file exist it: 3 ndt: 7 idate: 2019-10-30 12:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103012.swath.wind_6h.tif input file exist it: 4 ndt: 7 idate: 2019-10-30 18:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103018.swath.wind_6h.tif input file exist it: 5 ndt: 7 idate: 2019-10-31 00:00:00 /mnt/input/grib_HWRF/2019/1000624/matmo23w.2019103100.swath.wind_6h.tif input file exist it: 6 ndt: 7 idate: 2019-10-31 06:00:00 meteo files not yet existing home dir /mnt/output/SSCS/2019/GDACS/1000624/1_HWRF/ ret -4 no impact estimation ==============================================================