search

openwfm / wrf-fire-matlab

Contents of the directory other/Matlab filtered from the original wrf-fire repository when it moved to WRF-SFIRE
4 stars 0 forks source link

Reconcile fuel diagnostics #5

Open janmandel opened 3 years ago

janmandel commented 3 years ago

The last step of the diagnostics described at https://wiki.openwfm.org/wiki/How_to_diagnose_fuel_properties_in_WRF-Fire#Step-by-step_instructions gives an error. See also openwfm/WRF-SFIRE#40.

''Run WRF-SFIRE, you can kill it right after the the first time step. There will be file fuels.m created in the current directory. Clone [1] and copy the file fuels.m to the directory created. cd in MATLAB to the directory and type startup Type fuels. This will create variable fuel. Say you want to check fuel 3. Type fuel(3) to see what is there. Type plot_fuel(fuel(3)). Type big(check_ros(fuel)) to make sure the calculation in Matlab is the same as in WRF-SFIRE, which will compare the calculation with ROS on a number of data points stored in file fuels.m. It should return rounding error only (less than 1e-4)'' However, it returns 0.5696 even if the error checks in plot_fuel on individual fuels are fine.

janmandel commented 3 years ago

The difference is for fuel 4 only. After x=check_ros(fuel) large numbers are only in x(4,:,:).