******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-09-13 10:14:08.236802 UTC inp1= 20190912.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000599/09l submitting calc 2019-09-12 18:00:00 2019-09-13 00:00:00 delft3d -88.5 -60.48 9.55 37.55 4.0 HWRF 72 15 True GDACS/1000599/3_HWRF 6 1 False False 5 aa.txt 20190912.18 1000599/09l 6 wind False *************---------------------****************** GDACS/1000599/3_HWRF 1000599 ndt: 2 it: 0 ndt: 2 idate: 2019-09-12 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_6h.tif input file exist it: 1 ndt: 2 idate: 2019-09-13 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_6h.tif input file exist GDACS/1000599/3_HWRF 1000599 ndt: 2 it: 0 ndt: 2 idate: 2019-09-12 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_72h.tif input file exist it: 1 ndt: 2 idate: 2019-09-13 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_72h.tif input file exist home dir /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/ ret 0 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-09-12 18:00:00 currdate 2019-09-13 00:00:00 ndt: 6 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-09-12 18:00:00', '2019-09-13 00:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files processing curr bull... 2019-09-13 00:00:00 -88.5 -60.48 271.5 299.52 0.05 0.05 Creating output file that is 560P x 560L. Processing input file /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_6h.tif. Using internal nodata values (e.g. -3.40282e+38) for image /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_6h.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 560P x 560L. Processing input file /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_72h.tif. Using internal nodata values (e.g. -3.40282e+38) for image /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_72h.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 560P x 560L. Processing input file /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_6h.tif. Using internal nodata values (e.g. -3.40282e+38) for image /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_6h.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 560P x 560L. Processing input file /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_72h.tif. Using internal nodata values (e.g. -3.40282e+38) for image /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.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-09-13 10:24:13.741424 UTC inp1= 20190912.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000599/09l submitting calc 2019-09-12 18:00:00 2019-09-13 00:00:00 delft3d -88.5 -60.48 9.55 37.55 4.0 HWRF 72 15 True GDACS/1000599/3_HWRF 6 1 False False 5 aa.txt 20190912.18 1000599/09l 6 wind False *************---------------------****************** GDACS/1000599/3_HWRF 1000599 ndt: 2 it: 0 ndt: 2 idate: 2019-09-12 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_6h.tif input file exist it: 1 ndt: 2 idate: 2019-09-13 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_6h.tif input file exist GDACS/1000599/3_HWRF 1000599 ndt: 2 it: 0 ndt: 2 idate: 2019-09-12 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_72h.tif input file exist it: 1 ndt: 2 idate: 2019-09-13 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_72h.tif input file exist home dir /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/ ret 0 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-09-12 18:00:00 currdate 2019-09-13 00:00:00 ndt: 6 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-09-12 18:00:00', '2019-09-13 00:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files processing curr bull... 2019-09-13 00:00:00 -88.5 -60.48 271.5 299.52 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.00/ /mnt/input/grib_HWRF/2019/1000599/ input file /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_6h.tif /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.00/20190913.00_wind_stept0.tif -r bilinear -te 271.5 9.55 299.52 37.55 -dstnodata 0 Resample tif file /usr/bin/gdalwarp -r bilinear -te 271.5 9.55 299.52 37.55 -dstnodata 0 "/mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_6h.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.00/20190913.00_wind_stept0.tif" outFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.00/20190913.00_wind_stept0.tif dd: 0.0500357142857 >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091300/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091300/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091300/all_inpData.txt 271.525017857 299.494982143 9.55501785714 37.5249821429 lon 560 lat 560 read TC track /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091300/all_inpData.txt currDate 2019-09-13 00:00:00 2019-09-16 00:00:00 sdate 2019-09-16 00:00:00 -4 CAT ['TD', 'TD', 'TS', 'TS', 'TD', 'TD', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'H', 'H', 'H', 'H', 'MH', 'H'] savemap /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.00/20190913.00_wind_stept0.jpg -88.5 -60.48 271.5 299.52 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.00/ /mnt/input/grib_HWRF/2019/1000599/ input file /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_72h.tif /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.00/20190913.00_wind.tif -r bilinear -te 271.5 9.55 299.52 37.55 -dstnodata 0 Resample tif file /usr/bin/gdalwarp -r bilinear -te 271.5 9.55 299.52 37.55 -dstnodata 0 "/mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_72h.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.00/20190913.00_wind.tif" outFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.00/20190913.00_wind.tif dd: 0.0500357142857 >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091300/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091300/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091300/all_inpData.txt 271.525017857 299.494982143 9.55501785714 37.5249821429 lon 560 lat 560 read TC track /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091300/all_inpData.txt currDate 2019-09-13 00:00:00 2019-09-16 00:00:00 sdate 2019-09-16 00:00:00 -4 CAT ['TD', 'TD', 'TS', 'TS', 'TD', 'TD', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'H', 'H', 'H', 'H', 'MH', 'H'] savemap /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.00/20190913.00_wind.jpg processing all past bull only if Past=True... True itdate, istime 2019-09-12 18:00:00 20190912.18 -88.5 -60.48 271.5 299.52 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190912.18/ /mnt/input/grib_HWRF/2019/1000599/ input file /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_6h.tif /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190912.18/20190912.18_wind_stept0.tif -r bilinear -te 271.5 9.55 299.52 37.55 -dstnodata 0 Resample tif file /usr/bin/gdalwarp -r bilinear -te 271.5 9.55 299.52 37.55 -dstnodata 0 "/mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_6h.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190912.18/20190912.18_wind_stept0.tif" outFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190912.18/20190912.18_wind_stept0.tif dd: 0.0500357142857 >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091218/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091218/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091218/all_inpData.txt 271.525017857 299.494982143 9.55501785714 37.5249821429 lon 560 lat 560 read TC track /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091218/all_inpData.txt currDate 2019-09-12 18:00:00 2019-09-15 18:00:00 sdate 2019-09-15 18:00:00 -4 CAT ['TD', 'TD', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TD', 'TD', 'TD', 'TD'] savemap /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190912.18/20190912.18_wind_stept0.jpg -88.5 -60.48 271.5 299.52 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190912.18/ /mnt/input/grib_HWRF/2019/1000599/ input file /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_72h.tif /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190912.18/20190912.18_wind.tif -r bilinear -te 271.5 9.55 299.52 37.55 -dstnodata 0 Resample tif file /usr/bin/gdalwarp -r bilinear -te 271.5 9.55 299.52 37.55 -dstnodata 0 "/mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_72h.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190912.18/20190912.18_wind.tif" outFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190912.18/20190912.18_wind.tif dd: 0.0500357142857 >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091218/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091218/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091218/all_inpData.txt 271.525017857 299.494982143 9.55501785714 37.5249821429 lon 560 lat 560 read TC track /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091218/all_inpData.txt currDate 2019-09-12 18:00:00 2019-09-15 18:00:00 sdate 2019-09-15 18:00:00 -4 CAT ['TD', 'TD', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TD', 'TD', 'TD', 'TD'] savemap /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190912.18/20190912.18_wind.jpg >> 2.2. Process final data files compfile /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/20190913.00_Final_completed_wind.txt FINAL alldate: DatetimeIndex(['2019-09-12 18:00:00', '2019-09-13 00:00:00'], dtype='datetime64[ns]', freq='6H') alldate.shape 2 date: 2019-09-13 00:00:00 python /mnt/output/SSCS/scripts/gdal_calc.py -A /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190912.18/20190912.18_wind_stept0.tif -B /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.00/20190913.00_wind.tif --outfile=/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif --calc="maximum(A,B)" max file created /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif stormname 1000599/09l stormname 1000599/09l TCbullsource 09l >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091300/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091300/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091300/all_inpData.txt 271.550035714 299.52 9.52999999999 37.4999642857 lon 560 lat 560 read TC track /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091300/all_inpData.txt currDate 2019-09-13 00:00:00 2019-09-16 00:00:00 sdate 2019-09-16 00:00:00 -4 CAT ['TD', 'TD', 'TS', 'TS', 'TD', 'TD', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'H', 'H', 'H', 'H', 'MH', 'H'] savemap /mnt/output/SSCS/2019/GDACS/1000599/3_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 3362P x 3362L. Processing input file /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190912.18/20190912.18_wind_stept0.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190912.18/20190912.18_wind_stept0.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190912.18/20190912.18_wind_stept0.tif to destination /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190912.18//wind10m_res_stept0.tif. 0Creating output file that is 3362P x 3362L. Processing input file /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.00/20190913.00_wind_stept0.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.00/20190913.00_wind_stept0.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.00/20190913.00_wind_stept0.tif to destination /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.00//wind10m_res_stept0.tif. 0.....10..10....20...20..30.....4030.....50..40....60...50..70....60.80......7090......80...90...100 - done. 100 - done. Creating output file that is 3362P x 3362L. 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/1000599/3_HWRF/class/20190912.18//wind10m_popfile_stept0_clipped.tif. 0..Creating output file that is 3362P x 3362L. 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/1000599/3_HWRF/class/20190913.00//wind10m_popfile_stept0_clipped.tif. 0..10....10..20....20..30....30..40....40..50..50...60.....70....6080.....90...70....80...90...100 - done. 100 - done. Creating output file that is 3362P x 3362L. 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/1000599/3_HWRF/class/20190913.00//wind10m_countryfile_stept0_clipped.tif. 0..Creating output file that is 3362P x 3362L. 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/1000599/3_HWRF/class/20190912.18//wind10m_countryfile_stept0_clipped.tif. 0..10...10....20..20.....3030......4040......5050.....60..60....70...70..80....80.90......90...100 - done. 100 - done. input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190912.18/20190912.18_wind_stept0.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190912.18// PopFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190912.18//wind10m_popfile_stept0_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190912.18//wind10m_countryfile_stept0_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190912.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 271.525017857 9.5083154762 299.541684524 37.5249821429 HWRF /usr/bin/gdalwarp -te -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190912.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 -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190912.18//wind10m_countryfile_stept0_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.00833333333334 cellsize 0.00833333333333 >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190912.18/wind_popDensValues_stept0.xml >> 7. remove files done input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.00/20190913.00_wind_stept0.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.00// PopFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.00//wind10m_popfile_stept0_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.00//wind10m_countryfile_stept0_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.00/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 271.525017857 9.5083154762 299.541684524 37.5249821429 HWRF /usr/bin/gdalwarp -te -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.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 -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.00//wind10m_countryfile_stept0_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.00833333333334 cellsize 0.00833333333333 >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.00/wind_popDensValues_stept0.xml >> 7. remove files done t0 completed Creating output file that is 3362P x 3362L. Processing input file /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.00/20190913.00_wind.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.00/20190913.00_wind.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.00/20190913.00_wind.tif to destination /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.00//wind10m_res_all.tif. 0.Creating output file that is 3362P x 3362L. Processing input file /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190912.18/20190912.18_wind.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190912.18/20190912.18_wind.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190912.18/20190912.18_wind.tif to destination /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190912.18//wind10m_res_all.tif. 0...10....10..20.....2030......3040......5040.....60..50....70..60...80....70..90.....80...90...100 - done. 100 - done. Creating output file that is 3362P x 3362L. 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/1000599/3_HWRF/class/20190913.00//wind10m_popfile_all_clipped.tif. 0Creating output file that is 3362P x 3362L. 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/1000599/3_HWRF/class/20190912.18//wind10m_popfile_all_clipped.tif. 0......1010......2020......3030.....40...40.50...50...60....70....60..80....70..90.....80...90...100 - done. 100 - done. Creating output file that is 3362P x 3362L. 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/1000599/3_HWRF/class/20190913.00//wind10m_countryfile_all_clipped.tif. 0Creating output file that is 3362P x 3362L. 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/1000599/3_HWRF/class/20190912.18//wind10m_countryfile_all_clipped.tif. 0......1010......2020......3030......4040......5050......6060......7070......8080......9090......100 - done. 100 - done. input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.00/20190913.00_wind.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.00// PopFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.00//wind10m_popfile_all_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.00//wind10m_countryfile_all_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.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 271.525017857 9.5083154762 299.541684524 37.5249821429 HWRF /usr/bin/gdalwarp -te -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.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 -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.00//wind10m_countryfile_all_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.00833333333334 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/1000599/3_HWRF/class/20190913.00/wind_popDensValues_all.xml >> 7. remove files done input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190912.18/20190912.18_wind.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190912.18// PopFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190912.18//wind10m_popfile_all_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190912.18//wind10m_countryfile_all_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190912.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 271.525017857 9.5083154762 299.541684524 37.5249821429 HWRF /usr/bin/gdalwarp -te -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190912.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 -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190912.18//wind10m_countryfile_all_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.00833333333334 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/1000599/3_HWRF/class/20190912.18/wind_popDensValues_all.xml >> 7. remove files done Creating output file that is 3362P x 3362L. Processing input file /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif. Using internal nodata values (e.g. 3.40282e+38) for image /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif to destination /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/final//wind10m_res_final.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 3362P x 3362L. 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/1000599/3_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 3362P x 3362L. 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/1000599/3_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/1000599/3_HWRF/class/20190912.18/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190912.18/wind_popDensValues_t0.xml copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190912.18/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190912.18/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190912.18/20190912.18_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190912.18/wind.jpg copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.00/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190913.00/wind_popDensValues_t0.xml copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.00/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190913.00/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.00/20190913.00_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190913.00/wind.jpg >> 3.2. Classify final folder /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/final/20190913.00_final_completed_wind.txt input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/final// PopFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/final//wind10m_popfile_final_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/final//wind10m_countryfile_final_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000599/3_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 271.525017857 9.5083154762 299.541684524 37.5249821429 HWRF /usr/bin/gdalwarp -te -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_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 -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/final//wind10m_countryfile_final_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.00833333333334 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/1000599/3_HWRF/class/final/wind_popDensValues_final.xml >> 7. remove files done copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/final/wind_popDensValues_final.xml in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/final/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_FINAL.jpg in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/final/wind.jpg ******************************************************* * Case Completed * ******************************************************* Case starterd at: 2019-09-13 10:24:13.741424 UTC Case completed at: 2019-09-13 10:27:15.535718 UTC **********E N D O F J O B*********************** ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-09-13 10:34:06.594705 UTC inp1= 20190912.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000599/09l submitting calc 2019-09-12 18:00:00 2019-09-13 00:00:00 delft3d -88.5 -60.48 9.55 37.55 4.0 HWRF 72 15 True GDACS/1000599/3_HWRF 6 1 False False 5 aa.txt 20190912.18 1000599/09l 6 wind False *************---------------------****************** GDACS/1000599/3_HWRF 1000599 ndt: 2 it: 0 ndt: 2 idate: 2019-09-12 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_6h.tif input file exist it: 1 ndt: 2 idate: 2019-09-13 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_6h.tif input file exist GDACS/1000599/3_HWRF 1000599 ndt: 2 it: 0 ndt: 2 idate: 2019-09-12 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_72h.tif input file exist it: 1 ndt: 2 idate: 2019-09-13 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_72h.tif input file exist home dir /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/ ret 0 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-09-12 18:00:00 currdate 2019-09-13 00:00:00 ndt: 6 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-09-12 18:00:00', '2019-09-13 00:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files processing curr bull... 2019-09-13 00:00:00 ********** no new data no new data ******************************************************* * Case Completed * ******************************************************* Case starterd at: 2019-09-13 10:34:06.594705 UTC Case completed at: 2019-09-13 10:34:07.573867 UTC **********E N D O F J O B*********************** ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-09-13 10:44:11.343939 UTC inp1= 20190912.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000599/09l submitting calc 2019-09-12 18:00:00 2019-09-13 00:00:00 delft3d -88.5 -60.48 9.55 37.55 4.0 HWRF 72 15 True GDACS/1000599/3_HWRF 6 1 False False 5 aa.txt 20190912.18 1000599/09l 6 wind False *************---------------------****************** GDACS/1000599/3_HWRF 1000599 ndt: 2 it: 0 ndt: 2 idate: 2019-09-12 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_6h.tif input file exist it: 1 ndt: 2 idate: 2019-09-13 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_6h.tif input file exist GDACS/1000599/3_HWRF 1000599 ndt: 2 it: 0 ndt: 2 idate: 2019-09-12 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_72h.tif input file exist it: 1 ndt: 2 idate: 2019-09-13 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_72h.tif input file exist home dir /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/ ret 0 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-09-12 18:00:00 currdate 2019-09-13 00:00:00 ndt: 6 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-09-12 18:00:00', '2019-09-13 00:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files processing curr bull... 2019-09-13 00:00:00 ********** no new data no new data ******************************************************* * Case Completed * ******************************************************* Case starterd at: 2019-09-13 10:44:11.343939 UTC Case completed at: 2019-09-13 10:44:11.564845 UTC **********E N D O F J O B*********************** ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-09-13 10:54:10.955297 UTC inp1= 20190912.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000599/09l submitting calc 2019-09-12 18:00:00 2019-09-13 00:00:00 delft3d -88.5 -60.48 9.55 37.55 4.0 HWRF 72 15 True GDACS/1000599/3_HWRF 6 1 False False 5 aa.txt 20190912.18 1000599/09l 6 wind False *************---------------------****************** GDACS/1000599/3_HWRF 1000599 ndt: 2 it: 0 ndt: 2 idate: 2019-09-12 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_6h.tif input file exist it: 1 ndt: 2 idate: 2019-09-13 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_6h.tif input file exist GDACS/1000599/3_HWRF 1000599 ndt: 2 it: 0 ndt: 2 idate: 2019-09-12 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_72h.tif input file exist it: 1 ndt: 2 idate: 2019-09-13 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_72h.tif input file exist home dir /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/ ret 0 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-09-12 18:00:00 currdate 2019-09-13 00:00:00 ndt: 6 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-09-12 18:00:00', '2019-09-13 00:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files processing curr bull... 2019-09-13 00:00:00 ********** no new data no new data ******************************************************* * Case Completed * ******************************************************* Case starterd at: 2019-09-13 10:54:10.955297 UTC Case completed at: 2019-09-13 10:54:11.109510 UTC **********E N D O F J O B*********************** Creating output file that is 560P x 560L. Processing input file /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091312.swath.wind_6h.tif. Using internal nodata values (e.g. -3.40282e+38) for image /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091312.swath.wind_6h.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 560P x 560L. Processing input file /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091312.swath.wind_72h.tif. Using internal nodata values (e.g. -3.40282e+38) for image /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091312.swath.wind_72h.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 560P x 560L. Processing input file /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091306.swath.wind_6h.tif. Using internal nodata values (e.g. -3.40282e+38) for image /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091306.swath.wind_6h.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 560P x 560L. Processing input file /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091306.swath.wind_72h.tif. Using internal nodata values (e.g. -3.40282e+38) for image /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091306.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-09-13 22:34:05.601023 UTC inp1= 20190912.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000599/09l submitting calc 2019-09-12 18:00:00 2019-09-13 12:00:00 delft3d -88.5 -60.48 9.55 37.55 4.0 HWRF 72 15 True GDACS/1000599/3_HWRF 6 1 False False 5 aa.txt 20190912.18 1000599/09l 6 wind False *************---------------------****************** GDACS/1000599/3_HWRF 1000599 ndt: 4 it: 0 ndt: 4 idate: 2019-09-12 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_6h.tif input file exist it: 1 ndt: 4 idate: 2019-09-13 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_6h.tif input file exist it: 2 ndt: 4 idate: 2019-09-13 06:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091306.swath.wind_6h.tif input file exist it: 3 ndt: 4 idate: 2019-09-13 12:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091312.swath.wind_6h.tif input file exist GDACS/1000599/3_HWRF 1000599 ndt: 4 it: 0 ndt: 4 idate: 2019-09-12 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_72h.tif input file exist it: 1 ndt: 4 idate: 2019-09-13 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_72h.tif input file exist it: 2 ndt: 4 idate: 2019-09-13 06:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091306.swath.wind_72h.tif input file exist it: 3 ndt: 4 idate: 2019-09-13 12:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091312.swath.wind_72h.tif input file exist home dir /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/ ret 0 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-09-12 18:00:00 currdate 2019-09-13 12:00:00 ndt: 18 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-09-12 18:00:00', '2019-09-13 00:00:00', '2019-09-13 06:00:00', '2019-09-13 12:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files processing curr bull... 2019-09-13 12:00:00 -88.5 -60.48 271.5 299.52 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.12/ /mnt/input/grib_HWRF/2019/1000599/ input file /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091312.swath.wind_6h.tif /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.12/20190913.12_wind_stept0.tif -r bilinear -te 271.5 9.55 299.52 37.55 -dstnodata 0 Resample tif file /usr/bin/gdalwarp -r bilinear -te 271.5 9.55 299.52 37.55 -dstnodata 0 "/mnt/input/grib_HWRF/2019/1000599/nine09l.2019091312.swath.wind_6h.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.12/20190913.12_wind_stept0.tif" outFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.12/20190913.12_wind_stept0.tif dd: 0.0500357142857 >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091312/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091312/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091312/all_inpData.txt 271.525017857 299.494982143 9.55501785714 37.5249821429 lon 560 lat 560 read TC track /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091312/all_inpData.txt currDate 2019-09-13 12:00:00 2019-09-16 12:00:00 sdate 2019-09-16 12:00:00 -4 CAT ['TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'H', 'H', 'H', 'H', 'H', 'H', 'H'] savemap /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.12/20190913.12_wind_stept0.jpg -88.5 -60.48 271.5 299.52 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.12/ /mnt/input/grib_HWRF/2019/1000599/ input file /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091312.swath.wind_72h.tif /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.12/20190913.12_wind.tif -r bilinear -te 271.5 9.55 299.52 37.55 -dstnodata 0 Resample tif file /usr/bin/gdalwarp -r bilinear -te 271.5 9.55 299.52 37.55 -dstnodata 0 "/mnt/input/grib_HWRF/2019/1000599/nine09l.2019091312.swath.wind_72h.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.12/20190913.12_wind.tif" outFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.12/20190913.12_wind.tif dd: 0.0500357142857 >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091312/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091312/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091312/all_inpData.txt 271.525017857 299.494982143 9.55501785714 37.5249821429 lon 560 lat 560 read TC track /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091312/all_inpData.txt currDate 2019-09-13 12:00:00 2019-09-16 12:00:00 sdate 2019-09-16 12:00:00 -4 CAT ['TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'H', 'H', 'H', 'H', 'H', 'H', 'H'] savemap /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.12/20190913.12_wind.jpg processing all past bull only if Past=True... True itdate, istime 2019-09-12 18:00:00 20190912.18 meteo-processing past forecast already completed itdate, istime 2019-09-13 00:00:00 20190913.00 meteo-processing past forecast already completed itdate, istime 2019-09-13 06:00:00 20190913.06 -88.5 -60.48 271.5 299.52 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.06/ /mnt/input/grib_HWRF/2019/1000599/ input file /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091306.swath.wind_6h.tif /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.06/20190913.06_wind_stept0.tif -r bilinear -te 271.5 9.55 299.52 37.55 -dstnodata 0 Resample tif file /usr/bin/gdalwarp -r bilinear -te 271.5 9.55 299.52 37.55 -dstnodata 0 "/mnt/input/grib_HWRF/2019/1000599/nine09l.2019091306.swath.wind_6h.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.06/20190913.06_wind_stept0.tif" outFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.06/20190913.06_wind_stept0.tif dd: 0.0500357142857 >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091306/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091306/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091306/all_inpData.txt 271.525017857 299.494982143 9.55501785714 37.5249821429 lon 560 lat 560 read TC track /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091306/all_inpData.txt currDate 2019-09-13 06:00:00 2019-09-16 06:00:00 sdate 2019-09-16 06:00:00 -4 CAT ['TD', 'TD', 'TD', 'TS', 'TD', 'TS', 'TS', 'TD', 'TS', 'TS', 'TS', 'TS', 'H', 'TS', 'H', 'MH', 'H', 'H', 'H', 'H', 'H'] savemap /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.06/20190913.06_wind_stept0.jpg -88.5 -60.48 271.5 299.52 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.06/ /mnt/input/grib_HWRF/2019/1000599/ input file /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091306.swath.wind_72h.tif /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.06/20190913.06_wind.tif -r bilinear -te 271.5 9.55 299.52 37.55 -dstnodata 0 Resample tif file /usr/bin/gdalwarp -r bilinear -te 271.5 9.55 299.52 37.55 -dstnodata 0 "/mnt/input/grib_HWRF/2019/1000599/nine09l.2019091306.swath.wind_72h.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.06/20190913.06_wind.tif" outFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.06/20190913.06_wind.tif dd: 0.0500357142857 >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091306/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091306/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091306/all_inpData.txt 271.525017857 299.494982143 9.55501785714 37.5249821429 lon 560 lat 560 read TC track /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091306/all_inpData.txt currDate 2019-09-13 06:00:00 2019-09-16 06:00:00 sdate 2019-09-16 06:00:00 -4 CAT ['TD', 'TD', 'TD', 'TS', 'TD', 'TS', 'TS', 'TD', 'TS', 'TS', 'TS', 'TS', 'H', 'TS', 'H', 'MH', 'H', 'H', 'H', 'H', 'H'] savemap /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.06/20190913.06_wind.jpg >> 2.2. Process final data files compfile /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/20190913.12_Final_completed_wind.txt FINAL remove maxtif alldate: DatetimeIndex(['2019-09-12 18:00:00', '2019-09-13 00:00:00', '2019-09-13 06:00:00', '2019-09-13 12:00:00'], dtype='datetime64[ns]', freq='6H') alldate.shape 4 date: 2019-09-13 00:00:00 python /mnt/output/SSCS/scripts/gdal_calc.py -A /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190912.18/20190912.18_wind_stept0.tif -B /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.00/20190913.00_wind_stept0.tif --outfile=/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif --calc="maximum(A,B)" date: 2019-09-13 06:00:00 python /mnt/output/SSCS/scripts/gdal_calc.py -A /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif -B /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.06/20190913.06_wind_stept0.tif --outfile=/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif --calc="maximum(A,B)" date: 2019-09-13 12:00:00 python /mnt/output/SSCS/scripts/gdal_calc.py -A /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif -B /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.12/20190913.12_wind.tif --outfile=/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif --calc="maximum(A,B)" max file created /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif stormname 1000599/09l stormname 1000599/09l TCbullsource 09l >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091312/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091312/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091312/all_inpData.txt 271.550035714 299.52 9.52999999999 37.4999642857 lon 560 lat 560 read TC track /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091312/all_inpData.txt currDate 2019-09-13 12:00:00 2019-09-16 12:00:00 sdate 2019-09-16 12:00:00 -4 CAT ['TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'H', 'H', 'H', 'H', 'H', 'H', 'H'] savemap /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_FINAL.jpg ret: 0 ============================================ 3. Classify meteo + GDACS index score ============================================ >> 3.1. Classify curr + past forecast input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.12/20190913.12_wind_stept0.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.12// PopFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.12//wind10m_popfile_stept0_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.12//wind10m_countryfile_stept0_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.12/wind_popDensValues_stept0.xml input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.06/20190913.06_wind_stept0.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.06// PopFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.06//wind10m_popfile_stept0_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.06//wind10m_countryfile_stept0_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.06/wind_popDensValues_stept0.xml t0 completed Creating output file that is 3362P x 3362L. Processing input file /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.06/20190913.06_wind.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.06/20190913.06_wind.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.06/20190913.06_wind.tif to destination /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.06//wind10m_res_all.tif. 0Creating output file that is 3362P x 3362L. Processing input file /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.12/20190913.12_wind.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.12/20190913.12_wind.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.12/20190913.12_wind.tif to destination /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.12//wind10m_res_all.tif. 0......1010......2020......3030......4040......5050......6060......7070......8080.....90..90.....100 - done. 100 - done. Creating output file that is 3362P x 3362L. 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/1000599/3_HWRF/class/20190913.06//wind10m_popfile_all_clipped.tif. 0Creating output file that is 3362P x 3362L. 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/1000599/3_HWRF/class/20190913.12//wind10m_popfile_all_clipped.tif. 0......1010......2020......30.30....40..40....50.50.....60..60....70..70....80..80....90..90.....100 - done. 100 - done. Creating output file that is 3362P x 3362L. 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/1000599/3_HWRF/class/20190913.12//wind10m_countryfile_all_clipped.tif. 0Creating output file that is 3362P x 3362L. 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/1000599/3_HWRF/class/20190913.06//wind10m_countryfile_all_clipped.tif. 0......1010......2020......3030......4040......5050......6060......7070......8080......9090......100 - done. 100 - done. input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.06/20190913.06_wind.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.06// PopFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.06//wind10m_popfile_all_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.06//wind10m_countryfile_all_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.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 271.525017857 9.5083154762 299.541684524 37.5249821429 HWRF /usr/bin/gdalwarp -te -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.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 -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.06//wind10m_countryfile_all_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.00833333333334 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/1000599/3_HWRF/class/20190913.06/wind_popDensValues_all.xml >> 7. remove files done input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.12/20190913.12_wind.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.12// PopFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.12//wind10m_popfile_all_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.12//wind10m_countryfile_all_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.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 271.525017857 9.5083154762 299.541684524 37.5249821429 HWRF /usr/bin/gdalwarp -te -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.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 -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.12//wind10m_countryfile_all_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.00833333333334 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/1000599/3_HWRF/class/20190913.12/wind_popDensValues_all.xml >> 7. remove files done Creating output file that is 3362P x 3362L. Processing input file /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif. Using internal nodata values (e.g. 3.40282e+38) for image /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif to destination /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/final//wind10m_res_final.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 3362P x 3362L. 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/1000599/3_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 3362P x 3362L. 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/1000599/3_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/1000599/3_HWRF/class/20190912.18/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190912.18/wind_popDensValues_t0.xml copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190912.18/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190912.18/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190912.18/20190912.18_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190912.18/wind.jpg copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.00/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190913.00/wind_popDensValues_t0.xml copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.00/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190913.00/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.00/20190913.00_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190913.00/wind.jpg copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.06/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190913.06/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.06/20190913.06_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190913.06/wind.jpg copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.12/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190913.12/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.12/20190913.12_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190913.12/wind.jpg >> 3.2. Classify final folder /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/final/20190913.12_final_completed_wind.txt input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/final// PopFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/final//wind10m_popfile_final_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/final//wind10m_countryfile_final_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000599/3_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 271.525017857 9.5083154762 299.541684524 37.5249821429 HWRF /usr/bin/gdalwarp -te -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_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 -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/final//wind10m_countryfile_final_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.00833333333334 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/1000599/3_HWRF/class/final/wind_popDensValues_final.xml >> 7. remove files done copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/final/wind_popDensValues_final.xml in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/final/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_FINAL.jpg in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/final/wind.jpg ******************************************************* * Case Completed * ******************************************************* Case starterd at: 2019-09-13 22:34:05.601023 UTC Case completed at: 2019-09-13 22:36:29.063934 UTC **********E N D O F J O B*********************** ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-09-13 22:44:05.840667 UTC inp1= 20190912.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000599/09l submitting calc 2019-09-12 18:00:00 2019-09-13 12:00:00 delft3d -88.5 -60.48 9.55 37.55 4.0 HWRF 72 15 True GDACS/1000599/3_HWRF 6 1 False False 5 aa.txt 20190912.18 1000599/09l 6 wind False *************---------------------****************** GDACS/1000599/3_HWRF 1000599 ndt: 4 it: 0 ndt: 4 idate: 2019-09-12 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_6h.tif input file exist it: 1 ndt: 4 idate: 2019-09-13 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_6h.tif input file exist it: 2 ndt: 4 idate: 2019-09-13 06:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091306.swath.wind_6h.tif input file exist it: 3 ndt: 4 idate: 2019-09-13 12:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091312.swath.wind_6h.tif input file exist GDACS/1000599/3_HWRF 1000599 ndt: 4 it: 0 ndt: 4 idate: 2019-09-12 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_72h.tif input file exist it: 1 ndt: 4 idate: 2019-09-13 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_72h.tif input file exist it: 2 ndt: 4 idate: 2019-09-13 06:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091306.swath.wind_72h.tif input file exist it: 3 ndt: 4 idate: 2019-09-13 12:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091312.swath.wind_72h.tif input file exist home dir /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/ ret 0 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-09-12 18:00:00 currdate 2019-09-13 12:00:00 ndt: 18 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-09-12 18:00:00', '2019-09-13 00:00:00', '2019-09-13 06:00:00', '2019-09-13 12:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files processing curr bull... 2019-09-13 12:00:00 ********** no new data no new data ******************************************************* * Case Completed * ******************************************************* Case starterd at: 2019-09-13 22:44:05.840667 UTC Case completed at: 2019-09-13 22:44:06.079228 UTC **********E N D O F J O B*********************** ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-09-13 22:54:11.290050 UTC inp1= 20190912.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000599/09l submitting calc 2019-09-12 18:00:00 2019-09-13 12:00:00 delft3d -88.5 -60.48 9.55 37.55 4.0 HWRF 72 15 True GDACS/1000599/3_HWRF 6 1 False False 5 aa.txt 20190912.18 1000599/09l 6 wind False *************---------------------****************** GDACS/1000599/3_HWRF 1000599 ndt: 4 it: 0 ndt: 4 idate: 2019-09-12 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_6h.tif input file exist it: 1 ndt: 4 idate: 2019-09-13 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_6h.tif input file exist it: 2 ndt: 4 idate: 2019-09-13 06:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091306.swath.wind_6h.tif input file exist it: 3 ndt: 4 idate: 2019-09-13 12:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091312.swath.wind_6h.tif input file exist GDACS/1000599/3_HWRF 1000599 ndt: 4 it: 0 ndt: 4 idate: 2019-09-12 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_72h.tif input file exist it: 1 ndt: 4 idate: 2019-09-13 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_72h.tif input file exist it: 2 ndt: 4 idate: 2019-09-13 06:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091306.swath.wind_72h.tif input file exist it: 3 ndt: 4 idate: 2019-09-13 12:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091312.swath.wind_72h.tif input file exist home dir /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/ ret 0 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-09-12 18:00:00 currdate 2019-09-13 12:00:00 ndt: 18 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-09-12 18:00:00', '2019-09-13 00:00:00', '2019-09-13 06:00:00', '2019-09-13 12:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files processing curr bull... 2019-09-13 12:00:00 ********** no new data no new data ******************************************************* * Case Completed * ******************************************************* Case starterd at: 2019-09-13 22:54:11.290050 UTC Case completed at: 2019-09-13 22:54:11.527487 UTC **********E N D O F J O B*********************** ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-09-13 23:04:06.857052 UTC inp1= 20190912.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000599/09l submitting calc 2019-09-12 18:00:00 2019-09-13 12:00:00 delft3d -88.5 -60.48 9.55 37.55 4.0 HWRF 72 15 True GDACS/1000599/3_HWRF 6 1 False False 5 aa.txt 20190912.18 1000599/09l 6 wind False *************---------------------****************** GDACS/1000599/3_HWRF 1000599 ndt: 4 it: 0 ndt: 4 idate: 2019-09-12 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_6h.tif input file exist it: 1 ndt: 4 idate: 2019-09-13 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_6h.tif input file exist it: 2 ndt: 4 idate: 2019-09-13 06:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091306.swath.wind_6h.tif input file exist it: 3 ndt: 4 idate: 2019-09-13 12:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091312.swath.wind_6h.tif input file exist GDACS/1000599/3_HWRF 1000599 ndt: 4 it: 0 ndt: 4 idate: 2019-09-12 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_72h.tif input file exist it: 1 ndt: 4 idate: 2019-09-13 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_72h.tif input file exist it: 2 ndt: 4 idate: 2019-09-13 06:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091306.swath.wind_72h.tif input file exist it: 3 ndt: 4 idate: 2019-09-13 12:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091312.swath.wind_72h.tif input file exist home dir /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/ ret 0 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-09-12 18:00:00 currdate 2019-09-13 12:00:00 ndt: 18 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-09-12 18:00:00', '2019-09-13 00:00:00', '2019-09-13 06:00:00', '2019-09-13 12:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files processing curr bull... 2019-09-13 12:00:00 ********** no new data no new data ******************************************************* * Case Completed * ******************************************************* Case starterd at: 2019-09-13 23:04:06.857052 UTC Case completed at: 2019-09-13 23:04:07.072523 UTC **********E N D O F J O B*********************** Creating output file that is 560P x 560L. Processing input file /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091400.swath.wind_6h.tif. Using internal nodata values (e.g. -3.40282e+38) for image /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091400.swath.wind_6h.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 560P x 560L. Processing input file /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091400.swath.wind_72h.tif. Using internal nodata values (e.g. -3.40282e+38) for image /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091400.swath.wind_72h.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 560P x 560L. Processing input file /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091318.swath.wind_6h.tif. Using internal nodata values (e.g. -3.40282e+38) for image /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091318.swath.wind_6h.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 560P x 560L. Processing input file /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091318.swath.wind_72h.tif. Using internal nodata values (e.g. -3.40282e+38) for image /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091318.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-09-14 11:14:05.150286 UTC inp1= 20190912.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000599/09l submitting calc 2019-09-12 18:00:00 2019-09-14 00:00:00 delft3d -88.5 -60.48 9.55 37.55 4.0 HWRF 72 15 True GDACS/1000599/3_HWRF 6 1 False False 5 aa.txt 20190912.18 1000599/09l 6 wind False *************---------------------****************** GDACS/1000599/3_HWRF 1000599 ndt: 6 it: 0 ndt: 6 idate: 2019-09-12 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_6h.tif input file exist it: 1 ndt: 6 idate: 2019-09-13 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_6h.tif input file exist it: 2 ndt: 6 idate: 2019-09-13 06:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091306.swath.wind_6h.tif input file exist it: 3 ndt: 6 idate: 2019-09-13 12:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091312.swath.wind_6h.tif input file exist it: 4 ndt: 6 idate: 2019-09-13 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091318.swath.wind_6h.tif input file exist it: 5 ndt: 6 idate: 2019-09-14 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091400.swath.wind_6h.tif input file exist GDACS/1000599/3_HWRF 1000599 ndt: 6 it: 0 ndt: 6 idate: 2019-09-12 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_72h.tif input file exist it: 1 ndt: 6 idate: 2019-09-13 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_72h.tif input file exist it: 2 ndt: 6 idate: 2019-09-13 06:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091306.swath.wind_72h.tif input file exist it: 3 ndt: 6 idate: 2019-09-13 12:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091312.swath.wind_72h.tif input file exist it: 4 ndt: 6 idate: 2019-09-13 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091318.swath.wind_72h.tif input file exist it: 5 ndt: 6 idate: 2019-09-14 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091400.swath.wind_72h.tif input file exist home dir /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/ ret 0 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-09-12 18:00:00 currdate 2019-09-14 00:00:00 ndt: 30 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-09-12 18:00:00', '2019-09-13 00:00:00', '2019-09-13 06:00:00', '2019-09-13 12:00:00', '2019-09-13 18:00:00', '2019-09-14 00:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files processing curr bull... 2019-09-14 00:00:00 -88.5 -60.48 271.5 299.52 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190914.00/ /mnt/input/grib_HWRF/2019/1000599/ input file /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091400.swath.wind_6h.tif /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190914.00/20190914.00_wind_stept0.tif -r bilinear -te 271.5 9.55 299.52 37.55 -dstnodata 0 Resample tif file /usr/bin/gdalwarp -r bilinear -te 271.5 9.55 299.52 37.55 -dstnodata 0 "/mnt/input/grib_HWRF/2019/1000599/nine09l.2019091400.swath.wind_6h.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190914.00/20190914.00_wind_stept0.tif" outFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190914.00/20190914.00_wind_stept0.tif dd: 0.0500357142857 >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091400/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091400/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091400/all_inpData.txt 271.525017857 299.494982143 9.55501785714 37.5249821429 lon 560 lat 560 read TC track /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091400/all_inpData.txt currDate 2019-09-14 00:00:00 2019-09-17 00:00:00 sdate 2019-09-17 00:00:00 -4 CAT ['TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'H', 'H', 'H', 'H', 'H', 'H', 'H', 'H', 'MH', 'H', 'MH'] savemap /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190914.00/20190914.00_wind_stept0.jpg -88.5 -60.48 271.5 299.52 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190914.00/ /mnt/input/grib_HWRF/2019/1000599/ input file /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091400.swath.wind_72h.tif /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190914.00/20190914.00_wind.tif -r bilinear -te 271.5 9.55 299.52 37.55 -dstnodata 0 Resample tif file /usr/bin/gdalwarp -r bilinear -te 271.5 9.55 299.52 37.55 -dstnodata 0 "/mnt/input/grib_HWRF/2019/1000599/nine09l.2019091400.swath.wind_72h.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190914.00/20190914.00_wind.tif" outFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190914.00/20190914.00_wind.tif dd: 0.0500357142857 >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091400/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091400/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091400/all_inpData.txt 271.525017857 299.494982143 9.55501785714 37.5249821429 lon 560 lat 560 read TC track /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091400/all_inpData.txt currDate 2019-09-14 00:00:00 2019-09-17 00:00:00 sdate 2019-09-17 00:00:00 -4 CAT ['TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'H', 'H', 'H', 'H', 'H', 'H', 'H', 'H', 'MH', 'H', 'MH'] savemap /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190914.00/20190914.00_wind.jpg processing all past bull only if Past=True... True itdate, istime 2019-09-12 18:00:00 20190912.18 meteo-processing past forecast already completed itdate, istime 2019-09-13 00:00:00 20190913.00 meteo-processing past forecast already completed itdate, istime 2019-09-13 06:00:00 20190913.06 meteo-processing past forecast already completed itdate, istime 2019-09-13 12:00:00 20190913.12 meteo-processing past forecast already completed itdate, istime 2019-09-13 18:00:00 20190913.18 -88.5 -60.48 271.5 299.52 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.18/ /mnt/input/grib_HWRF/2019/1000599/ input file /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091318.swath.wind_6h.tif /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.18/20190913.18_wind_stept0.tif -r bilinear -te 271.5 9.55 299.52 37.55 -dstnodata 0 Resample tif file /usr/bin/gdalwarp -r bilinear -te 271.5 9.55 299.52 37.55 -dstnodata 0 "/mnt/input/grib_HWRF/2019/1000599/nine09l.2019091318.swath.wind_6h.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.18/20190913.18_wind_stept0.tif" outFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.18/20190913.18_wind_stept0.tif dd: 0.0500357142857 >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091318/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091318/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091318/all_inpData.txt 271.525017857 299.494982143 9.55501785714 37.5249821429 lon 560 lat 560 read TC track /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091318/all_inpData.txt currDate 2019-09-13 18:00:00 2019-09-16 18:00:00 sdate 2019-09-16 18:00:00 -4 CAT ['TD', 'TD', 'TD', 'TD', 'TD', 'TS', 'TS', 'TS', 'TS', 'TS', 'H', 'TS', 'TS', 'H', 'H', 'H', 'H', 'H', 'H', 'MH', 'MH', 'H', 'H'] savemap /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.18/20190913.18_wind_stept0.jpg -88.5 -60.48 271.5 299.52 0.05 0.05 /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.18/ /mnt/input/grib_HWRF/2019/1000599/ input file /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091318.swath.wind_72h.tif /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.18/20190913.18_wind.tif -r bilinear -te 271.5 9.55 299.52 37.55 -dstnodata 0 Resample tif file /usr/bin/gdalwarp -r bilinear -te 271.5 9.55 299.52 37.55 -dstnodata 0 "/mnt/input/grib_HWRF/2019/1000599/nine09l.2019091318.swath.wind_72h.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.18/20190913.18_wind.tif" outFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.18/20190913.18_wind.tif dd: 0.0500357142857 >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091318/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091318/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091318/all_inpData.txt 271.525017857 299.494982143 9.55501785714 37.5249821429 lon 560 lat 560 read TC track /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091318/all_inpData.txt currDate 2019-09-13 18:00:00 2019-09-16 18:00:00 sdate 2019-09-16 18:00:00 -4 CAT ['TD', 'TD', 'TD', 'TD', 'TD', 'TS', 'TS', 'TS', 'TS', 'TS', 'H', 'TS', 'TS', 'H', 'H', 'H', 'H', 'H', 'H', 'MH', 'MH', 'H', 'H'] savemap /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.18/20190913.18_wind.jpg >> 2.2. Process final data files compfile /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/20190914.00_Final_completed_wind.txt FINAL remove maxtif alldate: DatetimeIndex(['2019-09-12 18:00:00', '2019-09-13 00:00:00', '2019-09-13 06:00:00', '2019-09-13 12:00:00', '2019-09-13 18:00:00', '2019-09-14 00:00:00'], dtype='datetime64[ns]', freq='6H') alldate.shape 6 date: 2019-09-13 00:00:00 python /mnt/output/SSCS/scripts/gdal_calc.py -A /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190912.18/20190912.18_wind_stept0.tif -B /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.00/20190913.00_wind_stept0.tif --outfile=/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif --calc="maximum(A,B)" date: 2019-09-13 06:00:00 python /mnt/output/SSCS/scripts/gdal_calc.py -A /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif -B /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.06/20190913.06_wind_stept0.tif --outfile=/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif --calc="maximum(A,B)" date: 2019-09-13 12:00:00 python /mnt/output/SSCS/scripts/gdal_calc.py -A /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif -B /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.12/20190913.12_wind_stept0.tif --outfile=/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif --calc="maximum(A,B)" date: 2019-09-13 18:00:00 python /mnt/output/SSCS/scripts/gdal_calc.py -A /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif -B /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.18/20190913.18_wind_stept0.tif --outfile=/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif --calc="maximum(A,B)" date: 2019-09-14 00:00:00 python /mnt/output/SSCS/scripts/gdal_calc.py -A /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif -B /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190914.00/20190914.00_wind.tif --outfile=/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif --calc="maximum(A,B)" max file created /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif stormname 1000599/09l stormname 1000599/09l TCbullsource 09l >> 0. SET INPUT/OUTPUT HWRF /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091400/all_inpData.txt /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091400/all_inpData.txt trackfile /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091400/all_inpData.txt 271.550035714 299.52 9.52999999999 37.4999642857 lon 560 lat 560 read TC track /mnt/output/GDACS/TC/TCTrack/1000599/HWRF/2019091400/all_inpData.txt currDate 2019-09-14 00:00:00 2019-09-17 00:00:00 sdate 2019-09-17 00:00:00 -4 CAT ['TD', 'TD', 'TD', 'TD', 'TD', 'TD', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'TS', 'H', 'H', 'H', 'H', 'H', 'H', 'H', 'H', 'MH', 'H', 'MH'] savemap /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_FINAL.jpg ret: 0 ============================================ 3. Classify meteo + GDACS index score ============================================ >> 3.1. Classify curr + past forecast input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190914.00/20190914.00_wind_stept0.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190914.00// PopFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190914.00//wind10m_popfile_stept0_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190914.00//wind10m_countryfile_stept0_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190914.00/wind_popDensValues_stept0.xml Creating output file that is 3362P x 3362L. Processing input file /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.12/20190913.12_wind_stept0.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.12/20190913.12_wind_stept0.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.12/20190913.12_wind_stept0.tif to destination /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.12//wind10m_res_stept0.tif. 0.Creating output file that is 3362P x 3362L. Processing input file /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.06/20190913.06_wind_stept0.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.06/20190913.06_wind_stept0.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.06/20190913.06_wind_stept0.tif to destination /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.06//wind10m_res_stept0.tif. 0Creating output file that is 3362P x 3362L. Processing input file /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.18/20190913.18_wind_stept0.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.18/20190913.18_wind_stept0.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.18/20190913.18_wind_stept0.tif to destination /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.18//wind10m_res_stept0.tif. 0......10...1010......20...2020......30...3030......40...4040......50...5050......60...6060......70...7070......80...8080......90...9090........100 - done. 100 - done. 100 - done. Creating output file that is 3362P x 3362L. 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/1000599/3_HWRF/class/20190913.18//wind10m_popfile_stept0_clipped.tif. 0Creating output file that is 3362P x 3362L. 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/1000599/3_HWRF/class/20190913.06//wind10m_popfile_stept0_clipped.tif. 0Creating output file that is 3362P x 3362L. 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/1000599/3_HWRF/class/20190913.12//wind10m_popfile_stept0_clipped.tif. 0.........101010.........202020.........3030.30.......40..4040......50..5050.........606060.........7070.70.......80.80.80.......90.90.90........100 - done. 100 - done. 100 - done. Creating output file that is 3362P x 3362L. 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/1000599/3_HWRF/class/20190913.06//wind10m_countryfile_stept0_clipped.tif. 0Creating output file that is 3362P x 3362L. 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/1000599/3_HWRF/class/20190913.12//wind10m_countryfile_stept0_clipped.tif. 0..Creating output file that is 3362P x 3362L. 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/1000599/3_HWRF/class/20190913.18//wind10m_countryfile_stept0_clipped.tif. 0......1010...10......2020...20......3030...30......4040...40......5050...50......60.60.60......70...7070......80...8080.....90....90.90.......100 - done. 100 - done. 100 - done. input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.12/20190913.12_wind_stept0.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.12// PopFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.12//wind10m_popfile_stept0_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.12//wind10m_countryfile_stept0_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.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 271.525017857 9.5083154762 299.541684524 37.5249821429 HWRF /usr/bin/gdalwarp -te -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.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 -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.12//wind10m_countryfile_stept0_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.00833333333334 cellsize 0.00833333333333 >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.12/wind_popDensValues_stept0.xml >> 7. remove files done input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.06/20190913.06_wind_stept0.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.06// PopFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.06//wind10m_popfile_stept0_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.06//wind10m_countryfile_stept0_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.06/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 271.525017857 9.5083154762 299.541684524 37.5249821429 HWRF /usr/bin/gdalwarp -te -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.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 -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.06//wind10m_countryfile_stept0_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.00833333333334 cellsize 0.00833333333333 >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.06/wind_popDensValues_stept0.xml >> 7. remove files done input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.18/20190913.18_wind_stept0.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.18// PopFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.18//wind10m_popfile_stept0_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.18//wind10m_countryfile_stept0_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.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 271.525017857 9.5083154762 299.541684524 37.5249821429 HWRF /usr/bin/gdalwarp -te -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.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 -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.18//wind10m_countryfile_stept0_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.00833333333334 cellsize 0.00833333333333 >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.18/wind_popDensValues_stept0.xml >> 7. remove files done t0 completed Creating output file that is 3362P x 3362L. Processing input file /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190914.00/20190914.00_wind.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190914.00/20190914.00_wind.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190914.00/20190914.00_wind.tif to destination /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190914.00//wind10m_res_all.tif. 0Creating output file that is 3362P x 3362L. Processing input file /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.18/20190913.18_wind.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.18/20190913.18_wind.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.18/20190913.18_wind.tif to destination /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.18//wind10m_res_all.tif. 0......1010......2020......3030......4040......5050......6060......7070......8080......9090......100 - done. 100 - done. Creating output file that is 3362P x 3362L. 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/1000599/3_HWRF/class/20190914.00//wind10m_popfile_all_clipped.tif. 0Creating output file that is 3362P x 3362L. 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/1000599/3_HWRF/class/20190913.18//wind10m_popfile_all_clipped.tif. 0......1010......2020......3030......4040......5050......6060......7070......8080......9090......100 - done. 100 - done. Creating output file that is 3362P x 3362L. 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/1000599/3_HWRF/class/20190914.00//wind10m_countryfile_all_clipped.tif. 0Creating output file that is 3362P x 3362L. 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/1000599/3_HWRF/class/20190913.18//wind10m_countryfile_all_clipped.tif. 0......1010......2020......3030......4040......5050......6060......7070......8080......9090......100 - done. 100 - done. input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190914.00/20190914.00_wind.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190914.00// PopFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190914.00//wind10m_popfile_all_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190914.00//wind10m_countryfile_all_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190914.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 271.525017857 9.5083154762 299.541684524 37.5249821429 HWRF /usr/bin/gdalwarp -te -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190914.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 -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190914.00//wind10m_countryfile_all_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.00833333333334 cellsize 0.00833333333333 >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190914.00/wind_popDensValues_all.xml >> 7. remove files done input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.18/20190913.18_wind.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.18// PopFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.18//wind10m_popfile_all_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.18//wind10m_countryfile_all_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.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 271.525017857 9.5083154762 299.541684524 37.5249821429 HWRF /usr/bin/gdalwarp -te -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.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 -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.18//wind10m_countryfile_all_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.00833333333334 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/1000599/3_HWRF/class/20190913.18/wind_popDensValues_all.xml >> 7. remove files done Creating output file that is 3362P x 3362L. Processing input file /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif. Using internal nodata values (e.g. 3.40282e+38) for image /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif to destination /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/final//wind10m_res_final.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 3362P x 3362L. 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/1000599/3_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 3362P x 3362L. 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/1000599/3_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/1000599/3_HWRF/class/20190912.18/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190912.18/wind_popDensValues_t0.xml copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190912.18/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190912.18/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190912.18/20190912.18_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190912.18/wind.jpg copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.00/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190913.00/wind_popDensValues_t0.xml copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.00/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190913.00/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.00/20190913.00_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190913.00/wind.jpg copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.06/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190913.06/wind_popDensValues_t0.xml copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.06/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190913.06/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.06/20190913.06_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190913.06/wind.jpg copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.12/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190913.12/wind_popDensValues_t0.xml copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.12/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190913.12/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.12/20190913.12_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190913.12/wind.jpg copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.18/wind_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190913.18/wind_popDensValues_t0.xml copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190913.18/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190913.18/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190913.18/20190913.18_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190913.18/wind.jpg copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/20190914.00/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190914.00/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/20190914.00/20190914.00_wind.jpg in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/20190914.00/wind.jpg >> 3.2. Classify final folder /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/final/20190914.00_final_completed_wind.txt input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_final.tif hurName: hdate: var: wind10m description: wind10m: _ (HWRF) OutDir: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/final// PopFile: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/final//wind10m_popfile_final_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/final//wind10m_countryfile_final_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000599/3_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 271.525017857 9.5083154762 299.541684524 37.5249821429 HWRF /usr/bin/gdalwarp -te -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_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 -88.4749821429 9.5083154762 -60.4583154762 37.5249821429 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/final//wind10m_countryfile_final_clipped.tif" >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.00833333333334 cellsize 0.00833333333333 >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/final/wind_popDensValues_final.xml >> 7. remove files done copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/class/final/wind_popDensValues_final.xml in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/final/wind_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/tif/final/wind_FINAL.jpg in: /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/delft3d/final/wind.jpg ******************************************************* * Case Completed * ******************************************************* Case starterd at: 2019-09-14 11:14:05.150286 UTC Case completed at: 2019-09-14 11:17:29.529906 UTC **********E N D O F J O B*********************** ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-09-14 11:24:12.357656 UTC inp1= 20190912.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000599/09l submitting calc 2019-09-12 18:00:00 2019-09-14 00:00:00 delft3d -88.5 -60.48 9.55 37.55 4.0 HWRF 72 15 True GDACS/1000599/3_HWRF 6 1 False False 5 aa.txt 20190912.18 1000599/09l 6 wind False *************---------------------****************** GDACS/1000599/3_HWRF 1000599 ndt: 6 it: 0 ndt: 6 idate: 2019-09-12 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_6h.tif input file exist it: 1 ndt: 6 idate: 2019-09-13 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_6h.tif input file exist it: 2 ndt: 6 idate: 2019-09-13 06:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091306.swath.wind_6h.tif input file exist it: 3 ndt: 6 idate: 2019-09-13 12:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091312.swath.wind_6h.tif input file exist it: 4 ndt: 6 idate: 2019-09-13 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091318.swath.wind_6h.tif input file exist it: 5 ndt: 6 idate: 2019-09-14 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091400.swath.wind_6h.tif input file exist GDACS/1000599/3_HWRF 1000599 ndt: 6 it: 0 ndt: 6 idate: 2019-09-12 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_72h.tif input file exist it: 1 ndt: 6 idate: 2019-09-13 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_72h.tif input file exist it: 2 ndt: 6 idate: 2019-09-13 06:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091306.swath.wind_72h.tif input file exist it: 3 ndt: 6 idate: 2019-09-13 12:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091312.swath.wind_72h.tif input file exist it: 4 ndt: 6 idate: 2019-09-13 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091318.swath.wind_72h.tif input file exist it: 5 ndt: 6 idate: 2019-09-14 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091400.swath.wind_72h.tif input file exist home dir /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/ ret 0 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-09-12 18:00:00 currdate 2019-09-14 00:00:00 ndt: 30 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-09-12 18:00:00', '2019-09-13 00:00:00', '2019-09-13 06:00:00', '2019-09-13 12:00:00', '2019-09-13 18:00:00', '2019-09-14 00:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files processing curr bull... 2019-09-14 00:00:00 ********** no new data no new data ******************************************************* * Case Completed * ******************************************************* Case starterd at: 2019-09-14 11:24:12.357656 UTC Case completed at: 2019-09-14 11:24:12.817303 UTC **********E N D O F J O B*********************** ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-09-14 11:34:37.892835 UTC inp1= 20190912.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000599/09l submitting calc 2019-09-12 18:00:00 2019-09-14 00:00:00 delft3d -88.5 -60.48 9.55 37.55 4.0 HWRF 72 15 True GDACS/1000599/3_HWRF 6 1 False False 5 aa.txt 20190912.18 1000599/09l 6 wind False *************---------------------****************** GDACS/1000599/3_HWRF 1000599 ndt: 6 it: 0 ndt: 6 idate: 2019-09-12 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_6h.tif input file exist it: 1 ndt: 6 idate: 2019-09-13 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_6h.tif input file exist it: 2 ndt: 6 idate: 2019-09-13 06:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091306.swath.wind_6h.tif input file exist it: 3 ndt: 6 idate: 2019-09-13 12:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091312.swath.wind_6h.tif input file exist it: 4 ndt: 6 idate: 2019-09-13 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091318.swath.wind_6h.tif input file exist it: 5 ndt: 6 idate: 2019-09-14 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091400.swath.wind_6h.tif input file exist GDACS/1000599/3_HWRF 1000599 ndt: 6 it: 0 ndt: 6 idate: 2019-09-12 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_72h.tif input file exist it: 1 ndt: 6 idate: 2019-09-13 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_72h.tif input file exist it: 2 ndt: 6 idate: 2019-09-13 06:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091306.swath.wind_72h.tif input file exist it: 3 ndt: 6 idate: 2019-09-13 12:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091312.swath.wind_72h.tif input file exist it: 4 ndt: 6 idate: 2019-09-13 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091318.swath.wind_72h.tif input file exist it: 5 ndt: 6 idate: 2019-09-14 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091400.swath.wind_72h.tif input file exist home dir /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/ ret 0 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-09-12 18:00:00 currdate 2019-09-14 00:00:00 ndt: 30 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-09-12 18:00:00', '2019-09-13 00:00:00', '2019-09-13 06:00:00', '2019-09-13 12:00:00', '2019-09-13 18:00:00', '2019-09-14 00:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files processing curr bull... 2019-09-14 00:00:00 ********** no new data no new data ******************************************************* * Case Completed * ******************************************************* Case starterd at: 2019-09-14 11:34:37.892835 UTC Case completed at: 2019-09-14 11:34:38.405042 UTC **********E N D O F J O B*********************** ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-09-14 11:44:07.877383 UTC inp1= 20190912.18 ncores= 5 var= wind lw= aa.txt stormname!!!!!!!!!!!!!!!= 1000599/09l submitting calc 2019-09-12 18:00:00 2019-09-14 00:00:00 delft3d -88.5 -60.48 9.55 37.55 4.0 HWRF 72 15 True GDACS/1000599/3_HWRF 6 1 False False 5 aa.txt 20190912.18 1000599/09l 6 wind False *************---------------------****************** GDACS/1000599/3_HWRF 1000599 ndt: 6 it: 0 ndt: 6 idate: 2019-09-12 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_6h.tif input file exist it: 1 ndt: 6 idate: 2019-09-13 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_6h.tif input file exist it: 2 ndt: 6 idate: 2019-09-13 06:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091306.swath.wind_6h.tif input file exist it: 3 ndt: 6 idate: 2019-09-13 12:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091312.swath.wind_6h.tif input file exist it: 4 ndt: 6 idate: 2019-09-13 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091318.swath.wind_6h.tif input file exist it: 5 ndt: 6 idate: 2019-09-14 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091400.swath.wind_6h.tif input file exist GDACS/1000599/3_HWRF 1000599 ndt: 6 it: 0 ndt: 6 idate: 2019-09-12 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091218.swath.wind_72h.tif input file exist it: 1 ndt: 6 idate: 2019-09-13 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091300.swath.wind_72h.tif input file exist it: 2 ndt: 6 idate: 2019-09-13 06:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091306.swath.wind_72h.tif input file exist it: 3 ndt: 6 idate: 2019-09-13 12:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091312.swath.wind_72h.tif input file exist it: 4 ndt: 6 idate: 2019-09-13 18:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091318.swath.wind_72h.tif input file exist it: 5 ndt: 6 idate: 2019-09-14 00:00:00 /mnt/input/grib_HWRF/2019/1000599/nine09l.2019091400.swath.wind_72h.tif input file exist home dir /mnt/output/SSCS/2019/GDACS/1000599/3_HWRF/ ret 0 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-09-12 18:00:00 currdate 2019-09-14 00:00:00 ndt: 30 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-09-12 18:00:00', '2019-09-13 00:00:00', '2019-09-13 06:00:00', '2019-09-13 12:00:00', '2019-09-13 18:00:00', '2019-09-14 00:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files processing curr bull... 2019-09-14 00:00:00 ********** no new data no new data ******************************************************* * Case Completed * ******************************************************* Case starterd at: 2019-09-14 11:44:07.877383 UTC Case completed at: 2019-09-14 11:44:08.431265 UTC **********E N D O F J O B***********************