请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 8074|回复: 6

请教:关于GRAPES三维同化中Jo的值总为NaN可能的原因

[复制链接]

新浪微博达人勋

发表于 2013-5-11 11:58:57 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

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]



VWJ$4%M[S~RUG%00T]4ITUT.jpg
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-5-11 12:00:16 | 显示全部楼层
可要有人看到哟
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-5-11 20:05:23 | 显示全部楼层
我没做过模式,不过个人感觉计算过程中值为nan,很有可能是除法中有的除数为零了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-5-11 21:21:43 | 显示全部楼层

谢谢提示,我刚上网查了NaN出现的可能原因,觉得有可能是过程中出现了不合数学运算的步骤,排查原因ing
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-5-26 18:07:13 | 显示全部楼层
我也是用GRAPES的  可以加好友一起交流下
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2020-11-30 08:14:54 | 显示全部楼层
飞王子 发表于 2013-5-11 21:21
谢谢提示,我刚上网查了NaN出现的可能原因,觉得有可能是过程中出现了不合数学运算的步骤,排查原因ing

请问lz是怎么解决的?我遇到了同样的问题。“不合数学运算的步骤”是啥?非常感谢
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2020-11-30 16:37:13 | 显示全部楼层
clarawho 发表于 2020-11-30 08:14
请问lz是怎么解决的?我遇到了同样的问题。“不合数学运算的步骤”是啥?非常感谢

比如,零除以零,这是很早之前的问题了,现在只记得运算中出现零除以零时就会出现NAN,或者记成缺测
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表