Context Navigation

source: branches/publications/ORCHIDEE-PEAT_r5488/src_stomate/lpj_cover.f90 @ 6890

Last change on this file since 6890 was 6890, checked in by bertrand.guenet, 4 years ago
update of forcesoil and lpj_cover
Property svn:keywords set to `HeadURL Date Author Revision`
File size: 23.1 KB

Rev	Line
[947]	1	! =================================================================================================================================
	2	! MODULE : lpj_cover
[8]	3	!
[947]	4	! CONTACT : orchidee-help _at_ ipsl.jussieu.fr
[8]	5	!
[947]	6	! LICENCE : IPSL (2006)
	7	! This software is governed by the CeCILL licence see ORCHIDEE/ORCHIDEE_CeCILL.LIC
	8	!
	9	!>\BRIEF Recalculate vegetation cover and LAI
	10	!!
	11	!!\n DESCRIPTION : None
	12	!!
[2268]	13	!! RECENT CHANGE(S) : Including permafrost carbon
[947]	14	!!
	15	!! REFERENCE(S) :
	16	!! Sitch, S., B. Smith, et al. (2003), Evaluation of ecosystem dynamics,
	17	!! plant geography and terrestrial carbon cycling in the LPJ dynamic
	18	!! global vegetation model, Global Change Biology, 9, 161-185.\n
	19	!! Smith, B., I. C. Prentice, et al. (2001), Representation of vegetation
	20	!! dynamics in the modelling of terrestrial ecosystems: comparing two
	21	!! contrasting approaches within European climate space,
	22	!! Global Ecology and Biogeography, 10, 621-637.\n
	23	!!
	24	!! SVN :
	25	!! $HeadURL$
	26	!! $Date$
	27	!! $Revision$
	28	!! \n
	29	!_ ================================================================================================================================
	30
[8]	31	MODULE lpj_cover
	32
	33	! modules used:
	34
[1536]	35	USE ioipsl_para
[511]	36	USE stomate_data
	37	USE pft_parameters
[2268]	38	USE constantes_soil_var
[8]	39
	40	IMPLICIT NONE
	41
	42	! private & public routines
	43
	44	PRIVATE
	45	PUBLIC cover
	46
	47	CONTAINS
	48
[947]	49	!! ================================================================================================================================
	50	!! SUBROUTINE : lpj_cover
	51	!!
	52	!>\BRIEF Recalculate vegetation cover and LAI
	53	!!
	54	!!\n DESCRIPTION : Veget_max is first renewed here according to newly calculated foliage biomass in this calculation step
	55	!! Then, litter, soil carbon, and biomass are also recalcuted with taking into account the changes in Veget_max (i.e. delta_veg)
	56	!! Grid-scale fpc (foliage projected coverage) is calculated to obtain the shadede ground area by leaf's light capture
	57	!! Finally, grid-scale fpc is adjusted not to exceed 1.0
	58	!!
	59	!! RECENT CHANGE(S) : None
	60	!!
	61	!! MAIN OUTPUT VARIABLE(S) : ::lai (leaf area index, @tex $(m^2 m^{-2})$ @endtex),
	62	!! :: veget (fractional vegetation cover, unitless)
	63	!!
	64	!! REFERENCE(S) : None
	65	!!
	66	!! FLOWCHART :
	67	!! \latexonly
	68	!! \includegraphics[scale=0.5]{lpj_cover_flowchart.png}
	69	!! \endlatexonly
	70	!! \n
	71	!_ ================================================================================================================================
	72
[8]	73	SUBROUTINE cover (npts, cn_ind, ind, biomass, &
[2655]	74	veget_max, veget_max_old, lai, &
	75	litter, litter_avail, litter_not_avail, carbon, &
	76	fuel_1hr, fuel_10hr, fuel_100hr, fuel_1000hr, &
	77	turnover_daily, bm_to_litter, &
[2590]	78	co2_to_bm, co2_fire, resp_hetero, resp_maint, resp_growth, gpp_daily, &
[2268]	79	deepC_a, deepC_s, deepC_p)
[8]	80
[947]	81	!! 0. Variable and parameter declaration
[8]	82
[947]	83	!! 0.1 Input variables
[8]	84
[947]	85	INTEGER(i_std), INTENT(in) :: npts !! Domain size (unitless)
	86	REAL(r_std), DIMENSION(npts,nvm), INTENT(in) :: cn_ind !! Crown area
	87	!! @tex $(m^2)$ @endtex per individual
	88	REAL(r_std), DIMENSION(npts,nvm), INTENT(in) :: ind !! Number of individuals
	89	!! @tex $(m^{-2})$ @endtex
	90	REAL(r_std), DIMENSION(npts,nvm), INTENT(in) :: veget_max_old !! "Maximal" coverage fraction of a PFT (LAI->
	91	!! infinity) on ground at beginning of time
	92	!! 0.2 Output variables
[8]	93
[947]	94	!! 0.3 Modified variables
[8]	95
[1536]	96	REAL(r_std), DIMENSION(npts,nvm), INTENT(inout) :: lai !! Leaf area index OF AN INDIVIDUAL PLANT
	97	!! @tex $(m^2 m^{-2})$ @endtex
	98	REAL(r_std), DIMENSION(npts,nlitt,nvm,nlevs,nelements), INTENT(inout) :: litter !! Metabolic and structural litter, above and
[2655]	99	REAL(r_std), DIMENSION(npts,nvm,nlitt,nelements), INTENT(inout) :: fuel_1hr
	100	REAL(r_std), DIMENSION(npts,nvm,nlitt,nelements), INTENT(inout) :: fuel_10hr
	101	REAL(r_std), DIMENSION(npts,nvm,nlitt,nelements), INTENT(inout) :: fuel_100hr
	102	REAL(r_std), DIMENSION(npts,nvm,nlitt,nelements), INTENT(inout) :: fuel_1000hr
[1536]	103	!! below ground @tex $(gC m^{-2})$ @endtex
[2615]	104	REAL(r_std), DIMENSION(npts,nlitt,nvm), INTENT(inout):: litter_avail
	105	REAL(r_std), DIMENSION(npts,nlitt,nvm) , INTENT(inout):: litter_not_avail
[1536]	106	REAL(r_std), DIMENSION(npts,ncarb,nvm), INTENT(inout) :: carbon !! Carbon pool: active, slow, or passive
	107	REAL(r_std), DIMENSION(npts,nvm,nparts,nelements), INTENT(inout) :: biomass !! Biomass @tex $(gC m^{-2})$ @endtex
	108	REAL(r_std), DIMENSION(npts,nvm), INTENT(inout) :: veget_max !! "Maximal" coverage fraction of a PFT (LAI->
	109	!! infinity) on ground (unitless)
	110	REAL(r_std), DIMENSION(npts,nvm,nparts,nelements), INTENT(inout) :: turnover_daily !! Turnover rates (gC m^{-2} day^{-1})
	111	REAL(r_std), DIMENSION(npts,nvm,nparts,nelements), INTENT(inout) :: bm_to_litter !! Conversion of biomass to litter
	112	!! @tex $(gC m^{-2} day^{-1})$ @endtex
[2590]	113	REAL(r_std), DIMENSION(npts,nvm), INTENT(inout) :: co2_to_bm !! biomass up take for establishment
	114	REAL(r_std), DIMENSION(npts,nvm), INTENT(inout) :: co2_fire
	115	REAL(r_std), DIMENSION(npts,nvm), INTENT(inout) :: resp_hetero
	116	REAL(r_std), DIMENSION(npts,nvm), INTENT(inout) :: resp_maint
	117	REAL(r_std), DIMENSION(npts,nvm), INTENT(inout) :: resp_growth
	118	REAL(r_std), DIMENSION(npts,nvm), INTENT(inout) :: gpp_daily
	119
[2268]	120	REAL(r_std), DIMENSION(npts,ndeep,nvm), INTENT(inout) :: deepC_a !! Permafrost soil carbon (g/m**3) active
	121	REAL(r_std), DIMENSION(npts,ndeep,nvm), INTENT(inout) :: deepC_s !! Permafrost soil carbon (g/m**3) slow
	122	REAL(r_std), DIMENSION(npts,ndeep,nvm), INTENT(inout) :: deepC_p !! Permafrost soil carbon (g/m**3) passive
	123
[947]	124	!! 0.4 Local variables
[8]	125
[3026]	126	INTEGER(i_std) :: i,j,k,m !! Index (unitless)
	127	REAL(r_std), DIMENSION(npts,nlitt,nlevs,nelements) :: dilu_lit !! Litter dilution @tex $(gC m^{-2})$ @endtex
	128	REAL(r_std), DIMENSION(npts,ncarb) :: dilu_soil_carbon !! Soil Carbondilution
	129	!! @tex $(gC m^{-2})$ @endtex
	130	REAL(r_std), DIMENSION(npts,nparts,nelements) :: dilu_bio !! Biomass dilution
	131	REAL(r_std), DIMENSION(npts) :: dilu_TCarbon
	132	REAL(r_std), DIMENSION(npts,nparts,nelements) :: dilu_turnover_daily
	133	REAL(r_std), DIMENSION(npts,nparts,nelements) :: dilu_bm_to_litter
	134	REAL(r_std), DIMENSION(npts) :: dilu_co2flux_new
	135	REAL(r_std), DIMENSION(npts) :: dilu_gpp_daily
	136	REAL(r_std), DIMENSION(npts) :: dilu_resp_growth
	137	REAL(r_std), DIMENSION(npts) :: dilu_resp_maint
	138	REAL(r_std), DIMENSION(npts) :: dilu_resp_hetero
	139	REAL(r_std), DIMENSION(npts) :: dilu_co2_to_bm
	140	REAL(r_std), DIMENSION(npts) :: dilu_co2_fire
	141	REAL(r_std), DIMENSION(npts,nvm) :: TCarbon
	142	REAL(r_std), DIMENSION(npts,nvm) :: co2flux_new
	143	REAL(r_std), DIMENSION(npts,nvm) :: co2flux_old
[2615]	144	REAL(r_std), DIMENSION(npts,ncarb,nvm) :: carbon_old
[2268]	145	REAL(r_std),DIMENSION(npts,ndeep,ncarb) :: dilu_soil_carbon_vertres !!vertically-resolved Soil Carbondilution (gC/mÂ²)
	146
[947]	147	REAL(r_std), DIMENSION(nvm) :: delta_veg !! Conversion factors (unitless)
	148	REAL(r_std), DIMENSION(nvm) :: reduct !! Conversion factors (unitless)
	149	REAL(r_std) :: delta_veg_sum !! Conversion factors (unitless)
	150	REAL(r_std) :: diff !! Conversion factors (unitless)
	151	REAL(r_std) :: sr !! Conversion factors (unitless)
	152	REAL(r_std), DIMENSION(npts) :: frac_nat !! Conversion factors (unitless)
	153	REAL(r_std), DIMENSION(npts) :: sum_vegettree !! Conversion factors (unitless)
	154	REAL(r_std), DIMENSION(npts) :: sum_vegetgrass !! Conversion factors (unitless)
	155	REAL(r_std), DIMENSION(npts) :: sum_veget_natveg !! Conversion factors (unitless)
[2590]	156	REAL(r_std), DIMENSION(npts) :: vartmp !! Temporary variable used to add history
[3026]	157	REAL(r_std), DIMENSION(npts,nlitt,nelements) :: dilu_f1hr !! Litter dilution @tex $(gC m^{-2})$ @endtex
	158	REAL(r_std), DIMENSION(npts,nlitt,nelements) :: dilu_f10hr !! Litter dilution @tex $(gC m^{-2})$ @endtex
	159	REAL(r_std), DIMENSION(npts,nlitt,nelements) :: dilu_f100hr !! Litter dilution @tex $(gC m^{-2})$ @endtex
	160	REAL(r_std), DIMENSION(npts,nlitt,nelements) :: dilu_f1000hr !! Litter dilution @tex $(gC m^{-2})$ @endtex
[947]	161	!_ ================================================================================================================================
[8]	162
[947]	163	!! 1. If the vegetation is dynamic, calculate new maximum vegetation cover for natural plants
	164
[3026]	165	IF ( ok_dgvm ) THEN
[8]	166
[947]	167	!! 1.1 Calculate initial values of vegetation cover
[186]	168	frac_nat(:) = un
	169	sum_veget_natveg(:) = zero
	170	veget_max(:,ibare_sechiba) = un
[3026]	171	co2flux_new = undef
	172	co2flux_old = undef
	173	TCarbon = undef
[186]	174
[2615]	175	carbon_old(:,:,:)=carbon(:,:,:)
	176
[947]	177	DO j = 2,nvm ! loop over PFTs
[6890]	178	IF ( natural(j) .AND. .NOT. pasture(j) .AND. .NOT. is_peat(j)) THEN
[3932]	179
[947]	180	! Summation of individual tree crown area to get total foliar projected coverage
[8]	181	veget_max(:,j) = ind(:,j) * cn_ind(:,j)
[186]	182	sum_veget_natveg(:) = sum_veget_natveg(:) + veget_max(:,j)
	183	ELSE
[947]	184
[186]	185	!fraction occupied by agriculture needs to be substracted for the DGVM
	186	!this is used below to constrain veget for natural vegetation, see below
	187	frac_nat(:) = frac_nat(:) - veget_max(:,j)
	188
[8]	189	ENDIF
	190
[947]	191	ENDDO ! loop over PFTs
[8]	192
[947]	193	DO i = 1, npts ! loop over grid points
[3026]	194
[947]	195	! Recalculation of vegetation projected coverage when ::frac_nat was below ::sum_veget_natveg
	196	! It means that non-natural vegetation will recover ::veget_max as natural vegetation
[186]	197	IF (sum_veget_natveg(i) .GT. frac_nat(i) .AND. frac_nat(i) .GT. min_stomate) THEN
	198
[947]	199	DO j = 2,nvm ! loop over PFTs
[6890]	200	IF( natural(j) .AND. .NOT. pasture(j) .AND. .NOT. is_peat(j)) THEN
[186]	201	veget_max(i,j) = veget_max(i,j) * frac_nat(i) / sum_veget_natveg(i)
	202	ENDIF
[947]	203	ENDDO ! loop over PFTs
[186]	204
	205	ENDIF
[947]	206	ENDDO ! loop over grid points
[3026]	207
[947]	208	! Renew veget_max of bare soil as 0 to difference of veget_max (ibare_sechiba)
	209	! to current veget_max
	210	DO j = 2,nvm ! loop over PFTs
[186]	211	veget_max(:,ibare_sechiba) = veget_max(:,ibare_sechiba) - veget_max(:,j)
[947]	212	ENDDO ! loop over PFTs
[8]	213	veget_max(:,ibare_sechiba) = MAX( veget_max(:,ibare_sechiba), zero )
	214
[947]	215	!! 1.2 Calculate carbon fluxes between PFTs to maintain mass balance
[3275]	216	!! Assure carbon closure when veget_max changes(delta_veg): if veget_max of some PFTs decrease, we use "dilu" to
	217	!! record the corresponding lost in carbon (biomass, litter, soil carbon, gpp, respiration etc.) for
	218	!! these PFTs, and re-allocate "dilu" to those PFTs with increasing veget_max.
[947]	219	DO i = 1, npts ! loop over grid points
	220
[3275]	221	! Calculate the change in veget_max between previous time step and current time step
[186]	222	delta_veg(:) = veget_max(i,:)-veget_max_old(i,:)
	223	delta_veg_sum = SUM(delta_veg,MASK=delta_veg.LT.zero)
[8]	224
[1536]	225	dilu_lit(i,:,:,:) = zero
[2655]	226	dilu_f1hr(i,:,:) = zero
	227	dilu_f10hr(i,:,:) = zero
	228	dilu_f100hr(i,:,:) = zero
	229	dilu_f1000hr(i,:,:) = zero
[186]	230	dilu_soil_carbon(i,:) = zero
[2590]	231	dilu_soil_carbon_vertres(i,:,:) = zero
	232
	233	dilu_bio(i,:,:) = zero
	234	dilu_TCarbon(i)=zero
	235
	236	dilu_turnover_daily(i,:,:)=zero
	237	dilu_bm_to_litter(i,:,:)=zero
	238	dilu_co2flux_new(i)=zero
	239	dilu_gpp_daily(i)=zero
	240	dilu_resp_growth(i)=zero
	241	dilu_resp_maint(i)=zero
	242	dilu_resp_hetero(i)=zero
	243	dilu_co2_to_bm(i)=zero
	244	dilu_co2_fire(i)=zero
	245
[3275]	246	! Calculate TCarbon: total carbon including biomass, litter and soil carbon, as well as "today's" turnover and
	247	! bm_to_litter due to mortality, because today's turnover and bm_to_litter are not yet added into "litter" until tomorrow.
	248	DO j=1, nvm
[2655]	249	TCarbon(i,j)=SUM(biomass(i,j,:,icarbon))+SUM(carbon(i,:,j))+SUM(litter(i,:,j,:,icarbon))+SUM(turnover_daily(i,j,:,icarbon))+SUM(bm_to_litter(i,j,:,icarbon))
[2590]	250	co2flux_old(i,j)=resp_maint(i,j)+resp_growth(i,j)+resp_hetero(i,j)+co2_fire(i,j)-co2_to_bm(i,j)-gpp_daily(i,j)
	251	co2flux_new(i,j)=resp_maint(i,j)+resp_growth(i,j)+resp_hetero(i,j)+co2_fire(i,j)-co2_to_bm(i,j)-gpp_daily(i,j)
	252	ENDDO
	253
[947]	254	DO j=1, nvm ! loop over PFTs
[2268]	255	IF ( delta_veg(j) < -min_stomate ) THEN
[1536]	256	dilu_lit(i,:,:,:) = dilu_lit(i,:,:,:) + delta_veg(j) * litter(i,:,j,:,:) / delta_veg_sum
[2655]	257	dilu_f1hr(i,:,:) = dilu_f1hr(i,:,:) + delta_veg(j) * fuel_1hr(i,j,:,:) / delta_veg_sum
	258	dilu_f10hr(i,:,:) = dilu_f10hr(i,:,:) + delta_veg(j) * fuel_10hr(i,j,:,:) / delta_veg_sum
	259	dilu_f100hr(i,:,:) = dilu_f100hr(i,:,:) + delta_veg(j) * fuel_100hr(i,j,:,:) / delta_veg_sum
	260	dilu_f1000hr(i,:,:) = dilu_f1000hr(i,:,:) + delta_veg(j) * fuel_1000hr(i,j,:,:) / delta_veg_sum
[2590]	261	dilu_soil_carbon(i,:) = dilu_soil_carbon(i,:) + delta_veg(j) * carbon(i,:,j) / delta_veg_sum
	262	dilu_TCarbon(i)= dilu_TCarbon(i) + delta_veg(j) * TCarbon(i,j) / delta_veg_sum
	263	dilu_turnover_daily(i,:,:)=dilu_turnover_daily(i,:,:)+delta_veg(j)*turnover_daily(i,j,:,:)/delta_veg_sum
	264	dilu_bm_to_litter(i,:,:)=dilu_bm_to_litter(i,:,:)+delta_veg(j)*bm_to_litter(i,j,:,:)/delta_veg_sum
	265	dilu_co2flux_new(i)=dilu_co2flux_new(i)+delta_veg(j)*co2flux_old(i,j)/delta_veg_sum
	266	dilu_gpp_daily(i)=dilu_gpp_daily(i)+delta_veg(j)*gpp_daily(i,j)/delta_veg_sum
	267	dilu_resp_growth(i)=dilu_resp_growth(i)+delta_veg(j)*resp_growth(i,j)/delta_veg_sum
	268	dilu_resp_maint(i)=dilu_resp_maint(i)+delta_veg(j)*resp_maint(i,j)/delta_veg_sum
	269	dilu_resp_hetero(i)=dilu_resp_hetero(i)+delta_veg(j)*resp_hetero(i,j)/delta_veg_sum
	270	dilu_co2_to_bm(i)=dilu_co2_to_bm(i)+delta_veg(j)*co2_to_bm(i,j)/delta_veg_sum
	271	dilu_co2_fire(i)=dilu_co2_fire(i)+delta_veg(j)*co2_fire(i,j)/delta_veg_sum
[186]	272	ENDIF
[947]	273	ENDDO ! loop over PFTs
[8]	274
[947]	275	DO j=1, nvm ! loop over PFTs
[186]	276	IF ( delta_veg(j) > min_stomate) THEN
[8]	277
[186]	278	! Dilution of reservoirs
[947]	279	! Recalculate the litter and soil carbon with taking into accout the change in
	280	! veget_max (delta_veg)
[186]	281	! Litter
[1536]	282	litter(i,:,j,:,:)=(litter(i,:,j,:,:) * veget_max_old(i,j) + dilu_lit(i,:,:,:) * delta_veg(j)) / veget_max(i,j)
[2655]	283	fuel_1hr(i,j,:,:)=(fuel_1hr(i,j,:,:) * veget_max_old(i,j) + dilu_f1hr(i,:,:) * delta_veg(j)) / veget_max(i,j)
	284	fuel_10hr(i,j,:,:)=(fuel_10hr(i,j,:,:) * veget_max_old(i,j) + dilu_f10hr(i,:,:) * delta_veg(j)) / veget_max(i,j)
	285	fuel_100hr(i,j,:,:)=(fuel_100hr(i,j,:,:) * veget_max_old(i,j) + dilu_f100hr(i,:,:) * delta_veg(j)) / veget_max(i,j)
	286	fuel_1000hr(i,j,:,:)=(fuel_1000hr(i,j,:,:) * veget_max_old(i,j) + dilu_f1000hr(i,:,:) * delta_veg(j)) / veget_max(i,j)
[2615]	287	!JCADD available and not available litter for grazing
	288	! only not available litter change, available litter will not
	289	! change, because tree litter can not be eaten
	290	IF (is_grassland_manag(j) .AND. is_grassland_grazed(j)) THEN
	291	litter_avail(i,:,j) = litter_avail(i,:,j) * veget_max_old(i,j) / veget_max(i,j)
	292	litter_not_avail(i,:,j) = litter(i,:,j,iabove,icarbon) - litter_avail(i,:,j)
	293	ENDIF
	294	!ENDJCADD
[2655]	295	!IF ( ok_pc ) THEN
	296	! deepC_a(i,:,j)=(deepC_a(i,:,j) * veget_max_old(i,j) + &
	297	! dilu_soil_carbon_vertres(i,:,iactive) * delta_veg(j)) / veget_max(i,j)
	298	! deepC_s(i,:,j)=(deepC_s(i,:,j) * veget_max_old(i,j) + &
	299	! dilu_soil_carbon_vertres(i,:,islow) * delta_veg(j)) / veget_max(i,j)
	300	! deepC_p(i,:,j)=(deepC_p(i,:,j) * veget_max_old(i,j) + &
	301	! dilu_soil_carbon_vertres(i,:,ipassive) * delta_veg(j)) / veget_max(i,j)
	302	!ENDIF
[2590]	303	! Soil carbon
	304	carbon(i,:,j)=(carbon(i,:,j) * veget_max_old(i,j) + dilu_soil_carbon(i,:) * delta_veg(j)) / veget_max(i,j)
[2615]	305	IF ( ok_pc ) THEN
	306	IF (carbon_old(i,iactive,j) .GT. min_stomate) THEN
	307	deepC_a(i,:,j)=deepC_a(i,:,j)*carbon(i,iactive,j)/carbon_old(i,iactive,j)
	308	ENDIF
	309	IF (carbon_old(i,islow,j) .GT. min_stomate) THEN
	310	deepC_s(i,:,j)=deepC_s(i,:,j)*carbon(i,islow,j)/carbon_old(i,islow,j)
	311	ENDIF
	312	IF (carbon_old(i,ipassive,j) .GT. min_stomate) THEN
	313	deepC_p(i,:,j)=deepC_p(i,:,j)*carbon(i,ipassive,j)/carbon_old(i,ipassive,j)
	314	ENDIF
	315	ENDIF
	316
[2655]	317	!biomass(i,j,:,:)=(biomass(i,j,:,:) * veget_max_old(i,j) + dilu_bio(i,:,:) * delta_veg(j)) / veget_max(i,j)
[2590]	318	TCarbon(i,j)=(TCarbon(i,j) * veget_max_old(i,j) + dilu_TCarbon(i) * delta_veg(j)) / veget_max(i,j)
[8]	319
[2590]	320	turnover_daily(i,j,:,:)=(turnover_daily(i,j,:,:)veget_max_old(i,j)+dilu_turnover_daily(i,:,:)delta_veg(j))/veget_max(i,j)
	321	bm_to_litter(i,j,:,:)=(bm_to_litter(i,j,:,:)veget_max_old(i,j)+dilu_bm_to_litter(i,:,:)delta_veg(j))/veget_max(i,j)
	322	co2flux_new(i,j)=(co2flux_old(i,j)veget_max_old(i,j)+dilu_co2flux_new(i)delta_veg(j))/veget_max(i,j)
	323	gpp_daily(i,j)=(gpp_daily(i,j)veget_max_old(i,j)+dilu_gpp_daily(i)delta_veg(j))/veget_max(i,j)
	324	resp_growth(i,j)=(resp_growth(i,j)veget_max_old(i,j)+dilu_resp_growth(i)delta_veg(j))/veget_max(i,j)
	325	resp_maint(i,j)=(resp_maint(i,j)veget_max_old(i,j)+dilu_resp_maint(i)delta_veg(j))/veget_max(i,j)
	326	resp_hetero(i,j)=(resp_hetero(i,j)veget_max_old(i,j)+dilu_resp_hetero(i)delta_veg(j))/veget_max(i,j)
	327	co2_to_bm(i,j)=(co2_to_bm(i,j)veget_max_old(i,j)+dilu_co2_to_bm(i)delta_veg(j))/veget_max(i,j)
	328	co2_fire(i,j)=(co2_fire(i,j)veget_max_old(i,j)+dilu_co2_fire(i)delta_veg(j))/veget_max(i,j)
[186]	329
	330	ENDIF
	331
[2615]	332	IF(veget_max(i,j).GT.min_stomate) THEN
	333
[3026]	334	! Correct biomass densities to conserve mass
	335	! since it's defined on veget_max
[2615]	336	biomass(i,j,:,:) = biomass(i,j,:,:) * veget_max_old(i,j) / veget_max(i,j)
	337
	338	ENDIF
	339
[947]	340	ENDDO ! loop over PFTs
[3026]	341	ENDDO ! loop over grid points
	342
	343	vartmp(:)=SUM(co2flux_new*veget_max,dim=2)
	344	CALL histwrite_p (hist_id_stomate, "tCO2FLUX", itime, vartmp, npts, hori_index)
	345	vartmp(:)=SUM(co2flux_old*veget_max_old,dim=2)
	346	CALL histwrite_p (hist_id_stomate, "tCO2FLUX_OLD", itime, vartmp, npts, hori_index)
	347	vartmp(:)=SUM(TCarbon*veget_max,dim=2)
	348	CALL histwrite_p (hist_id_stomate, "tCARBON", itime, vartmp, npts, hori_index)
[3057]	349	vartmp(:)=SUM(gpp_daily*veget_max,dim=2)
	350	CALL histwrite_p (hist_id_stomate, "tGPP", itime, vartmp, npts, hori_index)
	351	vartmp(:)=SUM(resp_growth*veget_max,dim=2)
	352	CALL histwrite_p (hist_id_stomate, "tRESP_GROWTH", itime, vartmp, npts, hori_index)
	353	vartmp(:)=SUM(resp_maint*veget_max,dim=2)
	354	CALL histwrite_p (hist_id_stomate, "tRESP_MAINT", itime, vartmp, npts, hori_index)
	355	vartmp(:)=SUM(resp_hetero*veget_max,dim=2)
	356	CALL histwrite_p (hist_id_stomate, "tRESP_HETERO", itime, vartmp, npts, hori_index)
	357	vartmp(:)=SUM(co2_to_bm*veget_max,dim=2)
	358	CALL histwrite_p (hist_id_stomate, "tCO2_TAKEN", itime, vartmp, npts, hori_index)
	359	vartmp(:)=SUM(co2_fire*veget_max,dim=2)
	360	CALL histwrite_p (hist_id_stomate, "tCO2_FIRE", itime, vartmp, npts, hori_index)
	361	vartmp(:)=SUM(SUM(biomass(:,:,:,icarbon),dim=3)*veget_max,dim=2)
	362	CALL histwrite_p (hist_id_stomate, "tBIOMASS", itime, vartmp, npts, hori_index)
	363	vartmp(:)=SUM(SUM(SUM(litter(:,:,:,:,icarbon),dim=4),dim=2)*veget_max,dim=2)
	364	CALL histwrite_p (hist_id_stomate, "tLITTER", itime, vartmp, npts, hori_index)
	365	vartmp(:)=SUM(SUM(fuel_1hr(:,:,:,icarbon),dim=3)*veget_max,dim=2)
	366	CALL histwrite_p (hist_id_stomate, "tFUEL1HR", itime, vartmp, npts, hori_index)
	367	vartmp(:)=SUM(SUM(fuel_10hr(:,:,:,icarbon),dim=3)*veget_max,dim=2)
	368	CALL histwrite_p (hist_id_stomate, "tFUEL10HR", itime, vartmp, npts, hori_index)
	369	vartmp(:)=SUM(SUM(fuel_100hr(:,:,:,icarbon),dim=3)*veget_max,dim=2)
	370	CALL histwrite_p (hist_id_stomate, "tFUEL100HR", itime, vartmp, npts, hori_index)
	371	vartmp(:)=SUM(SUM(fuel_1000hr(:,:,:,icarbon),dim=3)*veget_max,dim=2)
	372	CALL histwrite_p (hist_id_stomate, "tFUEL1000HR", itime, vartmp, npts, hori_index)
	373	vartmp(:)=SUM(SUM(carbon,dim=2)*veget_max,dim=2)
	374	CALL histwrite_p (hist_id_stomate, "tSOILC", itime, vartmp, npts, hori_index)
[3026]	375
[3057]	376	IF ( ok_pc ) THEN
	377	vartmp(:)=SUM(SUM(deepC_a,dim=2)*veget_max,dim=2)
	378	CALL histwrite_p (hist_id_stomate, "tDEEPCa", itime, vartmp, npts, hori_index)
	379	vartmp(:)=SUM(SUM(deepC_s,dim=2)*veget_max,dim=2)
	380	CALL histwrite_p (hist_id_stomate, "tDEEPCs", itime, vartmp, npts, hori_index)
	381	vartmp(:)=SUM(SUM(deepC_p,dim=2)*veget_max,dim=2)
	382	CALL histwrite_p (hist_id_stomate, "tDEEPCp", itime, vartmp, npts, hori_index)
	383	ENDIF
	384
[3026]	385	ENDIF
[8]	386
	387	END SUBROUTINE cover
	388
	389	END MODULE lpj_cover

Note: See TracBrowser for help on using the repository browser.

Download in other formats: