Difference between revisions of "CABLE"
| Line 3: | Line 3: | ||
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. | 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 | + | 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 points. |
| − | 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. |
| + | |||
| + | The mapping of the 62433 land points to lat/lon is provided in the file latlon_cable.csv, in the format (point lat lon). | ||
In the GEOS-Chem code, I am modifying co2_mod.f, and have added an option "LBIOCABLE" and subroutine "READ_BIO_CABLE" | 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. | 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:49, 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 points.
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.
The mapping of the 62433 land points to lat/lon is provided in the file latlon_cable.csv, in the format (point lat lon).
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.
