Gen_be problem

Issues with running 3DVAR.

Gen_be problem

Postby soleiman » Tue Jun 01, 2010 1:54 pm


I'm new WRF user and now i want to assimilated GPSPW data into WRF
i prepared ascii file with GPSPW file and convrted to little_r format
the i made a link to wrfout_d01_* and also link to ob_ascii
but i could not create my own be.dat file
for this purpose i used gen_be.ksh and modified these fields :


Code: Select all
export REL_DIR=${REL_DIR:-$HOME/WRF}
export WRFVAR_DIR=${WRFVAR_DIR:-$REL_DIR/WRFDA_V3.1.1}
export SCRIPTS_DIR=${SCRIPTS_DIR:-$WRFVAR_DIR/var/scripts}
export REG_DIR=${REL_DIR}DATA/DOMAINS/florida
export EXP_DIR=${REL_DIR}/DATA/DOMAINS/florida/OI
export FC_DIR=$EXP_DIR/FC
export START_DATE=2010050900
export END_DATE=2010051000
export INTERVAL=24
export NUM_LEVELS=35
export STRIDE=1
export N_SMTH_SL=0


After running gen_be.ksh i have got this error :
In gen_be_diags.log
Code: Select all
At line 2276 of file da_gen_be.f (unit = 23, file = 'gen_be_stage3.chi_u.dat')
Fortran runtime error: End of file

 [1] Gather regression coefficients.

 [2] Gather vertical error eigenvectors, eigenvalues.



Also in In gen_be_stage3.t_u.log, gen_be_stage3.chi_u.log
Code: Select all
 ---------------------------- FATAL ERROR -----------------------
 Fatal error in file:  da_eof_decomposition.inc  LINE:      40
 ----------------------------------------------------------------
 [1] Initialize namelist variables and other scalars.
 Computing vertical error statistics for dates 2010050900 to 2010051000
 Interval between dates =       24hours.
 Number of ensemble members at each time =        1
 [2] Process fields for variable t_u
    Processing data for date 2010050900, variable t_u and member    1
    i, j, k dimensions are      149     169      34
    Processing data for date 2010051000, variable t_u and member    1
 [3] Calculate eigenvectors and eigenvalues for variable t_u
 Calculate eigenvectors and eigenvalues for bin      1 of      1
 Error in decomposition, info =   33


here also the gen_be.ksh log
Code: Select all
+ REL_DIR=/home/wrf/WRF
+ export REL_DIR
+ WRFVAR_DIR=/home/wrf/WRF/WRFDA_V3.1.1
+ export WRFVAR_DIR
+ SCRIPTS_DIR=/home/wrf/WRF/WRFDA_V3.1.1/var/scripts
+ export SCRIPTS_DIR
+ REG_DIR=/home/wrf/WRFDATA/DOMAINS/florida
+ export REG_DIR
+ EXP_DIR=/home/wrf/WRF/DATA/DOMAINS/florida/OI
+ export EXP_DIR
+ FC_DIR=/home/wrf/WRF/DATA/DOMAINS/florida/OI/FC
+ export FC_DIR
+ START_DATE=2010050900
+ export START_DATE
+ END_DATE=2010051000
+ export END_DATE
+ INTERVAL=24
+ export INTERVAL
+ NUM_LEVELS=35
+ export NUM_LEVELS
+ STRIDE=1
+ export STRIDE
+ N_SMTH_SL=0
+ export N_SMTH_SL
+ . /home/wrf/WRF/DATA/DOMAINS/florida/OI/gen_be_set_defaults.ksh
+ set -x
+ RUN_GEN_BE_STAGE0=true
+ export RUN_GEN_BE_STAGE0
+ RUN_GEN_BE_STAGE1=true
+ export RUN_GEN_BE_STAGE1
+ RUN_GEN_BE_STAGE2=true
+ export RUN_GEN_BE_STAGE2
+ RUN_GEN_BE_STAGE2A=true
+ export RUN_GEN_BE_STAGE2A
+ RUN_GEN_BE_STAGE3=true
+ export RUN_GEN_BE_STAGE3
+ RUN_GEN_BE_STAGE4=true
+ export RUN_GEN_BE_STAGE4
+ RUN_GEN_BE_DIAGS=true
+ export RUN_GEN_BE_DIAGS
+ RUN_GEN_BE_DIAGS_READ=true
+ export RUN_GEN_BE_DIAGS_READ
+ RUN_GEN_BE_MULTICOV=true
+ export RUN_GEN_BE_MULTICOV
+ DOMAIN=01
+ export DOMAIN
+ START_DATE=2010050900
+ export START_DATE
+ END_DATE=2010051000
+ export END_DATE
+ FCST_RANGE=24
+ export FCST_RANGE
+ FCST_RANGE1=24
+ export FCST_RANGE1
+ FCST_RANGE2=12
+ export FCST_RANGE2
+ INTERVAL=24
+ export INTERVAL
+ BE_METHOD=NMC
+ export BE_METHOD
+ NE=1
+ export NE
+ BIN_TYPE=5
+ export BIN_TYPE
+ LAT_MIN=19.0
+ export LAT_MIN
+ LAT_MAX=32.0
+ export LAT_MAX
+ BINWIDTH_LAT=25.0
+ export BINWIDTH_LAT
+ BINWIDTH_HGT=1000.0
+ export BINWIDTH_HGT
+ HGT_MIN=0.0
+ export HGT_MIN
+ HGT_MAX=20000.0
+ export HGT_MAX
+ REMOVE_MEAN=.true.
+ export REMOVE_MEAN
+ GAUSSIAN_LATS=.false.
+ export GAUSSIAN_LATS
+ TESTING_EOFS=.true.
+ export TESTING_EOFS
+ NUM_PASSES=0
+ export NUM_PASSES
+ RF_SCALE=1.0
+ export RF_SCALE
+ USE_GLOBAL_EOFS=.true.
+ export USE_GLOBAL_EOFS
+ DATA_ON_LEVELS=.false.
+ export DATA_ON_LEVELS
+ GLOBAL=false
+ export GLOBAL
+ NUM_LEVELS=35
+ export NUM_LEVELS
+ N_SMTH_SL=0
+ export N_SMTH_SL
+ STRIDE=1
+ export STRIDE
+ NBINS=1
+ export NBINS
+ IBIN=1
+ export IBIN
+ TESTING_SPECTRAL=.false.
+ export TESTING_SPECTRAL
+ LOCAL=true
+ export LOCAL
+ NUM_JOBS=1
+ export NUM_JOBS
+ MACHINES=' node1 node1 node2 node2 node3 node3 node4 node4                              node5 node5 node6 node6 node7 node7 node8 node8'
+ export MACHINES
+ REGION=con200
+ export REGION
+ EXPT=noobs
+ export EXPT
+ ID=gen_be
+ export ID
+ ID1=NMC.bin_type5
+ export ID1
+ VARIABLE1=chi_u
+ export VARIABLE1
+ VARIABLE2=chi
+ export VARIABLE2
+ CLEAN=false
+ export CLEAN
+ FILE_TYPE=wrfout
+ export FILE_TYPE
+ REL_DIR=/home/wrf/WRF
+ export REL_DIR
+ WRFVAR_DIR=/home/wrf/WRF/WRFDA_V3.1.1
+ export WRFVAR_DIR
+ BUILD_DIR=/home/wrf/WRF/WRFDA_V3.1.1/var/da
+ export BUILD_DIR
+ DAT_DIR=/home/wrf/WRF/VAR_DATA
+ export DAT_DIR
+ REG_DIR=/home/wrf/WRFDATA/DOMAINS/florida
+ export REG_DIR
+ EXP_DIR=/home/wrf/WRF/DATA/DOMAINS/florida/OI
+ export EXP_DIR
+ FC_DIR=/home/wrf/WRF/DATA/DOMAINS/florida/OI/FC
+ export FC_DIR
+ RUN_DIR=/home/wrf/WRF/DATA/DOMAINS/florida/OI/gen_be5
+ export RUN_DIR
+ WORK_DIR=/home/wrf/WRF/DATA/DOMAINS/florida/OI/gen_be5/working
+ export WORK_DIR
+ STAGE0_DIR=/home/wrf/WRF/DATA/DOMAINS/florida/OI/gen_be5/working/stage0
+ export STAGE0_DIR
+ false
+ UH_METHOD=scale
+ export UH_METHOD
+ CONTROL_VARIABLES=' psi chi_u t_u rh ps_u '
+ export CONTROL_VARIABLES
+ DELETE_DIRS=' '
+ export DELETE_DIRS
+ [[ ! -d /home/wrf/WRF/DATA/DOMAINS/florida/OI/gen_be5 ]]
+ mkdir -p /home/wrf/WRF/DATA/DOMAINS/florida/OI/gen_be5
+ [[ ! -d /home/wrf/WRF/DATA/DOMAINS/florida/OI/gen_be5/working/stage0 ]]
+ mkdir -p /home/wrf/WRF/DATA/DOMAINS/florida/OI/gen_be5/working/stage0
+ mkdir -p /home/wrf/WRF/DATA/DOMAINS/florida/OI/gen_be5/working
+ cd /home/wrf/WRF/DATA/DOMAINS/florida/OI/gen_be5/working
+ mkdir -p fullflds
+ mkdir -p psi
+ mkdir -p chi
+ mkdir -p t
+ mkdir -p rh
+ mkdir -p ps
+ mkdir -p psi
+ mkdir -p chi_u
+ mkdir -p t_u
+ mkdir -p rh
+ mkdir -p ps_u
+ echo

+ date
+ echo Tue Jun 1 19:38:41 MSD 2010 Start
Tue Jun 1 19:38:41 MSD 2010 Start
+ svnversion /home/wrf/WRF/WRFDA_V3.1.1
+ echo 'WRFVAR_DIR is' /home/wrf/WRF/WRFDA_V3.1.1 exported
WRFVAR_DIR is /home/wrf/WRF/WRFDA_V3.1.1 exported
+ true
+ echo ---------------------------------------------------------------
---------------------------------------------------------------
+ echo 'Run Stage 0: Calculate ensemble perturbations from model forecasts.'
Run Stage 0: Calculate ensemble perturbations from model forecasts.
+ echo ---------------------------------------------------------------
---------------------------------------------------------------
+ date
+ BEGIN_CPU='Tue Jun  1 19:38:41 MSD 2010'
+ export BEGIN_CPU
+ echo 'Beginning CPU time: Tue Jun  1 19:38:41 MSD 2010'
Beginning CPU time: Tue Jun  1 19:38:41 MSD 2010
+ /home/wrf/WRF/WRFDA_V3.1.1/var/scripts/gen_be/gen_be_stage0_wrf.ksh
gen_be_stage0_wrf: Calculating standard perturbation fields valid at time  2010050900
2010050800 /home/wrf/WRF/DATA/DOMAINS/florida/OI/FC/2010050800/wrfout_d01_2010-05-09_00:00:00 /home/wrf/WRF/DATA/DOMAINS/florida/OI/FC/2010050900/wrfout_d01_2010-05-09_00:00:00
gen_be_stage0_wrf: Calculating standard perturbation fields valid at time  2010051000
2010050900 /home/wrf/WRF/DATA/DOMAINS/florida/OI/FC/2010050900/wrfout_d01_2010-05-10_00:00:00 /home/wrf/WRF/DATA/DOMAINS/florida/OI/FC/2010051000/wrfout_d01_2010-05-10_00:00:00
+ RC=0
+ [[ 0 != 0 ]]
+ date
+ END_CPU='Tue Jun  1 19:39:02 MSD 2010'
+ export END_CPU
+ echo 'Ending CPU time: Tue Jun  1 19:39:02 MSD 2010'
Ending CPU time: Tue Jun  1 19:39:02 MSD 2010
+ true
+ echo ---------------------------------------------------------------
---------------------------------------------------------------
+ echo 'Run Stage 1: Read standard' 'fields, and remove time/ensemble/area mean.'
Run Stage 1: Read standard fields, and remove time/ensemble/area mean.
+ echo ---------------------------------------------------------------
---------------------------------------------------------------
+ date
+ BEGIN_CPU='Tue Jun  1 19:39:02 MSD 2010'
+ export BEGIN_CPU
+ echo 'Beginning CPU time: Tue Jun  1 19:39:02 MSD 2010'
Beginning CPU time: Tue Jun  1 19:39:02 MSD 2010
+ ln -sf /home/wrf/WRF/WRFDA_V3.1.1/var/da/gen_be_stage1.exe .
+ cat
+ 1> gen_be_stage1_nl.nl 0<< \EOF
&gen_be_stage1_nl
    start_date = '2010050900',
    end_date = '2010051000',
    interval = 24,
    be_method = 'NMC',
    ne = 1,
    bin_type = 5,
    lat_min = 19.0,
    lat_max = 32.0,
    binwidth_lat = 25.0,
    hgt_min = 0.0,
    hgt_max = 20000.0,
    binwidth_hgt = 1000.0,
    dat_dir = '/home/wrf/WRF/DATA/DOMAINS/florida/OI/gen_be5/working/stage0' /
EOF
+ ./gen_be_stage1.exe
+ 1> gen_be_stage1.log 2>& 1
+ RC=0
+ [[ 0 != 0 ]]
+ date
+ END_CPU='Tue Jun  1 19:39:04 MSD 2010'
+ export END_CPU
+ echo 'Ending CPU time: Tue Jun  1 19:39:04 MSD 2010'
Ending CPU time: Tue Jun  1 19:39:04 MSD 2010
+ true
+ echo ---------------------------------------------------------------
---------------------------------------------------------------
+ echo 'Run Stage 2: Calculate regression coefficients.'
Run Stage 2: Calculate regression coefficients.
+ echo ---------------------------------------------------------------
---------------------------------------------------------------
+ date
+ BEGIN_CPU='Tue Jun  1 19:39:04 MSD 2010'
+ export BEGIN_CPU
+ echo 'Beginning CPU time: Tue Jun  1 19:39:04 MSD 2010'
Beginning CPU time: Tue Jun  1 19:39:04 MSD 2010
+ ln -sf /home/wrf/WRF/WRFDA_V3.1.1/var/da/gen_be_stage2.exe .
+ cat
+ 1> gen_be_stage2_nl.nl 0<< \EOF
&gen_be_stage2_nl
    start_date = '2010050900',
    end_date = '2010051000',
    interval = 24,
    ne = 1,
    testing_eofs = .true. /
EOF
+ ./gen_be_stage2.exe
+ 1> gen_be_stage2.log 2>& 1
+ RC=0
+ [[ 0 != 0 ]]
+ date
+ END_CPU='Tue Jun  1 19:39:05 MSD 2010'
+ export END_CPU
+ echo 'Ending CPU time: Tue Jun  1 19:39:05 MSD 2010'
Ending CPU time: Tue Jun  1 19:39:05 MSD 2010
+ true
+ echo ---------------------------------------------------------------
---------------------------------------------------------------
+ echo 'Run Stage 2a: Calculate control variable fields.'
Run Stage 2a: Calculate control variable fields.
+ echo ---------------------------------------------------------------
---------------------------------------------------------------
+ date
+ BEGIN_CPU='Tue Jun  1 19:39:05 MSD 2010'
+ export BEGIN_CPU
+ echo 'Beginning CPU time: Tue Jun  1 19:39:05 MSD 2010'
Beginning CPU time: Tue Jun  1 19:39:05 MSD 2010
+ ln -sf /home/wrf/WRF/WRFDA_V3.1.1/var/da/gen_be_stage2a.exe .
+ cat
+ 1> gen_be_stage2a_nl.nl 0<< \EOF
&gen_be_stage2a_nl
    start_date = '2010050900',
    end_date = '2010051000',
    interval = 24,
    ne = 1,
    num_passes = 0,
    rf_scale = 1.0 /
EOF
+ ./gen_be_stage2a.exe
+ 1> gen_be_stage2a.log 2>& 1
+ RC=0
+ [[ 0 != 0 ]]
+ rm -rf
+ 2> /dev/null
+ date
+ END_CPU='Tue Jun  1 19:39:06 MSD 2010'
+ export END_CPU
+ echo 'Ending CPU time: Tue Jun  1 19:39:06 MSD 2010'
Ending CPU time: Tue Jun  1 19:39:06 MSD 2010
+ true
+ echo ---------------------------------------------------------------
---------------------------------------------------------------
+ echo 'Run Stage 3: Read 3D control variable fields, and calculate vertical covariances.'
Run Stage 3: Read 3D control variable fields, and calculate vertical covariances.
+ echo ---------------------------------------------------------------
---------------------------------------------------------------
+ date
+ BEGIN_CPU='Tue Jun  1 19:39:06 MSD 2010'
+ export BEGIN_CPU
+ echo 'Beginning CPU time: Tue Jun  1 19:39:06 MSD 2010'
Beginning CPU time: Tue Jun  1 19:39:06 MSD 2010
+ ln -sf /home/wrf/WRF/WRFDA_V3.1.1/var/da/gen_be_stage3.exe .
+ cat
+ 1> gen_be_stage3_nl.nl 0<< \EOF
&gen_be_stage3_nl
    start_date = '2010050900',
    end_date = '2010051000',
    interval = 24,
    variable = 'psi',
    ne = 1,
    bin_type = 5,
    lat_min = 19.0,
    lat_max = 32.0,
    binwidth_lat = 25.0,
    hgt_min = 0.0,
    hgt_max = 20000.0,
    binwidth_hgt = 1000.0,
    testing_eofs = .true.,
    use_global_eofs = .true.,
    data_on_levels = .false. /
EOF
+ ./gen_be_stage3.exe
+ 1> gen_be_stage3.psi.log 2>& 1
+ RC=0
+ [[ 0 != 0 ]]
+ cat
+ 1> gen_be_stage3_nl.nl 0<< \EOF
&gen_be_stage3_nl
    start_date = '2010050900',
    end_date = '2010051000',
    interval = 24,
    variable = 'chi_u',
    ne = 1,
    bin_type = 5,
     lat_min = 19.0,
    lat_max = 32.0,
    binwidth_lat = 25.0,
    hgt_min = 0.0,
    hgt_max = 20000.0,
    binwidth_hgt = 1000.0,
    testing_eofs = .true.,
    use_global_eofs = .true.,
    data_on_levels = .false. /
EOF
+ ./gen_be_stage3.exe
+ 1> gen_be_stage3.chi_u.log 2>& 1
+ RC=0
+ [[ 0 != 0 ]]
+ cat
+ 1> gen_be_stage3_nl.nl 0<< \EOF
&gen_be_stage3_nl
    start_date = '2010050900',
    end_date = '2010051000',
    interval = 24,
    variable = 't_u',
    ne = 1,
    bin_type = 5,
    lat_min = 19.0,
    lat_max = 32.0,
    binwidth_lat = 25.0,
    hgt_min = 0.0,
    hgt_max = 20000.0,
    binwidth_hgt = 1000.0,
    testing_eofs = .true.,
    use_global_eofs = .true.,
    data_on_levels = .false. /
EOF
+ ./gen_be_stage3.exe
+ 1> gen_be_stage3.t_u.log 2>& 1
+ RC=0
+ [[ 0 != 0 ]]
+ cat
+ 1> gen_be_stage3_nl.nl 0<< \EOF
&gen_be_stage3_nl
    start_date = '2010050900',
    end_date = '2010051000',
    interval = 24,
    variable = 'rh',
    ne = 1,
    bin_type = 5,
    lat_min = 19.0,
    lat_max = 32.0,
    binwidth_lat = 25.0,
    hgt_min = 0.0,
    hgt_max = 20000.0,
    binwidth_hgt = 1000.0,
    testing_eofs = .true.,
    use_global_eofs = .true.,
    data_on_levels = .false. /
EOF
+ ./gen_be_stage3.exe
+ 1> gen_be_stage3.rh.log 2>& 1
+ RC=0
+ [[ 0 != 0 ]]
+ cat
+ 1> gen_be_stage3_nl.nl 0<< \EOF
&gen_be_stage3_nl
    start_date = '2010050900',
    end_date = '2010051000',
    interval = 24,
    variable = 'ps_u',
    ne = 1,
    bin_type = 5,
    lat_min = 19.0,
    lat_max = 32.0,
    binwidth_lat = 25.0,
    hgt_min = 0.0,
    hgt_max = 20000.0,
    binwidth_hgt = 1000.0,
    testing_eofs = .true.,
    use_global_eofs = .true.,
    data_on_levels = .false. /
EOF
+ ./gen_be_stage3.exe
+ 1> gen_be_stage3.ps_u.log 2>& 1
+ RC=0
+ [[ 0 != 0 ]]
+ date
+ END_CPU='Tue Jun  1 19:39:10 MSD 2010'
+ export END_CPU
+ echo 'Ending CPU time: Tue Jun  1 19:39:10 MSD 2010'
Ending CPU time: Tue Jun  1 19:39:10 MSD 2010
+ true
+ date
+ BEGIN_CPU='Tue Jun  1 19:39:10 MSD 2010'
+ export BEGIN_CPU
+ echo 'Beginning CPU time: Tue Jun  1 19:39:10 MSD 2010'
Beginning CPU time: Tue Jun  1 19:39:10 MSD 2010
+ false
+ echo ---------------------------------------------------------------
---------------------------------------------------------------
+ echo 'Run Stage 4: Calculate horizontal covariances (regional lengthscales).'
Run Stage 4: Calculate horizontal covariances (regional lengthscales).
+ echo ---------------------------------------------------------------
---------------------------------------------------------------
+ /home/wrf/WRF/WRFDA_V3.1.1/var/scripts/gen_be/gen_be_stage4_regional.ksh
+ 1> gen_be_stage4_regional.log 2>& 1
+ RC=0
+ [[ 0 != 0 ]]
+ date
+ END_CPU='Tue Jun  1 19:41:49 MSD 2010'
+ export END_CPU
+ echo 'Ending CPU time: Tue Jun  1 19:41:49 MSD 2010'
Ending CPU time: Tue Jun  1 19:41:49 MSD 2010
+ true
+ ln -sf /home/wrf/WRF/WRFDA_V3.1.1/var/da/gen_be_diags.exe .
+ cat
+ 1> gen_be_diags_nl.nl 0<< \EOF
&gen_be_diags_nl
   uh_method = 'scale',
   n_smth_sl = 0, /
EOF
+ ./gen_be_diags.exe
+ 1> gen_be_diags.log 2>& 1
+ RC=2
+ [[ 2 != 0 ]]
+ echo 'Stage gen_be_diags failed with error' 2
Stage gen_be_diags failed with error 2
+ echo gen_be_diags
+ 1> /home/wrf/WRF/DATA/DOMAINS/florida/OI/gen_be5/FAIL
+ exit 1
               



Please help us
Thanks in advance for your help
soleiman
 
Posts: 12
Joined: Sat Oct 25, 2008 11:02 am

Return to Runtime Problems

Who is online

Users browsing this forum: No registered users and 2 guests