******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-04-23 09:55:43.837941 UTC inp1= 20190423.00 ncores= 5 var= wind submitting calc 2019-04-23 00:00:00 2019-04-23 00:00:00 delft3d 29.3 49.62 -26.1 -3.0 4.0 GFS 72 15 True GDACS/1000559/1_GFS 6 1 False False 5 20190423.00 wind False *************---------------------****************** ndt: 1 it: 0 ndt: 1 idate: 2019-04-23 00:00:00 running case from 2019-04-23 00:00:00 for 72 h start= 1 var wind **** gometeo: 72 listWindows rundate:20190423.00 RUNNING 2019-04-23 00:00:00 for 72 hours prevCalcDate 2019-04-22 18:00:00 Nothing to do, case already completed ret= -3 newcase= False forceFinal= False forceBulletin= False home dir /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/ wind and rainfall classification ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-04-23 00:00:00 currdate 2019-04-23 00:00:00 ndt: 0 delta: 6 nt 72 alldate: DatetimeIndex(['2019-04-23'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files processing curr bull... 2019-04-23 00:00:00 forcing GFS verifying that input file is present start reading nc... wind ntNC: 73 ntmax 73 use all data in nc file ...create velAll 73 93 81 ...start calculating velAll ...end calculating velAll 73 93 81 nt,nx,ny, ntmax 73 93 81 73 [29.5, 0.25, 0, -3.0, 0, -0.25] *********** 6 1 6 varMAX.shape (93, 81) 29.5 49.5 -26.0 -3.0 lon 81 lat 93 savemap /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.00/20190423.00_wind.jpg dtk,nt,ntmax 6 73 73 67 k1, k2, ht 0 6 0 k1, k2, ht 6 12 6 k1, k2, ht 12 18 12 k1, k2, ht 18 24 18 k1, k2, ht 24 30 24 k1, k2, ht 30 36 30 k1, k2, ht 36 42 36 k1, k2, ht 42 48 42 k1, k2, ht 48 54 48 k1, k2, ht 54 60 54 k1, k2, ht 60 66 60 k1, k2, ht 66 72 66 processing all past bull only if Past=True... False no past >> 2.2. Process final data files compfile /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/20190423.00_Final_completed_wind.txt FINAL alldate: DatetimeIndex(['2019-04-23'], dtype='datetime64[ns]', freq='6H') 1 **FIRST cp /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.00/20190423.00_wind.tif /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/wind_final.tif max file created /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/wind_final.tif 29.625 49.625 -26.125 -3.125 lon 81 lat 93 savemap /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/wind_FINAL.jpg ret: 0 ============================================ 3. Classify meteo + GDACS index score ============================================ >> 3.1. Classify curr + past forecast Creating output file that is 2430P x 2790L. Processing input file /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.00/20190423.00_wind_t0.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.00/20190423.00_wind_t0.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.00/20190423.00_wind_t0.tif to destination /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.00//wind10m_res_t0.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2430P x 2790L. 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/1000559/1_GFS/class/20190423.00//wind10m_popfile_t0_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2430P x 2790L. 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/1000559/1_GFS/class/20190423.00//wind10m_countryfile_t0_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.00/20190423.00_wind_t0.tif hurName: hdate: var: wind10m description: wind10m: _ (GFS) OutDir: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.00// PopFile: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.00//wind10m_popfile_t0_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.00//wind10m_countryfile_t0_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.00/wind_popDensValues_t0.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 >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.00833333333333 cellsize 0.00833333333333 >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.00/wind_popDensValues_t0.xml >> 7. remove files done t0 completed Creating output file that is 2430P x 2790L. Processing input file /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.00/20190423.00_wind.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.00/20190423.00_wind.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.00/20190423.00_wind.tif to destination /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.00//wind10m_res_all.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2430P x 2790L. 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/1000559/1_GFS/class/20190423.00//wind10m_popfile_all_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2430P x 2790L. 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/1000559/1_GFS/class/20190423.00//wind10m_countryfile_all_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.00/20190423.00_wind.tif hurName: hdate: var: wind10m description: wind10m: _ (GFS) OutDir: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.00// PopFile: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.00//wind10m_popfile_all_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.00//wind10m_countryfile_all_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.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 >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.00833333333333 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/1000559/1_GFS/class/20190423.00/wind_popDensValues_all.xml >> 7. remove files done Creating output file that is 2430P x 2790L. Processing input file /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/wind_final.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/wind_final.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/wind_final.tif to destination /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/final//wind10m_res_final.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2430P x 2790L. 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/1000559/1_GFS/class/final//wind10m_popfile_final_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2430P x 2790L. 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/1000559/1_GFS/class/final//wind10m_countryfile_final_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. t0 completed copy: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.00/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/delft3d/20190423.00/wind_popDensValues_all.xml >> 3.2. Classify final folder /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/final/20190423.00_final_completed_wind.txt input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/wind_final.tif hurName: hdate: var: wind10m description: wind10m: _ (GFS) OutDir: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/final// PopFile: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/final//wind10m_popfile_final_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/final//wind10m_countryfile_final_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/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 >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.00833333333333 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/1000559/1_GFS/class/final/wind_popDensValues_final.xml >> 7. remove files done copy: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/final/wind_popDensValues_final.xml in: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/delft3d/final/wind_popDensValues_final.xml ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-04-23 09:56:59.767082 UTC inp1= 20190423.00 ncores= 5 var= rain submitting calc 2019-04-23 00:00:00 2019-04-23 00:00:00 delft3d 29.3 49.62 -26.1 -3.0 4.0 GFS 72 15 True GDACS/1000559/1_GFS 6 1 False False 5 20190423.00 rain False *************---------------------****************** ndt: 1 it: 0 ndt: 1 idate: 2019-04-23 00:00:00 running case from 2019-04-23 00:00:00 for 72 h start= 1 var rain **** gometeo: 72 listWindows rundate:20190423.00 RUNNING 2019-04-23 00:00:00 for 72 hours prevCalcDate 2019-04-22 18:00:00 Nothing to do, case already completed ret= -3 newcase= False forceFinal= False forceBulletin= False home dir /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/ wind and rainfall classification ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-04-23 00:00:00 currdate 2019-04-23 00:00:00 ndt: 0 delta: 6 nt 72 alldate: DatetimeIndex(['2019-04-23'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files processing curr bull... 2019-04-23 00:00:00 processing all past bull only if Past=True... False no past >> 2.2. Process final data files compfile /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/20190423.00_Final_completed_rain.txt ============================================ 3. Classify meteo + GDACS index score ============================================ >> 3.1. Classify curr + past forecast t0 completed t0 completed >> 3.2. Classify final folder /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/final/20190423.00_final_completed_rain.txt ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-04-23 09:58:39.549906 UTC inp1= 20190423.00 ncores= 5 var= wind submitting calc 2019-04-23 00:00:00 2019-04-23 00:00:00 delft3d 29.3 49.62 -26.1 -3.0 4.0 GFS 72 15 True GDACS/1000559/1_GFS 6 1 False False 5 20190423.00 wind False *************---------------------****************** ndt: 1 it: 0 ndt: 1 idate: 2019-04-23 00:00:00 running case from 2019-04-23 00:00:00 for 72 h start= 1 var wind **** gometeo: 72 listWindows rundate:20190423.00 RUNNING 2019-04-23 00:00:00 for 72 hours prevCalcDate 2019-04-22 18:00:00 Nothing to do, case already completed ret= -3 newcase= False forceFinal= False forceBulletin= False home dir /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/ wind and rainfall classification ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-04-23 00:00:00 currdate 2019-04-23 00:00:00 ndt: 0 delta: 6 nt 72 alldate: DatetimeIndex(['2019-04-23'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files processing curr bull... 2019-04-23 00:00:00 processing all past bull only if Past=True... False no past >> 2.2. Process final data files compfile /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/20190423.00_Final_completed_wind.txt ============================================ 3. Classify meteo + GDACS index score ============================================ >> 3.1. Classify curr + past forecast t0 completed t0 completed >> 3.2. Classify final folder /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/final/20190423.00_final_completed_wind.txt ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-04-23 09:58:49.580425 UTC inp1= 20190423.00 ncores= 5 var= rain submitting calc 2019-04-23 00:00:00 2019-04-23 00:00:00 delft3d 29.3 49.62 -26.1 -3.0 4.0 GFS 72 15 True GDACS/1000559/1_GFS 6 1 False False 5 20190423.00 rain False *************---------------------****************** ndt: 1 it: 0 ndt: 1 idate: 2019-04-23 00:00:00 running case from 2019-04-23 00:00:00 for 72 h start= 1 var rain **** gometeo: 72 listWindows rundate:20190423.00 RUNNING 2019-04-23 00:00:00 for 72 hours prevCalcDate 2019-04-22 18:00:00 Nothing to do, case already completed ret= -3 newcase= False forceFinal= False forceBulletin= False home dir /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/ wind and rainfall classification ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-04-23 00:00:00 currdate 2019-04-23 00:00:00 ndt: 0 delta: 6 nt 72 alldate: DatetimeIndex(['2019-04-23'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files processing curr bull... 2019-04-23 00:00:00 processing all past bull only if Past=True... False no past >> 2.2. Process final data files compfile /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/20190423.00_Final_completed_rain.txt ============================================ 3. Classify meteo + GDACS index score ============================================ >> 3.1. Classify curr + past forecast t0 completed t0 completed >> 3.2. Classify final folder /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/final/20190423.00_final_completed_rain.txt ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-04-23 11:16:12.575308 UTC inp1= 20190423.00 ncores= 5 var= wind submitting calc 2019-04-23 00:00:00 2019-04-23 06:00:00 delft3d 29.3 49.62 -26.1 -3.0 4.0 GFS 72 15 True GDACS/1000559/1_GFS 6 1 False False 5 20190423.00 wind False *************---------------------****************** ndt: 2 it: 0 ndt: 2 idate: 2019-04-23 00:00:00 running case from 2019-04-23 00:00:00 for 72 h start= 1 var wind **** gometeo: 72 listWindows rundate:20190423.00 RUNNING 2019-04-23 00:00:00 for 72 hours prevCalcDate 2019-04-22 18:00:00 Nothing to do, case already completed ret= -3 newcase= False forceFinal= False forceBulletin= False it: 1 ndt: 2 idate: 2019-04-23 06:00:00 running case from 2019-04-23 06:00:00 for 72 h start= 0 var wind **** gometeo: 72 listWindows rundate:20190423.06 RUNNING 2019-04-23 06:00:00 for 72 hours prevCalcDate 2019-04-23 00:00:00 0 .. 20 .. 50 .. 70 .. 100 - Done Nothing to do, case already completed ret= -3 newcase= False forceFinal= False forceBulletin= False home dir /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/ wind and rainfall classification ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-04-23 00:00:00 currdate 2019-04-23 06:00:00 ndt: 6 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-04-23 00:00:00', '2019-04-23 06:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files processing curr bull... 2019-04-23 06:00:00 forcing GFS verifying that input file is present start reading nc... wind ntNC: 73 ntmax 73 use all data in nc file ...create velAll 73 93 81 ...start calculating velAll ...end calculating velAll 73 93 81 nt,nx,ny, ntmax 73 93 81 73 [29.5, 0.25, 0, -3.0, 0, -0.25] *********** 6 1 6 varMAX.shape (93, 81) 29.5 49.5 -26.0 -3.0 lon 81 lat 93 savemap /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.06/20190423.06_wind.jpg dtk,nt,ntmax 6 73 73 67 k1, k2, ht 0 6 0 k1, k2, ht 6 12 6 k1, k2, ht 12 18 12 k1, k2, ht 18 24 18 k1, k2, ht 24 30 24 k1, k2, ht 30 36 30 k1, k2, ht 36 42 36 k1, k2, ht 42 48 42 k1, k2, ht 48 54 48 k1, k2, ht 54 60 54 k1, k2, ht 60 66 60 k1, k2, ht 66 72 66 processing all past bull only if Past=True... True itdate, istime 2019-04-23 00:00:00 20190423.00 meteo-processing past forecast already completed >> 2.2. Process final data files compfile /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/20190423.06_Final_completed_wind.txt FINAL remove maxtif alldate: DatetimeIndex(['2019-04-23 00:00:00', '2019-04-23 06:00:00'], dtype='datetime64[ns]', freq='6H') 2 date: 2019-04-23 06:00:00 python /mnt/output/SSCS/scripts_dev/gdal_calc.py -A /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.00/20190423.00_wind_stept0.tif -B /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.06/20190423.06_wind.tif --outfile=/mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/wind_final.tif --calc="maximum(A,B)" max file created /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/wind_final.tif 29.625 49.625 -26.125 -3.125 lon 81 lat 93 savemap /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/wind_FINAL.jpg ret: 0 ============================================ 3. Classify meteo + GDACS index score ============================================ >> 3.1. Classify curr + past forecast Creating output file that is 2430P x 2790L. Processing input file /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.06/20190423.06_wind_t0.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.06/20190423.06_wind_t0.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.06/20190423.06_wind_t0.tif to destination /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.06//wind10m_res_t0.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2430P x 2790L. 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/1000559/1_GFS/class/20190423.06//wind10m_popfile_t0_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2430P x 2790L. 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/1000559/1_GFS/class/20190423.06//wind10m_countryfile_t0_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.06/20190423.06_wind_t0.tif hurName: hdate: var: wind10m description: wind10m: _ (GFS) OutDir: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.06// PopFile: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.06//wind10m_popfile_t0_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.06//wind10m_countryfile_t0_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.06/wind_popDensValues_t0.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 >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.00833333333333 cellsize 0.00833333333333 >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.06/wind_popDensValues_t0.xml >> 7. remove files done t0 completed Creating output file that is 2430P x 2790L. Processing input file /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.06/20190423.06_wind.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.06/20190423.06_wind.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.06/20190423.06_wind.tif to destination /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.06//wind10m_res_all.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2430P x 2790L. 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/1000559/1_GFS/class/20190423.06//wind10m_popfile_all_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2430P x 2790L. 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/1000559/1_GFS/class/20190423.06//wind10m_countryfile_all_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.06/20190423.06_wind.tif hurName: hdate: var: wind10m description: wind10m: _ (GFS) OutDir: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.06// PopFile: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.06//wind10m_popfile_all_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.06//wind10m_countryfile_all_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.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 >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.00833333333333 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/1000559/1_GFS/class/20190423.06/wind_popDensValues_all.xml >> 7. remove files done Creating output file that is 2430P x 2790L. Processing input file /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/wind_final.tif. Using internal nodata values (e.g. 3.40282e+38) for image /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/wind_final.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/wind_final.tif to destination /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/final//wind10m_res_final.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2430P x 2790L. 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/1000559/1_GFS/class/final//wind10m_popfile_final_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2430P x 2790L. 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/1000559/1_GFS/class/final//wind10m_countryfile_final_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. t0 completed copy: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.06/wind_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/delft3d/20190423.06/wind_popDensValues_all.xml >> 3.2. Classify final folder /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/final/20190423.06_final_completed_wind.txt input var: wind10m Input File: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/wind_final.tif hurName: hdate: var: wind10m description: wind10m: _ (GFS) OutDir: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/final// PopFile: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/final//wind10m_popfile_final_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/final//wind10m_countryfile_final_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/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 >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.00833333333333 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/1000559/1_GFS/class/final/wind_popDensValues_final.xml >> 7. remove files done ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-04-23 11:18:12.737794 UTC inp1= 20190423.00 ncores= 5 var= rain submitting calc 2019-04-23 00:00:00 2019-04-23 06:00:00 delft3d 29.3 49.62 -26.1 -3.0 4.0 GFS 72 15 True GDACS/1000559/1_GFS 6 1 False False 5 20190423.00 rain False *************---------------------****************** ndt: 2 it: 0 ndt: 2 idate: 2019-04-23 00:00:00 running case from 2019-04-23 00:00:00 for 72 h start= 1 var rain **** gometeo: 72 listWindows rundate:20190423.00 RUNNING 2019-04-23 00:00:00 for 72 hours prevCalcDate 2019-04-22 18:00:00 Nothing to do, case already completed ret= -3 newcase= False forceFinal= False forceBulletin= False it: 1 ndt: 2 idate: 2019-04-23 06:00:00 running case from 2019-04-23 06:00:00 for 72 h start= 0 var rain **** gometeo: 72 listWindows rundate:20190423.06 RUNNING 2019-04-23 06:00:00 for 72 hours prevCalcDate 2019-04-23 00:00:00 0 .. 20 .. 50 .. 70 .. 100 - Done Nothing to do, case already completed ret= -3 newcase= False forceFinal= False forceBulletin= False home dir /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/ wind and rainfall classification ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-04-23 00:00:00 currdate 2019-04-23 06:00:00 ndt: 6 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-04-23 00:00:00', '2019-04-23 06:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files processing curr bull... 2019-04-23 06:00:00 forcing GFS verifying that input file is present start reading nc... rain use all data in nc file nt,nx,ny, ntmax 72 93 81 73 nt,nx,ny, ntmax 72 93 81 73 [29.5, 0.25, 0, -3.0, 0, -0.25] *********** 6 1 6 varMAX.shape (93, 81) 29.5 49.5 -26.0 -3.0 lon 81 lat 93 savemap /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.06/20190423.06_rain.jpg dtk,nt,ntmax 6 72 73 67 k1, k2, ht 0 6 0 k1, k2, ht 6 12 6 k1, k2, ht 12 18 12 k1, k2, ht 18 24 18 k1, k2, ht 24 30 24 k1, k2, ht 30 36 30 k1, k2, ht 36 42 36 k1, k2, ht 42 48 42 k1, k2, ht 48 54 48 k1, k2, ht 54 60 54 k1, k2, ht 60 66 60 k1, k2, ht 66 72 66 processing all past bull only if Past=True... True itdate, istime 2019-04-23 00:00:00 20190423.00 meteo-processing past forecast already completed >> 2.2. Process final data files compfile /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/20190423.06_Final_completed_rain.txt FINAL remove maxtif alldate: DatetimeIndex(['2019-04-23 00:00:00', '2019-04-23 06:00:00'], dtype='datetime64[ns]', freq='6H') 2 date: 2019-04-23 06:00:00 python /mnt/output/SSCS/scripts_dev/gdal_calc.py -A /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.00/20190423.00_rain_stept0.tif -B /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.06/20190423.06_rain.tif --outfile=/mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/rain_final.tif --calc="A+B" max file created /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/rain_final.tif 29.625 49.625 -26.125 -3.125 lon 81 lat 93 savemap /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/rain_FINAL.jpg ret: 0 ============================================ 3. Classify meteo + GDACS index score ============================================ >> 3.1. Classify curr + past forecast Creating output file that is 2430P x 2790L. Processing input file /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.06/20190423.06_rain_t0.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.06/20190423.06_rain_t0.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.06/20190423.06_rain_t0.tif to destination /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.06//rain_res_t0.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2430P x 2790L. 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/1000559/1_GFS/class/20190423.06//rain_popfile_t0_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2430P x 2790L. 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/1000559/1_GFS/class/20190423.06//rain_countryfile_t0_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. input var: rain Input File: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.06/20190423.06_rain_t0.tif hurName: hdate: var: rain description: rain: _ (GFS) OutDir: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.06// PopFile: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.06//rain_popfile_t0_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.06//rain_countryfile_t0_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.06/rain_popDensValues_t0.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 >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.00833333333333 cellsize 0.00833333333333 >> 6. print / save output in /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.06/rain_popDensValues_t0.xml >> 7. remove files done t0 completed Creating output file that is 2430P x 2790L. Processing input file /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.06/20190423.06_rain.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.06/20190423.06_rain.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.06/20190423.06_rain.tif to destination /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.06//rain_res_all.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2430P x 2790L. 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/1000559/1_GFS/class/20190423.06//rain_popfile_all_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2430P x 2790L. 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/1000559/1_GFS/class/20190423.06//rain_countryfile_all_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. input var: rain Input File: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/20190423.06/20190423.06_rain.tif hurName: hdate: var: rain description: rain: _ (GFS) OutDir: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.06// PopFile: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.06//rain_popfile_all_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.06//rain_countryfile_all_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.06/rain_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 >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.00833333333333 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/1000559/1_GFS/class/20190423.06/rain_popDensValues_all.xml >> 7. remove files done Creating output file that is 2430P x 2790L. Processing input file /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/rain_final.tif. Using internal nodata values (e.g. 3.40282e+38) for image /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/rain_final.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/rain_final.tif to destination /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/final//rain_res_final.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2430P x 2790L. 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/1000559/1_GFS/class/final//rain_popfile_final_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 2430P x 2790L. 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/1000559/1_GFS/class/final//rain_countryfile_final_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. t0 completed copy: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/20190423.06/rain_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/delft3d/20190423.06/rain_popDensValues_all.xml >> 3.2. Classify final folder /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/final/20190423.06_final_completed_rain.txt input var: rain Input File: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/rain_final.tif hurName: hdate: var: rain description: rain: _ (GFS) OutDir: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/final// PopFile: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/final//rain_popfile_final_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/final//rain_countryfile_final_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/final/rain_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 >> 3b. extract a piece of countries corresponding to the required bounding box and resolution/proj of pop density >> 4. classify the vmax file creating another array of values classified cellsize 0.00833333333333 cellsize 0.00833333333333 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/1000559/1_GFS/class/final/rain_popDensValues_final.xml >> 7. remove files done ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-04-23 11:29:23.516321 UTC inp1= 20190423.00 ncores= 5 var= wind submitting calc 2019-04-23 00:00:00 2019-04-23 06:00:00 delft3d 29.3 49.62 -26.1 -3.0 4.0 GFS 72 15 True GDACS/1000559/1_GFS 6 1 False False 5 20190423.00 wind False *************---------------------****************** ndt: 2 it: 0 ndt: 2 idate: 2019-04-23 00:00:00 running case from 2019-04-23 00:00:00 for 72 h start= 1 var wind **** gometeo: 72 listWindows rundate:20190423.00 RUNNING 2019-04-23 00:00:00 for 72 hours prevCalcDate 2019-04-22 18:00:00 Nothing to do, case already completed ret= -3 newcase= False forceFinal= False forceBulletin= False it: 1 ndt: 2 idate: 2019-04-23 06:00:00 running case from 2019-04-23 06:00:00 for 72 h start= 0 var wind **** gometeo: 72 listWindows rundate:20190423.06 RUNNING 2019-04-23 06:00:00 for 72 hours prevCalcDate 2019-04-23 00:00:00 Nothing to do, case already completed ret= -3 newcase= False forceFinal= False forceBulletin= False home dir /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/ wind and rainfall classification ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-04-23 00:00:00 currdate 2019-04-23 06:00:00 ndt: 6 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-04-23 00:00:00', '2019-04-23 06:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files processing curr bull... 2019-04-23 06:00:00 processing all past bull only if Past=True... True itdate, istime 2019-04-23 00:00:00 20190423.00 meteo-processing past forecast already completed >> 2.2. Process final data files compfile /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/20190423.06_Final_completed_wind.txt ============================================ 3. Classify meteo + GDACS index score ============================================ >> 3.1. Classify curr + past forecast t0 completed t0 completed >> 3.2. Classify final folder /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/final/20190423.06_final_completed_wind.txt ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-04-23 11:29:27.817040 UTC inp1= 20190423.00 ncores= 5 var= rain submitting calc 2019-04-23 00:00:00 2019-04-23 06:00:00 delft3d 29.3 49.62 -26.1 -3.0 4.0 GFS 72 15 True GDACS/1000559/1_GFS 6 1 False False 5 20190423.00 rain False *************---------------------****************** ndt: 2 it: 0 ndt: 2 idate: 2019-04-23 00:00:00 running case from 2019-04-23 00:00:00 for 72 h start= 1 var rain **** gometeo: 72 listWindows rundate:20190423.00 RUNNING 2019-04-23 00:00:00 for 72 hours prevCalcDate 2019-04-22 18:00:00 Nothing to do, case already completed ret= -3 newcase= False forceFinal= False forceBulletin= False it: 1 ndt: 2 idate: 2019-04-23 06:00:00 running case from 2019-04-23 06:00:00 for 72 h start= 0 var rain **** gometeo: 72 listWindows rundate:20190423.06 RUNNING 2019-04-23 06:00:00 for 72 hours prevCalcDate 2019-04-23 00:00:00 Nothing to do, case already completed ret= -3 newcase= False forceFinal= False forceBulletin= False home dir /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/ wind and rainfall classification ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-04-23 00:00:00 currdate 2019-04-23 06:00:00 ndt: 6 delta: 6 nt1=delta 6 nt 72 alldate: DatetimeIndex(['2019-04-23 00:00:00', '2019-04-23 06:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files processing curr bull... 2019-04-23 06:00:00 processing all past bull only if Past=True... True itdate, istime 2019-04-23 00:00:00 20190423.00 meteo-processing past forecast already completed >> 2.2. Process final data files compfile /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/tif/final/20190423.06_Final_completed_rain.txt ============================================ 3. Classify meteo + GDACS index score ============================================ >> 3.1. Classify curr + past forecast t0 completed t0 completed >> 3.2. Classify final folder /mnt/output/SSCS/2019/GDACS/1000559/1_GFS/class/final/20190423.06_final_completed_rain.txt ==============================================================