******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-11-25 23:05:06.505488 UTC inp1= 20191125.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000636/29w submitting calc 2019-11-25 18:00:00 2019-11-25 18:00:00 delft3d 127.25 154.63 5.0 16.45 4.0 HWRF 72 15 True GDACS/1000636/1_HWRF 6 1 False False 5 aa.txt 20191125.18 1000636/29w 6 wind False *************---------------------****************** GDACS/1000636/1_HWRF 1000636 ndt: 1 it: 0 ndt: 1 idate: 2019-11-25 18:00:00 meteo files not yet existing /mnt/input/grib_HWRF/2019/1000636/.2019112518.swath.wind_6h.tif ********* ------ >input file doesn't exist /mnt/input/grib_HWRF/2019/1000636/.2019112518.swath.wind_6h.tif home dir /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-11-25 23:15:06.488647 UTC inp1= 20191125.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000636/29w submitting calc 2019-11-25 18:00:00 2019-11-25 18:00:00 delft3d 127.25 154.63 5.0 16.45 4.0 HWRF 72 15 True GDACS/1000636/1_HWRF 6 1 False False 5 aa.txt 20191125.18 1000636/29w 6 wind False *************---------------------****************** GDACS/1000636/1_HWRF 1000636 ndt: 1 it: 0 ndt: 1 idate: 2019-11-25 18:00:00 meteo files not yet existing /mnt/input/grib_HWRF/2019/1000636/.2019112518.swath.wind_6h.tif ********* ------ >input file doesn't exist /mnt/input/grib_HWRF/2019/1000636/.2019112518.swath.wind_6h.tif home dir /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-11-25 23:25:10.471669 UTC inp1= 20191125.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000636/29w submitting calc 2019-11-25 18:00:00 2019-11-25 18:00:00 delft3d 127.25 154.63 5.0 16.45 4.0 HWRF 72 15 True GDACS/1000636/1_HWRF 6 1 False False 5 aa.txt 20191125.18 1000636/29w 6 wind False *************---------------------****************** GDACS/1000636/1_HWRF 1000636 ndt: 1 it: 0 ndt: 1 idate: 2019-11-25 18:00:00 meteo files not yet existing /mnt/input/grib_HWRF/2019/1000636/.2019112518.swath.wind_6h.tif ********* ------ >input file doesn't exist /mnt/input/grib_HWRF/2019/1000636/.2019112518.swath.wind_6h.tif home dir /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-11-25 23:35:06.778522 UTC inp1= 20191125.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000636/29w submitting calc 2019-11-25 18:00:00 2019-11-25 18:00:00 delft3d 127.25 154.63 5.0 16.45 4.0 HWRF 72 15 True GDACS/1000636/1_HWRF 6 1 False False 5 aa.txt 20191125.18 1000636/29w 6 wind False *************---------------------****************** GDACS/1000636/1_HWRF 1000636 ndt: 1 it: 0 ndt: 1 idate: 2019-11-25 18:00:00 meteo files not yet existing /mnt/input/grib_HWRF/2019/1000636/.2019112518.swath.wind_6h.tif ********* ------ >input file doesn't exist /mnt/input/grib_HWRF/2019/1000636/.2019112518.swath.wind_6h.tif home dir /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-11-25 23:45:06.868814 UTC inp1= 20191125.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000636/29w submitting calc 2019-11-25 18:00:00 2019-11-25 18:00:00 delft3d 127.25 154.63 5.0 16.45 4.0 HWRF 72 15 True GDACS/1000636/1_HWRF 6 1 False False 5 aa.txt 20191125.18 1000636/29w 6 wind False *************---------------------****************** GDACS/1000636/1_HWRF 1000636 ndt: 1 it: 0 ndt: 1 idate: 2019-11-25 18:00:00 meteo files not yet existing /mnt/input/grib_HWRF/2019/1000636/.2019112518.swath.wind_6h.tif ********* ------ >input file doesn't exist /mnt/input/grib_HWRF/2019/1000636/.2019112518.swath.wind_6h.tif home dir /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-11-25 23:55:10.010125 UTC inp1= 20191125.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000636/29w submitting calc 2019-11-25 18:00:00 2019-11-25 18:00:00 delft3d 127.25 154.63 5.0 16.45 4.0 HWRF 72 15 True GDACS/1000636/1_HWRF 6 1 False False 5 aa.txt 20191125.18 1000636/29w 6 wind False *************---------------------****************** GDACS/1000636/1_HWRF 1000636 ndt: 1 it: 0 ndt: 1 idate: 2019-11-25 18:00:00 meteo files not yet existing /mnt/input/grib_HWRF/2019/1000636/.2019112518.swath.wind_6h.tif ********* ------ >input file doesn't exist /mnt/input/grib_HWRF/2019/1000636/.2019112518.swath.wind_6h.tif home dir /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-11-26 00:05:07.260347 UTC inp1= 20191125.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000636/29w submitting calc 2019-11-25 18:00:00 2019-11-25 18:00:00 delft3d 127.25 154.63 5.0 16.45 4.0 HWRF 72 15 True GDACS/1000636/1_HWRF 6 1 False False 5 aa.txt 20191125.18 1000636/29w 6 wind False *************---------------------****************** GDACS/1000636/1_HWRF 1000636 ndt: 1 it: 0 ndt: 1 idate: 2019-11-25 18:00:00 meteo files not yet existing /mnt/input/grib_HWRF/2019/1000636/.2019112518.swath.wind_6h.tif ********* ------ >input file doesn't exist /mnt/input/grib_HWRF/2019/1000636/.2019112518.swath.wind_6h.tif home dir /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-11-26 00:15:07.079194 UTC inp1= 20191125.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000636/29w submitting calc 2019-11-25 18:00:00 2019-11-25 18:00:00 delft3d 127.25 154.63 5.0 16.45 4.0 HWRF 72 15 True GDACS/1000636/1_HWRF 6 1 False False 5 aa.txt 20191125.18 1000636/29w 6 wind False *************---------------------****************** GDACS/1000636/1_HWRF 1000636 ndt: 1 it: 0 ndt: 1 idate: 2019-11-25 18:00:00 meteo files not yet existing /mnt/input/grib_HWRF/2019/1000636/.2019112518.swath.wind_6h.tif ********* ------ >input file doesn't exist /mnt/input/grib_HWRF/2019/1000636/.2019112518.swath.wind_6h.tif home dir /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/ ret -4 no impact estimation ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-11-26 00:25:10.808986 UTC inp1= 20191125.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000636/29w submitting calc 2019-11-25 18:00:00 2019-11-25 18:00:00 delft3d 127.25 154.63 5.0 16.45 4.0 HWRF 72 15 True GDACS/1000636/1_HWRF 6 1 False False 5 aa.txt 20191125.18 1000636/29w 6 wind False *************---------------------****************** GDACS/1000636/1_HWRF 1000636 ndt: 1 it: 0 ndt: 1 idate: 2019-11-25 18:00:00 /mnt/input/grib_HWRF/2019/1000636/twenty-ni29w.2019112518.swath.wind_6h.tif input file exist GDACS/1000636/1_HWRF 1000636 ndt: 1 it: 0 ndt: 1 idate: 2019-11-25 18:00:00 /mnt/input/grib_HWRF/2019/1000636/twenty-ni29w.2019112518.swath.wind_72h.tif input file exist home dir /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/ ret 0 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-11-25 18:00:00 currdate 2019-11-25 18:00:00 ndt: 0 delta: 6 nt 72 alldate: DatetimeIndex(['2019-11-25 18:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files lon0,lon1,lat0,lat1 127.25 154.63 5.0 16.45 processing curr bull... 2019-11-25 18:00:00 127.25 154.63 127.25 154.63 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/tif/20191125.18/ /mnt/input/grib_HWRF/2019/1000636/ Creating output file that is 548P x 229L. Processing input file /mnt/input/grib_HWRF/2019/1000636/twenty-ni29w.2019112518.swath.wind_6h.tif. Using internal nodata values (e.g. -3.40282e+38) for image /mnt/input/grib_HWRF/2019/1000636/twenty-ni29w.2019112518.swath.wind_6h.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 548P x 229L. Processing input file /mnt/input/grib_HWRF/2019/1000636/twenty-ni29w.2019112518.swath.wind_72h.tif. Using internal nodata values (e.g. -3.40282e+38) for image /mnt/input/grib_HWRF/2019/1000636/twenty-ni29w.2019112518.swath.wind_72h.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-11-26 00:35:06.225298 UTC inp1= 20191125.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000636/29w submitting calc 2019-11-25 18:00:00 2019-11-25 18:00:00 delft3d 127.25 154.63 5.0 16.45 4.0 HWRF 72 15 True GDACS/1000636/1_HWRF 6 1 False False 5 aa.txt 20191125.18 1000636/29w 6 wind False *************---------------------****************** GDACS/1000636/1_HWRF 1000636 ndt: 1 it: 0 ndt: 1 idate: 2019-11-25 18:00:00 /mnt/input/grib_HWRF/2019/1000636/twenty-ni29w.2019112518.swath.wind_6h.tif input file exist GDACS/1000636/1_HWRF 1000636 ndt: 1 it: 0 ndt: 1 idate: 2019-11-25 18:00:00 /mnt/input/grib_HWRF/2019/1000636/twenty-ni29w.2019112518.swath.wind_72h.tif input file exist home dir /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/ ret 0 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-11-25 18:00:00 currdate 2019-11-25 18:00:00 ndt: 0 delta: 6 nt 72 alldate: DatetimeIndex(['2019-11-25 18:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files lon0,lon1,lat0,lat1 127.25 154.63 5.0 16.45 processing curr bull... 2019-11-25 18:00:00 127.25 154.63 127.25 154.63 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/tif/20191125.18/ /mnt/input/grib_HWRF/2019/1000636/ input file /mnt/input/grib_HWRF/2019/1000636/twenty-ni29w.2019112518.swath.wind_6h.tif /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/tif/20191125.18/20191125.18_wind_stept0.tif -r bilinear -te 127.25 5.0 154.63 16.45 -dstnodata 0 Resample tif file /usr/bin/gdalwarp -r bilinear -te 127.25 5.0 154.63 16.45 -dstnodata 0 "/mnt/input/grib_HWRF/2019/1000636/twenty-ni29w.2019112518.swath.wind_6h.tif" "/mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/tif/20191125.18/20191125.18_wind_stept0.tif" outFile: /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/tif/20191125.18/20191125.18_wind_stept0.tif dd: 0.0499635036496 >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000636/HWRF/2019112518/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000636/HWRF/2019112518/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000636/HWRF/2019112518/all_inpData.txt 127.274981752 154.605018248 5.03333941606 16.4250182482 lon 548 lat 229 read TC track /mnt/output/GDACS/TC/TCTrack/1000636/HWRF/2019112518/all_inpData.txt currDate 2019-11-25 18:00:00 f 1 nt-1 5 nt 6 currDate1 2019-11-25 18:00:00 2019-11-26 00:00:00 sdate 2019-11-26 00:00:00 2019-11-26 00:00:00 1 CAT ['TD', 'TS'] savemap /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/tif/20191125.18/20191125.18_wind_stept0.jpg 127.25 154.63 127.25 154.63 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/tif/20191125.18/ /mnt/input/grib_HWRF/2019/1000636/ input file /mnt/input/grib_HWRF/2019/1000636/twenty-ni29w.2019112518.swath.wind_72h.tif /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/tif/20191125.18/20191125.18_wind.tif -r bilinear -te 127.25 5.0 154.63 16.45 -dstnodata 0 Resample tif file /usr/bin/gdalwarp -r bilinear -te 127.25 5.0 154.63 16.45 -dstnodata 0 "/mnt/input/grib_HWRF/2019/1000636/twenty-ni29w.2019112518.swath.wind_72h.tif" "/mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/tif/20191125.18/20191125.18_wind.tif" outFile: /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/tif/20191125.18/20191125.18_wind.tif dd: 0.0499635036496 >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000636/HWRF/2019112518/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000636/HWRF/2019112518/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000636/HWRF/2019112518/all_inpData.txt 127.274981752 154.605018248 5.03333941606 16.4250182482 lon 548 lat 229 read TC track /mnt/output/GDACS/TC/TCTrack/1000636/HWRF/2019112518/all_inpData.txt currDate 2019-11-25 18:00:00 f 1 nt-1 71 nt 72 currDate1 2019-11-25 18:00:00 2019-11-28 18:00:00 sdate 2019-11-28 18:00:00 2019-11-28 18:00:00 12 CAT ['TD', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'H', 'TS', 'H', 'H', 'H', 'H'] savemap /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/tif/20191125.18/20191125.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/1000636/1_HWRF/tif/final/20191125.18_Final_completed_wind.txt FINAL **FIRST alldate: DatetimeIndex(['2019-11-25 18:00:00'], dtype='datetime64[ns]', freq='6H') alldate.shape 1 **FIRST cp /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/tif/20191125.18/20191125.18_wind.tif /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/tif/final/wind_final.tif max file created /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/tif/final/wind_final.tif stormname 1000636/29w stormname 1000636/29w TCbullsource 29w >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000636/HWRF/2019112518/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000636/HWRF/2019112518/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000636/HWRF/2019112518/all_inpData.txt 127.299963504 154.63 5.00835766423 16.4000364964 lon 548 lat 229 read TC track /mnt/output/GDACS/TC/TCTrack/1000636/HWRF/2019112518/all_inpData.txt currDate 2019-11-25 18:00:00 f 1 nt-1 71 nt 72 currDate1 2019-11-25 18:00:00 2019-11-28 18:00:00 sdate 2019-11-28 18:00:00 2019-11-28 18:00:00 12 CAT ['TD', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'H', 'TS', 'H', 'H', 'H', 'H'] savemap /mnt/output/SSCS/2019/GDACS/1000636/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 3286P x 1373L. Processing input file /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/tif/20191125.18/20191125.18_wind_stept0.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/tif/20191125.18/20191125.18_wind_stept0.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/tif/20191125.18/20191125.18_wind_stept0.tif to destination /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/class/20191125.18//wind10m_res_stept0.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 3286P x 1373L. 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/1000636/1_HWRF/class/20191125.18//wind10m_popfile_stept0_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 3286P x 1373L. 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/1000636/1_HWRF/class/20191125.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/1000636/1_HWRF/tif/20191125.18/20191125.18_wind_stept0.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/class/20191125.18// PopFile: /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/class/20191125.18//wind10m_popfile_stept0_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/class/20191125.18//wind10m_countryfile_stept0_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/class/20191125.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 127.274981752 4.98335158151 154.658315085 16.4250182482 HWRF /usr/bin/gdalwarp -te 127.274981752 4.98335158151 154.658315085 16.4250182482 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/class/20191125.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 127.274981752 4.98335158151 154.658315085 16.4250182482 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/class/20191125.18//wind10m_countryfile_stept0_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.00833333333323 cellsize 0.00833333333333 >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/class/20191125.18/wind_popDensValues_stept0.xml >> 7. remove files done t0 completed Creating output file that is 3286P x 1373L. Processing input file /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/tif/20191125.18/20191125.18_wind.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/tif/20191125.18/20191125.18_wind.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/tif/20191125.18/20191125.18_wind.tif to destination /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/class/20191125.18//wind10m_res_all.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 3286P x 1373L. 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/1000636/1_HWRF/class/20191125.18//wind10m_popfile_all_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 3286P x 1373L. 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/1000636/1_HWRF/class/20191125.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/1000636/1_HWRF/tif/20191125.18/20191125.18_wind.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/class/20191125.18// PopFile: /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/class/20191125.18//wind10m_popfile_all_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/class/20191125.18//wind10m_countryfile_all_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/class/20191125.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 127.274981752 4.98335158151 154.658315085 16.4250182482 HWRF /usr/bin/gdalwarp -te 127.274981752 4.98335158151 154.658315085 16.4250182482 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/class/20191125.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 127.274981752 4.98335158151 154.658315085 16.4250182482 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/class/20191125.18//wind10m_countryfile_all_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.00833333333323 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/1000636/1_HWRF/class/20191125.18/wind_popDensValues_all.xml >> 7. remove files done Creating output file that is 3286P x 1373L. Processing input file /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/tif/final/wind_final.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/tif/final/wind_final.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/tif/final/wind_final.tif to destination /mnt/output/SSCS/2019/GDACS/1000636/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 3286P x 1373L. 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/1000636/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 3286P x 1373L. 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/1000636/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/1000636/1_HWRF/class/20191125.18/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/delft3d/20191125.18/wind_popDensValues_t0.xml check all xml copy: /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/class/20191125.18/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/delft3d/20191125.18/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/tif/20191125.18/20191125.18_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/delft3d/20191125.18/wind.jpg >> !!!! Check Copy !!!! , listfiles_copy , /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/delft3d/20191125.18/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/delft3d/20191125.18/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/delft3d/20191125.18/wind_popDensValues.xml , listfiles_all , /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/delft3d/20191125.18/wind_popDensValues_t0.xml , /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/delft3d/20191125.18/wind_popDensValues.xml , /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/delft3d/20191125.18/wind_popDensValues.xml all files copied copycheck: /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/delft3d/20191125.18/20191125.18_copy_completed_wind.txt >> 3.2. Classify final folder /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/class/final/20191125.18_final_completed_wind.txt input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/tif/final/wind_final.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/class/final// PopFile: /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/class/final//wind10m_popfile_final_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/class/final//wind10m_countryfile_final_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/class/final/wind_popDensValues_final.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 127.274981752 4.98335158151 154.658315085 16.4250182482 HWRF /usr/bin/gdalwarp -te 127.274981752 4.98335158151 154.658315085 16.4250182482 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000636/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 127.274981752 4.98335158151 154.658315085 16.4250182482 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000636/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.00833333333323 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/1000636/1_HWRF/class/final/wind_popDensValues_final.xml >> 7. remove files done copy: /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/class/final/wind_popDensValues_final.xml in: /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/delft3d/final/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/tif/final/wind_FINAL.jpg in: /mnt/output/SSCS/2019/GDACS/1000636/1_HWRF/delft3d/final/wind.jpg ******************************************************* * Case Completed * ******************************************************* Case starterd at: 2019-11-26 00:35:06.225298 UTC Case completed at: 2019-11-26 00:36:15.901422 UTC **********E N D O F J O B***********************