cesm单点模拟

开心，进入新part

步骤

#— 这个前提是已经建立好地表数据和大气强迫数据

#— 环境设置

export GRIDNAME=1x1_xx #这个是自己命名的
export CDATE='21xxxx' #按照自己的数据的情况改
export GRIDFILE=.../cesm_2.1.3/components/clm/tools/mkmapgrids/SCRIPgrid_1x1_xx_nomask_c210527.nc 
export MAPFILE=.../cesm_2.1.3/components/clm/tools/mkmapdata/map_1x1_mq_noocean_to_1x1_mq_nomask_aave_da_210527.nc 
export OCNDOM=domain.ocn_noocean.nc #这个可以改成自己的domain文件的名称
export ATMDOM=domain.lnd.${GRIDNAME}_noocean.nc

#— 创建案例

cd .../cesm_2.1.3/cime/scripts
./create_newcase -case casename -res CLM_USRDAT -compset I1PtClm50SpGs -mach bscca3 --run-unsupported

#— 更改相应设置

./xmlchange ATM_DOMAIN_PATH=.../cesm_2.1.3/cime/tools/mapping/gen_domain_files,LND_DOMAIN_PATH=.../cesm_2.1.3/cime/tools/mapping/gen_domain_files
./xmlchange LND_DOMAIN_PATH=.../cesm_2.1.3/cime/tools/mapping/gen_domain_files,LND_DOMAIN_PATH=.../cesm_2.1.3/cime/tools/mapping/gen_domain_files
./xmlchange ATM_DOMAIN_FILE=$ATMDOM
./xmlchange LND_DOMAIN_FILE=$ATMDOM
./xmlchange CLM_USRDAT_NAME=$GRIDNAME
./xmlchange STOP_OPTION=ndays
./xmlchange STOP_N=31
#./xmlchange REST_OPTION=$REST_OPTION
#./xmlchange REST_N=$REST_N
./xmlchange RUN_STARTDATE=2017-06-01 #这里因为强迫场数据的连续性
#./xmlchange DATM_MODE=CLM1PT
./case.setup

#— 设置user_nl_clm

echo "fsurdat='.../cesm_2.1.3/inputdata/lnd/clm2/surfdata_map/surfdata_1x1_xxx_hist_16pfts_Irrig_CMIP6_simyr2000.nc'" >> user_nl_clm
echo "finidat='.../cesm_2.1.3/inputdata/lnd/clm2/initdata_map/clmi.I2000Clm50BgcCrop.2011-01-01.1.9x2.5_gx1v7_gl4_simyr2000_c190312.nc'" >> user_nl_clm
echo "use_init_interp = .true." >> user_nl_clm
echo "hist_empty_htapes = .true." >> user_nl_clm
echo "hist_fincl1=" " >> user_nl_clm #设置需要变量
echo "hist_nhtfrq = -24" >> user_nl_clm
echo "hist_mfilt = 30" >> user_nl_clm

#— 设置user_nl_datm

echo "streams = "datm.streams.txt.CLM1PT.CLM_USRDAT 2017 2017 2017","datm.streams.txt.presaero.clim_2000 2010 2010 2010"" >> user_nl_datm

#— 修改流文件

cd CaseDocs/
cp datm.streams.txt.CLM1PT.CLM_USRDAT \/public1/home/sc91286/cesm_2.1.3/cime/scripts/mqtext/user_datm.streams.txt.CLM1PT.CLM_USRDAT
cd ..
vim user_datm.streams.txt.CLM1PT.CLM_USRDAT #这个自己打开改一下路径和大气强迫资料的名称

#— build
./case.build
./case.submit

存在问题

！cp的时候出现 cannot create regular file
是路径的问题。在“/”不能加空格
另外，这一步要删去cdate

! domainfile 没有找到
命名错误了，改一下名称

！流文件不知道在哪里
在Case/CaseDocs/

! 长时间计算，没有文件输出（这里注意了，不然就是金钱的大量损失，对不起我老板！）
这个主要的问题居然是再user_nl_clm中没有设置hist_fincl1就算是你都要，也要设为空值！！！金钱警告！
变量在这里哦：clm的main文件夹下有个histFldsMod的脚本（大佬的）

！大气强迫中太阳辐射不能小于等于0

！显示错误：LVD not found, all data is before yearfirst
这个是user_nl_datm 中流文件时间设置不对

！显示，大气返回的xxx数据有问题，算不下去的
拜托检查一下自己制作的大气强迫数据是不是有问题
检查单位等
不然就会像我一样卡n天！

！大气强迫数据模拟出来的时间

峰值提前，猜测是由于大气强迫资料的时间设置有问题

感谢老板、万能师姐、大腿、大神等！！
有人罩是多么快乐的一件事情啊
还在坑底，希望有爬出来的一天！！