Difference between revisions of "CABLE"
Line 6: | Line 6: | ||
Of course, this creates the headache of mapping the flux to the correct grid point, especially as the lat/lon coordinates are not given for each netCDF file, except for the original 0.5x0.5° resolution. So I will regrid to 0.5x0.5° (basically map the flux to match the corresponding points if it is a land point, or set to 0 if not). Then I can worry about aggregating to other resolutions. | Of course, this creates the headache of mapping the flux to the correct grid point, especially as the lat/lon coordinates are not given for each netCDF file, except for the original 0.5x0.5° resolution. So I will regrid to 0.5x0.5° (basically map the flux to match the corresponding points if it is a land point, or set to 0 if not). Then I can worry about aggregating to other resolutions. | ||
+ | |||
+ | In the GEOS-Chem code, I am modifying co2_mod.f, and have added an option "LBIOCABLE" and subroutine "READ_BIO_CABLE" | ||
+ | |||
+ | At present I have made it is far as where I need to build a filename to read the fluxes. Until I set the re-gridded files up, that's all I can probably do. |
Revision as of 06:04, 24 March 2015
CABLE - CSIRO Atmosphere Biosphere Land Exchange model
I have fluxes provided by Chris Lu and Rachel Law that I want to use to investigate differences in NEE over Australia, particularly in comparison to CarbonTracker.
These fluxes are currently on a land only grid, with hourly resolution. The grid is 0.5x0.5°, so potentially 720 x 360 (=259200 points). File size is saved by only including land point fluxes (it is a biosphere model, after all), leaving 62433 point.
Of course, this creates the headache of mapping the flux to the correct grid point, especially as the lat/lon coordinates are not given for each netCDF file, except for the original 0.5x0.5° resolution. So I will regrid to 0.5x0.5° (basically map the flux to match the corresponding points if it is a land point, or set to 0 if not). Then I can worry about aggregating to other resolutions.
In the GEOS-Chem code, I am modifying co2_mod.f, and have added an option "LBIOCABLE" and subroutine "READ_BIO_CABLE"
At present I have made it is far as where I need to build a filename to read the fluxes. Until I set the re-gridded files up, that's all I can probably do.