| 170 | |
| 171 | == merge de slowproc == |
| 172 | * à vérifier avec les autres driver : déplacement/redéfinition de diaglev + validation stomate activé |
| 173 | * ancien bug toujours présent dans la version LMD dans la définition des soiltiles : |
| 174 | {{{ |
| 175 | ! Soiltiles are only used in hydrol, but we fix them in here because some time it might depend |
| 176 | ! on a changing vegetation (but then some adaptation should be made to hydrol) and be also used |
| 177 | ! in the other modules to perform separated energy balances |
| 178 | njsc(:) = 0 |
| 179 | soiltile(:,:) = zero |
| 180 | DO ji = 1, kjpindex |
| 181 | njsc(ji) = MAXLOC(soilclass(ji,:),1) |
| 182 | soiltile(:,1) = SUM(frac_nobio(ji,:)) |
| 183 | ENDDO |
| 184 | DO jv = 1, nvm |
| 185 | jst = pref_soil_veg(jv) |
| 186 | DO ji = 1, kjpindex |
| 187 | soiltile(ji,jst) = soiltile(ji,jst) + veget(ji,jv) |
| 188 | ENDDO |
| 189 | ENDDO |
| 190 | }}} |
| 191 | On prend alors ici toujours le frac_nobio du dernier point. Le code corrigé est : |
| 192 | {{{ |
| 193 | ! Soiltiles are only used in hydrol, but we fix them in here because some time it might depend |
| 194 | ! on a changing vegetation (but then some adaptation should be made to hydrol) and be also used |
| 195 | ! in the other modules to perform separated energy balances |
| 196 | njsc(:) = 0 |
| 197 | soiltile(:,:) = zero |
| 198 | DO ji = 1, kjpindex |
| 199 | njsc(ji) = MAXLOC(soilclass(ji,:),1) |
| 200 | soiltile(ji,1) = SUM(frac_nobio(ji,:)) |
| 201 | ENDDO |
| 202 | DO jv = 1, nvm |
| 203 | jst = pref_soil_veg(jv) |
| 204 | DO ji = 1, kjpindex |
| 205 | soiltile(ji,jst) = soiltile(ji,jst) + veget(ji,jv) |
| 206 | ENDDO |
| 207 | ENDDO |
| 208 | }}} |
| 209 | }}} |