Error with 4DVar single observation experiment

Setting up the 3DVAR package within WRF.

Error with 4DVar single observation experiment

Postby deejeek » Thu Feb 22, 2018 5:42 am

Hi all,

I am trying to run single observation experiment using WRFDA-4DVar v.3.8.1. Attached herewith the namelist used for the test and the rsl.out file. The model is not taking the pseudo observation and the run is getting finished without any assimilation. However, similar run was successful with 3DVar. It would be a great help, if anyone can give some clarification in this regard.

Thanks,
Deepak.

namelist.input

&wrfvar1
var4d=true,
var4d_lbc=false,
var4d_bin=3600,
var4d_bin_rain=21600,
print_detail_outerloop=true,
print_detail_grad=true,
/
&wrfvar2
/
&wrfvar3
ob_format=2,
/
&wrfvar4
use_synopobs=true,
use_shipsobs=true,
use_metarobs=true,
use_soundobs=true,
use_pilotobs=true,
use_airepobs=true,
use_geoamvobs=true,
use_polaramvobs=true,
use_bogusobs=true,
use_buoyobs=true,
use_profilerobs=true,
use_satemobs=true,
use_gpspwobs=false,
use_gpsrefobs=false,
use_qscatobs=true,
use_rainobs=false,
/
&wrfvar5
check_max_iv=true,
/
&wrfvar6
max_ext_its=1,
ntmax=200,
orthonorm_gradient=true,
/
&wrfvar7
cv_options=5,
/
&wrfvar8
/
&wrfvar9
/
&wrfvar10
test_transforms=false,
test_gradient=false,
/
&wrfvar11
cv_options_hum=1,
check_rh=0,
sfc_assi_options=1,
calculate_cg_cost_fn=true,
/
&wrfvar12
/
&wrfvar13
/
&wrfvar14
/
&wrfvar15
num_pseudo = 1,
pseudo_x = 23.0,
pseudo_y = 23.0,
pseudo_z = 14.0,
pseudo_err = 1.0,
pseudo_val = 1.0,
/
&wrfvar16
/
&wrfvar17
/
&wrfvar18
analysis_date=2013-06-16_00:00:00
/
&wrfvar19
pseudo_var = “u”,
/
&wrfvar20
/
&wrfvar21
time_window_min=2013-06-16_00:00:00
/
&wrfvar22
time_window_max=2013-06-16_06:00:00
/
&time_control
run_hours=6,
start_year=2013
start_month=06
start_day=16
start_hour=00
end_year=2013
end_month=06
end_day=16
end_hour=06
interval_seconds=21600,
debug_level=0,
/
&domains
time_step=120,
e_we=375,
e_sn=250,
e_vert=36,
p_top_requested=5000,
interp_type=1,
t_extrap_type=1,
dx=27000,
dy=27000,
smooth_option=0,
nproc_x=0,
hypsometric_opt=2,
/
&fdda
/
&dfi_control
/
&tc
/
&physics
mp_physics=4,
ra_lw_physics=1,
ra_sw_physics=1,
radt=30,
sf_sfclay_physics=1,
sf_surface_physics=2,
bl_pbl_physics=1,
cu_physics=1,
cudt=5,
num_soil_layers=4,
mp_zero_out=2,
co2tf=0,
/
&scm
/
&dynamics
w_damping=1,
diff_opt=1,
km_opt=4,
dampcoef=0.01,
time_step_sound=4,
base_temp=290.0,
use_baseparam_fr_nml=true,
/
&bdy_control
specified=true,
real_data_init_type=3,
/
&grib2
/
&namelist_quilt
/
&perturbation
trajectory_io=true,
enable_identity=false,
jcdfi_use=false,
jcdfi_diag=1,
jcdfi_penalty=1000.0,
/

rsl.out.0000

Namelist logging not found in namelist.input. Using registry defaults for varia
bles in logging.

*** VARIATIONAL ANALYSIS ***
WRFDA V3.8.1

Ntasks in X 8, ntasks in Y 10
*************************************
Parent domain
ids,ide,jds,jde 1 375 1 250
ims,ime,jms,jme -4 54 -4 32
ips,ipe,jps,jpe 1 47 1 25
*************************************
DYNAMICS OPTION: Eulerian Mass Coordinate
alloc_space_field: domain 1, 132558304 bytes allocated
Timing for processing lateral boundary for domain 1: 0.92859 elapsed seconds
num_fgat_time is: 7

Tile Strategy is not specified. Assuming 1D-Y
WRF TILE 1 IS 1 IE 47 JS 1 JE 25
WRF NUMBER OF TILES = 1
Domain mapping info:
map_proj = 1
cen_lat = 0.300000E+02
cen_lon = 0.750000E+02
truelat1 = 0.500000E+01
truelat2 = 0.250000E+02
start_lat = -0.542975E+01
start_lon = 0.320338E+02
pole_lat = 0.900000E+02
dsm = 0.270000E+02

Set up observations (ob)

Single OBS Test:: Turn off all the OBS switches ***

Using ASCII format observation input

--------------------------- WARNING ---------------------------
WARNING FROM FILE: <A HREF="http://www.mmm.ucar.edu/people/wrfhelp/wrfvar/code/trunk/da_scan_obs_ascii.html">da_scan_obs_ascii.inc</a> LINE: 74
Error -1reading gts obs header on unit 25
---------------------------------------------------------------
--------------------------- WARNING ---------------------------
WARNING FROM FILE: <A HREF="http://www.mmm.ucar.edu/people/wrfhelp/wrfvar/code/trunk/da_scan_obs_ascii.html">da_scan_obs_ascii.inc</a> LINE: 74
Error -1reading gts obs header on unit 26
---------------------------------------------------------------
--------------------------- WARNING ---------------------------
WARNING FROM FILE: <A HREF="http://www.mmm.ucar.edu/people/wrfhelp/wrfvar/code/trunk/da_scan_obs_ascii.html">da_scan_obs_ascii.inc</a> LINE: 61
Error 29 opening gts obs file ob03.ascii
---------------------------------------------------------------
--------------------------- WARNING ---------------------------
WARNING FROM FILE: <A HREF="http://www.mmm.ucar.edu/people/wrfhelp/wrfvar/code/trunk/da_scan_obs_ascii.html">da_scan_obs_ascii.inc</a> LINE: 61
Error 29 opening gts obs file ob04.ascii
---------------------------------------------------------------
--------------------------- WARNING ---------------------------
WARNING FROM FILE: <A HREF="http://www.mmm.ucar.edu/people/wrfhelp/wrfvar/code/trunk/da_scan_obs_ascii.html">da_scan_obs_ascii.inc</a> LINE: 61
Error 29 opening gts obs file ob05.ascii
---------------------------------------------------------------
--------------------------- WARNING ---------------------------
WARNING FROM FILE: <A HREF="http://www.mmm.ucar.edu/people/wrfhelp/wrfvar/code/trunk/da_scan_obs_ascii.html">da_scan_obs_ascii.inc</a> LINE: 61
Error 29 opening gts obs file ob06.ascii
---------------------------------------------------------------
--------------------------- WARNING ---------------------------
WARNING FROM FILE: <A HREF="http://www.mmm.ucar.edu/people/wrfhelp/wrfvar/code/trunk/da_scan_obs_ascii.html">da_scan_obs_ascii.inc</a> LINE: 61
Error 29 opening gts obs file ob07.ascii
---------------------------------------------------------------
--------------------------- WARNING ---------------------------
WARNING FROM FILE: <A HREF="http://www.mmm.ucar.edu/people/wrfhelp/wrfvar/code/trunk/da_read_obs_ascii.html">da_read_obs_ascii.inc</a> LINE: 112
Problem reading gts obs header, skipping file
---------------------------------------------------------------
--------------------------- WARNING ---------------------------
WARNING FROM FILE: <A HREF="http://www.mmm.ucar.edu/people/wrfhelp/wrfvar/code/trunk/da_read_obs_ascii.html">da_read_obs_ascii.inc</a> LINE: 112
Problem reading gts obs header, skipping file
---------------------------------------------------------------
--------------------------- WARNING ---------------------------
WARNING FROM FILE: <A HREF="http://www.mmm.ucar.edu/people/wrfhelp/wrfvar/code/trunk/da_read_obs_ascii.html">da_read_obs_ascii.inc</a> LINE: 100
Error 29 opening gts obs file ob03.ascii
---------------------------------------------------------------
--------------------------- WARNING ---------------------------
WARNING FROM FILE: <A HREF="http://www.mmm.ucar.edu/people/wrfhelp/wrfvar/code/trunk/da_read_obs_ascii.html">da_read_obs_ascii.inc</a> LINE: 100
Error 29 opening gts obs file ob04.ascii
---------------------------------------------------------------
--------------------------- WARNING ---------------------------
WARNING FROM FILE: <A HREF="http://www.mmm.ucar.edu/people/wrfhelp/wrfvar/code/trunk/da_read_obs_ascii.html">da_read_obs_ascii.inc</a> LINE: 100
Error 29 opening gts obs file ob05.ascii
---------------------------------------------------------------
--------------------------- WARNING ---------------------------
WARNING FROM FILE: <A HREF="http://www.mmm.ucar.edu/people/wrfhelp/wrfvar/code/trunk/da_read_obs_ascii.html">da_read_obs_ascii.inc</a> LINE: 100
Error 29 opening gts obs file ob06.ascii
---------------------------------------------------------------
--------------------------- WARNING ---------------------------
WARNING FROM FILE: <A HREF="http://www.mmm.ucar.edu/people/wrfhelp/wrfvar/code/trunk/da_read_obs_ascii.html">da_read_obs_ascii.inc</a> LINE: 100
Error 29 opening gts obs file ob07.ascii
---------------------------------------------------------------
Observation summary
ob time 1
pseudo 1 global, 1 local
ob time 2
ob time 3
ob time 4
ob time 5
ob time 6
ob time 7

Set up background errors for regional application for cv_options = 5
:
:
:
wrf: back from adjoint integrate

Diagnostics
Final cost function J = 0.00

Total number of obs. = 0
Final value of J = 0.00000
Final value of Jo = 0.00000
Final value of Jd = 0.00000
Final value of Jb = 0.00000
Final value of Jc = 0.00000
Final value of Je = 0.00000
Final value of Jp = 0.00000
Final value of Jl = 0.00000
Jb factor used(1) = 1.00000
Jb factor used(2) = 1.00000
Jb factor used(3) = 1.00000
Jb factor used(4) = 1.00000
Jb factor used(5) = 1.00000
Jb factor used = 1.00000
Je factor used = 1.00000
VarBC factor used = 1.00000

d01 2013-06-16_00:00:00 xtraj linked list is initialized
d01 2013-06-16_00:00:00 ad_forcing linked list is initialized
d01 2013-06-16_00:00:00 tl_pert linked list is initialized
d01 2013-06-16_00:00:00 wrf: SUCCESS COMPLETE WRF
*** WRF-Var completed successfully ***
deejeek
 
Posts: 2
Joined: Sat Sep 13, 2014 2:11 pm

Return to Installation

Who is online

Users browsing this forum: No registered users and 1 guest