登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
先来说一下我的过程吧——
我是使用的脚本grapes_single.sh来运行GRAPES,前几步准备数据和生成背景场都没有显示出现问题
这是前几步的过程:
#--- [0.0] Set directory ------------------------------------
GRAPESROOT=/home2/xx/x/grapes_3hr GEODATA_DIR=$GRAPESROOT/geodata #directory of geodata T213_DIR=$GRAPESROOT/../../dataset/T213 #directory of T213 files forGRAPES OBS_DIR=$GRAPESROOT/../../dataset/GTS #directory of obs files(aob) GRAPES_MODEL_DIR=$GRAPESROOT/grapes_model #directory of GRAPES model GRAPES_M3DV_DIR=$GRAPESROOT/grapes_m3dv #directory of GRAPES 3dvar GRAPES_GPS_DIR=$GRAPESROOT/grapes_gps #directory of GRAPESpreprocess PATH_NEWDATE=$GRAPESROOT/grapes_sh #directory of script files(newdate) if [ !-e ${GRAPESROOT} ]; then echo 'FATAL ERROR : GRAPES DIR ${GRAPESROOT} DOES NOT EXIST !!!!!' exit 1 fi if [ !-e ${GEODATA_DIR} ]; then echo 'WARNING : GEODATA DIR ${GEODATA_DIR} DOES NOT EXIST !!!!!' fi #--- [0.1] Set parameters ------------------------------------ YYYYMMDDHH=$1 mfcast_len=$2 #forecast time length in hour do_static=$3 interval_seconds=21600 output_interval=6 #modeloutput inteval in hour model_dt=180 #modeltime step in second begintime=$YYYYMMDDHH # --- [0.2] calculate parameter of time ----------- hour=`echo $begintime | cut -c 9-10` if [$hour -eq 06 -o $hour -eq 18 ]; then YYYYMMDDHH=$( ${PATH_NEWDATE}/newdate $begintime -6) else YYYYMMDDHH=$( ${PATH_NEWDATE}/newdate $begintime -12) fi YYYYMMDDHH_xb=$( ${PATH_NEWDATE}/newdate $begintime -12) #--- [1.0] prepare data and make xb for m3dvar------------ echo---------------------------------- echoStep 1: Preprocess T213 data for 3dvar endtime=$( ${PATH_NEWDATE}/newdate $begintime +$mfcast_len) step_max=`expr $mfcast_len \* 3600 \/ $model_dt` step_output=`expr $output_interval \* 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` # --- link geodata directory --------------- cd ${GRAPESROOT} if [! -e ${GRAPESROOT}/geodata ]; then ln-sf $GEODATA_DIR ${GRAPESROOT} fi # --- link T213 data ----------------------- cd ${GRAPES_MODEL_DIR}/run if [! -e T213 ]; then mkdir T213 fi rm -f./T213/T213* ln-sf $T213_DIR/T213${YYYYMMDDHH}* ./T213 #- --- make namelist.input file and run sisystem---------------- rm -fgrapesbdy grapesinput make_surface_t=.false. ./namelist_grapes.sh $model_dt $step_max $step_output \ $YYYY1 $MM1 $DD1 $HH1 $YYYY2 $MM2 $DD2 $HH2 $do_static \ $YYYYMMDDHH ${make_surface_t} ${interval_seconds} ./si.exe > ./worklog/si.log.$YYYY1$MM1$DD1$HH1 2>&1 if [-e grapesinput ]; then echo grapesinput created! mvgrapesinput ${GRAPES_M3DV_DIR}/data/input/backgrd/grapesinput${begintime} else echo grapesinput NOT created! exit 1 fi #--- [2.0] preprocess Obs data for GRAPESm3dvar --------------------- echo---------------------------------- echoStep 2: Preprocess Obs data for 3dvar rm -f${GRAPES_M3DV_DIR}/data/input/GTS/* cd${GRAPES_M3DV_DIR}/preproc/ObsGts if [! -e ObsProc.exe ]; then make fi ./namelist.sh $OBS_DIR $begintime ./ObsProc.exe sleep 1到这里都还是正常的,在往下就出现了问题,机器显示:
char=TEMP
END of Setting up TEMP wind & mass field
max_obs= 5000
Num of valid wind obs. = 213
Num of valid mass obs. = 188
max_level= 250
nwnd>0 213
ndata_u ndata_phi ndata_t ndata_hum 3789 0 0
0
nmss>0 188
ndata_u ndata_phi ndata_t ndata_hum 0 5928 6001
4099
3 : both wind and mass field data
Put wind & mass fields together and write out
Ndata= 216
Num of valid obs. as one file output = 216
char=SYNOP
END of Setting up SYNOP wind & mass field
max_obs= 15000
Num of valid wind obs. = 1379
Num of valid mass obs. = 1382
max_level= 2
nwnd>0 1379
ndata_u ndata_phi ndata_t ndata_hum 1379 0 0
0
nmss>0 1382
ndata_u ndata_phi ndata_t ndata_hum 0 2428 1362
1369
3 : both wind and mass field data
Put wind & mass fields together and write out
Ndata= 1382
Num of valid obs. as one file output = 1382
char=SHIPS
END of Setting up SHIPS wind & mass field
max_obs= 5000
Num of valid wind obs. = 35
Num of valid mass obs. = 45
Num of valid obs. as one file output = 45
char=AIREP
END of Setting up AIREP wind & mass field
max_obs= 20000
Num of valid wind obs. = 484
Num of valid mass obs. = 484
Num of valid obs. as one file output = 484
char=SATOB
END of Setting up SATOB wind & mass field
max_obs= 20000
Num of valid wind obs. = 0
Num of valid mass obs. = 0
There are NO DATA forSATOB
Num of valid obs. as one file output = 484
Buffs Preprocessing is finished
Obs Pre-process is finished
++ sleep 10 以上为第二步运行结束的结果
++ echo ----------------------------------
----------------------------------
++ echo Step 3: Run Grapes m3dvar
Step 3: Run Grapes m3dvar
++ cd /home2/xx/x/grapes_3hr/grapes_m3dv/m3dv/rundir
++ ./namelist_3dv.sh 2008010300
++ ./Grapes3dvar.exe 卡在这里了
查看动态文件,机器上显示的是
10 nallfun=101 J= NaN gnorm= NaN step= 8.738D+04 finish= F
jo_GTS_u,v,mss,hum= NaN : NaN NaN NaN NaN
jo_GTS_u,v,mss,hum= NaN : NaN NaN NaN NaN
J = Jb + Jo : NaN : NaN NaN
jo_GTS_u,v,mss,hum= NaN : NaN NaN NaN NaN
jo_GTS_u,v,mss,hum= NaN : NaN NaN NaN NaN
J = Jb + Jo : NaN : NaN NaN
jo_GTS_u,v,mss,hum= NaN : NaN NaN NaN NaN
jo_GTS_u,v,mss,hum= NaN : NaN NaN NaN NaN
我是刚刚接触模式,水平有限,于是在家园里发个帖子试试运气,恳请高
人指点一二,不胜感激。以上如有表述不清的,还请大家多多指正
俺的邮箱是 kingqueen1201@163.com
[img]file:///C:\Documents and Settings\Administrator\Application Data\Tencent\Users\1162069562\QQ\WinTemp\RichOle\VWJ$4%M[S~RUG%00T]4ITUT.jpg[/img]
|