参数化方案对耦合模式影响试验

今天建立几个实验

HIST_model_test_lthf_tbf

HIST_model_test_lthf_phiq

HIST_model_test_sshf_tbf

HIST_model_test_pblh_zli

HIST_model_test_nocouple

十天过去了，这次运行的速度显然慢了好多，还没运行到2000年，没办法只能就地取材了。

我准备分析1992年到1996年的数据

将之前的文件夹复制一份

cp -r 10yearmean 5yearmean

然后对CTRL模式输出资料进行平均，1992-1996年平均

ls B20TR_C35_y100.iap.h0.199[2-6]*|xargs -I{} cdo yearmean {} CASESM.atm.1992-1996ym

对EXP1模式输出资料进行平均，1992-1996年平均

ls B20TR_C35_y100.iap.h0.199[2-6]*|xargs -I{} cdo yearmean {} CASESM.atm.1992-1996ym

再准备观测资料

cd /sharedata/dataset/historical/SD013-NCEP_GODAScp dbss_obml.199[2-6]* /data/chengxl/obs_duibiyanjiu/10yearmean/

obml GODAS

salt GODAS

skint 皮层温度 ERA-in

slp NCEP reanalysis

sp ERA——interim

t2m ERA——interim

u10m ERA——interim

ucur NCEP_GODAS

vcur NCEP_GODAS

v10mERA——interim

uwnd NCEP reanalysis

vwnd NCEP reanalysis

将这些资料都复制进来，然后进行多年平均

ls dbss_obml.199[2-6].*|xargs -I{} cdo -b f32 yearmean {} obml.1992-1996ym.nc
ls salt.199[2-6]*|xargs -I{} cdo -b f32  yearmean  {} salt.1992-1996ym.nc 
ls skint_199[2-6]*|xargs -I{} cdo -b f32  yearmean  {} skint.1992-1996ym.nc
ls slp.199[2-6]*|xargs -I{} cdo -b f32  yearmean  {} slp.1992-1996ym.nc
ls sst_199[2-6]*|xargs -I{} cdo -b f32  yearmean  {} sst.1992-1996ym.nc 
ls t2m_199[2-6]*|xargs -I{} cdo -b f32  yearmean  {} t2m.1992-1996ym.nc
ls u10m_199[2-6]*|xargs -I{} cdo -b f32  yearmean  {} u10m.1992-1996ym.nc 
ls ucur.199[2-6]*|xargs -I{} cdo -b f32  yearmean  {} ucur.1992-1996ym.nc 
ls vcur.199[2-6]*|xargs -I{} cdo -b f32  yearmean  {} vcur.1992-1996ym.nc 
ls uwnd.*199[2-6]*|xargs -I{} cdo -b f32 yearmean {} uwnd.1992-1996ym.nc 
ls vwnd.*199[2-6]*|xargs -I{} cdo -b f32 yearmean {} vwnd.1992-1996ym.nc 
ls v10m_*199[2-6]*|xargs -I{} cdo -b f32 yearmean {} v10.1992-1996ym.nc

准备好了资料可以对脚本进行处理了

(base) [chengxl@login02 5yearmean]$ ls
autorun.sh                      plot_2000_casesm_psl_atm_7.ncl   plot_2000_casesm_u10_atm.ncl
panel_two_sets.ncl              plot_2000_casesm_psl_atm.ncl     plot_2000_casesm_uv_atm.ncl
plot_2000_casesm_lthf_atm.ncl   plot_2000_casesm_psl_uv_atm.ncl
plot_2000_casesm_psl_atm_1.ncl  plot_2000_casesm_ts_atm.ncl

先查看需要批量修改哪几个路径

首先将2000-2009转为1992-1996

sed -i "s/2000-2009/1992-1996/g" *.ncl

sed -i "s/HIST_model_test_finial/HIST_model_test_nocouple/g" *.

./autorun.sh

sed -i "s/HIST_model_test_nocouple/HIST_model_test_lthf_tbf/g" *.ncl

./autorun.sh

sed -i "s/HIST_model_test_nocouple/HIST_model_test_lthf_phiq/g" *.ncl

./autorun.sh

cd pblh_zil/

sed -i "s/HIST_model_test_nocouple/HIST_model_test_pblh_zli/g" *.ncl

./autorun.sh

大气变量已经输出好了，现在对海洋变量进行输出

首先对模式输出CTRL和EXP1 计算年平均，观测资料已经处理好了

分别对四个文件夹文件使用命令

ls MMEAN199[2-6]* |xargs -I{} cdo -b f32 yearmean {} MMEAN1992-1996ym.nc

然后新建一个文件夹将之前写的画图脚本批量修改2000-2009变为1992-1996

sed -i "s/2000-2009/1992-1996/g" *.ncl1127 sed -i "s/2000_/1992-1996_/g" *.ncl1128 sed -i "s/HIST_model_test_finial/HIST_model_test_nocouple/g" *.ncl1129 rm *.png1130 ./autorun.sh1131 vim plot_2000_casesm_sst.ncl1132 ./autorun.sh1136 cp -r nocouple/ lthf_tbf1137 cd lthf_tbf/1139 sed -i "s//HIST_model_test_nocouple/HIST_model_test_lthf_tbf/g" *.ncl1140 sed -i "s/HIST_model_test_nocouple/HIST_model_test_lthf_tbf/g" *.ncl1141 ./autorun.sh1142 cd ..1144 cp -r lthf_tbf/ lthf_phiq1145 cd lthf_phiq/1147 rm *.png1148 history1150 sed -i "s/HIST_model_test_lthf_tbf/HIST_model_test_lthf_phiq/g" *.ncl1151 ./autorun.sh1152 cd ..1154 cp -r lthf_tbf pblh_zil1155 cd pblh_zil/1157 rm *.png1158 sed -i "s/HIST_model_test_lthf_tbf/HIST_model_test_pblh_zli/g" *.ncl1159 ./autorun.sh

然后重复替换变量和画图