| 19 | == New input file with USDA soil orders == |
| 20 | We use the USDA soil order map of Sun et al. (in revision). The gap filling of the missing values was done like this |
| 21 | (see script /home/users/dgoll/ORC_data/USDA_soilorders/try_2_fill/fill_deserts.ksh): |
| 22 | |
| 23 | We use the soiltext class 6 for cold desert from the input file "soils_param.nc" to fill missing values in these region with Gelisols (5). |
| 24 | After we set cold desert, we use the level 1 of the maxvegetfrac in the input file "PFTmap.nc" to fill missing values of the remaining desert to Aridisols (2). |
| 25 | |
| 26 | In case of the remaining missing points, mostly coast line points, ORCHIDEE will assign the value of the nearest neighbour in the routine "get_soil_orders" in src_stomate_io.f90. |
| 27 | |