WRFHELP Please~ Radiance Assimilation

Issues with running 3DVAR.

WRFHELP Please~ Radiance Assimilation

Postby chantingfai » Mon Aug 03, 2009 4:28 am

Dear all,

Everything was fine when running da_wrfvar.exe until......Cold-starting.....it stopped.

taskid: 0 hostname: node0001
*** VARIATIONAL ANALYSIS ***
Ntasks in X 2 , ntasks in Y 3
*************************************
Parent domain
ids,ide,jds,jde 1 401 1 321
ims,ime,jms,jme -4 207 -4 114
ips,ipe,jps,jpe 1 200 1 107
*************************************
DYNAMICS OPTION: Eulerian Mass Coordinate
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; ^M
SpcCoeff RELEASE.VERSION: 7.03 N_CHANNELS=15^M
amsua_n15 AntCorr RELEASE.VERSION: 1.04 N_FOVS=30 N_CHANNELS=15
Read_SpcCoeff_Binary(INFORMATION) : FILE: crtm_coeffs/amsua_n16.SpcCoeff.bin; ^M
SpcCoeff RELEASE.VERSION: 7.03 N_CHANNELS=15^M
amsua_n16 AntCorr RELEASE.VERSION: 1.04 N_FOVS=30 N_CHANNELS=15
Read_SpcCoeff_Binary(INFORMATION) : FILE: crtm_coeffs/amsua_n18.SpcCoeff.bin; ^M
SpcCoeff RELEASE.VERSION: 7.03 N_CHANNELS=15^M
amsua_n18 AntCorr RELEASE.VERSION: 1.04 N_FOVS=30 N_CHANNELS=15
Read_SpcCoeff_Binary(INFORMATION) : FILE: crtm_coeffs/amsua_aqua.SpcCoeff.bin; ^M
SpcCoeff RELEASE.VERSION: 7.01 N_CHANNELS=15
Read_SpcCoeff_Binary(INFORMATION) : FILE: crtm_coeffs/amsua_metop-a.SpcCoeff.bin; ^M
SpcCoeff RELEASE.VERSION: 7.03 N_CHANNELS=15^M
amsua_metop-a AntCorr RELEASE.VERSION: 1.04 N_FOVS=30 N_CHANNELS=15
Read_TauCoeff_Binary(INFORMATION) : FILE: crtm_coeffs/amsua_n15.TauCoeff.bin; ^M
TauCoeff RELEASE.VERSION: 5.06 N_ORDERS=10 N_PREDICTORS= 6 N_ABSORBERS= 3 N_CHANNELS= 15 N_SENSORS= 1
Read_TauCoeff_Binary(INFORMATION) : FILE: crtm_coeffs/amsua_n16.TauCoeff.bin; ^M
TauCoeff RELEASE.VERSION: 5.06 N_ORDERS=10 N_PREDICTORS= 6 N_ABSORBERS= 3 N_CHANNELS= 15 N_SENSORS= 1
Read_TauCoeff_Binary(INFORMATION) : FILE: crtm_coeffs/amsua_n18.TauCoeff.bin; ^M
TauCoeff RELEASE.VERSION: 5.06 N_ORDERS=10 N_PREDICTORS= 6 N_ABSORBERS= 3 N_CHANNELS= 15 N_SENSORS= 1
Read_TauCoeff_Binary(INFORMATION) : FILE: crtm_coeffs/amsua_aqua.TauCoeff.bin; ^M
TauCoeff RELEASE.VERSION: 5.06 N_ORDERS=10 N_PREDICTORS= 6 N_ABSORBERS= 3 N_CHANNELS= 15 N_SENSORS= 1
Read_TauCoeff_Binary(INFORMATION) : FILE: crtm_coeffs/amsua_metop-a.TauCoeff.bin; ^M
TauCoeff RELEASE.VERSION: 5.06 N_ORDERS=10 N_PREDICTORS= 6 N_ABSORBERS= 3 N_CHANNELS= 15 N_SENSORS= 1
Read_CloudCoeff_Binary(INFORMATION) : FILE: crtm_coeffs/CloudCoeff.bin; ^M
CloudCoeff RELEASE.VERSION: 2.03 N_FREQUENCIES(MW)= 31 N_FREQUENCIES(IR)= 601 N_RADII(MW)=10 N_RADII(IR)=10 N_TEMPERATURES= 5 N_DENSITIES= 3 N_LEGENDRE_TERMS=38 N_PHASE_ELEMENTS= 1
Read_AerosolCoeff_Binary(INFORMATION) : FILE: crtm_coeffs/AerosolCoeff.bin; ^M
AerosolCoeff RELEASE.VERSION: 2.02 N_WAVELENGTHS= 91 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; ^M
EmisCoeff RELEASE.VERSION: 2.01 N_ANGLES= 67 N_FREQUENCIES= 153 N_WIND_SPEEDS= 23
Reading radiance 1b data from amsua.bufr
Bufr file date is 2008 6 23 0
amsua.bufr
num_tovs_file num_tovs_global num_tovs_local num_tovs_used num_tovs_thinned
324329 446 0 0 0
Reading radiance 1b data from amsub.bufr
Bufr file date is 2008 6 23 0
amsub.bufr
num_tovs_file num_tovs_global num_tovs_local num_tovs_used num_tovs_thinned
2191830 0 0 0 0
Observation summary
ob time 1
sound 63 global, 6 local
synop 0 global, 0 local
pilot 0 global, 0 local
satem 0 global, 0 local
geoamv 3214 global, 455 local
polaramv 0 global, 0 local
airep 269 global, 3 local
gpspw 0 global, 0 local
gpsrf 0 global, 0 local
metar 0 global, 0 local
ships 0 global, 0 local
ssmi_rv 0 global, 0 local
ssmi_tb 0 global, 0 local
ssmt1 0 global, 0 local
ssmt2 0 global, 0 local
qscat 0 global, 0 local
profiler 0 global, 0 local
buoy 0 global, 0 local
bogus 0 global, 0 local
pseudo 0 global, 0 local
radar 0 global, 0 local
radiance 446 global, 0 local
airs retrieval 0 global, 0 local
sonde_sfc 63 global, 6 local
mtgirs 0 global, 0 local
tamdar 0 global, 0 local

----------------------------------------------------------
[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= 49
VARBC: Reading VARBC.in file
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 10 2 15
VARBC: no matching for platform/satid/sensor 2 16 10
VARBC: no matching for platform/satid/sensor 9 2 11
Calculate innovation vector(iv)

CRTM Reference profile is: U.S. Standard Atmosphere
VARBC: Cold-starting noaa-15-amsua 1 --> 0.00
VARBC: Cold-starting noaa-15-amsua 2 --> 0.00
VARBC: Cold-starting noaa-15-amsua 3 --> 0.00
VARBC: Cold-starting noaa-15-amsua 4 --> 0.00
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-15-amsua 10 --> 0.00
VARBC: Cold-starting noaa-15-amsua 11 --> 0.00
VARBC: Cold-starting noaa-15-amsua 12 --> 0.00
VARBC: Cold-starting noaa-15-amsua 13 --> 0.00
VARBC: Cold-starting noaa-15-amsua 14 --> 0.00
VARBC: Cold-starting noaa-15-amsua 15 --> 0.00


I checked from rsl.error.0000 file. It said that....

DYNAMICS OPTION: Eulerian Mass Coordinate
WRF NUMBER OF TILES = 1
0: SUM: mask array must be default logical kind


What does it mean?
Does anyone can help me to solve the problem?

Thanks!
Last edited by chantingfai on Mon Aug 03, 2009 5:29 am, edited 1 time in total.
chantingfai
 
Posts: 26
Joined: Tue Apr 14, 2009 10:10 am

Re: HELP Please~ Radiance Assimilation

Postby chantingfai » Mon Aug 03, 2009 5:29 am

Or are there something wrong with SUM(array,[MASK]) in da_varbc.f ?
(Remark: I used PGI-f90 as $FC to compile WRF-Var. Does my compiler relate to this issue?)


I quoted the part which may be related in da_varbc.f
! Accumulate statistics for predictor mean/std
! ---------------------------------------------
do i = 1, npredmax
mean(i) = SUM(iv%instid(inst)%varbc_info%pred(i,1:num_rad), &
MASK=iv%instid(inst)%info%proc_domain(1,1:num_rad)) ! do not count HALO
rms(i) = SUM(iv%instid(inst)%varbc_info%pred(i,1:num_rad)**2, &
MASK=iv%instid(inst)%info%proc_domain(1,1:num_rad)) ! do not count HALO
end do
Last edited by chantingfai on Wed Aug 05, 2009 10:28 pm, edited 1 time in total.
chantingfai
 
Posts: 26
Joined: Tue Apr 14, 2009 10:10 am

Re: WRFHELP Please~ Radiance Assimilation

Postby chantingfai » Tue Aug 04, 2009 9:55 pm

One more finding!!

I found that if I just ONLY use air.bufr (gdas1.t00z.airsev.tm00.bufr_d), WRFVar can be run successfully!

However, it seems that no observation on my domain..... Does any observations of air.bufr cover south China domain?

For other files (e.g. gdas1.t00z.1bamua.tm00.bufr_d), WRFVar stop with previous error.
0: SUM: mask array must be default logical kind


Looking forward to your help~
chantingfai
 
Posts: 26
Joined: Tue Apr 14, 2009 10:10 am

Re: WRFHELP Please~ Radiance Assimilation

Postby chantingfai » Thu Aug 13, 2009 12:51 am

Problem was solved after installing WRF-Var Version 3.1.1 :)
chantingfai
 
Posts: 26
Joined: Tue Apr 14, 2009 10:10 am

Re: WRFHELP Please~ Radiance Assimilation

Postby taoyufeng » Thu Apr 08, 2010 10:10 am

Can you tell me your E-mail, I will turn to you for some questions.
My e-mail is tao1985yf@126.com.
Thank you !
taoyufeng
 
Posts: 4
Joined: Thu Apr 08, 2010 9:30 am


Return to Runtime Problems

Who is online

Users browsing this forum: No registered users and 2 guests