Getting met forcing from ACCESS-ESM1.5
Timeseries output from grid-cells is needed. The list of grid-cells is supplied through the domain field defined in the STASHC file for an ACCESS-ESM1.5 run. Examples are DGTS_PT (21 locations, probably mostly flux towers), DCO2TS (locations of atmospheric CO2 monitoring sites), DCO2TSL (similar to DCO2TS but for variables defined on model levels rather than just at the surface).
The required time field in STASHC is TALLTS.
The variables to be output are listed in the STASHC file by section and number. This is what I have currently used to match up with what CABLE requires
| Section | Number | StashMaster | CABLE |
|---|---|---|---|
| 0 | 409 | SURFACE PRESSURE AFTER TIMESTEP | Psurf |
| 1 | 235 | TOTAL DOWNWARD SURFACE SW FLUX | SWdown |
| 2 | 207 | DOWNWARD LW RAD FLUX: SURFACE | LWdown |
| 3 | 230 | 10 METRE WIND SPEED ON C GRID | Wind |
| 3 | 236 | TEMPERATURE AT 1.5M | Tair |
| 3 | 237 | SPECIFIC HUMIDITY AT 1.5M | Qair |
| 5 | 216 | TOTAL PRECIPITATION RATE KG/M2/S | Rainf |
In an ESM1.5 set-up, the output is requested in the STASHC file using syntax like this:
&STREQ IMOD= 1, ISEC=0, ITEM=409, IDOM=33, ITIM=19, IUSE=12 /
In this case the DGTS_PT domain is the 33rd in the file, and the TALLTS time field is the 19th in the file. IUSE determines which output file the data is written to. In this case it goes to *.pg files.
If adding output requests to the STASHC file, remember to increase the number (NUM_REQ) at the top of the file
&STASHNUM NUM_REQ= 280, NUM_DOM=37, NUM_TIM=20, NUM_USE=12 /
Note that the radiation output will be at the 3 hourly radiation timestep. Other output will be half-hourly.
Update June 2024
CABLE standalone needs to be run with the temperature and specific humidity from the lowest ACCESS model level not the 1.5 m values. These stash outputs were used.
| Section | Number | StashMaster | CABLE |
|---|---|---|---|
| 16 | 4 | TEMPERATURE ON THETA LEVELS | Tair |
| 0 | 10 | SPECIFIC HUMIDITY AFTER TIMESTEP | Qair |
Note that since these are multi-level fields, a new output domain is required to specify the level. This has been set up for the pft-focussed 13-site set of grid-cells described in a later post.