80 likes | 120 Views
Recent STILT work at Jena by Christoph Gerbig and Stefan Körner from Max-Planck-Institute for Biogeochemistry includes modifications to address uncertainties in vertical mixing, CO2 residuals, and more. Planned future STILT developments involve coupling with TM3 global model and nested inversions for improved resolution. Keeping STILT stable during development phases is vital, with suggested methods and a controlled facility setup. Developers are encouraged to follow specific procedures for compiling and verifying code on the Harvard Grid computer.
E N D
Recent STILT work at Jena Christoph Gerbig and Stefan Körner Max-Planck-Institute for Biogeochemistry STILT users at Jena:K. Dhanyalekshmi, Kristina Trusilova, Ronald Macatangay (also at U. Bremen) STILT telecon April 4th 2007
Modifications to STILT • Uncertainties in vertical mixing propagated • Approach similar to Lin and Gerbig ’05 • Uncertainties comparable to uncertainties due to wind errors • Paper in prep. CO2 residuals (model-truth) CO2 biosph (GSB) Day in June 05 CO2 transport error (propagated)
Modifications to STILT • VPRM preprocessor • input: configuration file “namelist” • Grid specifications (projection, limits) • Time frame • output: • EVI, LSWI & Vegetation map at desired resolution for one year (min) • netcdf format • Uses: Modis reflectances MOD09A1 from local database • Tool tells what tiles are required for downloading
Modifications to STILT • Coupling hourly emission fluxes within Trajecvprm() • Emissions in netcdf format • “regridding” to coarser grids online • EDGAR emission converter • get.edgar() to create hourly fluxes at desired resolution (needs upgrade to create ncdf files)
STILT plans at Jena • coupling to TM3 global model as Boundary Condition • TM3 runs using fluxes from inversions, coupled to Biome BGC • Check if nesting improves in forward runs • Nested inversions coupling STILT + TM3 adjoint • European nest for high resolution (~ 20 km) (Kristina Trusilova) • Tall tower data plus global network as input • STILT-WRF-VPRM for mesoscale inversions • At specific tall tower locations in Europe + aircraft data CERES
Future: keeping STILT stable • Keeping STILT stable throughout development phases • Method: set of hymodelc test runs with a result control facility • How (suggestion): • @ Harvard Grid computer • each group provides a (set of) test runs with all required inputs • SK provides a fuzzy checker tool to compare PARTICLE.DATs as core of semi-automated verification (random off, can only check when results are not supposed to change) • developer should run one batch file only which is part of the STILT CVS repository itself • it takes f90 files from user (scp to local machine) and compiles, runs and compares on GRID