WAM2layers / Moisture_tracking_intercomparison

Apache License 2.0
4 stars 0 forks source link

Update combine_data.py to read WRF-WVT results from netCDF #85

Closed crespoalfredo closed 3 months ago

crespoalfredo commented 4 months ago

Now WRF-WVT results can are readed from netCDF files. This allows the same functions to be used in Functions.py as with the other methods, but requires the netCDF files which are in the "results WRF-WVT" folders since the beginning of this month (before, there were only CSV files). Observe that you may want to omit WRF-WVT results when representing gridded results.

crespoalfredo commented 2 months ago

Hi, sorry for the late replay, but I completely forgot to answer after my holidays.

You are right, these numbers should be the same, and the correct one is 119.66 mm. There was a problem when generating the files Australia_WRF-WVTs.nc (the same for Pakistan and Scotland), because I used the area of the sink region computed in the WRF domain, which is slightly smaller. Because of this, the values of the variable "moisture_sources_grid" in these netCDF files are a bit lower than they should be (around a 5 %).

I will update these files in the drive. This only affects the absolute results, all other results calculated from fractional moisture sources are correct (for example, the barplots for the fraction of moisture uptake (%) will not change).

crespoalfredo commented 2 months ago

Hi, sorry for the late replay, but I completely forgot to answer after my holidays.

You are right, these numbers should be the same, and the correct one is 119.66 mm. There was a problem when generating the files Australia_WRF-WVTs.nc (the same for Pakistan and Scotland), because I used the area of the sink region computed in the WRF domain, which is slightly smaller. Because of this, the values of the variable "moisture_sources_grid" in these netCDF files are a bit lower than they should be (around a 5 %).

I will update these files in the drive. This only affects the absolute results, all other results calculated from fractional moisture sources are correct (for example, the barplots for the fraction of moisture uptake (%) will not change).

Sorry, the name of the file is AustraliaCase_Final.nc, I confused it with the name of the same file in my computer.