wiki:ticket/350/TicketSummary

Version 3 (modified by Martin Dix, 11 months ago) (diff)

--

I've done a one year CABLE run starting from a Jan 1 initial state from the au-aa655 CABLE AMIP run and compared to a similar one year JULES run also starting from an AMIP initial state.

My CABLE run is basically a copy of Jhan's u-ap748 with a couple of tweaks to the reconfiguration to enable it to work with the AMIP file.

Time series of various global means at https://accessdev.nci.org.au/~mrd599/cmm_cable/cmm.html (log in with your NCI account).

Maps of annual mean differences attached. NetCDF files with the monthly data are /short/p66/mrd599/am567a.pm0001.nc (JULES) and ao854a.pm0001.nc (CABLE).

In the annual mean the TOA net flux is 1.2 W/m2 with JULES and -2.3 with CABLE, consistent with the cooling.

Net heat flux in bare soil regions is about 5 W/m2 in CABLE, even though the depth averaged soil temperature doesn't change much over the year.

To further investigate the surface energy balance I set up JULES and CABLE versions of the single column model, suites u-ar806 and u-ar807 respectively.

Used a configuration with bare soil, zero initial soil moisture and relaxation to initial atmospheric state to allow long runs with an extreme diurnal cycle.

In routine stempv, added a calculation of the integrated soil heat content and the initial and corrected ground heat flux.

   ! Integrated heat content
   tmp = 0.
   DO k = 1, ms
      tmp(:) = tmp(:) + ssnow%tgg(:,k)*ssnow%gammzz(:,k)
   END DO
   write(6,*) "COLUMN HEAT", tmp
   write(6,*) "DT", ssnow%tgg(:,1) - tgg_save(:)
   write(6,*) "GA", canopy%ga, canopy%dgdtg, canopy%ga + canopy%dgdtg*(ssnow%tgg(:,1) - tgg_save(:))

This can be compared with the net surface heat flux calculated from the SCM diagnostics surfsw + surf_lw - sens_ht - lat_ht.

The old access 1.3 version of the SCM doesn't have quite the same options for atmospheric relaxation but it's possible to set up a case that has a similarly large diurnal cycle and soil warming.

Averaged over the last 40 days of the runs

10.6 7.3
SCM net 10.87 6.77
Soil change 5.18 6.51
GA 5.87 8.69
GA corrected 5.22 9.34

In access 1.3 there was a good agreement between the SCM net flux and the soil temperature change, though for some reason the agreement with GA is poorer (is the soil temperature changing again after stempv?). In 10.6 CABLE seems to be losing about 5 W/m2.

Attachments (27)