WRF 运行
WPS
设置链接
ln -sf /home/ahri/WRFV4/WPS/geogrid.exe
ln -sf /home/ahri/WRFV4/WPS/ungrib.exe
ln -sf /home/ahri/WRFV4/WPS/metgrid.exe
ln -sf /home/ahri/WRFV4/WPS/link_grib.csh
geogrid.exe
定义模型域并将静态地理数据插值到网格
配置文件 namelist.wps (doc)
查看静态地图文件来确定区域是否正确
geogrid.exe
ungrib.exe
从 GRIB 格式的文件中提取气象字段（单核运行）
设置 Vtable
链接数据
ungrib.exe
查看内容（可选）
metgrid.exe
将 ungrib.exe 提取的气象数据水平插值到 geogrid.exe 定义的模型网格中
WRF(&GSI)
real.exe
气象场的垂直插值、生成初始条件和边界条件、一致性检查
配置文件 namelist.input
ra_lw_physics Longwave Radiation 长波辐射
ra_sw-physics Shortwave Radiation 短波辐射
sf_surface_physics land-surface option 路面过程
real.exe
# in: met_em.d0*.YYYY-MM-DD_HH:MM:SS.nc
mpirun -np 4 ./real.exe
# out: wrfbdy_d0*
# wrfinput_d0*
# ...
wrfda/gsi(可选)
观测数据同化
3dvar
官方示例
export WRFDA_DIR=/home/ahri/WRFV4/WRFDA
export DATA_DIR=/home/ahri/WRFV4/test_wrfda/data
ln -sf $WRFDA_DIR/run/LANDUSE.TBL
ln -sf $WRFDA_DIR/var/build/da_wrfvar.exe
# 链接背景误差文件
ln -sf $DATA_DIR/be/be.dat
# 适用于 NCEP 的 背景误差文件
ln -sf /home/ahri/WRFV4/WRFDA/var/run/be.dat.cv3 be.dat
# 链接 ascii 观测文件
ln -sf $DATA_DIR/ob/2008020512/ob.ascii
# 或者链接 prepbufr 常规观测文件
ln -fs $DATA_DIR/ob/2008020512/gdas1.t12z.prepbufr.nr ob.bufr
# 链接第一次猜测文件
ln -sf $DATA_DIR/rc/2008020512/wrfinput_d01 ./fg
# 复制 namelist.input
cp $DATA_DIR/namelist.input.3dvar namelist.inpb ut
ulimit -s unlimited
mpirun -np 16 ./da_wrfvar.exe
生成 be.dat
复制 gen_be_wrapper.ksh
参考：
- https://zhuanlan.zhihu.com/p/565561915
4dvar
官方示例(运行失败)
export WRFDA_DIR=/home/ahri/WRFV4/WRFDA
export DATA_DIR=/home/ahri/WRFV4/test_wrfda/data
cd /home/ahri/WRFV4/test_wrfda
ln -fs $DATA_DIR/ob/2008020512/ob01.ascii ob01.ascii
ln -fs $DATA_DIR/ob/2008020513/ob02.ascii ob02.ascii
ln -fs $DATA_DIR/ob/2008020514/ob03.ascii ob03.ascii
ln -fs $DATA_DIR/ob/2008020515/ob04.ascii ob04.ascii
ln -fs $DATA_DIR/ob/2008020516/ob05.ascii ob05.ascii
ln -fs $DATA_DIR/ob/2008020517/ob06.ascii ob06.ascii
ln -fs $DATA_DIR/ob/2008020518/ob07.ascii ob07.ascii
ln -fs $DATA_DIR/ob/2008020512/gdas1.t12z.prepbufr.nr ob01.bufr
ln -fs $DATA_DIR/ob/2008020518/gdas1.t18z.prepbufr.nr ob02.bufr
ln -fs $DATA_DIR/rc/2008020512/wrfinput_d01
ln -fs $DATA_DIR/rc/2008020512/wrfbdy_d01
ln -fs wrfinput_d01 fg
ln -fs $DATA_DIR/be/be.dat
ln -fs $WRFDA_DIR/run/LANDUSE.TBL
ln -fs $WRFDA_DIR/run/GENPARM.TBL
ln -fs $WRFDA_DIR/run/SOILPARM.TBL
ln -fs $WRFDA_DIR/run/VEGPARM.TBL
ln -fs $WRFDA_DIR/run/RRTM_DATA_DBL RRTM_DATA
ln -sf $WRFDA_DIR/run/CAMtr_volume_mixing_ratio
ln -sf $WRFDA_DIR/var/build/da_wrfvar.exe
cp $DATA_DIR/namelist.input.4dvar namelist.input
生成 be.dat 背景误差文件
gsi.x
观测数据同化
DTC || Community GSI Users Page
CRTM (Community Radiative Transfer Model 社区辐射转移模型)
调用 ksh 脚本运行 gsi.x
comgsi_run_regional.ksh
#!/bin/ksh
#####################################################
# machine set up (users should change this part)
#####################################################
set -x
#
# GSIPROC = processor number used for GSI analysis
#------------------------------------------------
GSIPROC=1
ARCH='LINUX_LSF'
# Supported configurations:
# IBM_LSF,
# LINUX, LINUX_LSF, LINUX_PBS,
# DARWIN_PGI
#
#####################################################
# case set up (users should change this part)
#####################################################
#
# ANAL_TIME= analysis time (YYYYMMDDHH)
# WORK_ROOT= working directory, where GSI runs
# PREPBURF = path of PreBUFR conventional obs
# BK_FILE = path and name of background file
# OBS_ROOT = path of observations files
# FIX_ROOT = path of fix files
# GSI_EXE = path and name of the gsi executable
# ENS_ROOT = path where ensemble background files exist
ANAL_TIME=2017051318
JOB_DIR=the_job_directory
#normally you put run scripts here and submit jobs form here, require a copy of gsi.x at this directory
RUN_NAME=a_descriptive_run_name_such_as_case05_3denvar_etc
OBS_ROOT=the_directory_where_observation_files_are_located
BK_ROOT=the_directory_where_background_files_are_located
GSI_ROOT=the_comgsi_main directory where src/ ush/ fix/ etc are located
CRTM_ROOT=the_CRTM_directory
ENS_ROOT=the_directory_where_ensemble_backgrounds_are_located
#ENS_ROOT is not required if not running hybrid EnVAR
HH=`echo $ANAL_TIME | cut -c9-10`
GSI_EXE=${JOB_DIR}/gsi.x #assume you have a copy of gsi.x here
WORK_ROOT=${JOB_DIR}/${RUN_NAME}
FIX_ROOT=${GSI_ROOT}/fix
GSI_NAMELIST=${GSI_ROOT}/ush/comgsi_namelist.sh
PREPBUFR=${OBS_ROOT}/nam.t${HH}z.prepbufr.tm00
BK_FILE=${BK_ROOT}/wrfinput_d01.${ANAL_TIME}
#
#------------------------------------------------
# bk_core= which WRF core is used as background (NMM or ARW or NMMB)
# bkcv_option= which background error covariance and parameter will be used
# (GLOBAL or NAM)
# if_clean = clean : delete temperal files in working directory (default)
# no : leave running directory as is (this is for debug only)
# if_observer = Yes : only used as observation operater for enkf
# if_hybrid = Yes : Run GSI as 3D/4D EnVar
# if_4DEnVar = Yes : Run GSI as 4D EnVar
# if_nemsio = Yes : The GFS background files are in NEMSIO format
# if_oneob = Yes : Do single observation test
if_hybrid=No # Yes, or, No -- case sensitive !
if_4DEnVar=No # Yes, or, No -- case sensitive (set if_hybrid=Yes first)!
if_observer=No # Yes, or, No -- case sensitive !
if_nemsio=No # Yes, or, No -- case sensitive !
if_oneob=No # Yes, or, No -- case sensitive !
bk_core=ARW
bkcv_option=NAM
if_clean=clean
#
# setup whether to do single obs test
if [ ${if_oneob} = Yes ]; then
if_oneobtest='.true.'
else
if_oneobtest='.false.'
fi
#
# setup for GSI 3D/4D EnVar hybrid
if [ ${if_hybrid} = Yes ] ; then
PDYa=`echo $ANAL_TIME | cut -c1-8`
cyca=`echo $ANAL_TIME | cut -c9-10`
gdate=`date -u -d "$PDYa $cyca -6 hour" +%Y%m%d%H` #guess date is 6hr ago
gHH=`echo $gdate |cut -c9-10`
datem1=`date -u -d "$PDYa $cyca -1 hour" +%Y-%m-%d_%H:%M:%S` #1hr ago
datep1=`date -u -d "$PDYa $cyca 1 hour" +%Y-%m-%d_%H:%M:%S` #1hr later
if [ ${if_nemsio} = Yes ]; then
if_gfs_nemsio='.true.'
ENSEMBLE_FILE_mem=${ENS_ROOT}/gdas.t${gHH}z.atmf006s.mem
else
if_gfs_nemsio='.false.'
ENSEMBLE_FILE_mem=${ENS_ROOT}/sfg_${gdate}_fhr06s_mem
fi
if [ ${if_4DEnVar} = Yes ] ; then
BK_FILE_P1=${BK_ROOT}/wrfout_d01_${datep1}
BK_FILE_M1=${BK_ROOT}/wrfout_d01_${datem1}
if [ ${if_nemsio} = Yes ]; then
ENSEMBLE_FILE_mem_p1=${ENS_ROOT}/gdas.t${gHH}z.atmf009s.mem
ENSEMBLE_FILE_mem_m1=${ENS_ROOT}/gdas.t${gHH}z.atmf003s.mem
else
ENSEMBLE_FILE_mem_p1=${ENS_ROOT}/sfg_${gdate}_fhr09s_mem
ENSEMBLE_FILE_mem_m1=${ENS_ROOT}/sfg_${gdate}_fhr03s_mem
fi
fi
fi
# The following two only apply when if_observer = Yes, i.e. run observation operator for EnKF
# no_member number of ensemble members
# BK_FILE_mem path and base for ensemble members
no_member=20
BK_FILE_mem=${BK_ROOT}/wrfarw.mem
#
#
#####################################################
# Users should NOT make changes after this point
#####################################################
#
BYTE_ORDER=Big_Endian
# BYTE_ORDER=Little_Endian
case $ARCH in
'IBM_LSF')
###### IBM LSF (Load Sharing Facility)
RUN_COMMAND="mpirun.lsf " ;;
'LINUX')
if [ $GSIPROC = 1 ]; then
#### Linux workstation - single processor
RUN_COMMAND=""
else
###### Linux workstation - mpi run
RUN_COMMAND="mpirun -np ${GSIPROC} "
fi ;;
'LINUX_LSF')
###### LINUX LSF (Load Sharing Facility)
RUN_COMMAND="mpirun.lsf " ;;
'LINUX_PBS')
#### Linux cluster PBS (Portable Batch System)
RUN_COMMAND="mpirun -np ${GSIPROC} " ;;
'DARWIN_PGI')
### Mac - mpi run
if [ $GSIPROC = 1 ]; then
#### Mac workstation - single processor
RUN_COMMAND=""
else
###### Mac workstation - mpi run
RUN_COMMAND="mpirun -np ${GSIPROC} -machinefile ~/mach "
fi ;;
* )
print "error: $ARCH is not a supported platform configuration."
exit 1 ;;
esac
##################################################################################
# Check GSI needed environment variables are defined and exist
#
# Make sure ANAL_TIME is defined and in the correct format
if [ ! "${ANAL_TIME}" ]; then
echo "ERROR: \$ANAL_TIME is not defined!"
exit 1
fi
# Make sure WORK_ROOT is defined and exists
if [ ! "${WORK_ROOT}" ]; then
echo "ERROR: \$WORK_ROOT is not defined!"
exit 1
fi
# Make sure the background file exists
if [ ! -r "${BK_FILE}" ]; then
echo "ERROR: ${BK_FILE} does not exist!"
exit 1
fi
# Make sure OBS_ROOT is defined and exists
if [ ! "${OBS_ROOT}" ]; then
echo "ERROR: \$OBS_ROOT is not defined!"
exit 1
fi
if [ ! -d "${OBS_ROOT}" ]; then
echo "ERROR: OBS_ROOT directory '${OBS_ROOT}' does not exist!"
exit 1
fi
# Set the path to the GSI static files
if [ ! "${FIX_ROOT}" ]; then
echo "ERROR: \$FIX_ROOT is not defined!"
exit 1
fi
if [ ! -d "${FIX_ROOT}" ]; then
echo "ERROR: fix directory '${FIX_ROOT}' does not exist!"
exit 1
fi
# Set the path to the CRTM coefficients
if [ ! "${CRTM_ROOT}" ]; then
echo "ERROR: \$CRTM_ROOT is not defined!"
exit 1
fi
if [ ! -d "${CRTM_ROOT}" ]; then
echo "ERROR: fix directory '${CRTM_ROOT}' does not exist!"
exit 1
fi
# Make sure the GSI executable exists
if [ ! -x "${GSI_EXE}" ]; then
echo "ERROR: ${GSI_EXE} does not exist!"
exit 1
fi
# Check to make sure the number of processors for running GSI was specified
if [ -z "${GSIPROC}" ]; then
echo "ERROR: The variable $GSIPROC must be set to contain the number of processors to run GSI"
exit 1
fi
#
##################################################################################
# Create the ram work directory and cd into it
workdir=${WORK_ROOT}
echo " Create working directory:" ${workdir}
if [ -d "${workdir}" ]; then
rm -rf ${workdir}
fi
mkdir -p ${workdir}
cd ${workdir}
#
##################################################################################
echo " Copy GSI executable, background file, and link observation bufr to working directory"
# Save a copy of the GSI executable in the workdir
cp ${GSI_EXE} gsi.x
# Bring over background field (it's modified by GSI so we can't link to it)
cp ${BK_FILE} ./wrf_inout
if [ ${if_4DEnVar} = Yes ] ; then
cp ${BK_FILE_P1} ./wrf_inou3
cp ${BK_FILE_M1} ./wrf_inou1
fi
# Link to the prepbufr data
ln -s ${PREPBUFR} ./prepbufr
# ln -s ${OBS_ROOT}/gdas1.t${HH}z.sptrmm.tm00.bufr_d tmirrbufr
# Link to the radiance data
srcobsfile[1]=${OBS_ROOT}/gdas1.t${HH}z.satwnd.tm00.bufr_d
gsiobsfile[1]=satwnd
srcobsfile[2]=${OBS_ROOT}/gdas1.t${HH}z.1bamua.tm00.bufr_d
gsiobsfile[2]=amsuabufr
srcobsfile[3]=${OBS_ROOT}/gdas1.t${HH}z.1bhrs4.tm00.bufr_d
gsiobsfile[3]=hirs4bufr
srcobsfile[4]=${OBS_ROOT}/gdas1.t${HH}z.1bmhs.tm00.bufr_d
gsiobsfile[4]=mhsbufr
srcobsfile[5]=${OBS_ROOT}/gdas1.t${HH}z.1bamub.tm00.bufr_d
gsiobsfile[5]=amsubbufr
srcobsfile[6]=${OBS_ROOT}/gdas1.t${HH}z.ssmisu.tm00.bufr_d
gsiobsfile[6]=ssmirrbufr
# srcobsfile[7]=${OBS_ROOT}/gdas1.t${HH}z.airsev.tm00.bufr_d
gsiobsfile[7]=airsbufr
srcobsfile[8]=${OBS_ROOT}/gdas1.t${HH}z.sevcsr.tm00.bufr_d
gsiobsfile[8]=seviribufr
srcobsfile[9]=${OBS_ROOT}/gdas1.t${HH}z.iasidb.tm00.bufr_d
gsiobsfile[9]=iasibufr
srcobsfile[10]=${OBS_ROOT}/gdas1.t${HH}z.gpsro.tm00.bufr_d
gsiobsfile[10]=gpsrobufr
srcobsfile[11]=${OBS_ROOT}/gdas1.t${HH}z.amsr2.tm00.bufr_d
gsiobsfile[11]=amsrebufr
srcobsfile[12]=${OBS_ROOT}/gdas1.t${HH}z.atms.tm00.bufr_d
gsiobsfile[12]=atmsbufr
srcobsfile[13]=${OBS_ROOT}/gdas1.t${HH}z.geoimr.tm00.bufr_d
gsiobsfile[13]=gimgrbufr
srcobsfile[14]=${OBS_ROOT}/gdas1.t${HH}z.gome.tm00.bufr_d
gsiobsfile[14]=gomebufr
srcobsfile[15]=${OBS_ROOT}/gdas1.t${HH}z.omi.tm00.bufr_d
gsiobsfile[15]=omibufr
srcobsfile[16]=${OBS_ROOT}/gdas1.t${HH}z.osbuv8.tm00.bufr_d
gsiobsfile[16]=sbuvbufr
srcobsfile[17]=${OBS_ROOT}/gdas1.t${HH}z.eshrs3.tm00.bufr_d
gsiobsfile[17]=hirs3bufrears
srcobsfile[18]=${OBS_ROOT}/gdas1.t${HH}z.esamua.tm00.bufr_d
gsiobsfile[18]=amsuabufrears
srcobsfile[19]=${OBS_ROOT}/gdas1.t${HH}z.esmhs.tm00.bufr_d
gsiobsfile[19]=mhsbufrears
srcobsfile[20]=${OBS_ROOT}/rap.t${HH}z.nexrad.tm00.bufr_d
gsiobsfile[20]=l2rwbufr
srcobsfile[21]=${OBS_ROOT}/rap.t${HH}z.lgycld.tm00.bufr_d
gsiobsfile[21]=larcglb
srcobsfile[22]=${OBS_ROOT}/gdas1.t${HH}z.glm.tm00.bufr_d
gsiobsfile[22]=
ii=1
while [[ $ii -le 21 ]]; do
if [ -r "${srcobsfile[$ii]}" ]; then
# ln -s ${srcobsfile[$ii]} ${gsiobsfile[$ii]}
echo "link source obs file ${srcobsfile[$ii]}"
fi
(( ii = $ii + 1 ))
done
#
##################################################################################
ifhyb=.false.
if [ ${if_hybrid} = Yes ] ; then
ls ${ENSEMBLE_FILE_mem}* > filelist02
if [ ${if_4DEnVar} = Yes ] ; then
ls ${ENSEMBLE_FILE_mem_p1}* > filelist03
ls ${ENSEMBLE_FILE_mem_m1}* > filelist01
fi
nummem=`more filelist02 | wc -l`
nummem=$((nummem -3 ))
if [[ ${nummem} -ge 5 ]]; then
ifhyb=.true.
${ECHO} " GSI hybrid uses ${ENSEMBLE_FILE_mem} with n_ens=${nummem}"
fi
fi
if4d=.false.
if [[ ${ifhyb} = .true. && ${if_4DEnVar} = Yes ]] ; then
if4d=.true.
fi
#
##################################################################################
echo " Copy fixed files and link CRTM coefficient files to working directory"
# Set fixed files
# berror = forecast model background error statistics
# specoef = CRTM spectral coefficients
# trncoef = CRTM transmittance coefficients
# emiscoef = CRTM coefficients for IR sea surface emissivity model
# aerocoef = CRTM coefficients for aerosol effects
# cldcoef = CRTM coefficients for cloud effects
# satinfo = text file with information about assimilation of brightness temperatures
# satangl = angle dependent bias correction file (fixed in time)
# pcpinfo = text file with information about assimilation of prepcipitation rates
# ozinfo = text file with information about assimilation of ozone data
# errtable = text file with obs error for conventional data (regional only)
# convinfo = text file with information about assimilation of conventional data
# lightinfo= text file with information about assimilation of GLM lightning data
# bufrtable= text file ONLY needed for single obs test (oneobstest=.true.)
# bftab_sst= bufr table for sst ONLY needed for sst retrieval (retrieval=.true.)
if [ ${bkcv_option} = GLOBAL ] ; then
echo ' Use global background error covariance'
BERROR=${FIX_ROOT}/${BYTE_ORDER}/nam_glb_berror.f77.gcv
OBERROR=${FIX_ROOT}/prepobs_errtable.global
if [ ${bk_core} = NMM ] ; then
ANAVINFO=${FIX_ROOT}/anavinfo_ndas_netcdf_glbe
fi
if [ ${bk_core} = ARW ] ; then
ANAVINFO=${FIX_ROOT}/anavinfo_arw_netcdf_glbe
fi
if [ ${bk_core} = NMMB ] ; then
ANAVINFO=${FIX_ROOT}/anavinfo_nems_nmmb_glb
fi
else
echo ' Use NAM background error covariance'
BERROR=${FIX_ROOT}/${BYTE_ORDER}/nam_nmmstat_na.gcv
OBERROR=${FIX_ROOT}/nam_errtable.r3dv
if [ ${bk_core} = NMM ] ; then
ANAVINFO=${FIX_ROOT}/anavinfo_ndas_netcdf
fi
if [ ${bk_core} = ARW ] ; then
ANAVINFO=${FIX_ROOT}/anavinfo_arw_netcdf
fi
if [ ${bk_core} = NMMB ] ; then
ANAVINFO=${FIX_ROOT}/anavinfo_nems_nmmb
fi
fi
SATANGL=${FIX_ROOT}/global_satangbias.txt
SATINFO=${FIX_ROOT}/global_satinfo.txt
CONVINFO=${FIX_ROOT}/global_convinfo.txt
OZINFO=${FIX_ROOT}/global_ozinfo.txt
PCPINFO=${FIX_ROOT}/global_pcpinfo.txt
LIGHTINFO=${FIX_ROOT}/global_lightinfo.txt
# copy Fixed fields to working directory
cp $ANAVINFO anavinfo
cp $BERROR berror_stats
cp $SATANGL satbias_angle
cp $SATINFO satinfo
cp $CONVINFO convinfo
cp $OZINFO ozinfo
cp $PCPINFO pcpinfo
cp $LIGHTINFO lightinfo
cp $OBERROR errtable
#
# # CRTM Spectral and Transmittance coefficients
CRTM_ROOT_ORDER=${CRTM_ROOT}/${BYTE_ORDER}
emiscoef_IRwater=${CRTM_ROOT_ORDER}/Nalli.IRwater.EmisCoeff.bin
emiscoef_IRice=${CRTM_ROOT_ORDER}/NPOESS.IRice.EmisCoeff.bin
emiscoef_IRland=${CRTM_ROOT_ORDER}/NPOESS.IRland.EmisCoeff.bin
emiscoef_IRsnow=${CRTM_ROOT_ORDER}/NPOESS.IRsnow.EmisCoeff.bin
emiscoef_VISice=${CRTM_ROOT_ORDER}/NPOESS.VISice.EmisCoeff.bin
emiscoef_VISland=${CRTM_ROOT_ORDER}/NPOESS.VISland.EmisCoeff.bin
emiscoef_VISsnow=${CRTM_ROOT_ORDER}/NPOESS.VISsnow.EmisCoeff.bin
emiscoef_VISwater=${CRTM_ROOT_ORDER}/NPOESS.VISwater.EmisCoeff.bin
emiscoef_MWwater=${CRTM_ROOT_ORDER}/FASTEM6.MWwater.EmisCoeff.bin
aercoef=${CRTM_ROOT_ORDER}/AerosolCoeff.bin
cldcoef=${CRTM_ROOT_ORDER}/CloudCoeff.bin
ln -s $emiscoef_IRwater ./Nalli.IRwater.EmisCoeff.bin
ln -s $emiscoef_IRice ./NPOESS.IRice.EmisCoeff.bin
ln -s $emiscoef_IRsnow ./NPOESS.IRsnow.EmisCoeff.bin
ln -s $emiscoef_IRland ./NPOESS.IRland.EmisCoeff.bin
ln -s $emiscoef_VISice ./NPOESS.VISice.EmisCoeff.bin
ln -s $emiscoef_VISland ./NPOESS.VISland.EmisCoeff.bin
ln -s $emiscoef_VISsnow ./NPOESS.VISsnow.EmisCoeff.bin
ln -s $emiscoef_VISwater ./NPOESS.VISwater.EmisCoeff.bin
ln -s $emiscoef_MWwater ./FASTEM6.MWwater.EmisCoeff.bin
ln -s $aercoef ./AerosolCoeff.bin
ln -s $cldcoef ./CloudCoeff.bin
# Copy CRTM coefficient files based on entries in satinfo file
for file in `awk '{if($1!~"!"){print $1}}' ./satinfo | sort | uniq` ;do
ln -s ${CRTM_ROOT_ORDER}/${file}.SpcCoeff.bin ./
ln -s ${CRTM_ROOT_ORDER}/${file}.TauCoeff.bin ./
done
# Only need this file for single obs test
bufrtable=${FIX_ROOT}/prepobs_prep.bufrtable
cp $bufrtable ./prepobs_prep.bufrtable
# for satellite bias correction
# Users may need to use their own satbias files for correct bias correction
cp ${GSI_ROOT}/fix/comgsi_satbias_in ./satbias_in
cp ${GSI_ROOT}/fix/comgsi_satbias_pc_in ./satbias_pc_in
#
##################################################################################
# Set some parameters for use by the GSI executable and to build the namelist
echo " Build the namelist "
# default is NAM
# as_op='1.0,1.0,0.5 ,0.7,0.7,0.5,1.0,1.0,'
vs_op='1.0,'
hzscl_op='0.373,0.746,1.50,'
if [ ${bkcv_option} = GLOBAL ] ; then
# as_op='0.6,0.6,0.75,0.75,0.75,0.75,1.0,1.0'
vs_op='0.7,'
hzscl_op='1.7,0.8,0.5,'
fi
if [ ${bk_core} = NMMB ] ; then
vs_op='0.6,'
fi
# default is NMM
bk_core_arw='.false.'
bk_core_nmm='.true.'
bk_core_nmmb='.false.'
bk_if_netcdf='.true.'
if [ ${bk_core} = ARW ] ; then
bk_core_arw='.true.'
bk_core_nmm='.false.'
bk_core_nmmb='.false.'
bk_if_netcdf='.true.'
fi
if [ ${bk_core} = NMMB ] ; then
bk_core_arw='.false.'
bk_core_nmm='.false.'
bk_core_nmmb='.true.'
bk_if_netcdf='.false.'
fi
if [ ${if_observer} = Yes ] ; then
nummiter=0
if_read_obs_save='.true.'
if_read_obs_skip='.false.'
else
nummiter=2
if_read_obs_save='.false.'
if_read_obs_skip='.false.'
fi
# Build the GSI namelist on-the-fly
. $GSI_NAMELIST
# modify the anavinfo vertical levels based on wrf_inout for WRF ARW and NMM
if [ ${bk_core} = ARW ] || [ ${bk_core} = NMM ] ; then
bklevels=`ncdump -h wrf_inout | grep "bottom_top =" | awk '{print $3}' `
bklevels_stag=`ncdump -h wrf_inout | grep "bottom_top_stag =" | awk '{print $3}' `
anavlevels=`cat anavinfo | grep ' sf ' | tail -1 | awk '{print $2}' ` # levels of sf, vp, u, v, t, etc
anavlevels_stag=`cat anavinfo | grep ' prse ' | tail -1 | awk '{print $2}' ` # levels of prse
sed -i 's/ '$anavlevels'/ '$bklevels'/g' anavinfo
sed -i 's/ '$anavlevels_stag'/ '$bklevels_stag'/g' anavinfo
fi
#
###################################################
# run GSI
###################################################
echo ' Run GSI with' ${bk_core} 'background'
case $ARCH in
'IBM_LSF')
${RUN_COMMAND} ./gsi.x < gsiparm.anl > stdout 2>&1 ;;
* )
${RUN_COMMAND} ./gsi.x > stdout 2>&1 ;;
esac
##################################################################
# run time error check
##################################################################
error=$?
if [ ${error} -ne 0 ]; then
echo "ERROR: ${GSI} crashed Exit status=${error}"
exit ${error}
fi
#
##################################################################
#
# GSI updating satbias_in
#
# GSI updating satbias_in (only for cycling assimilation)
# Copy the output to more understandable names
ln -s stdout stdout.anl.${ANAL_TIME}
ln -s wrf_inout wrfanl.${ANAL_TIME}
ln -s fort.201 fit_p1.${ANAL_TIME}
ln -s fort.202 fit_w1.${ANAL_TIME}
ln -s fort.203 fit_t1.${ANAL_TIME}
ln -s fort.204 fit_q1.${ANAL_TIME}
ln -s fort.207 fit_rad1.${ANAL_TIME}
# Loop over first and last outer loops to generate innovation
# diagnostic files for indicated observation types (groups)
#
# NOTE: Since we set miter=2 in GSI namelist SETUP, outer
# loop 03 will contain innovations with respect to
# the analysis. Creation of o-a innovation files
# is triggered by write_diag(3)=.true. The setting
# write_diag(1)=.true. turns on creation of o-g
# innovation files.
#
loops="01 03"
for loop in $loops; do
case $loop in
01) string=ges;;
03) string=anl;;
*) string=$loop;;
esac
# Collect diagnostic files for obs types (groups) below
# listall="conv amsua_metop-a mhs_metop-a hirs4_metop-a hirs2_n14 msu_n14 \
# sndr_g08 sndr_g10 sndr_g12 sndr_g08_prep sndr_g10_prep sndr_g12_prep \
# sndrd1_g08 sndrd2_g08 sndrd3_g08 sndrd4_g08 sndrd1_g10 sndrd2_g10 \
# sndrd3_g10 sndrd4_g10 sndrd1_g12 sndrd2_g12 sndrd3_g12 sndrd4_g12 \
# hirs3_n15 hirs3_n16 hirs3_n17 amsua_n15 amsua_n16 amsua_n17 \
# amsub_n15 amsub_n16 amsub_n17 hsb_aqua airs_aqua amsua_aqua \
# goes_img_g08 goes_img_g10 goes_img_g11 goes_img_g12 \
# pcp_ssmi_dmsp pcp_tmi_trmm sbuv2_n16 sbuv2_n17 sbuv2_n18 \
# omi_aura ssmi_f13 ssmi_f14 ssmi_f15 hirs4_n18 amsua_n18 mhs_n18 \
# amsre_low_aqua amsre_mid_aqua amsre_hig_aqua ssmis_las_f16 \
# ssmis_uas_f16 ssmis_img_f16 ssmis_env_f16 mhs_metop_b \
# hirs4_metop_b hirs4_n19 amusa_n19 mhs_n19 goes_glm_16"
listall=`ls pe* | cut -f2 -d"." | awk '{print substr($0, 0, length($0)-3)}' | sort | uniq `
for type in $listall; do
count=`ls pe*${type}_${loop}* | wc -l`
if [[ $count -gt 0 ]]; then
cat pe*${type}_${loop}* > diag_${type}_${string}.${ANAL_TIME}
fi
done
done
# Clean working directory to save only important files
ls -l * > list_run_directory
if [[ ${if_clean} = clean && ${if_observer} != Yes ]]; then
echo ' Clean working directory after GSI run'
rm -f *Coeff.bin # all CRTM coefficient files
rm -f pe0* # diag files on each processor
rm -f obs_input.* # observation middle files
rm -f siganl sigf0? # background middle files
rm -f fsize_* # delete temperal file for bufr size
fi
#
#
#################################################
# start to calculate diag files for each member
#################################################
#
if [ ${if_observer} = Yes ] ; then
string=ges
for type in $listall; do
count=0
if [[ -f diag_${type}_${string}.${ANAL_TIME} ]]; then
mv diag_${type}_${string}.${ANAL_TIME} diag_${type}_${string}.ensmean
fi
done
mv wrf_inout wrf_inout_ensmean
# Build the GSI namelist on-the-fly for each member
nummiter=0
if_read_obs_save='.false.'
if_read_obs_skip='.true.'
. $GSI_NAMELIST
# Loop through each member
loop="01"
ensmem=1
while [[ $ensmem -le $no_member ]];do
rm pe0*
print "\$ensmem is $ensmem"
ensmemid=`printf %3.3i $ensmem`
# get new background for each member
if [[ -f wrf_inout ]]; then
rm wrf_inout
fi
BK_FILE=${BK_FILE_mem}${ensmemid}
echo $BK_FILE
ln -s $BK_FILE wrf_inout
# run GSI
echo ' Run GSI with' ${bk_core} 'for member ', ${ensmemid}
case $ARCH in
'IBM_LSF')
${RUN_COMMAND} ./gsi.x < gsiparm.anl > stdout_mem${ensmemid} 2>&1 ;;
* )
${RUN_COMMAND} ./gsi.x > stdout_mem${ensmemid} 2>&1 ;;
esac
# run time error check and save run time file status
error=$?
if [ ${error} -ne 0 ]; then
echo "ERROR: ${GSI} crashed for member ${ensmemid} Exit status=${error}"
exit ${error}
fi
ls -l * > list_run_directory_mem${ensmemid}
# generate diag files
for type in $listall; do
count=`ls pe*${type}_${loop}* | wc -l`
if [[ $count -gt 0 ]]; then
cat pe*${type}_${loop}* > diag_${type}_${string}.mem${ensmemid}
fi
done
# next member
(( ensmem += 1 ))
done
fi
exit 0
wrf.exe
模式运行，生成 wrfout 文件