run RVIC using sample rasm files throw unexpected errors

[iacopoff]

Hi, I have just installed RVIC on the linux server and I wanted to run the test cases that I have downloaded from RVIC/samples/. I am wondering why, provided that the data inputs are fine, I am getting errors that seems more related to the code as you can see from the print of the log file below, but I may be wrong also in that...

Could you tell where I am making some mistake?

many thanks

iacopo

INFO:init_logger>> -------------------- INITIALIZED RVIC LOG ------------------ INFO:init_logger>> LOG LEVEL: DEBUG INFO:init_logger>> Logging To Console: True INFO:init_logger>> LOG FILE: ./samples/cases/sample_rasm_parameters/logs/RVIC-20180912-123423.log INFO:init_logger>> ----------------------------------------------------------

INFO:gen_uh_init>> plots directory is ./samples/cases/sample_rasm_parameters/plots INFO:gen_uh_init>> logs directory is ./samples/cases/sample_rasm_parameters/logs INFO:gen_uh_init>> params directory is ./samples/cases/sample_rasm_parameters/params INFO:gen_uh_init>> inputs directory is ./samples/cases/sample_rasm_parameters/inputs INFO:gen_uh_init>> aggregated directory is ./samples/cases/sample_rasm_parameters/temp/aggregated INFO:gen_uh_init>> remapped directory is ./samples/cases/sample_rasm_parameters/temp/remapped INFO:gen_uh_init>> Opened Pour Points File: ./samples/cases/sample_rasm_parameters/inputs/rasm_sample_pour_points.csv INFO:gen_uh_init>> Opened UHbox File: ./samples/cases/sample_rasm_parameters/inputs/UH_RASM_hourly.csv DEBUG:read_netcdf>> Reading input data variables: ['Basin_ID', 'Flow_Direction', 'Flow_Distance', 'Land_Mask', 'Source_Area', 'crs', 'lat', 'lon'], from file: ./samples/cases/sample_rasm_parameters/inputs/Wu_routing_inputs_060313.nc DEBUG:gen_uh_init>> Flow Direction inputs came in upside down, flipping everything now. DEBUG:gen_uh_init>> flipping Basin_ID DEBUG:gen_uh_init>> flipping Flow_Direction DEBUG:gen_uh_init>> flipping Flow_Distance DEBUG:gen_uh_init>> flipping Land_Mask DEBUG:gen_uh_init>> flipping Source_Area DEBUG:gen_uh_init>> flipping crs DEBUG:gen_uh_init>> flipping lat DEBUG:gen_uh_init>> flipping lon INFO:gen_uh_init>> Added velocity grid to fdr_data INFO:gen_uh_init>> Added diffusion grid to fdr_data INFO:gen_uh_init>> Opened FDR File: ./samples/cases/sample_rasm_parameters/inputs/Wu_routing_inputs_060313.nc DEBUG:read_netcdf>> Reading input data variables: ['xc', 'yc', 'xv', 'yv', 'mask', 'area', 'frac'], from file: ./samples/cases/sample_rasm_parameters/inputs/domain.lnd.wr50a_ar9v4.100920.nc INFO:gen_uh_init>> Opened Domain File: ./samples/cases/sample_rasm_parameters/inputs/domain.lnd.wr50a_ar9v4.100920.nc INFO:make_agg_pairs>> Finding addresses now... INFO:latlon2yx>> adjusted grid lon to astronomical conventions INFO:make_agg_pairs>> ------------------ SUMMARY OF MakeAggPairs ------------------ INFO:make_agg_pairs>> NUMBER OF POUR POINTS IN INPUT LIST: 10 INFO:make_agg_pairs>> NUMBER OF POINTS TO AGGREGATE TO: 5 INFO:make_agg_pairs>> NUMBER OF POUR POINTS AGGREGATED: 10 INFO:make_agg_pairs>> EFFECIENCY OF: 100.00 % INFO:make_agg_pairs>> UNASSIGNED POUR POINTS: 0

INFO:make_agg_pairs>> -------------------------------------------------------------

INFO:gen_uh_init>> Finished making agg pairs of pour points and outlet grid cells DEBUG:gen_uh_init>> OrderedDict([(46050, -- Point --
name: outlet_58.4107_93.4308 lat: 58.4107 lon: 93.4308 domy: 167 domx: 125 routy: -9999 routx: -9999 ), (27396, -- Point --
name: outlet83.7711-42.9300 lat: 83.7711 lon: -42.9300 domy: 99 domx: 171 routy: -9999 routx: -9999 ), (27673, -- Point --
name: outlet83.0264-36.9934 lat: 83.0264 lon: -36.9934 domy: 100 domx: 173 routy: -9999 routx: -9999 ), (25467, -- Point --
name: outlet83.3996-75.7979 lat: 83.3996 lon: -75.7979 domy: 92 domx: 167 routy: -9999 routx: -9999 ), (25468, -- Point --
name: outlet83.1114-72.8403 lat: 83.1114 lon: -72.8403 domy: 92 domx: 168 routy: -9999 routx: -9999 )]) INFO:gen_uh_init>> Finished with gen_uh_init INFO:gen_uh_init>> -------------------------------------------------------------

INFO:gen_uh_run>> Running outlet cell id 46050 INFO:gen_uh_run>> On pour_point #1 of 1 INFO:rout>> Starting routing program for point: Point(outlet_58.3500_93.5500, lat:58.3500, lon:93.5500, y:167, x:125) INFO:rout>> Input Latitude: 58.350000 INFO:rout>> Input Longitude: 93.550000 INFO:rout>> Global Basid ID: 35684 DEBUG:rout>> Grid cells in subset: 204402 ERROR:write>> Traceback (most recent call last): ERROR:write>> File "./lib/python/rvic-1.1.1-py3.6-linux-x86_64.egg/EGG-INFO/scripts/rvic", line 61, in ERROR:write>> main() ERROR:write>> File "./lib/python/rvic-1.1.1-py3.6-linux-x86_64.egg/EGG-INFO/scripts/rvic", line 43, in main ERROR:write>> parameters(args.config_file, numofproc=args.numofproc) ERROR:write>> File "/projects/mas1261/wp3/RVIC/lib/python/rvic-1.1.1-py3.6-linux-x86_64.egg/rvic/parameters.py", line 90, in parameters ERROR:write>> config_dict, directories) ERROR:write>> File "/projects/mas1261/wp3/RVIC/lib/python/rvic-1.1.1-py3.6-linux-x86_64.egg/rvic/parameters.py", line 454, in gen_uh_run ERROR:write>> config_dict['ROUTING']) ERROR:write>> File "/projects/mas1261/wp3/RVIC/lib/python/rvic-1.1.1-py3.6-linux-x86_64.egg/rvic/core/make_uh.py", line 84, in rout ERROR:write>> glons=basin['lon']) ERROR:write>> File "/projects/mas1261/wp3/RVIC/lib/python/rvic-1.1.1-py3.6-linux-x86_64.egg/rvic/core/utilities.py", line 36, in latlon2yx ERROR:write>> if (plons.max() > 180): ERROR:write>> AttributeError ERROR:write>> : ERROR:write>> 'float' object has no attribute 'max'

[solomonsg]

Hi,

How did you solve this issue? I am getting the same error.

Thanks Solomon

[iacopoff]

Hi no I gave up, I am using this tool: https://github.com/nmizukami/mizuRoute