[Wrf-users] WRFDA negative height problem
Sean Crowell
Sean.Crowell at noaa.gov
Tue Nov 29 12:38:28 MST 2011
I am attempting to run WRFDA, and I have generated a background error
file, a set of observations, and a first guess from WRF using NAM data
as initial and boundary conditions. When I run da_wrfvar.exe, I get
the following error message in my rsl.error.*:
Ntasks in X 2, ntasks in Y 3
*************************************
Parent domain
ids,ide,jds,jde 1 266 1 172
ims,ime,jms,jme 127 271 -4 64
ips,ipe,jps,jpe 134 266 1 57
*************************************
DYNAMICS OPTION: Eulerian Mass Coordinate
alloc_space_field: domain 1, 852723024
bytes allocat
ed
WRF TILE 1 IS 134 IE 265 JS 1 JE 57
WRF NUMBER OF TILES = 1
---------------------------- FATAL ERROR -----------------------
Fatal error in file: <A HREF="http://www.mmm.ucar.edu/people/wrfhelp/wrfvar/cod
e/trunk/da_transfer_xatowrf.html">da_transfer_xatowrf.inc</a>
LINE: 442
Negative height found
222 22 ht = -1.19 terr = 7.11
----------------------------------------------------------------
Can you help me to understand what might be going wrong here? I'm
certain that I've set something wrong, but I don't know what. I
copied my WRF namelist over to my WRFDA namelist, so I know the
settings should match. For completeness, here's the rsl.out.* file:
*** VARIATIONAL ANALYSIS ***
Ntasks in X 2, ntasks in Y 3
*************************************
Parent domain
ids,ide,jds,jde 1 266 1 172
ims,ime,jms,jme 127 271 -4 64
ips,ipe,jps,jpe 134 266 1 57
*************************************
DYNAMICS OPTION: Eulerian Mass Coordinate
alloc_space_field: domain 1, 852723024
bytes allocat
ed
WRF TILE 1 IS 134 IE 265 JS 1 JE 57
WRF NUMBER OF TILES = 1
Set up observations (ob)
Using ASCII format observation input
scan obs ascii
end scan obs ascii
Observation summary
ob time 1
sound 79 global, 13 local
synop 961 global, 114 local
pilot 0 global, 0 local
satem 0 global, 0 local
geoamv 0 global, 0 local
polaramv 0 global, 0 local
airep 1603 global, 369 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 0 global, 0 local
airs retrieval 0 global, 0 local
sonde_sfc 79 global, 13 local
mtgirs 0 global, 0 local
tamdar 0 global, 0 local
Set up background errors for regional application for cv_options = 5
Using the averaged regression coefficients for unbalanced part
WRF-Var dry control variables are:psi, chi_u, t_u and ps_u
Humidity control variable is rh
Vertical truncation for psi = 13( 99.00%)
Vertical truncation for chi_u = 16( 99.00%)
Vertical truncation for t_u = 25( 99.00%)
Vertical truncation for rh = 23( 99.00%)
>>> Save the variances and scale-lengths in outer-loop 1
Scaling: var, len, ds: 0.100000E+01 0.100000E+01 0.180000E+05
Scaling: var, len, ds: 0.100000E+01 0.100000E+01 0.180000E+05
Scaling: var, len, ds: 0.100000E+01 0.100000E+01 0.180000E+05
Scaling: var, len, ds: 0.100000E+01 0.100000E+01 0.180000E+05
Scaling: var, len, ds: 0.100000E+01 0.100000E+01 0.180000E+05
Calculate innovation vector(iv)
Minimize cost function using CG method
jo_sound 0.262194959846E+04
jo%sound_u 0.149525238198E+04
jo%sound_v 0.756634251656E+03
jo%sound_t 0.253029644007E+03
jo%sound_q 0.117033320818E+03
jo_sonde_sfc 0.000000000000E+00
jo%sonde_sfc_u 0.000000000000E+00
jo%sonde_sfc_v 0.000000000000E+00
jo%sonde_sfc_t 0.000000000000E+00
jo%sonde_sfc_p 0.000000000000E+00
jo%sonde_sfc_q 0.000000000000E+00
jo_synop 0.636955314734E+08
jo%synop_u 0.491152971520E+03
jo%synop_v 0.304559122743E+04
jo%synop_t 0.829944173431E+07
jo%synop_p 0.553925320036E+08
jo%synop_q 0.209912537612E+02
jo%total 0.637064792883E+08
jo_sound 0.262194959846E+04
jo_sonde_sfc 0.000000000000E+00
jo_geoamv 0.000000000000E+00
jo_polaramv 0.000000000000E+00
jo_synop 0.636955314734E+08
jo_satem 0.000000000000E+00
jo_pilot 0.000000000000E+00
jo_airep 0.832586532811E+04
jo_metar 0.000000000000E+00
jo_ships 0.000000000000E+00
jo_gpspw 0.000000000000E+00
jo_ssmi_tb 0.000000000000E+00
jo_ssmi_rv 0.000000000000E+00
jo_ssmt1 0.000000000000E+00
jo_ssmt2 0.000000000000E+00
jo_pseudo 0.000000000000E+00
jo_qscat 0.000000000000E+00
jo_profiler 0.000000000000E+00
jo_buoy 0.000000000000E+00
jo_radar 0.000000000000E+00
jo_gpsref 0.000000000000E+00
jo_bogus 0.000000000000E+00
jo_radiance 0.000000000000E+00
jo_airsr 0.000000000000E+00
jo_mtgirs 0.000000000000E+00
jo_tamdar 0.000000000000E+00
jo_tamdar_sfc 0.000000000000E+00
Starting outer iteration : 1
Starting cost function: 4.25254199D+08, Gradient= 1.98261367D+05
For this outer iteration gradient target is: 1.98261367D+03
----------------------------------------------------------
Iter Cost Function Gradient Step
1 1.71179200D+08 8.41493011D+04 1.29275354D-02
2 9.54624377D+07 4.47172535D+04 2.13855787D-02
3 7.16577641D+07 3.08483006D+04 2.38090454D-02
4 5.95942218D+07 1.72014545D+04 2.53537591D-02
5 5.47628019D+07 1.56966639D+04 3.26568609D-02
6 5.16598630D+07 9.96024256D+03 2.51876992D-02
7 4.95550876D+07 8.22891171D+03 4.24322367D-02
8 4.82312735D+07 6.35948495D+03 3.90995895D-02
9 4.73919322D+07 5.09046293D+03 4.15073184D-02
10 4.69090645D+07 4.20757477D+03 3.72686462D-02
11 4.66196330D+07 3.79322360D+03 3.26973125D-02
12 4.64207492D+07 2.83758913D+03 2.76447402D-02
13 4.62677592D+07 2.03484654D+03 3.80009081D-02
14 4.61558173D+07 1.48366366D+03 5.40703801D-02
----------------------------------------------------------
Inner iteration stopped after 14 iterations
jo_sound 0.686012931808E+06
jo%sound_u 0.122287370327E+06
jo%sound_v 0.149678539580E+05
jo%sound_t 0.548725663831E+06
jo%sound_q 0.320436913736E+02
jo_sonde_sfc 0.000000000000E+00
jo%sonde_sfc_u 0.000000000000E+00
jo%sonde_sfc_v 0.000000000000E+00
jo%sonde_sfc_t 0.000000000000E+00
jo%sonde_sfc_p 0.000000000000E+00
jo%sonde_sfc_q 0.000000000000E+00
jo_synop 0.209394523833E+07
jo%synop_u 0.272392793587E+05
jo%synop_v 0.115704332054E+05
jo%synop_t 0.180881916304E+07
jo%synop_p 0.246310098594E+06
jo%synop_q 0.626412325805E+01
jo%total 0.329586223340E+07
jo_sound 0.686012931808E+06
jo_sonde_sfc 0.000000000000E+00
jo_geoamv 0.000000000000E+00
jo_polaramv 0.000000000000E+00
jo_synop 0.209394523833E+07
jo_satem 0.000000000000E+00
jo_pilot 0.000000000000E+00
jo_airep 0.515904063267E+06
jo_metar 0.000000000000E+00
jo_ships 0.000000000000E+00
jo_gpspw 0.000000000000E+00
jo_ssmi_tb 0.000000000000E+00
jo_ssmi_rv 0.000000000000E+00
jo_ssmt1 0.000000000000E+00
jo_ssmt2 0.000000000000E+00
jo_pseudo 0.000000000000E+00
jo_qscat 0.000000000000E+00
jo_profiler 0.000000000000E+00
jo_buoy 0.000000000000E+00
jo_radar 0.000000000000E+00
jo_gpsref 0.000000000000E+00
jo_bogus 0.000000000000E+00
jo_radiance 0.000000000000E+00
jo_airsr 0.000000000000E+00
jo_mtgirs 0.000000000000E+00
jo_tamdar 0.000000000000E+00
jo_tamdar_sfc 0.000000000000E+00
Final: 14 iter, J= 4.61558173D+07, g= 1.48366366D+03
----------------------------------------------------------
---------------------------- FATAL ERROR -----------------------
Fatal error in file: <A HREF="http://www.mmm.ucar.edu/people/wrfhelp/wrfvar/cod
e/trunk/da_transfer_xatowrf.html">da_transfer_xatowrf.inc</a>
LINE: 442
Negative height found
222 22 ht = -1.19 terr = 7.11
----------------------------------------------------------------
taskid: 1 hostname: 4351-crowell.winstorm.nssl
taskid: 1 hostname: 4351-crowell.winstorm.nssl
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/wrf-users/attachments/20111129/00f48b27/attachment-0001.html
More information about the Wrf-users
mailing list