Dear all:
The AD SPINUP model was used to preheat clm5.0 for 100 years using self-produced atmospheric forcing data. However, after checking the history file, I find that the values of NNP and NEP becomes 0 after 10 years. Is this normal? Here is my code for running clm. Is there any error? ZIP file is my NEE and NPP of history file.
I have one more question. I will simulate the carbon dynamics from 2015 to 2100. Can I set RUN_TYPE=branch
RUN_REFCASE, RUN_REFDATE and other related information to run?
#! /bin/bash
cd /home/user/clm5.0.35/clm5new/cime/scripts
./create_newcase --case ~/cases/preheat/ACCESS.board1.preheat1206 --res CLM_USRDAT --compset SSP585_DATM%1PT_CLM50%BGC_SICE_SOCN_SROF_SGLC_SWAV --run-unsupported --compiler intel --mach myintel
cd /home/user/cases/preheat/ACCESS.board1.preheat1206
./xmlchange CLM_ACCELERATED_SPINUP="on"
./xmlchange CLM_USRDAT_NAME='tibet'
./xmlchange ATM_DOMAIN_PATH='/home/user/tibet/myout_board1/'
./xmlchange LND_DOMAIN_PATH='/home/user/tibet/myout_board1/'
./xmlchange ATM_DOMAIN_FILE='domain.lnd.fv0.47x0.63_gx1v7.180521_myout_board1_board1.nc'
./xmlchange LND_DOMAIN_FILE='domain.lnd.fv0.47x0.63_gx1v7.180521_myout_board1_board1.nc'
./xmlchange NTASKS=1
./xmlchange STOP_N=1
./xmlchange STOP_OPTION=nyears
./xmlchange RESUBMIT=99
./xmlchange DATM_CLMNCEP_YR_START=2015
./xmlchange DATM_CLMNCEP_YR_END=2015
./xmlchange DATM_CLMNCEP_YR_ALIGN=2015
./xmlchange RUN_REFDATE=0001-01-01
./xmlchange RUN_STARTDATE=0001-01-01
./case.setup
echo "fsurdat = '/home/user/tibet/myout_board1/surfdata_0.47x0.63_78pfts_CMIP6_simyr2000_c180508_myout_board1_board1.nc'
flanduse_timeseries = '/home/user/tibet/myin/landuse.timeseries_0.9x1.25_SSP5-8.5_78pfts_CMIP6_simyr1850-2100_c181209.nc'
finidat = '/home/user/cesm/inputdata/lnd/clm2/initdata_map/clmi.I2000Clm50BgcCrop.2011-01-01.1.9x2.5_gx1v7_gl4_simyr2000_c190312.nc'
! hist_empty_htapes = .true.
hist_fincl1='GPP:A','NPP:A','NEP:A','NEE:A','NBP:A','TOTECOSYSC:A','TOTSOMC:A','TOTVEGC:A','TLAI:A','TWS:A','H2OSNO:A'
hist_fincl2='GPP:X','NPP:X','NEP:X','NEE:X','NBP:X','TOTECOSYSC:X','TOTSOMC:X','TOTVEGC:X','TLAI:X','TWS:X','H2OSNO:X'
hist_nhtfrq= 0,0,0
hist_mfilt = 3600,3600,3600
use_init_interp = .true.
do_harvest = .false.
do_transient_crops = .false.
do_transient_pfts = .false." >> ./user_nl_clm
echo "<?xml version='1.0'?>
<file id='stream' version='1.0'>
<dataSource>
GENERIC
</dataSource>
<domainInfo>
<variableNames>
time time
xc lon
yc lat
area area
mask mask
</variableNames>
<filePath>
/home/user/tibet/myout_board1/
</filePath>
<fileNames>
domain.lnd.fv0.47x0.63_gx1v7.180521_myout_board1_board1.nc
</fileNames>
</domainInfo>
<fieldInfo>
<variableNames>
ZBOT z
TBOT tbot
QBOT shum
WIND wind
PRECTmms precn
FSDS swdn
PSRF pbot
FLDS lwdn
</variableNames>
<filePath>
/home/user/tibet/15pt_forcing
</filePath>
<fileNames>
ACC.a86.board1.nc
</fileNames>
<offset>
0
</offset>
</fieldInfo>
</file>" >> ./user_datm.streams.txt.CLM1PT.CLM_USRDAT
echo "<?xml version='1.0'?>
<file id='stream' version='1.0'>
<dataSource>
GENERIC
</dataSource>
<domainInfo>
<variableNames>
time time
lon lon
lat lat
area area
mask mask
</variableNames>
<filePath>
/home/user/cesm/inputdata/atm/cam/chem/trop_mozart_aero/aero
</filePath>
<fileNames>
aerodep_clm_SSP585_b.e21.BSSP585cmip6.f09_g17.CMIP6-SSP5-8.5.001_2014-2101_monthly_0.9x1.25_c190419.nc
</fileNames>
</domainInfo>
<fieldInfo>
<variableNames>
BCDEPWET bcphiwet
BCPHODRY bcphodry
BCPHIDRY bcphidry
OCDEPWET ocphiwet
OCPHIDRY ocphidry
OCPHODRY ocphodry
DSTX01WD dstwet1
DSTX01DD dstdry1
DSTX02WD dstwet2
DSTX02DD dstdry2
DSTX03WD dstwet3
DSTX03DD dstdry3
DSTX04WD dstwet4
DSTX04DD dstdry4
</variableNames>
<filePath>
/home/user/cesm/inputdata/atm/cam/chem/trop_mozart_aero/aero
</filePath>
<fileNames>
aerodep_clm_SSP585_b.e21.BSSP585cmip6.f09_g17.CMIP6-SSP5-8.5.001_2014-2101_monthly_0.9x1.25_c190419.nc
</fileNames>
<offset>
0
</offset>
</fieldInfo>
</file>" >> ./user_datm.streams.txt.presaero.SSP5-8.5
./case.build
./case.submit
The AD SPINUP model was used to preheat clm5.0 for 100 years using self-produced atmospheric forcing data. However, after checking the history file, I find that the values of NNP and NEP becomes 0 after 10 years. Is this normal? Here is my code for running clm. Is there any error? ZIP file is my NEE and NPP of history file.
I have one more question. I will simulate the carbon dynamics from 2015 to 2100. Can I set RUN_TYPE=branch
RUN_REFCASE, RUN_REFDATE and other related information to run?
#! /bin/bash
cd /home/user/clm5.0.35/clm5new/cime/scripts
./create_newcase --case ~/cases/preheat/ACCESS.board1.preheat1206 --res CLM_USRDAT --compset SSP585_DATM%1PT_CLM50%BGC_SICE_SOCN_SROF_SGLC_SWAV --run-unsupported --compiler intel --mach myintel
cd /home/user/cases/preheat/ACCESS.board1.preheat1206
./xmlchange CLM_ACCELERATED_SPINUP="on"
./xmlchange CLM_USRDAT_NAME='tibet'
./xmlchange ATM_DOMAIN_PATH='/home/user/tibet/myout_board1/'
./xmlchange LND_DOMAIN_PATH='/home/user/tibet/myout_board1/'
./xmlchange ATM_DOMAIN_FILE='domain.lnd.fv0.47x0.63_gx1v7.180521_myout_board1_board1.nc'
./xmlchange LND_DOMAIN_FILE='domain.lnd.fv0.47x0.63_gx1v7.180521_myout_board1_board1.nc'
./xmlchange NTASKS=1
./xmlchange STOP_N=1
./xmlchange STOP_OPTION=nyears
./xmlchange RESUBMIT=99
./xmlchange DATM_CLMNCEP_YR_START=2015
./xmlchange DATM_CLMNCEP_YR_END=2015
./xmlchange DATM_CLMNCEP_YR_ALIGN=2015
./xmlchange RUN_REFDATE=0001-01-01
./xmlchange RUN_STARTDATE=0001-01-01
./case.setup
echo "fsurdat = '/home/user/tibet/myout_board1/surfdata_0.47x0.63_78pfts_CMIP6_simyr2000_c180508_myout_board1_board1.nc'
flanduse_timeseries = '/home/user/tibet/myin/landuse.timeseries_0.9x1.25_SSP5-8.5_78pfts_CMIP6_simyr1850-2100_c181209.nc'
finidat = '/home/user/cesm/inputdata/lnd/clm2/initdata_map/clmi.I2000Clm50BgcCrop.2011-01-01.1.9x2.5_gx1v7_gl4_simyr2000_c190312.nc'
! hist_empty_htapes = .true.
hist_fincl1='GPP:A','NPP:A','NEP:A','NEE:A','NBP:A','TOTECOSYSC:A','TOTSOMC:A','TOTVEGC:A','TLAI:A','TWS:A','H2OSNO:A'
hist_fincl2='GPP:X','NPP:X','NEP:X','NEE:X','NBP:X','TOTECOSYSC:X','TOTSOMC:X','TOTVEGC:X','TLAI:X','TWS:X','H2OSNO:X'
hist_nhtfrq= 0,0,0
hist_mfilt = 3600,3600,3600
use_init_interp = .true.
do_harvest = .false.
do_transient_crops = .false.
do_transient_pfts = .false." >> ./user_nl_clm
echo "<?xml version='1.0'?>
<file id='stream' version='1.0'>
<dataSource>
GENERIC
</dataSource>
<domainInfo>
<variableNames>
time time
xc lon
yc lat
area area
mask mask
</variableNames>
<filePath>
/home/user/tibet/myout_board1/
</filePath>
<fileNames>
domain.lnd.fv0.47x0.63_gx1v7.180521_myout_board1_board1.nc
</fileNames>
</domainInfo>
<fieldInfo>
<variableNames>
ZBOT z
TBOT tbot
QBOT shum
WIND wind
PRECTmms precn
FSDS swdn
PSRF pbot
FLDS lwdn
</variableNames>
<filePath>
/home/user/tibet/15pt_forcing
</filePath>
<fileNames>
ACC.a86.board1.nc
</fileNames>
<offset>
0
</offset>
</fieldInfo>
</file>" >> ./user_datm.streams.txt.CLM1PT.CLM_USRDAT
echo "<?xml version='1.0'?>
<file id='stream' version='1.0'>
<dataSource>
GENERIC
</dataSource>
<domainInfo>
<variableNames>
time time
lon lon
lat lat
area area
mask mask
</variableNames>
<filePath>
/home/user/cesm/inputdata/atm/cam/chem/trop_mozart_aero/aero
</filePath>
<fileNames>
aerodep_clm_SSP585_b.e21.BSSP585cmip6.f09_g17.CMIP6-SSP5-8.5.001_2014-2101_monthly_0.9x1.25_c190419.nc
</fileNames>
</domainInfo>
<fieldInfo>
<variableNames>
BCDEPWET bcphiwet
BCPHODRY bcphodry
BCPHIDRY bcphidry
OCDEPWET ocphiwet
OCPHIDRY ocphidry
OCPHODRY ocphodry
DSTX01WD dstwet1
DSTX01DD dstdry1
DSTX02WD dstwet2
DSTX02DD dstdry2
DSTX03WD dstwet3
DSTX03DD dstdry3
DSTX04WD dstwet4
DSTX04DD dstdry4
</variableNames>
<filePath>
/home/user/cesm/inputdata/atm/cam/chem/trop_mozart_aero/aero
</filePath>
<fileNames>
aerodep_clm_SSP585_b.e21.BSSP585cmip6.f09_g17.CMIP6-SSP5-8.5.001_2014-2101_monthly_0.9x1.25_c190419.nc
</fileNames>
<offset>
0
</offset>
</fieldInfo>
</file>" >> ./user_datm.streams.txt.presaero.SSP5-8.5
./case.build
./case.submit