Changeset 136
- Timestamp:
- 2011-03-08T10:04:07+01:00 (14 years ago)
- Location:
- branches/ORCHIDEE_EXT/ORCHIDEE/src_stomate
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/ORCHIDEE_EXT/ORCHIDEE/src_stomate/lpj_fire.f90
r64 r136 40 40 41 41 SUBROUTINE fire (npts, dt, litterpart, & 42 litterhum_daily, t2m_daily, lignin_struc, &42 litterhum_daily, t2m_daily, lignin_struc,veget_max, & 43 43 fireindex, firelitter, biomass, ind, & 44 44 litter, dead_leaves, bm_to_litter, black_carbon, & … … 64 64 ! (gC/(m**2 of ground)) 65 65 REAL(r_std), DIMENSION(npts,nvm,nlevs), INTENT(in) :: lignin_struc 66 ! NVui & DS Add veget_max for externalisation 02/03/2011 67 REAL(r_std), DIMENSION(npts,nvm), INTENT(in) :: veget_max 66 68 67 69 ! 0.2 modified fields … … 374 376 375 377 black_carbon(:) = & 376 black_carbon(:) + bcfrac(:) * residue(:) 378 black_carbon(:) + bcfrac(:) * residue(:)* veget_max(:,j) 379 377 380 !MM in SZ black_carbon(:) + bcfrac(:) * residue(:) * biomass(:,j,k) 381 !NVu & DS 02032011: multiply black_carbon by veget_max for the externalised version 382 !black_carbon(:) + bcfrac(:) * residue(:) 378 383 379 384 ! 4.3.2.3 The rest (largest part) of the residue becomes litter. … … 468 473 469 474 black_carbon(:) = & 470 black_carbon(:) + bcfrac(:) * residue(:) 475 black_carbon(:) + bcfrac(:) * residue(:)* veget_max(:,j) 471 476 !MM in SZ black_carbon(:) + bcfrac(:) * residue(:) * litter(:,iwoody,j,iabove) 477 !NVu & DS 02032011: multiply black_carbon by veget_max for the externalised version 478 !black_carbon(:) = black_carbon(:) + bcfrac(:) * residue(:) 472 479 473 480 ! 5.2.6 The rest (largest part) of the residue remains litter. Remaining litter -
branches/ORCHIDEE_EXT/ORCHIDEE/src_stomate/stomate_lpj.f90
r64 r136 502 502 503 503 CALL fire (npts, dt_days, litterpart, & 504 litterhum_daily, t2m_daily, lignin_struc, &504 litterhum_daily, t2m_daily, lignin_struc,veget_max, & 505 505 fireindex, firelitter, biomass, ind, & 506 506 litter, dead_leaves, bm_to_litter, black_carbon, &
Note: See TracChangeset
for help on using the changeset viewer.