Changes between Version 18 and Version 19 of DevelopmentActivities/ORCHIDEE-MICT-IMBALANCE-P/Evaluation


Ignore:
Timestamp:
2015-09-15T14:54:19+02:00 (9 years ago)
Author:
ajornet
Comment:

--

Legend:

Unmodified
Added
Removed
Modified

  • DevelopmentActivities/ORCHIDEE-MICT-IMBALANCE-P/Evaluation

    v18 v19  
    229229 
    230230Albert: done in r2912. 
     231 
     232== Chao == 
     233 
     234Please put these changes in the '''src_stomate/Grassland_Management.f90'''. They're used to protect the division-by-zero from happening. But I don't think this correction is related negative fuel problem. But it could lead to Nan value in fuel if you activate the GM (which is not your case as I checked the run.def you sent before). 
     235 
     236From  
     237{{{ 
     238!spitfire 
     239fuel_all_type(:,:,:) = fuel_1hr(:,:,:) + fuel_10hr(:,:,:) + & 
     240                            fuel_100hr(:,:,:) + fuel_1000hr(:,:,:) 
     241fuel_type_frac(:,:,:,1) = fuel_1hr(:,:,:)/fuel_all_type(:,:,:) 
     242fuel_type_frac(:,:,:,2) = fuel_10hr(:,:,:)/fuel_all_type(:,:,:) 
     243fuel_type_frac(:,:,:,3) = fuel_100hr(:,:,:)/fuel_all_type(:,:,:) 
     244fuel_type_frac(:,:,:,4) = fuel_1000hr(:,:,:)/fuel_all_type(:,:,:) 
     245!end spitfire    
     246}}} 
     247 
     248 
     249to 
     250{{{ 
     251!spitfire 
     252fuel_type_frac(:,:,:,:) = zero 
     253fuel_all_type(:,:,:) = fuel_1hr(:,:,:) + fuel_10hr(:,:,:) + & 
     254                            fuel_100hr(:,:,:) + fuel_1000hr(:,:,:) 
     255WHERE (fuel_all_type(:,:,:) .GT. min_stomate) 
     256   fuel_type_frac(:,:,:,1) = fuel_1hr(:,:,:)/fuel_all_type(:,:,:) 
     257   fuel_type_frac(:,:,:,2) = fuel_10hr(:,:,:)/fuel_all_type(:,:,:) 
     258   fuel_type_frac(:,:,:,3) = fuel_100hr(:,:,:)/fuel_all_type(:,:,:) 
     259   fuel_type_frac(:,:,:,4) = fuel_1000hr(:,:,:)/fuel_all_type(:,:,:) 
     260ENDWHERE 
     261!end spitfire    
     262}}}