FHIST_BGC方案

1. 任务：

VPD 控制实验：
默认实验：不作任何修改，FHIST_BGC案例
控制实验：循环输入默认实验输出的1901-1910年的VPD。

输出变量：hist_fincl1="RAIN","QOVER","QH2OSFC","H2OSOI","GPP","NPP","NEP","TLAI","QBOT","RH","PBOT","TBOT","Rnet","FSH","SWdown","SWup","LWup","LWdown","QVEGE","QVEGT","QSOIL","FCEV","FCTR","FGEV","VPD"

月值输出：径流量、土壤湿度、gpp、npp、nep、lai、总辐射、净辐射、植物蒸腾、土壤蒸发、冠层蒸发、感热、降水
半小时输出：VPD

case: FHIST-BGC
1890-1900：预热10年，输出年值
1901-1910：历史输出前10年，月值和半小时
1911-2000：历史输出后90年，月值

控制实验基本过程：

1. 1890-1901年预热：D5_nocrop_spinup
2. 默认实验先跑10年，用于输出VPD强迫文件：D5_ncrop
3. 默认实验继续跑剩余的90年：D5
4. 控制实验跑100年：C5

变量名	long name	standard name	unit	note
降水	atmospheric rain	RAIN	mm/s
径流量		QOVER总径流， QH2OSFC地表径流	mm/s
总径流	surface runoff	QOVER	mm/s
地表径流	surface water runoff	QH2OSFC	mm/s
渗透	infiltration	QINFL	mm/s
截留	interception	QINTR	mm/s
穿透	throughfall	QDRIP	mm/s
土壤湿度	volumetric soil water (vegetated landunits only)	H2OSOI	mm3/mm3
总初级生产力	gross primary production	GPP	gC/m^2/s
净初级生产力	net primary production	NPP	gC/m^2/s
地上NPP	aboveground NPP	AGNPP	gC/m^2/s
地下NPP	belowground NPP	BGNPP	gC/m^2/s
自养呼吸	autotrophic respiration (MR + GR)	AR	gC/m^2/s
总异养呼吸	total heterotrophic respiration	HR	gC/m^2/s
生态系统总呼吸	total ecosystem respiration, autotrophic + heterotrophic	ER	gC/m^2/s
净生态系统生产力	net ecosystem production, excludes fire, landuse, and harvest flux, positive for sink	NEP	gC/m^2/s
净生态系统碳交换	net ecosystem exchange of carbon, includes fire, landuse, harvest, and hrv_xsmrpool flux, positive for source	NEE	gC/m^2/s
TLAI	total projected leaf area index	TLAI	none
2M比湿	2m specific humidity	Q2M	kg/kg
大气比湿	atmospheric specific humidity	QBOT	kg/kg	QBOT=Qair
大气比湿	atmospheric specific humidity	Qair	kg/kg
大气相对湿度	atmospheric relative humidity	RH	%
2M相对湿度	2m relative humidity	RH2M	%
大气压	atmospheric pressure	PBOT	Pa	PBOT=PSurf
表面气压	surface pressure	PSurf	Pa
大气CO2分压	atmospheric partial pressure of CO2	PCO2	Pa
2M大气温度	2m air temperature	TSA	K
植被温度	vegetation temperature	TV	K
地表温度	ground temperature	TG	K
大气空气潜在温度	atmospheric air potential temperature	THBOT	K
大气空气温度	atmospheric air temperature	TBOT	K
净辐射	net radiation	Rnet	W/m^2
大气长波辐射	atmospheric longwave radiation	FLDS	W/m^2	FLDS=LWdown
大气入射短波辐射	atmospheric incident solar radiation	FSDS	W/m^2	FSDS=SWdown
大气入射太阳辐射	atmospheric incident solar radiation	SWdown	W/m^2
向上短波辐射	upwelling shortwave radiation	SWup	W/m^2
大气长波辐射	atmospheric longwave radiation	LWdown	W/m^2
上行长波辐射	upwelling longwave radiation	LWup	W/m^2
吸收的太阳辐射	absorbed solar radiation	FSA	W/m^2
感热	sensible heat	FSH	W/m^2	FSH=Qh
感热	sensible heat	Qh	W/m^2
人为热通量	anthropogenic heat flux	Qanth	W/m^2
冠层以下的下行长波辐射	downward longwave radiation below the canopy	DLRAD	W/m^2
总潜热通量至大气	total latent heat flux [+ to atm]	EFLX_LH_TOT	W/m^2	NO
冠层蒸发	canopy evaporation	QVEGE	mm/s
冠层蒸发	canopy evaporation	FCEV	W/m^2
冠层蒸腾	canopy transpiration	QVEGT	mm/s
冠层蒸腾	canopy transpiration	FCTR	W/m^2
土壤蒸发	Ground evaporation (soil/snow evaporation + soil/snow sublimation - dew)	QSOIL	mm/s
土壤蒸发	ground evaporation	FGEV	W/m^2
VPD			hPa
阴PAR	shade PAR	PAR_shade	umol/m2/s
阳PAR	sunlit PAR	PAR_sun	umol/m2/s
阴叶光合	shaded leaf photosynthesis	PSNSHA	umolCO2/m^2/s
阳叶光合	sunlit leaf photosynthesis	PSNSUN	umolCO2/m^2/s
光合	photosynthesis	FPSN	umol/m2s
阴叶气孔阻力	shaded leaf stomatal resistance	RSSHA	s/m
阳叶气孔阻力	sunlit leaf stomatal resistance	RSSUN	s/m

2. 默认实验

2.1 预热11年：1891-1901 D5_nocrop_spinup

案例：预热
1. 关掉作物模块
2. 默认案例，输出VPD

%% PCT_NAT_PFT，最大值设为100%，其它为0
%% PCT_CROP设为0
%% Default_1D_none    flanduse_timeseries

%% PCT_NAT_PFT，最大值设为100%，其它为0
%% PCT_NATVEG = PCT_CROP + PCT_NATVEG
%% PCT_CROP设为0
%% Default_1D_none    surfdata

新建案例
# FHIST_BGC实验，1度分辨率
./create_newcase --compset FHIST_BGC --res f09_f09_mg17 --case ~/CESM2/scratch/D5_nocrop_spinup
调用节点
./xmlchange NTASKS=-60
./xmlchange NTASKS_ESP=60
./xmlquery NTASKS,NTHRDS,ROOTPE

输出
./case.setup
# 年值输出
vim user_nl_clm
hist_empty_htapes = .true
hist_avgflag_pertape = 'A'
hist_nhtfrq = -8760
hist_mfilt = 1
hist_fincl1="RAIN","QOVER","QH2OSFC","H2OSOI","GPP","NPP","NEP","TLAI","QBOT","RH","PBOT","TBOT","Rnet","FSH","SWdown","SWup","LWup","LWdown","QVEGE","QVEGT","QSOIL","FCEV","FCTR","FGEV","VPD"
fsurdat = '/BIGDATA1/sysu_syli_3/CESM2/inputdata/lnd/clm2/surfdata_map/release-clm5.0.18/surfdata_0.9x1.25_hist_16pfts_Irrig_CMIP6_simyr1850_c190214_Default_1D_none.nc'
flanduse_timeseries = '/BIGDATA1/sysu_syli_3/CESM2/inputdata/lnd/clm2/surfdata_map/landuse.timeseries_0.9x1.25_hist_16pfts_Irrig_CMIP6_simyr1850-2015_c170824_Default_1D_none.nc'

关闭大气模块的默认输出
vim user_nl_cam
empty_htapes = .true. 
fincl1=' '
fincl2=' '
fincl3=' '
fincl4=' '
fincl5=' '
fincl6=' '

% 关闭作物模块
./xmlquery CLM_BLDNML_OPTS
./xmlchange CLM_BLDNML_OPTS="-bgc bgc"


% 运行方式，起始时段，运行时长设置
./xmlchange RUN_TYPE=startup,RUN_STARTDATE=1891-01-01,GET_REFCASE=FALSE
./xmlchange DOUT_S=TRUE
./xmlchange JOB_WALLCLOCK_TIME=24:00:00
./xmlchange STOP_N=2,STOP_OPTION=nyears
./xmlchange RESUBMIT=4         % 推荐一次运行任务时段减小，多提交几次任务

./xmlquery STOP_OPTION,STOP_N,RUN_TYPE,CONTINUE_RUN,RUN_REFCASE,RUN_REFDATE,RESUBMIT,RUN_STARTDATE,RUN_REFDIR,DOUT_S,DOUT_S_ROOT,JOB_WALLCLOCK_TIME,GET_REFCASE


复制修改后的源码文件放在case目录下的SourceMods/src.clm/目录下就可以了

./case.build --skip-provenance-check &> build.log
./case.submit

判断预热是否达到平衡态？
下图分别以净辐射Rnet的时间序列和1899年GPP年均值的空间分布说明，预热已经达到平衡状态，变量的空间分布正常。

2.2 默认实验生成VPD forcing文件：1901-1910 D5_nocrop

新建一个FHIST_BGC_branch案例，以branch run的方式运行，接着预热案例，从1901-01-01开始。

例：默认实验生成VPD forcing文件：1901-1910 
1. 关掉作物模块
2. 默认案例，输出VPD

%% PCT_NAT_PFT，最大值设为100%，其它为0
%% PCT_CROP设为0
%% Default_1D_none    flanduse_timeseries

%% PCT_NAT_PFT，最大值设为100%，其它为0
%% PCT_NATVEG = PCT_CROP + PCT_NATVEG
%% PCT_CROP设为0
%% Default_1D_none    surfdata


./create_newcase --compset FHIST_BGC --res f09_f09_mg17 --case ~/CESM2/scratch/D5_nocrop

./xmlchange NTASKS=-60
./xmlchange NTASKS_ESP=60
./xmlquery NTASKS,NTHRDS,ROOTPE

% 输出
./case.setup
vim user_nl_clm
hist_empty_htapes = .true.
hist_avgflag_pertape = 'A','A'
hist_nhtfrq = 0,1
hist_mfilt = 12,17520
hist_fincl1="RAIN","QOVER","QH2OSFC","H2OSOI","GPP","NPP","NEP","TLAI","QBOT","RH","PBOT","TBOT","Rnet","FSH","SWdown","SWup","LWup","LWdown","QVEGE","QVEGT","QSOIL","FCEV","FCTR","FGEV","VPD"
hist_fincl2='VPD'
fsurdat = '/BIGDATA1/sysu_syli_3/CESM2/inputdata/lnd/clm2/surfdata_map/release-clm5.0.18/surfdata_0.9x1.25_hist_16pfts_Irrig_CMIP6_simyr1850_c190214_Default_1D_none.nc'
flanduse_timeseries = '/BIGDATA1/sysu_syli_3/CESM2/inputdata/lnd/clm2/surfdata_map/landuse.timeseries_0.9x1.25_hist_16pfts_Irrig_CMIP6_simyr1850-2015_c170824_Default_1D_none.nc'

# 关闭大气模块的默认输出
vim user_nl_cam
empty_htapes = .true. 
fincl1=' '
fincl2=' '
fincl3=' '
fincl4=' '
fincl5=' '
fincl6=' '

% 关闭作物
./xmlquer