GRAPES MESO模式学习笔记04——脚本运行之资料同化
对资料同化我一无所知，只能修改已有的脚本运行。
摘自GRPAES手册的GRAPES VAR系统流程图
分外循环和内循环两个部分。外循环在高分辨率下计算背景场和观测之间的偏差，即观测增量。内循环在低分辨率下进行极小化迭代。
外循环
对应可执行程序4dvar_h.exe
输入
背景场资料：标准初始化得到的背景场信息（grapesinput）
观测资料：观测资料预处理得到的观测资料信息（GTS，常规观测资料目录）
背景误差估计：rmsdata
还需要一系列静态文件资料。包括TBL文件等。
参数
namelist_h.input和namelist.4dvar
输出
观测增量 (d_***YYYYMMDDHH00_nod)
初猜场信息（xb_h.input,xb_l.input）
内循环
对应可执行程序4dvar.exe
输入
观测增量 (d_***YYYYMMDDHH00_nod)
初猜场信息（xb_h.input,xb_l.input）
参数
namelist.input和namelist.4dvar
输出
xb_h.grd：GrADS 格式（顺序无格式）的外循环高分辨率的背景场
dxa_h.grd：GrADS 格式的外循环高分辨率的分析增量,dxa_h.ctl 可画图查看
xa_h.grd：GrADS 格式的外循环高分辨率的分析场
xb_l.grd：GrADS 格式的内循环低分辨率的背景场
dxa_l.grd：GrADS 格式的内循环低分辨率的分析增量，dxa_l.ctl 可画图查看
xa_l.grd： GrADS 格式的内循环低分辨率的分析场
grapesinput_h-4dv-YYYYMMDDHH：与外循环高分辨率背景场 grapesinput_h 数据结构一致的分析场
遇到的问题
运行4dvar.exe遇到同一个错误卡了两天，错误信息：
源代码
[code]
!————————————————————————-
! [5.0] Output for test
!————————————————————————-
IF(on_monitor() ) then
OPEN(9,file=’../../data/Eigen’//FILOUT(ivar)//’.dat’)
write(9,’(2x,a)’) ‘i,rmspsi(i), rmschi, rmspiu,rmsrh, rmsThu,rmsq,rmspu’
write(9,’(1x,i3,2E21.12, 2E21.12, 4E21.12)’) (i,rmspsi(i), rmschi(i), &
rmspiu(i),rmsrh(i), rmsThu(i),rmsq(i),rmspu(i),rmspu(i),i=1,N)
WRITE(9,’(A,A)’) ‘symetric matrix:’,FILOUT(ivar)
DO i=kms, kme
WRITE(9,30) i, (rr(i,j),j=1,31 ) !kms,kme)
END DO
WRITE(9,’(A,A)’) ‘the eigenvalues of :’,FILOUT(ivar)
DO i=1,N
WRITE(9,30) i, EVAL(i)
END DO
WRITE(9,’(A,A)’) ‘the eigenvectors of :’,FILOUT(ivar)
DO k=1,N
WRITE(9,30) k, (evec(k,m),m=1,N)
END DO
CLOSE(9)
30 FORMAT(i3,31E15.6)
ENDIF
[/code]
就是file=’../../data/Eigen’//FILOUT(ivar)//’.dat’这句害得我三天都没有结果。将路径写死在文件中，为啥我前两天没有发现这个问题呢？
修改为当前目录，不使用上级或上上级目录。可见，模式运行需要的文件远比手册中写的要多，好多静态文件隐藏在运行目录中。
m4dvar
这步骤更新边界条件，我不知道应该叫什么名字，用字母代替。对应可执行程序m4dvar.exe
输入
grapesinput-4dv-${begintime} grapesinput${begintime}
grapesbdy grapesbdy_old
输出
grapesbdy_new
grapesinput${begintime}
脚本
生成namelist文件
[shell collapse=”true”]
#!/bin/ksh
# generate namelists for 4dvar
# $14 => ${14} parameter positon in function
##################
# head
##################
set -x
set -e
set -u
date
###############
# directory
###############
BASE_DIR=./nwp/system/grapes_meso/my_grapes_meso
BASE_SOURCE_DIR=./nwp/system/grapes_meso/GRAPES_MESO3.3.2.4
RUN_DIR=${BASE_DIR}/run/an
BCKG_DIR=${BASE_DIR}/run/bckg
OBS_DIR=${BASE_DIR}/run/obs
GEODATA_DIR=/cma/g2/COMMDATA/static/rfs/geog/v3
CONDATA_DIR=${BASE_DIR}/condat
AN_CONDAT_DIR=${CONDATA_DIR}/an
NEWDATE_BIN=${BASE_SOURCE_DIR}/fcst/grapes_model/run/newdate
############
# parameter
############
YYYYMMDDHH=2013041012
mfcast_len=24
do_static=.true.
do_3dv=.true.
output_inteval=3 #model output inteval in hour
model_dt=90 #model time step in second
make_surface_t=.false.
while getopts “t:b:r:g:c:” arg
do
case $arg in
t)
YYYYMMDDHH=$OPTARG
;;
b)
BCKG_DIR=$OPTARG
;;
r)
RUN_DIR=$OPTARG
;;
esac
done
begintime=$YYYYMMDDHH
endtime=$( ${NEWDATE_BIN} $begintime +$mfcast_len)
step_max=`expr $mfcast_len * 3600 / $model_dt`
step_output=`expr $output_inteval * 3600 / $model_dt`
YYYY1=`echo $begintime | cut -c1-4`
MM1=`echo $begintime | cut -c5-6`
DD1=`echo $begintime | cut -c7-8`
HH1=`echo $begintime | cut -c9-10`
YYYY2=`echo $endtime | cut -c1-4`
MM2=`echo $endtime | cut -c5-6`
DD2=`echo $endtime | cut -c7-8`
HH2=`echo $endtime | cut -c9-10`
########################
# enter work dir
########################
test -d ${RUN_DIR} || mkdir -p ${RUN_DIR}
cd ${RUN_DIR}
test -d input || mkdir -p input
test -d output || mkdir -p output
test -d output/check || mkdir -p output/check
test -d grads || mkdir grads
#ln -sf ${GEODATA_DIR} geog_data
#######################
# prepare data
#######################
#################
# namelists
#################
# namelist_h.input
generate_namelist_h_input(){
set -x
test -f namelist_h.input || rm -f namelist_h.input
cat > namelist_h.input « END_OF_DATA &namelist_01 s_we = 1, e_we = 502, s_sn = 1, e_sn = 330, s_vert = 1, e_vert = 31, spec_bdy_width = 5, dt = $1, time_step_max = $2, time_step_count_output = $3, step_phytend_output = 0, maxit = 2, interp = 0, irec = 0, nh = 1, writetofile = .false., xs_we = 70., ys_sn = 15., cen_lat = 40., xd = 0.15, yd = 0.15, ztop = 35000.0, epson_depart = 1.0, i_parent_start = 1, j_parent_start = 1, shw = 1, parent_grid_ratio = 1 / &namelist_02 dyn_opt = 4, diff_opt = 1, km_opt = 1, damp_opt = 0, spec_zone = 1, relax_zone = 4, isfflx = 1, ifsnow = 0, icloud = 1, num_soil_layers = 4, julyr = 0, julday = 1, gmt = 0, radt = 60., bldt = 0, cudt = 5, time_step_begin_restart = 0 / &namelist_03 mp_lc_logical = .true., ! .true. for microphysics, .false. for largescale_rain mp_physics = 13, lc_physics = 1, ! now only 1 can selected(031212) ra_lw_physics = 1, ra_sw_physics = 1, bl_sfclay_physics = 1, bl_surface_physics = 1, bl_pbl_physics = 1, cu_physics = 2 / &namelist_05 init_date = ${13}, start_year = $4, start_month = $5, start_day = $6, start_hour = $7, start_minute = 00, start_second = 00, end_year = $8, end_month = $9, end_day = ${10}, end_hour = ${11}, end_minute = 00, end_second = 00, interval_seconds = 10800, real_data_init_type = 2, bdyfrq = 10800, nested = .false., specified = .false. / &namelist_bckg nz_bckg = 17, xs_we_bckg = 0., ys_sn_bckg = 0., xe_we_bckg = 180., ye_sn_bckg = 90., xd_bckg = 0.28125, yd_bckg = 0.28125, p_bckg = 1000,925, 850, 700, 600, 500, 400, 300, 250, 200, 150, 100, 70, 50, 30, 20, 10/ &namelist_si ideal_flags = 0, !ideal_flags=1 for idealtest and ideal_flags=0 for realdata do_static_data = ${12}, do_3dv = .false., do_9210 = .false., do_surface_t = ${14}, ! .true. not use surface temperature in T213 data set hinterp_method = 2, vinterp_method = 2, op_ver_lev = 2 / ! op_ver_lev=0 for uniform level and op_ver_lev=1 for uneven level &namelist_dfilter do_df = .false., df_period=10800.0 / ! 2*T=2*10800sec=6hr :(-T,T), unit:sec END_OF_DATA } generate_namelist_h_input $model_dt $step_max $step_output \ $YYYY1 $MM1 $DD1 $HH1 $YYYY2 $MM2 $DD2 $HH2 $do_static \ $YYYYMMDDHH ${make_surface_t} # namelist.input generate_namelist_input(){ set -x test -f namelist.input || rm -f namelist.input cat > namelist.input « END_OF_DATA &namelist_01 s_we = 1, e_we = 502, s_sn = 1, e_sn = 330, s_vert = 1, e_vert = 31, spec_bdy_width = 5, dt = $1, time_step_max = $2, time_step_count_output = $3, step_phytend_output = 0, maxit = 2, interp = 0, irec = 0, nh = 1, writetofile = .false., xs_we = 70., ys_sn = 15., cen_lat = 40., xd = 0.15, yd = 0.15, ztop = 35000.0, epson_depart = 1.0, i_parent_start = 1, j_parent_start = 1, shw = 1, parent_grid_ratio = 1 / &namelist_02 dyn_opt = 4, diff_opt = 1, km_opt = 1, damp_opt = 0, spec_zone = 1, relax_zone = 4, isfflx = 1, ifsnow = 0, icloud = 1, num_soil_layers = 4, julyr = 0, julday = 1, gmt = 0, radt = 60., bldt = 0, cudt = 5, time_step_begin_restart = 0 / &namelist_03 mp_lc_logical = .false., ! .true. for microphysics, .false. for largescale_rain mp_physics = 0, lc_physics = 0, ! now only 1 can selected(031212) ra_lw_physics = 0, ra_sw_physics = 0, bl_sfclay_physics = 0, bl_surface_physics = 0, bl_pbl_physics = 0, cu_physics = 0 / &namelist_05 init_date = ${13}, start_year = $4, start_month = $5, start_day = $6, start_hour = $7, start_minute = 00, start_second = 00, end_year = $8, end_month = $9, end_day = ${10}, end_hour = ${11}, end_minute = 00, end_second = 00, interval_seconds = 10800, real_data_init_type = 2, bdyfrq = 10800, nested = .false., specified = .true. / &namelist_bckg nz_bckg = 17, xs_we_bckg = 0., ys_sn_bckg = 0., xe_we_bckg = 180., ye_sn_bckg = 90., xd_bckg = 0.28125, yd_bckg = 0.28125, p_bckg = 1000,925, 850, 700, 600, 500, 400, 300, 250, 200, 150, 100, 70, 50, 30, 20, 10/ &namelist_si ideal_flags = 0, !ideal_flags=1 for idealtest and ideal_flags=0 for realdata do_static_data = ${12}, do_3dv = .false., do_9210 = .false., do_surface_t = ${14}, ! .true. not use surface temperature in bckg data set hinterp_method = 2, vinterp_method = 2, op_ver_lev = 2 / ! op_ver_lev=0 for uniform level and op_ver_lev=1 for uneven level &namelist_dfilter do_df = .false., df_period=10800.0 / ! 2*T=2*10800sec=6hr :(-T,T), unit:sec END_OF_DATA } generate_namelist_input $model_dt $step_max $step_output \ $YYYY1 $MM1 $DD1 $HH1 $YYYY1 $MM1 $DD1 $HH1 $do_static \ $YYYYMMDDHH ${make_surface_t} # namelist.4dvar generate_namelist_4dvar(){ set -x test -f namelist.4dvar || rm -f namelist.4dvar cat > namelist.4dvar « END_OF_DATA &record1 ANALYSIS_VERSION = ‘GRAPES_4DVAR_2.0’ / &record2 VERT_COOR = ‘Terrian-Following Height’ / &record3 ANALYSIS_DATE = ‘$1’ / &record4 miy_h= 330, mjx_h= 502, miy_l= 330, mjx_l= 502, mks = 1, mkp = 31 / &record5 lats = 15.0, lonw = 70.0, dlat_h = 0.15, dlon_h = 0.15, dlat_l = 0.15, dlon_l = 0.15 / &record6 tkminlat = 15.0, tkmaxlat = 64.3, tkminlon = 70.0, tkmaxlon = 145.0, ptop_sound= 10.0 / &record7 iobs_sound= 1, id_obsuv_sound = 1, id_obsmss_sound= 1, id_obshum_sound= 1 / &record8 iobs_synop= 1, id_obsuv_synop = 0, id_obsmss_synop= 1, id_obshum_synop= 0 / &record9 iobs_ships= 1, id_obsuv_ships = 1, id_obsmss_ships= 1, id_obshum_ships= 1 / &record10 iobs_airep= 1, id_obsuv_airep = 1, id_obsT_airep = 2 / ! just=2, only T for airep. &record11 iobs_satob= 0, id_obsuv_satob = 1 / &record12 iobs_bda= 0, id_obsVt_bda = 0, id_obsmss_bda= 0, id_obshum_bda= 0 / &record13 iobs_gpspw= 0 / &record14 iobs_noaa16_hirs = 0, iobs_noaa16_amsua= 0, iobs_noaa16_amsub= 0, iobs_noaa17_hirs = 0, iobs_noaa17_amsua= 0, iobs_noaa17_amsub= 0, iobs_noaa15_amsua= 0, iobs_noaa15_amsub= 0 / &record15 lbiasprep = 0, lbias = 0 / &record16 iobs_vad= 0, id_obsuv_vad = 1, id_obsmss_vad= 1 / &record20 id_modmss = 1, id_modhum = 3, id_bal = 2, id_output_ana = 1 / &record21 file_output = ‘./output/’ / &record22 file_sound = ‘./input/GTS/TEMP’, file_synop = ‘./input/GTS/SYNOP’, file_ships = ‘./input/GTS/SHIPS’, file_airep = ‘./input/GTS/AIREP’, file_satob = ‘./input/GTS/SATOB’, file_vad = ‘./input/VAD/VAD’ / &record23 file_bda = ‘./input/bogus/$1.BDA’ / &record24 file_gpspw = ‘./input/gpspw/GPSG’ / &record25 file_noaa16_hirs = ‘./input/rt7data/noaa16_hrs$1’, file_noaa16_amsua = ‘./input/rt7data/noaa16_ama$1’, file_noaa16_amsub = ‘./input/rt7data/noaa16_amb$1’, file_noaa17_hirs = ‘./input/rt7data/noaa17_hrs$1’, file_noaa17_amsua = ‘./input/rt7data/noaa17_ama$1’, file_noaa17_amsub = ‘./input/rt7data/noaa17_amb$1’, file_noaa15_amsua = ‘./input/rt7data/noaa15_ama$1’, file_noaa15_amsub = ‘./input/rt7data/noaa15_amb$1’, file_pftrop = ‘./input/rt7para/prof_trop.dat’, file_input_xxx= ‘./input/rt7para/input_hirs.dat’, ‘./input/rt7para/input_amsua.dat’, ‘./input/rt7para/input_amsub.dat’, ‘./input/rt7para/input_hirs.dat’, ‘./input/rt7para/input_amsua.dat’, ‘./input/rt7para/input_amsub.dat’, ‘./input/rt7para/input_hirs.dat’, ‘./input/rt7para/input_amsua.dat’, ‘./input/rt7para/input_amsub.dat’ / &record30 file_xbchk = ‘./output/check/xbchk.dat’, file_dxa_binary= ‘./output/dxa$1.dat’, file_xa_binary = ‘./output/xa$1.dat’, file_dxa_grads = ‘./output/grads/dxa_grads.dat’, file_xa_grads = ‘./output/grads/xa_grads.dat’, file_dxa_NetCDF= ‘./output/dxa_NetCDF.dat’, file_xa_NetCDF = ‘./output/xa_NetCDF.dat’ / &record31 file_gtschk = ‘./output/check/chk’, file_bdachk = ‘./output/check/chk’, file_gpspwchk = ‘./output/check/chkgps’, file_noaa16hirschk = ‘./output/check/noaa16hirschk.dat’, file_noaa16amsuachk = ‘./output/check/noaa16amsuachk.dat’, file_noaa16amsubchk = ‘./output/check/noaa16amsubchk.dat’, file_noaa17hirschk = ‘./output/check/noaa17hirschk.dat’, file_noaa17amsuachk = ‘./output/check/noaa17amsuachk.dat’, file_noaa17amsubchk = ‘./output/check/noaa17amsubchk.dat’, file_noaa15amsuachk = ‘./output/check/noaa15amsuachk.dat’, file_noaa15amsubchk = ‘./output/check/noaa15amsubchk.dat’ / &record32 file_dsound= ‘./output/check/dsound$1’, file_dsynop= ‘./output/check/dsynop$1’, file_dships= ‘./output/check/dships$1’, file_dairep= ‘./output/check/dairep$1’, file_dsatob= ‘./output/check/dsatob$1’, file_dbda = ‘./output/check/dbda$1’, file_dgpspw= ‘./output/check/dgpspw$1’, file_dvad= ‘./output/check/dvad’, file_innotbnoaa16hirs = ‘./output/check/innotbnoaa16hirs.dat’, file_innotbnoaa16amsua= ‘./output/check/innotbnoaa16amsua.dat’, file_innotbnoaa16amsub= ‘./output/check/innotbnoaa16amsub.dat’, file_innotbnoaa17hirs = ‘./output/check/innotbnoaa17hirs.dat’, file_innotbnoaa17amsua= ‘./output/check/innotbnoaa17amsua.dat’, file_innotbnoaa17amsub= ‘./output/check/innotbnoaa17amsub.dat’ file_innotbnoaa15amsua= ‘./output/check/innotbnoaa15amsua.dat’, file_innotbnoaa15amsub= ‘./output/check/innotbnoaa15amsub.dat’ / &record33 file_mininfo = ‘./output/check/mini_info$1.dat’, file_minerror= ‘./output/check/mini_error$1.dat’ / &record34 surfout = ‘./output/’ / &record35 dqcu_sound = 4.0, dqcv_sound = 4.0, dqcmss_sound= 2.5, dqchum_sound= 3.5, dqcu_synop = 4.0, dqcv_synop = 4.0, dqcmss_synop= 4.0, dqchum_synop= 4.0, dqcu_ships = 4.0, dqcv_ships = 4.0, dqcmss_ships= 4.0, dqchum_ships= 4.0, dqcu_airep = 4.0, dqcv_airep = 4.0, dqcT_airep = 4.0, dqcu_vad = 4.0, dqcv_vad = 4.0, dqcmss_vad = 2.5, dqcu_satob = 4.0, dqcv_satob = 4.0 / &record36 dqcVt_bda = 3.5, dqcmss_bda = 2.5, dqchum_bda = 3.5 / &record37 dqc_gpspw = 4.0 / &record40 mcut=31, Lkppsi=5, Lkpchi=5, Lkpmss=5, Lkphum=2, rf_passes=10, rf_lengthpsi=500000.0, rf_lengthchi=500000.0, rf_lengthmss=500000.0, rf_lengthhum=200000.0 / &record41 eps0 = 0.001, ntmax = 50 / &record42 inmc_mss = 1, dpath_mss=’./input/rmsdata/piu1112.dat’, inmc_psi= 1, dpath_psi=’./input/rmsdata/psi1112.dat’, inmc_chi= 1, dpath_chi=’./input/rmsdata/chi1112.dat’, inmc_hum= 1, dpath_rh=’./input/rmsdata/rh1112.dat’, dpath_q=’./input/rmsdata/q1112.dat’ / &record51 int_obs = -36, first_outloop = .true. / END_OF_DATA } generate_namelist_4dvar $begintime ################# # clean up ################# ################ # tail ################ date set +x set +e set +u [/shell]
外循环
[shell collapse=”true”]
#!/bin/ksh
# @ job_type = parallel
# @ comment = GRAPES
# @ input = /dev/null
# @ output = ./out/4dvar_h_$(jobid).out
# @ error = ./out/4dvar_h_$(jobid).err
# @ initialdir = ./
# @ notification = complete
# @ checkpoint = no
# @ node = 1
# @ tasks_per_node = 32
# @ comment = GRAPES
# @ class = normal
## @ class = smalljob
## @ class = minijob
## @ wall_clock_limit=10800
# @ job_name = GRAPES
# @ network.MPI = sn_all,,US
# @ queue
## must be set equal to (CPUs-per-node / tasks_per_node)
export OMP_NUM_THREADS=1
## suggestion from Jim Edwards to reintroduce XLSMPOPTS on 11/13/03
export XLSMPOPTS=”stack=256000000″
export AIXTHREAD_SCOPE=S
export MALLOCMULTIHEAP=TRUE
export OMP_DYNAMIC=FALSE
## Do our best to get sufficient stack memory
ulimit -s unlimited
ulimit -c unlimited
set -x
set -e
set -u
###############
# directory
###############
BASE_DIR=./nwp/system/grapes_meso/my_grapes_meso
BASE_SOURCE_DIR=./nwp/system/grapes_meso/GRAPES_MESO3.3.2.4
RUN_DIR=${BASE_DIR}/run/an
#RUN_DIR=${BASE_SOURCE_DIR}/4dvar/m4dv/rundir
OBS_DIR=${BASE_DIR}/run/obs
CONDAT_DIR=${BASE_DIR}/condat
AN_CONDAT_DIR=${CONDAT_DIR}/an
SI_RUN_DIR=${RUN_DIR}/si
M4DVAR_H_BIN=${BASE_SOURCE_DIR}/4dvar/m4dv/rundir/4dvar_h.exe
##########################
# enter work directory
##########################
test -d ${RUN_DIR} || mkdir -p ${RUN_DIR}
cd ${RUN_DIR}
cp -f ${SI_RUN_DIR}/grapesinput .
###############
# Prepare data
###############
rm -f rsl*
# OBS GTS
ln -sf ${OBS_DIR} ${RUN_DIR}/input/GTS
# rmsdate
ln -sf ${AN_CONDAT_DIR}/input/rmsdata ${RUN_DIR}/input/rmsdata
# tbl data
test -f GENPARM.TBL || cp ${CONDAT_DIR}/fcst/GENPARM.TBL .
test -f GEOGRID.TBL || cp ${CONDAT_DIR}/fcst/GEOGRID.TBL .
test -f LANDUSE.TBL || cp ${CONDAT_DIR}/fcst/LANDUSE.TBL .
test -f SOILPARM.TBL || cp ${CONDAT_DIR}/fcst/SOILPARM.TBL .
test -f VEGPARM.TBL || cp ${CONDAT_DIR}/fcst/VEGPARM.TBL .
test -f RRTM_DATA || cp ${CONDAT_DIR}/fcst/RRTM_DATA .
###############
# run
################
${M4DVAR_H_BIN}
[/shell]
内循环
[shell collapse=”true”]
#!/bin/ksh
# @ job_type = parallel
# @ comment = GRAPES
# @ input = /dev/null
# @ output = ./out/4dvar_$(jobid).out
# @ error = ./out/4dvar_$(jobid).err
# @ initialdir = ./
# @ notification = complete
# @ checkpoint = no
# @ node = 1
# @ tasks_per_node = 32
# @ comment = GRAPES
# @ class = normal
## @ class = smalljob
## @ class = minijob
## @ wall_clock_limit=10800
# @ job_name = GRAPES
# @ network.MPI = sn_all,,US
# @ queue
## must be set equal to (CPUs-per-node / tasks_per_node)
export OMP_NUM_THREADS=1
## suggestion from Jim Edwards to reintroduce XLSMPOPTS on 11/13/03
export XLSMPOPTS=”stack=256000000″
export AIXTHREAD_SCOPE=S
export MALLOCMULTIHEAP=TRUE
export OMP_DYNAMIC=FALSE
## Do our best to get sufficient stack memory
ulimit -s unlimited
ulimit -c unlimited
set -x
set -e
set -u
###############
# directory
###############
BASE_DIR=./nwp/system/grapes_meso/my_grapes_meso
BASE_SOURCE_DIR=./nwp/system/grapes_meso/GRAPES_MESO3.3.2.4
RUN_DIR=${BASE_DIR}/run/an
#RUN_DIR=${BASE_SOURCE_DIR}/4dvar/m4dv/rundir
CONDAT_DIR=${BASE_DIR}/condat
M4DVAR_BIN=${BASE_SOURCE_DIR}/4dvar/m4dv/rundir/4dvar.exe
##########################
# enter work directory
##########################
test -d ${RUN_DIR} || mkdir -p ${RUN_DIR}
cd ${RUN_DIR}
###############
# Prepare data
###############
rm -f rsl*
# *.dat
test -f Eigenchi.dat || cp ${CONDAT_DIR}/an/Eigenchi.dat .
test -f Eigenhum.dat || cp ${CONDAT_DIR}/an/Eigenhum.dat .
test -f Eigenmss.dat || cp ${CONDAT_DIR}/an/Eigenmss.dat .
test -f Eigenpsi.dat || cp ${CONDAT_DIR}/an/Eigenpsi.dat .
###############
# run
################
${M4DVAR_BIN}
[/shell]
更新边界
[shell collapse=”true”]
#!/bin/ksh
# @ job_type=serial
# @ output = ./out/m4dvar.out
# @ error = ./out/m4dvar.err
# @ notification = error
# @ checkpoint = no
# @ restart = no
# @ comment = GRAPES
# @ class= normal
set -x
set -e
set -u
date
######################
# directory param
######################
BASE_DIR=./nwp/system/grapes_meso/my_grapes_meso
BASE_SOURCE_DIR=./nwp/system/grapes_meso/GRAPES_MESO3.3.2.4
BASE_RUN_DIR=${BASE_DIR}/run
AN_RUN_DIR=${BASE_RUN_DIR}/an
SI_RUN_DIR=${BASE_RUN_DIR}/si
OBS_RUN_DIR=${BASE_RUN_DIR}/obs
FCST_RUN_DIR=${BASE_RUN_DIR}/fcst
RUN_DIR=${AN_RUN_DIR}/mdvar_proc
CONDAT_DIR=${BASE_DIR}/condat
AN_CONDAT_DIR=${CONDAT_DIR}/an
MDVAR_BIN=${BASE_SOURCE_DIR}/fcst/modelvar_proc/mdvar.exe
NEWDATE_BIN=${BASE_SOURCE_DIR}/fcst/grapes_model/run/newdate
############
# parameter
############
YYYYMMDDHH=2013041012
mfcast_len=24
do_static=.true.
do_3dv=.true.
output_inteval=3 #model output inteval in hour
model_dt=90 #model time step in second
make_surface_t=.false.
begintime=$YYYYMMDDHH
endtime=$( ${NEWDATE_BIN} $begintime +$mfcast_len)
step_max=`expr $mfcast_len * 3600 / $model_dt`
step_output=`expr $output_inteval * 3600 / $model_dt`
YYYY1=`echo $begintime | cut -c1-4`
MM1=`echo $begintime | cut -c5-6`
DD1=`echo $begintime | cut -c7-8`
HH1=`echo $begintime | cut -c9-10`
YYYY2=`echo $endtime | cut -c1-4`
MM2=`echo $endtime | cut -c5-6`
DD2=`echo $endtime | cut -c7-8`
HH2=`echo $endtime | cut -c9-10`
#####################
# enter work space
#####################
test -d ${RUN_DIR} || mkdir -p ${RUN_DIR}
cd ${RUN_DIR}
rm -f grapesinput
rm -f grapesbdy_old
#############
# prepare data
#############
cp ${AN_RUN_DIR}/output/grapesinput-4dv-${begintime} grapesinput${begintime}
cp ${SI_RUN_DIR}/grapesbdy grapesbdy_old
########################
# generate namelist
###########################
generate_namelist_mdvar(){
test -f namelist_mdvar || rm -r namelist_mdvar
cat > namelist_mdvar «eof &namelist_01 s_we = 1, e_we = 502, s_sn = 1, e_sn = 330, s_vert = 1, e_vert = 31, xs_we = 70., ys_sn = 15., xd = 0.15, yd = 0.15 / &namelist_02 start_year = $1, start_month = $2, start_day = $3, start_hour = $4, end_year = $5, end_month = $6, end_day = $7, end_hour = $8, interval_seconds = 10800 / &namelist_03 num_soil_layers = 4, spec_bdy_width = 5 / EOF } generate_namelist_mdvar $YYYY1 $MM1 $DD1 $HH1 $YYYY2 $MM2 $DD2 $HH2 ############### # run program ############## ${MDVAR_BIN} ############# # post ############# sleep 10 test -d ${FCST_RUN_DIR} || mkdir -p ${FCST_RUN_DIR} mv grapesbdy_new ${FCST_RUN_DIR}/grapesbdy mv grapesinput${begintime} ${FCST_RUN_DIR}/grapesinput ############# # END ############## date [/shell]
