******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-08-27 12:45:48.781724 UTC inp1= 20190824.12 ncores= 5 var= rain submitting calc 2019-08-24 12:00:00 2019-08-24 12:00:00 delft3d -84.21 -41.76 5.0 32.44 4.0 GFS 24 15 True GDACS/1000588/003_GFS 6 1 False False 5 20190824.12 rain False *************---------------------****************** stormname ndt: 1 it: 0 ndt: 1 idate: 2019-08-24 12:00:00 running case from 2019-08-24 12:00:00 for 24 h start= 1 var rain **** gometeo: 24 listWindows rundate:20190824.12 RUNNING 2019-08-24 12:00:00 for 24 hours Nothing to do, case already completed ret= -3 newcase= False forceFinal= False forceBulletin= False netcdf already created home dir /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/ ret -3 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-08-24 12:00:00 currdate 2019-08-24 12:00:00 ndt: 0 delta: 6 nt 24 alldate: DatetimeIndex(['2019-08-24 12:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files processing curr bull... 2019-08-24 12:00:00 forcing GFS verifying that input file is present start reading nc... rain use all data in nc file nt,nx,ny, ntmax 24 110 169 25 nt,nx,ny, ntmax 24 110 169 25 forcing GFS >> 0. SET INPUT/OUTPUT GFS /mnt/output/GDACS/TC/TCTrack/1000588/GFS/2019082412/all_inpData.txt error track file-->map without track trackfile 0 -84.0 -42.0 5.0 32.25 lon 169 lat 110 savemap /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/20190824.12/20190824.12_rain.jpg processing all past bull only if Past=True... False no past >> 2.2. Process final data files compfile /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/final/20190824.12_Final_completed_rain.txt FINAL alldate: DatetimeIndex(['2019-08-24 12:00:00'], dtype='datetime64[ns]', freq='6H') alldate.shape 1 **FIRST cp /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/20190824.12/20190824.12_rain.tif /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/final/rain_final.tif max file created /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/final/rain_final.tif stormname stormname TCbullsource >> 0. SET INPUT/OUTPUT GFS /mnt/output/GDACS/TC/TCTrack/1000588/GFS/2019082412/all_inpData.txt error track file-->map without track trackfile 0 -83.875 -41.875 4.875 32.125 lon 169 lat 110 savemap /mnt/output/SSCS/2019/GDACS/1000588/003_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 5070P x 3300L. Processing input file /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/20190824.12/20190824.12_rain_stept0.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/20190824.12/20190824.12_rain_stept0.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/20190824.12/20190824.12_rain_stept0.tif to destination /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12//rain_res_stept0.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 5070P x 3300L. 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/1000588/003_GFS/class/20190824.12//rain_popfile_stept0_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 5070P x 3300L. 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/1000588/003_GFS/class/20190824.12//rain_countryfile_stept0_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. input var: rain Input File: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/20190824.12/20190824.12_rain_stept0.tif hurName: hdate: var: rain description: rain: _ (GFS) OutDir: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12// PopFile: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12//rain_popfile_stept0_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12//rain_countryfile_stept0_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12/rain_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 -84.0 4.75000000001 -41.75 32.25 GFS /usr/bin/gdalwarp -te -84.0 4.75000000001 -41.75 32.25 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12//rain_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 -84.0 4.75000000001 -41.75 32.25 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12//rain_countryfile_stept0_clipped.tif" >> 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/1000588/003_GFS/class/20190824.12/rain_popDensValues_stept0.xml >> 7. remove files done t0 completed Creating output file that is 5070P x 3300L. Processing input file /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/20190824.12/20190824.12_rain.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/20190824.12/20190824.12_rain.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/20190824.12/20190824.12_rain.tif to destination /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12//rain_res_all.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 5070P x 3300L. 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/1000588/003_GFS/class/20190824.12//rain_popfile_all_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 5070P x 3300L. 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/1000588/003_GFS/class/20190824.12//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/1000588/003_GFS/tif/20190824.12/20190824.12_rain.tif hurName: hdate: var: rain description: rain: _ (GFS) OutDir: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12// PopFile: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12//rain_popfile_all_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12//rain_countryfile_all_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12/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 split,rastershift False -1 -84.0 4.75000000001 -41.75 32.25 GFS /usr/bin/gdalwarp -te -84.0 4.75000000001 -41.75 32.25 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12//rain_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 -84.0 4.75000000001 -41.75 32.25 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12//rain_countryfile_all_clipped.tif" >> 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/1000588/003_GFS/class/20190824.12/rain_popDensValues_all.xml >> 7. remove files done Creating output file that is 5070P x 3300L. Processing input file /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/final/rain_final.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/final/rain_final.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/final/rain_final.tif to destination /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/final//rain_res_final.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 5070P x 3300L. 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/1000588/003_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 5070P x 3300L. 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/1000588/003_GFS/class/final//rain_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/1000588/003_GFS/class/20190824.12/rain_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/delft3d/20190824.12/rain_popDensValues_t0.xml copy: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12/rain_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/delft3d/20190824.12/rain_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/20190824.12/20190824.12_rain.jpg in: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/delft3d/20190824.12/rain.jpg >> 3.2. Classify final folder /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/final/20190824.12_final_completed_rain.txt input var: rain Input File: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/final/rain_final.tif hurName: hdate: var: rain description: rain: _ (GFS) OutDir: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/final// PopFile: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/final//rain_popfile_final_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/final//rain_countryfile_final_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/final/rain_popDensValues_final.xml outDir created popfile: LandScan popCellSize= 0.00833333333333 projection= GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]] >> 1. resample the tif file to the resolution and proj of pop density 0.00833333333333 deg >> 2. read the charactristics of the input file >> 3a. extract a piece of pop. file corresponding to the required bounding box split,rastershift False -1 -84.0 4.75000000001 -41.75 32.25 GFS /usr/bin/gdalwarp -te -84.0 4.75000000001 -41.75 32.25 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/final//rain_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 -84.0 4.75000000001 -41.75 32.25 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/final//rain_countryfile_final_clipped.tif" >> 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/1000588/003_GFS/class/final/rain_popDensValues_final.xml >> 7. remove files done copy: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/final/rain_popDensValues_final.xml in: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/delft3d/final/rain_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/final/rain_FINAL.jpg in: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/delft3d/final/rain.jpg ============================================================== ******************************************************* * Storm Surge Calculation System (SSCS) * ******************************************************* Now is : 2019-08-27 13:22:13.375841 UTC inp1= 20190824.12 ncores= 5 var= rain submitting calc 2019-08-24 12:00:00 2019-08-24 12:00:00 delft3d -84.21 -41.76 5.0 32.44 4.0 GFS 24 15 True GDACS/1000588/003_GFS 6 1 False False 5 20190824.12 6 rain False *************---------------------****************** stormname ndt: 1 it: 0 ndt: 1 idate: 2019-08-24 12:00:00 running case from 2019-08-24 12:00:00 for 24 h start= 1 var rain **** gometeo: 24 listWindows rundate:20190824.12 RUNNING 2019-08-24 12:00:00 for 24 hours Nothing to do, case already completed ret= -3 newcase= False forceFinal= False forceBulletin= False netcdf already created home dir /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/ ret -3 classifications ============================================ 1. Set INPUT/OUTPUT ============================================ startdate 2019-08-24 12:00:00 currdate 2019-08-24 12:00:00 ndt: 0 delta: 6 nt 24 alldate: DatetimeIndex(['2019-08-24 12:00:00'], dtype='datetime64[ns]', freq='6H') ============================================ 2. Processing meteo files: nc2tif ============================================ >> 2.1. Process curr + past files processing curr bull... 2019-08-24 12:00:00 forcing GFS verifying that input file is present start reading nc... rain use all data in nc file nt,nx,ny, ntmax 24 110 169 25 nt,nx,ny, ntmax 24 110 169 25 forcing GFS >> 0. SET INPUT/OUTPUT GFS /mnt/output/GDACS/TC/TCTrack/1000588/GFS/2019082412/all_inpData.txt error track file-->map without track trackfile 0 -84.0 -42.0 5.0 32.25 lon 169 lat 110 savemap /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/20190824.12/20190824.12_rain.jpg processing all past bull only if Past=True... False no past >> 2.2. Process final data files compfile /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/final/20190824.12_Final_completed_rain.txt FINAL alldate: DatetimeIndex(['2019-08-24 12:00:00'], dtype='datetime64[ns]', freq='6H') alldate.shape 1 **FIRST cp /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/20190824.12/20190824.12_rain.tif /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/final/rain_final.tif max file created /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/final/rain_final.tif stormname stormname TCbullsource >> 0. SET INPUT/OUTPUT GFS /mnt/output/GDACS/TC/TCTrack/1000588/GFS/2019082412/all_inpData.txt error track file-->map without track trackfile 0 -83.875 -41.875 4.875 32.125 lon 169 lat 110 savemap /mnt/output/SSCS/2019/GDACS/1000588/003_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 5070P x 3300L. Processing input file /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/20190824.12/20190824.12_rain_stept0.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/20190824.12/20190824.12_rain_stept0.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/20190824.12/20190824.12_rain_stept0.tif to destination /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12//rain_res_stept0.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 5070P x 3300L. 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/1000588/003_GFS/class/20190824.12//rain_popfile_stept0_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 5070P x 3300L. 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/1000588/003_GFS/class/20190824.12//rain_countryfile_stept0_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. input var: rain Input File: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/20190824.12/20190824.12_rain_stept0.tif hurName: hdate: var: rain description: rain: _ (GFS) OutDir: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12// PopFile: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12//rain_popfile_stept0_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12//rain_countryfile_stept0_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12/rain_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 -84.0 4.75000000001 -41.75 32.25 GFS /usr/bin/gdalwarp -te -84.0 4.75000000001 -41.75 32.25 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12//rain_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 -84.0 4.75000000001 -41.75 32.25 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12//rain_countryfile_stept0_clipped.tif" >> 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/1000588/003_GFS/class/20190824.12/rain_popDensValues_stept0.xml >> 7. remove files done t0 completed Creating output file that is 5070P x 3300L. Processing input file /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/20190824.12/20190824.12_rain.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/20190824.12/20190824.12_rain.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/20190824.12/20190824.12_rain.tif to destination /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12//rain_res_all.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 5070P x 3300L. 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/1000588/003_GFS/class/20190824.12//rain_popfile_all_clipped.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 5070P x 3300L. 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/1000588/003_GFS/class/20190824.12//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/1000588/003_GFS/tif/20190824.12/20190824.12_rain.tif hurName: hdate: var: rain description: rain: _ (GFS) OutDir: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12// PopFile: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12//rain_popfile_all_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12//rain_countryfile_all_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12/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 split,rastershift False -1 -84.0 4.75000000001 -41.75 32.25 GFS /usr/bin/gdalwarp -te -84.0 4.75000000001 -41.75 32.25 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12//rain_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 -84.0 4.75000000001 -41.75 32.25 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12//rain_countryfile_all_clipped.tif" >> 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/1000588/003_GFS/class/20190824.12/rain_popDensValues_all.xml >> 7. remove files done Creating output file that is 5070P x 3300L. Processing input file /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/final/rain_final.tif. Using internal nodata values (e.g. -999) for image /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/final/rain_final.tif. Copying nodata values from source /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/final/rain_final.tif to destination /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/final//rain_res_final.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 5070P x 3300L. 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/1000588/003_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 5070P x 3300L. 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/1000588/003_GFS/class/final//rain_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/1000588/003_GFS/class/20190824.12/rain_popDensValues_stept0.xml in: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/delft3d/20190824.12/rain_popDensValues_t0.xml copy: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/20190824.12/rain_popDensValues_all.xml in: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/delft3d/20190824.12/rain_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/20190824.12/20190824.12_rain.jpg in: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/delft3d/20190824.12/rain.jpg >> 3.2. Classify final folder /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/final/20190824.12_final_completed_rain.txt input var: rain Input File: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/final/rain_final.tif hurName: hdate: var: rain description: rain: _ (GFS) OutDir: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/final// PopFile: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/final//rain_popfile_final_clipped.tif country: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/final//rain_countryfile_final_clipped.tif outxml file: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/final/rain_popDensValues_final.xml outDir created popfile: LandScan popCellSize= 0.00833333333333 projection= GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]] >> 1. resample the tif file to the resolution and proj of pop density 0.00833333333333 deg >> 2. read the charactristics of the input file >> 3a. extract a piece of pop. file corresponding to the required bounding box split,rastershift False -1 -84.0 4.75000000001 -41.75 32.25 GFS /usr/bin/gdalwarp -te -84.0 4.75000000001 -41.75 32.25 "/mnt/output/GDACS/DATA/lspop20141.tif" "/mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/final//rain_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 -84.0 4.75000000001 -41.75 32.25 "/mnt/output/GDACS/DATA/countries.tif" "/mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/final//rain_countryfile_final_clipped.tif" >> 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/1000588/003_GFS/class/final/rain_popDensValues_final.xml >> 7. remove files done copy: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/class/final/rain_popDensValues_final.xml in: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/delft3d/final/rain_popDensValues.xml copy: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/tif/final/rain_FINAL.jpg in: /mnt/output/SSCS/2019/GDACS/1000588/003_GFS/delft3d/final/rain.jpg ==============================================================