- Timestamp:
- 2011-06-17T14:02:17+02:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/ORCHIDEE_EXT/ORCHIDEE/src_stomate/stomate_turnover.f90
r64 r257 196 196 tl(:)*tl(:) * senescence_temp(j,3) 197 197 198 WHERE ( ( biomass(:,j,ileaf) .GT. 0.0) .AND. &198 WHERE ( ( biomass(:,j,ileaf) .GT. zero ) .AND. & 199 199 ( leaf_meanage(:,j) .GT. min_leaf_age_for_senescence(j) ) .AND. & 200 200 ( t2m_month(:) .LT. t_crit(:) ) .AND. ( t2m_week(:) .LT. t2m_month(:) ) ) … … 215 215 nosenescence_hum(j) ) 216 216 217 WHERE ( ( biomass(:,j,ileaf) .GT. 0.0) .AND. &217 WHERE ( ( biomass(:,j,ileaf) .GT. zero ) .AND. & 218 218 ( leaf_meanage(:,j) .GT. min_leaf_age_for_senescence(j) ) .AND. & 219 219 ( moiavail_week(:,j) .LT. moiavail_crit(:) ) ) … … 240 240 241 241 ! critical temperature for senescence may depend on long term annual mean temperature 242 WHERE ( ( biomass(:,j,ileaf) .GT. 0.0) .AND. &242 WHERE ( ( biomass(:,j,ileaf) .GT. zero ) .AND. & 243 243 ( leaf_meanage(:,j) .GT. min_leaf_age_for_senescence(j) ) .AND. & 244 244 ( ( moiavail_week(:,j) .LT. moiavail_crit(:) ) .OR. & … … 319 319 turnover(:,j,ifruit) = biomass(:,j,ifruit) * dt / turnover_time(:,j) 320 320 ELSEWHERE 321 turnover(:,j,ileaf)= 0.0322 turnover(:,j,isapabove) = 0.0323 turnover(:,j,iroot) = 0.0324 turnover(:,j,ifruit) = 0.0321 turnover(:,j,ileaf)= zero 322 turnover(:,j,isapabove) = zero 323 turnover(:,j,iroot) = zero 324 turnover(:,j,ifruit) = zero 325 325 ENDWHERE 326 326 biomass(:,j,ileaf) = biomass(:,j,ileaf) - turnover(:,j,ileaf) … … 364 364 365 365 DO m = 1, nleafages 366 turnover_rate(:) = 0366 turnover_rate(:) = zero 367 367 WHERE ( leaf_age(:,j,m) .GT. leaf_age_crit(:,j)/2. ) 368 368 … … 454 454 leaf_frac(:,j,m) = ( leaf_frac(:,j,m)*lm_old(:) + delta_lm(:,m) ) / biomass(:,j,ileaf) 455 455 ELSEWHERE 456 leaf_frac(:,j,m) = 0.0456 leaf_frac(:,j,m) = zero 457 457 ENDWHERE 458 458 … … 489 489 ! check whether we shed the remaining leaves 490 490 491 WHERE ( ( biomass(:,j,ileaf) .GT. 0.0) .AND. senescence(:,j) .AND. &491 WHERE ( ( biomass(:,j,ileaf) .GT. zero ) .AND. senescence(:,j) .AND. & 492 492 ( biomass(:,j,ileaf) .LT. (lai_initmin(j) / 2.)/sla(j) ) ) 493 493 … … 498 498 turnover(:,j,ifruit) = turnover(:,j,ifruit) + biomass(:,j,ifruit) 499 499 500 biomass(:,j,ileaf) = 0.0501 biomass(:,j,iroot) = 0.0502 biomass(:,j,ifruit) = 0.0500 biomass(:,j,ileaf) = zero 501 biomass(:,j,iroot) = zero 502 biomass(:,j,ifruit) = zero 503 503 504 504 505 505 506 506 ! reset leaf age 507 leaf_meanage(:,j) = 0.0507 leaf_meanage(:,j) = zero 508 508 509 509 ENDWHERE … … 519 519 ! Shed the remaining leaves if LAI very low. 520 520 521 WHERE ( ( biomass(:,j,ileaf) .GT. 0.0) .AND. senescence(:,j) .AND. &521 WHERE ( ( biomass(:,j,ileaf) .GT. zero ) .AND. senescence(:,j) .AND. & 522 522 ( biomass(:,j,ileaf) .LT. (lai_initmin(j) / 2.)/sla(j) )) 523 523 … … 529 529 turnover(:,j,ifruit) = turnover(:,j,ifruit) + biomass(:,j,ifruit) 530 530 531 biomass(:,j,ileaf) = 0.0532 biomass(:,j,isapabove) = 0.0533 biomass(:,j,iroot) = 0.0534 biomass(:,j,ifruit) = 0.0531 biomass(:,j,ileaf) = zero 532 biomass(:,j,isapabove) = zero 533 biomass(:,j,iroot) = zero 534 biomass(:,j,ifruit) = zero 535 535 536 536 537 537 538 538 ! reset leaf age 539 leaf_meanage(:,j) = 0.0539 leaf_meanage(:,j) = zero 540 540 541 541 ENDWHERE … … 551 551 WHERE ( shed_rest(:) ) 552 552 553 leaf_age(:,j,m) = 0.0554 leaf_frac(:,j,m) = 0.0553 leaf_age(:,j,m) = zero 554 leaf_frac(:,j,m) = zero 555 555 556 556 ENDWHERE … … 679 679 hw_new(:) = biomass(:,j,iheartabove) + biomass(:,j,iheartbelow) 680 680 681 WHERE ( hw_new(:) .GT. 0.0)681 WHERE ( hw_new(:) .GT. zero ) 682 682 683 683 age(:,j) = age(:,j) * hw_old(:)/hw_new(:)
Note: See TracChangeset
for help on using the changeset viewer.