Difference between revisions of "CABLE"

From Geos-chem
Jump to: navigation, search
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.