WRFDA Radiation Data Assimilation with AMSUA Satellite

Any issues with the actual running of the WRF.

WRFDA Radiation Data Assimilation with AMSUA Satellite

Postby Sabrina93 » Thu Jun 24, 2021 2:16 am

Dear All,
I am try WRFDA-3DVAR to assimilate radiation with AMSUA satellite data. The DA run is completed successfully but with two warnings to open 2 gts obs file ob.ascii files. Also I could not get the all the time steps. Only one time step is available. I am really facing trouble as I am new in WRFDA. For your reference I am attaching rsl.out file and namelist.input file here.
Anyone will be appreciated to give me some advice on how to find the reasons and solving the problem? Thanks!

*** VARIATIONAL ANALYSIS ***
WRFDA V4.3

Ntasks in X 1 , ntasks in Y 1
*************************************
Parent domain
ids,ide,jds,jde 1 215 1 215
ims,ime,jms,jme -4 220 -4 220
ips,ipe,jps,jpe 1 215 1 215
*************************************
DYNAMICS OPTION: Eulerian Mass Coordinate
alloc_space_field: domain 1 , 1421526244 bytes allocated
Input data is acceptable to use: fg
hybrid_opt = 2

use_theta_m = 1

Tile Strategy is not specified. Assuming 1D-Y
WRF TILE 1 IS 1 IE 214 JS 1 JE 214
WRF NUMBER OF TILES = 1
Domain mapping info:
map_proj = 6
cen_lat = 0.220000E+02
cen_lon = 0.900000E+02
truelat1 = 0.220000E+02
truelat2 = 0.220000E+02
start_lat = 0.133786E+02
start_lon = 0.813786E+02
pole_lat = 0.900000E+02
dsm = 0.900000E+01

Set up observations (ob)

Using ASCII format observation input

--------------------------- WARNING ---------------------------
WARNING FROM FILE: da_scan_obs_ascii.inc LINE: 68
Error 2 opening gts obs file ob.ascii
---------------------------------------------------------------
--------------------------- WARNING ---------------------------
WARNING FROM FILE: da_read_obs_ascii.inc LINE: 106
Error 2 opening gts obs file ob.ascii
---------------------------------------------------------------
Using NCEP BUFR radiance 1b input

ACCoeff_ReadFile(Binary)(INFORMATION) : FILE: ./crtm_coeffs/amsua_n15.SpcCoeff.bin;
ACCoeff RELEASE.VERSION: 1.04
N_FOVS=30 N_CHANNELS=15
SpcCoeff_ReadFile(Binary)(INFORMATION) : FILE: ./crtm_coeffs/amsua_n15.SpcCoeff.bin;
SpcCoeff RELEASE.VERSION: 8.01
N_CHANNELS=15
ACCoeff_ReadFile(Binary)(INFORMATION) : FILE: ./crtm_coeffs/amsua_n16.SpcCoeff.bin;
ACCoeff RELEASE.VERSION: 1.04
N_FOVS=30 N_CHANNELS=15
SpcCoeff_ReadFile(Binary)(INFORMATION) : FILE: ./crtm_coeffs/amsua_n16.SpcCoeff.bin;
SpcCoeff RELEASE.VERSION: 8.01
N_CHANNELS=15
ACCoeff_ReadFile(Binary)(INFORMATION) : FILE: ./crtm_coeffs/amsua_n18.SpcCoeff.bin;
ACCoeff RELEASE.VERSION: 1.04
N_FOVS=30 N_CHANNELS=15
SpcCoeff_ReadFile(Binary)(INFORMATION) : FILE: ./crtm_coeffs/amsua_n18.SpcCoeff.bin;
SpcCoeff RELEASE.VERSION: 8.01
N_CHANNELS=15
SpcCoeff_ReadFile(Binary)(INFORMATION) : FILE: ./crtm_coeffs/amsua_aqua.SpcCoeff.bin;
SpcCoeff RELEASE.VERSION: 8.01
N_CHANNELS=15
ACCoeff_ReadFile(Binary)(INFORMATION) : FILE: ./crtm_coeffs/amsua_metop-a.SpcCoeff.bin;
ACCoeff RELEASE.VERSION: 1.04
N_FOVS=30 N_CHANNELS=15
SpcCoeff_ReadFile(Binary)(INFORMATION) : FILE: ./crtm_coeffs/amsua_metop-a.SpcCoeff.bin;
SpcCoeff RELEASE.VERSION: 8.01
N_CHANNELS=15
Read_ODPS_Binary(INFORMATION) : FILE: ./crtm_coeffs/amsua_n15.TauCoeff.bin;
ODPS RELEASE.VERSION: 2.01 N_LAYERS=100 N_COMPONENTS=2 N_ABSORBERS=1 N_CHANNELS=15 N_COEFFS=21600
Read_ODPS_Binary(INFORMATION) : FILE: ./crtm_coeffs/amsua_n16.TauCoeff.bin;
ODPS RELEASE.VERSION: 2.01 N_LAYERS=100 N_COMPONENTS=2 N_ABSORBERS=1 N_CHANNELS=15 N_COEFFS=21600
Read_ODPS_Binary(INFORMATION) : FILE: ./crtm_coeffs/amsua_n18.TauCoeff.bin;
ODPS RELEASE.VERSION: 2.01 N_LAYERS=100 N_COMPONENTS=2 N_ABSORBERS=1 N_CHANNELS=15 N_COEFFS=21600
Read_ODPS_Binary(INFORMATION) : FILE: ./crtm_coeffs/amsua_aqua.TauCoeff.bin;
ODPS RELEASE.VERSION: 2.01 N_LAYERS=100 N_COMPONENTS=2 N_ABSORBERS=1 N_CHANNELS=15 N_COEFFS=21600
Read_ODPS_Binary(INFORMATION) : FILE: ./crtm_coeffs/amsua_metop-a.TauCoeff.bin;
ODPS RELEASE.VERSION: 2.01 N_LAYERS=100 N_COMPONENTS=2 N_ABSORBERS=1 N_CHANNELS=15 N_COEFFS=21600
CloudCoeff_ReadFile(Binary)(INFORMATION) : FILE: ./crtm_coeffs/CloudCoeff.bin;
CloudCoeff RELEASE.VERSION: 3.04 N_FREQUENCIES(MW)= 31 N_FREQUENCIES(IR)= 61 N_RADII(MW)=10 N_RADII(IR)=10 N_TEMPERATURES= 5 N_DENSITIES= 3 N_LEGENDRE_TERMS=38 N_PHASE_ELEMENTS= 1
AerosolCoeff_ReadFile(Binary)(INFORMATION) : FILE: ./crtm_coeffs/AerosolCoeff.bin;
AerosolCoeff RELEASE.VERSION: 3.04 N_WAVELENGTHS= 61 N_RADII= 36 N_TYPES= 8 N_RH= 36 N_LEGENDRE_TERMS=37 N_PHASE_ELEMENTS= 1
FitCoeff_ReadFile(INFORMATION) : FILE: ./crtm_coeffs/FASTEM5.MWwater.EmisCoeff.bin;
FitCoeff RELEASE.VERSION : 1.5; DIMENSIONS= 2
FitCoeff_ReadFile(INFORMATION) : FILE: ./crtm_coeffs/FASTEM5.MWwater.EmisCoeff.bin;
FitCoeff RELEASE.VERSION : 1.5; DIMENSIONS= 6
FitCoeff_ReadFile(INFORMATION) : FILE: ./crtm_coeffs/FASTEM5.MWwater.EmisCoeff.bin;
FitCoeff RELEASE.VERSION : 1.5; DIMENSIONS= 3, 7, 2
FitCoeff_ReadFile(INFORMATION) : FILE: ./crtm_coeffs/FASTEM5.MWwater.EmisCoeff.bin;
FitCoeff RELEASE.VERSION : 1.5; DIMENSIONS= 10, 4, 3
FitCoeff_ReadFile(INFORMATION) : FILE: ./crtm_coeffs/FASTEM5.MWwater.EmisCoeff.bin;
FitCoeff RELEASE.VERSION : 1.5; DIMENSIONS= 8
FitCoeff_ReadFile(INFORMATION) : FILE: ./crtm_coeffs/FASTEM5.MWwater.EmisCoeff.bin;
FitCoeff RELEASE.VERSION : 1.5; DIMENSIONS= 3, 6, 2
MWwaterLUT_ReadFile(INFORMATION) : FILE: ./crtm_coeffs/FASTEM5.MWwater.EmisCoeff.bin;
MWwaterLUT RELEASE.VERSION: 1.01
N_ANGLES=15 N_FREQUENCIES=18 N_TEMPERATURES=12 N_WIND_SPEEDS=16
MWwaterCoeff_ReadFile(INFORMATION) : FILE: ./crtm_coeffs/FASTEM5.MWwater.EmisCoeff.bin;
MWwaterCoeff RELEASE.VERSION: 1.5
Reading radiance 1b data from amsua.bufr
Bufr file date is 2021 6 20 0 amsua
num_tovs_file num_tovs_global num_tovs_local num_tovs_used num_tovs_thinned
481524 0 0 0 0
Observation summary
ob time 1

----------------------------------------------------------
[3.0] Set up background errors (be) for cv_option: 3
----------------------------------------------------------

DA_Setup_Background_Errors: 3DVAR dry control variables are:
psi, chi_u, t_u and psfc
DA_Setup_Background_Errors: 3DVAR humidity control variable is q/qsg

Number of vertical level for stats = 42
Number of latitude nlath = 96

Number of vertical level for WRFVar= 32

VARBC: Reading VARBC.in file
VARBC: no matching for platform/satid/sensor 1 19 3
VARBC: no matching for platform/satid/sensor 23 1 40
VARBC: no matching for platform/satid/sensor 23 2 40
VARBC: no matching for platform/satid/sensor 10 1 3
VARBC: no matching for platform/satid/sensor 1 15 4
VARBC: no matching for platform/satid/sensor 1 16 4
VARBC: no matching for platform/satid/sensor 1 17 4
VARBC: no matching for platform/satid/sensor 1 18 15
VARBC: no matching for platform/satid/sensor 1 19 15
VARBC: no matching for platform/satid/sensor 23 1 41
VARBC: no matching for platform/satid/sensor 23 2 41
VARBC: no matching for platform/satid/sensor 23 3 73
VARBC: no matching for platform/satid/sensor 10 1 15
VARBC: no matching for platform/satid/sensor 10 2 15
VARBC: no matching for platform/satid/sensor 17 0 19
VARBC: no matching for platform/satid/sensor 2 16 10
VARBC: no matching for platform/satid/sensor 9 2 11
VARBC: no matching for platform/satid/sensor 19 0 0
VARBC: no matching for platform/satid/sensor 1 11 0
VARBC: no matching for platform/satid/sensor 1 12 0
VARBC: no matching for platform/satid/sensor 1 14 0
VARBC: no matching for platform/satid/sensor 19 0 1
VARBC: no matching for platform/satid/sensor 1 11 1
VARBC: no matching for platform/satid/sensor 1 12 1
VARBC: no matching for platform/satid/sensor 1 14 1
VARBC: no matching for platform/satid/sensor 10 1 16
VARBC: no matching for platform/satid/sensor 10 2 16
VARBC: no matching for platform/satid/sensor 12 2 21
VARBC: no matching for platform/satid/sensor 12 3 21
VARBC: no matching for platform/satid/sensor 29 1 63
VARBC: no matching for platform/satid/sensor 31 8 56
VARBC: no matching for platform/satid/sensor 4 13 22
VARBC: no matching for platform/satid/sensor 4 14 22
Calculate innovation vector(iv)

VARBC: Cold-starting noaa-15-amsua 5 --> 0.00
VARBC: Cold-starting noaa-15-amsua 6 --> 0.00
VARBC: Cold-starting noaa-15-amsua 7 --> 0.00
VARBC: Cold-starting noaa-15-amsua 8 --> 0.00
VARBC: Cold-starting noaa-15-amsua 9 --> 0.00
VARBC: Cold-starting noaa-18-amsua 5 --> 0.00
VARBC: Cold-starting noaa-18-amsua 6 --> 0.00
VARBC: Cold-starting noaa-18-amsua 7 --> 0.00
VARBC: Cold-starting noaa-18-amsua 8 --> 0.00
VARBC: Cold-starting eos-2-amsua 6 --> 0.00
VARBC: Cold-starting eos-2-amsua 8 --> 0.00
VARBC: Cold-starting eos-2-amsua 9 --> 0.00
VARBC: Cold-starting metop-2-amsua 5 --> 0.00
VARBC: Cold-starting metop-2-amsua 6 --> 0.00
VARBC: Cold-starting metop-2-amsua 8 --> 0.00
VARBC: Cold-starting metop-2-amsua 9 --> 0.00
VARBC: Applying bias correction for noaa-15-amsua
VARBC: Applying bias correction for noaa-16-amsua
VARBC: Applying bias correction for noaa-18-amsua
VARBC: Applying bias correction for eos-2-amsua
VARBC: Applying bias correction for metop-2-amsua
VARBC: Estimate Hessian for preconditioning
VARBC: 0 active observations for noaa-15-amsua channel 5
VARBC: 0 active observations for noaa-15-amsua channel 6
VARBC: 0 active observations for noaa-15-amsua channel 7
VARBC: 0 active observations for noaa-15-amsua channel 8
VARBC: 0 active observations for noaa-15-amsua channel 9
VARBC: 0 active observations for noaa-18-amsua channel 5
VARBC: 0 active observations for noaa-18-amsua channel 6
VARBC: 0 active observations for noaa-18-amsua channel 7
VARBC: 0 active observations for noaa-18-amsua channel 8
VARBC: 0 active observations for eos-2-amsua channel 6
VARBC: 0 active observations for eos-2-amsua channel 8
VARBC: 0 active observations for eos-2-amsua channel 9
VARBC: 0 active observations for metop-2-amsua channel 5
VARBC: 0 active observations for metop-2-amsua channel 6
VARBC: 0 active observations for metop-2-amsua channel 8
VARBC: 0 active observations for metop-2-amsua channel 9
Minimize cost function using CG method


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
Final value of Jm = 0.00000
Final value of Jt = 0.00000
Jb factor used = 1.00000
Je factor used = 1.00000
VarBC factor used = 1.00000

Total number of radiances = 0
Cost function for radiances = 0.00000

Writing radiance OMA ascii file

VARBC: Updating bias parameters
VARBC: Not enough data to keep statistics for noaa-15-amsua 5
VARBC: Not enough data to keep statistics for noaa-15-amsua 6
VARBC: Not enough data to keep statistics for noaa-15-amsua 7
VARBC: Not enough data to keep statistics for noaa-15-amsua 8
VARBC: Not enough data to keep statistics for noaa-15-amsua 9
VARBC: Not enough data to keep statistics for noaa-18-amsua 5
VARBC: Not enough data to keep statistics for noaa-18-amsua 6
VARBC: Not enough data to keep statistics for noaa-18-amsua 7
VARBC: Not enough data to keep statistics for noaa-18-amsua 8
VARBC: Not enough data to keep statistics for eos-2-amsua 6
VARBC: Not enough data to keep statistics for eos-2-amsua 8
VARBC: Not enough data to keep statistics for eos-2-amsua 9
VARBC: Not enough data to keep statistics for metop-2-amsua 5
VARBC: Not enough data to keep statistics for metop-2-amsua 6
VARBC: Not enough data to keep statistics for metop-2-amsua 8
VARBC: Not enough data to keep statistics for metop-2-amsua 9
VARBC: Writing information in VARBC.out file
convert T to THM when use_theta_m = 1

*** WRF-Var completed successfully ***


******************namelist.input********************
&wrfvar1
var4d=.false.
print_detail_grad=.false.,
/
&wrfvar2
/
&wrfvar3
fg_format=1,
ob_format=1,
ob_format_gpsro=1,
/
&wrfvar4
thin_conv=.true.
use_hirs4obs=.false.
use_mhsobs=.false.
use_amsuaobs=.true.
use_amsubobs=.false.
use_airsobs=.false.
use_atmsobs=.false.
use_iasiobs=.false.
use_seviriobs=.false.
/
&wrfvar5
/
&wrfvar6
/
&wrfvar7
cv_options=3,
/
&wrfvar8
/
&wrfvar9
/
&wrfvar10
/
&wrfvar11
/
&wrfvar12
/
&wrfvar13
/
&wrfvar14
rtminit_nsensor=5
rtminit_platform=1,1,1,9,10
rtminit_satid=15,16,18,2,2
rtminit_sensor=3,3,3,3,3
thinning_mesh=120.0,120.0,120.0,120.0,120.0
thinning=true,
qc_rad=true,
write_iv_rad_ascii=false,
write_oa_rad_ascii=true,
rtm_option=2,
only_sea_rad=false,
use_varbc=true,
use_crtm_kmatrix=true,
use_rttov_kmatrix=true,
crtm_coef_path='./crtm_coeffs'
crtm_mwwater_coef='FASTEM5.MWwater.EmisCoeff.bin'
crtm_irwater_coef='Nalli.IRwater.EmisCoeff.bin'
crtm_irland_coef='USGS.IRland.EmisCoeff.bin'
rttov_emis_atlas_ir=0,
rttov_emis_atlas_mw=0,
/
&wrfvar15
/
&wrfvar16
/
&wrfvar17
/
&wrfvar18
analysis_date="2021-06-23_00:00:00.0000",
/
&wrfvar19
/
&wrfvar20
/
&wrfvar21
time_window_min="2021-06-21_00:00:00.0000",
/
&wrfvar22
time_window_max="2021-06-21_18:00:00.0000",
/
&time_control
start_year=2021,
start_month=06,
start_day=21,
start_hour=00,
end_year=2021,
end_month=06,
end_day=30,
end_hour=00,
/
&fdda
/
&domains
e_we=215,
e_sn=215,
e_vert=33,
dx=9000,
dy=9000,
/
&dfi_control
/
&tc
/
&physics
mp_physics=3,
ra_lw_physics=1,
ra_sw_physics=1,
radt=60,
sf_sfclay_physics=1,
sf_surface_physics=1,
bl_pbl_physics=1,
cu_physics=1,
cudt=5,
num_soil_layers=5,
mp_zero_out=2,
co2tf=0,
/
&scm
/
&dynamics
/
&bdy_control
/
&grib2
/
&fire
/
&namelist_quilt
/
&perturbation
/
Sabrina93
 
Posts: 2
Joined: Sun May 16, 2021 2:17 am

Return to Runtime Problems

Who is online

Users browsing this forum: No registered users and 9 guests

cron