Reading radiance 1b data from amsua.bufr.......

[bellanchor]

Hi, everyone
I get a problem.
Today I tried to run WRFDA with Radiance data, and I just followed the User's Guide.
Downloaded two bufr files, converted into little endian, and tried to use them.
But the program da_wrfvar.exe was stuck in the ' Reading radiance 1b data from amsua.bufr ' line for more than one hour, and still, but no error occured. Like this:

*** VARIATIONAL ANALYSIS ***
DYNAMICS OPTION: Eulerian Mass Coordinate
alloc_space_field: domain 1 , 1632696704 bytes allocated
WRF TILE 1 IS 1 IE 159 JS 1 JE 139
WRF NUMBER OF TILES = 1
Set up observations (ob)

Using ASCII format observation input

scan obs ascii
end scan obs ascii
Using NCEP BUFR radiance 1b input

Read_SpcCoeff_Binary(INFORMATION) : FILE: crtm_coeffs/amsua_n15.SpcCoeff.bin;
SpcCoeff RELEASE.VERSION: 7.03 N_CHANNELS=15
amsua_n15 AntCorr RELEASE.VERSION: 1.04 N_FOVS=30 N_CHANNELS=15
Read_SpcCoeff_Binary(INFORMATION) : FILE: crtm_coeffs/amsua_n16.SpcCoeff.bin;
SpcCoeff RELEASE.VERSION: 7.03 N_CHANNELS=15
amsua_n16 AntCorr RELEASE.VERSION: 1.04 N_FOVS=30 N_CHANNELS=15
Read_SpcCoeff_Binary(INFORMATION) : FILE: crtm_coeffs/amsua_n18.SpcCoeff.bin;
SpcCoeff RELEASE.VERSION: 7.03 N_CHANNELS=15
amsua_n18 AntCorr RELEASE.VERSION: 1.04 N_FOVS=30 N_CHANNELS=15
Read_SpcCoeff_Binary(INFORMATION) : FILE: crtm_coeffs/amsua_n19.SpcCoeff.bin;
SpcCoeff RELEASE.VERSION: 7.03 N_CHANNELS=15
amsua_n19 AntCorr RELEASE.VERSION: 1.05 N_FOVS=30 N_CHANNELS=15
Read_SpcCoeff_Binary(INFORMATION) : FILE: crtm_coeffs/amsua_metop-a.SpcCoeff.bin;
SpcCoeff RELEASE.VERSION: 7.02 N_CHANNELS=15
amsua_metop-a AntCorr RELEASE.VERSION: 1.01 N_FOVS=30 N_CHANNELS=15
Read_SpcCoeff_Binary(INFORMATION) : FILE: crtm_coeffs/amsub_n15.SpcCoeff.bin;
SpcCoeff RELEASE.VERSION: 7.03 N_CHANNELS=5
amsub_n15 AntCorr RELEASE.VERSION: 1.04 N_FOVS=90 N_CHANNELS=5
Read_SpcCoeff_Binary(INFORMATION) : FILE: crtm_coeffs/amsub_n16.SpcCoeff.bin;
SpcCoeff RELEASE.VERSION: 7.03 N_CHANNELS=5
amsub_n16 AntCorr RELEASE.VERSION: 1.04 N_FOVS=90 N_CHANNELS=5
Read_SpcCoeff_Binary(INFORMATION) : FILE: crtm_coeffs/amsub_n17.SpcCoeff.bin;
SpcCoeff RELEASE.VERSION: 7.03 N_CHANNELS=5
amsub_n17 AntCorr RELEASE.VERSION: 1.04 N_FOVS=90 N_CHANNELS=5
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_n19.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
Read_ODPS_Binary(INFORMATION) : FILE: crtm_coeffs/amsub_n15.TauCoeff.bin;
ODPS RELEASE.VERSION: 2.01 N_LAYERS=100 N_COMPONENTS=2 N_ABSORBERS=1 N_CHANNELS=5 N_COEFFS=10300
Read_ODPS_Binary(INFORMATION) : FILE: crtm_coeffs/amsub_n16.TauCoeff.bin;
ODPS RELEASE.VERSION: 2.01 N_LAYERS=100 N_COMPONENTS=2 N_ABSORBERS=1 N_CHANNELS=5 N_COEFFS=10300
Read_ODPS_Binary(INFORMATION) : FILE: crtm_coeffs/amsub_n17.TauCoeff.bin;
ODPS RELEASE.VERSION: 2.01 N_LAYERS=100 N_COMPONENTS=2 N_ABSORBERS=1 N_CHANNELS=5 N_COEFFS=10300
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
Read_EmisCoeff_Binary(INFORMATION) : FILE: crtm_coeffs/EmisCoeff.bin;
EmisCoeff RELEASE.VERSION: 2.01 N_ANGLES= 67 N_FREQUENCIES= 153 N_WIND_SPEEDS= 23
Reading radiance 1b data from amsua.bufr

I don't know what the meaning of this is.
Anyone could give me a hand?
Thanks!

Best regards!

[bellanchor]

I'm using Gfortran , under Ubuntu 10.04.
And the WRF Var Known Problems said:

WRF Var Version 3.1 : Known Problems and Fixes
unblock and reblock radiance bufr data with gfortran (Posted 4/9/09)

Problem:
On linux with gfortran compiler, the radiance bufr data unblocked and reblocked by cwordsh might be unable to be read by WRF-Var.

Solution:
Peding

I don't know what's the meaning of " the radiance bufr data unblocked and reblocked by cwordsh".
And I installed the Version 3.3 not 3.1.
Anyone could help me?

[hnlim]

Download the cwordsh package from the following link

http://www.nco.ncep.noaa.gov/sib/decode ... c/cwordsh/

Then run your amsu data through that package.

Reason for doing so is explained in the website.

Agnes

[bellanchor]

Hi,Agnes.
Thank you for your reply!
But I'm not so clear about the relationship between the ssrc.exe utility and cwordsh utility.
Both do the same thing? Should I use ssrc.exe to convert the big endian to little endian and use cwordsh to reblock? Or just use cwordsh?
As I tried in some ways, it wasn't going any further, still.

[bellanchor]

I just took a look at the User's Guide of WRFV3.1, it mentioned the cwordsh utility.
And in the WRFV3.3, the Guide uses ssrc instead.
So I'm thinking one of these two utilities is enough.
But neither works well.

[hnlim]

I have not used the ssrc utility.

The cwordsh utility should work fine if you download your AMSU data from NCAR.

Agnes

[mkudsy]

You are not alone.
I let it run for one night but nothing out.
I have a friend that was successfully run it using version 3.2.1. Strange enough that I could not run it using either vesions 3.1 3.2 3.2.1 and 3.3
I use Fedora 14, Gfortran and gcc
The testing data of 3.3 has gdas data probably made on IBM computer, so it must be in big endian, and fortran blocked
To use it under Linux it must be unblocked using cwordsh, change it to little-endian using ssrc.exe , and then reblocked using cwordsh...This is my idea.
So let's try it and let informed.
Regards
Mahally

[bellanchor]

Hi Mahally,
I'm wondering whether you have solved this problem or not.
Following your idea, "standard input: unexpected short record 1" showed out while using ssrc.exe.

Regards.
Liu

[xinhe]

Hi, I use WRFDA_V3.2.1 to run the radiance data assimilation,and I can successfully complete it.
I use ifort ,when I set 'USE_CRTM_KMATRIX=true' in my namelist.input,the the run time is about 20 minutes,but when I set 'USE_CRTM_KMATRIX=false',the run time will reached about 40 minutes.I don't know if this tips will hlep you. Good luck!

[bellanchor]

Hi xinhe,
thanks for your suggestion,but unfortunately the problem still exists.
Maybe I should change a fortran compiler and try again.

Liu