Following commands executed to create 0.05deg soil, snowbands and lai data files. At the end of command execution, it is only created 0.05deg soil and snowbands without lai file. There is no error message is generated. So, pls help me to double check the commands.
Note: lancover.tiff and dem.tiff are in RHEAS/scripts
Following commands executed to create 0.05deg soil, snowbands and lai data files. At the end of command execution, it is only created 0.05deg soil and snowbands without lai file. There is no error message is generated. So, pls help me to double check the commands.
Note: lancover.tiff and dem.tiff are in RHEAS/scripts
in ipython
from resample_vic import * resample_soil(11, 13, 108, 110, 0.05) resample_vegetation("landcover.tif", "../data/vic/global_soil_0.05deg.txt", "../data/vic/vic_veglib.txt", "../data/global_lai_0.05.veg", 0.05) resample_snowbands("dem.tif", "../data/vic/global_soil_0.05deg.txt", "../data/vic/global_snowbands_0.05deg.txt", 0.05, 5)
outside ipython
echo "begin;" > vicsoils5km.sql awk '{printf("insert into vic.soils (id,line,depths,elev,geom,resolution) values (%d, '\''%s'\'', ARRAY[%f,%f,%f],%.2f,st_geomfromtext('\''POINT(%f %f)'\'',4326),0.05);\n",$2,$0,$23,$24,$25,$22,$4,$3)}' data/vic/global_soil_0.05deg.txt >> vicsouls5km.sql echo "end;" >> vicsoils5km.sql
./bin/psql -d rheas -f vicsoils5km.sql
inside psql
insert into vic.input values (0.05, 'vic/global_snowbands_0.05deg.txt', 'vic/global_lai_0.05deg.txt', 'vic/vic_veglib.txt', 'vic/global_soil_0.05deg.txt', 2, 'vic/dssat.inp.base');