| 107 | = Comparison of the old and new drivers = |
| 108 | |
| 109 | ORCHIDEE has now 2 drivers which are quite different in their conception. As discussed above the temporal interpolation is a little of a challenge so it is worth verifying how well it is done bu both driver and how important the differences are. This effort is undertaken here. |
| 110 | |
| 111 | The issue is particularly critical with the global forcing which are available at 3 or 6 hour interval and any shift in the date attributed to them impacts the interpolation in time and the variables given to ORCHIDEE. This issue is largely attenuated for point forcing where atmospheric conditions are available at a frequency higher than hourly. Thus the test here are carried out over the 5 global forcing data sets which are currently available to the community : |
| 112 | |
| 113 | - CRU-NCEP : This is the NCEP re-analysis corrected by Nicolas Viovy using CRU data. The version 5.4 which is used here only has data all the 6 hours. |
| 114 | - CSWP : is the forcing prepared by U Tokyo for the international GSWP excise and is based on C20C re-analysis of ECMWF. |
| 115 | - PGF : Is the Princeton forcing prepared which combines the NCEP re-analysis, CRU and GPCC data. |
| 116 | - WFD : Based on ERA-40 an corrected in the same way as WFDEI (This is the WATCH data set). |
| 117 | - WFDEI : Based on ERA-I and corrected by the UKMO using the independent CRU information (An evolution of the WATCH data set) |
| 118 | |
| 119 | == Modelling set-up == |
| 120 | |
| 121 | The ORCHIDEE trunk version of May 2016 is used for these tests. All forcing are used at a 900s time-step and the model is configured to output through XIOS at this same time step. This ensures that all the post-processing done XIOS on the forcing variables ORCHIDEE outputs is negligible compared to the assumptions on the time axis in the forcing. |
| 122 | |
| 123 | From the forcing files a small region (lon=-10.3:5.3 lat=35.3:45.3) was extracted to reduce computing time and ensure that the original data set could easily be compared to the model output. |
| 124 | |
| 125 | CRU-NCEP, WFD and WFDEI came already with the attributes describing the time averaging of the fluxes set. Two cases are distinguished there and they are : |
| 126 | - Fluxes : time: mean(end) (SWdown LWdown Rainf Snowf) |
| 127 | - Scalar : time: instantaneous (Tair Qair PSurf Wind) |
| 128 | |
| 129 | For GSWP and PGF this information was added after extracting the region of interest with the following commands : |
| 130 | - GSWP Fluxes : |
| 131 | {{{ |
| 132 | for v in $FLUX ; do |
| 133 | ncatted -h -a cell_methods,$v,o,c,"time: mean(center)" ES_GSWP3_${i}.nc |
| 134 | done |
| 135 | }}} |
| 136 | |
| 137 | - GSWP Scalar : |
| 138 | |
| 139 | {{{ |
| 140 | for v in $SCALAR ; do |
| 141 | ncatted -h -a cell_methods,$v,o,c,"time: instantaneous" ES_GSWP3_${i}.nc |
| 142 | done |
| 143 | }}} |
| 144 | |
| 145 | - PGF Fluxes : |
| 146 | |
| 147 | {{{ |
| 148 | for v in $FLUX ; do |
| 149 | ncatted -h -a cell_methods,$v,o,c,"time: mean(start)" ES_PGF_${i}.nc |
| 150 | done |
| 151 | }}} |
| 152 | |
| 153 | - PGF Scalar : |
| 154 | |
| 155 | {{{ |
| 156 | for v in $SCALAR ; do |
| 157 | ncatted -h -a cell_methods,$v,o,c,"time: instantaneous" ES_PGF_${i}.nc |
| 158 | done |
| 159 | }}} |
| 160 | |
| 161 | Furthermore for the fluxes in GSWP we tested all three possible hypothesis : |
| 162 | - GSWPe : "time: mean(end)" |
| 163 | - GSWPc : "time: mean(center)" |
| 164 | - GSWMs : "time: mean(start)" |
| 165 | |
| 166 | == Analysis == |
| 167 | |
| 168 | == Results == |
| 169 | |
| 170 | - WFDEI |
| 171 | - WFD |
| 172 | - CRU-NCEP |
| 173 | - GSWP |
| 174 | - PGF |
| 175 | |
| 176 | == Conclusions == |