taskid: 0 hostname: master *** VARIATIONAL ANALYSIS *** Ntasks in X 2, ntasks in Y 2 ************************************* Parent domain ids,ide,jds,jde 1 90 1 60 ims,ime,jms,jme -4 52 -4 37 ips,ipe,jps,jpe 1 45 1 30 ************************************* DYNAMICS OPTION: Eulerian Mass Coordinate alloc_space_field: domain 1, 105610432 bytes allocated WRF TILE 1 IS 1 IE 45 JS 1 JE 30 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/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_TauCoeff_Binary(INFORMATION) : FILE: crtm_coeffs/amsua_n15.TauCoeff.bin; 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; 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; 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/amsub_n15.TauCoeff.bin; TauCoeff RELEASE.VERSION: 5.06 N_ORDERS=10 N_PREDICTORS= 6 N_ABSORBERS= 3 N_CHANNELS= 5 N_SENSORS= 1 Read_TauCoeff_Binary(INFORMATION) : FILE: crtm_coeffs/amsub_n16.TauCoeff.bin; TauCoeff RELEASE.VERSION: 5.06 N_ORDERS=10 N_PREDICTORS= 6 N_ABSORBERS= 3 N_CHANNELS= 5 N_SENSORS= 1 Read_TauCoeff_Binary(INFORMATION) : FILE: crtm_coeffs/amsub_n17.TauCoeff.bin; TauCoeff RELEASE.VERSION: 5.06 N_ORDERS=10 N_PREDICTORS= 6 N_ABSORBERS= 3 N_CHANNELS= 5 N_SENSORS= 1 Read_CloudCoeff_Binary(INFORMATION) : FILE: crtm_coeffs/CloudCoeff.bin; 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; 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; 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 2 5 12 amsua.bufr num_tovs_file num_tovs_global num_tovs_local num_tovs_used num_tovs_thinned 320937 4384 1401 406 995 Allocating space for radiance innov structure 1 noaa-15-amsua 325 Allocating space for radiance innov structure 2 noaa-16-amsua 81 Reading radiance 1b data from amsub.bufr Bufr file date is 2008 2 5 12 amsub.bufr num_tovs_file num_tovs_global num_tovs_local num_tovs_used num_tovs_thinned 2166160 39863 12564 428 12136 Allocating space for radiance innov structure 4 noaa-15-amsub 332 Allocating space for radiance innov structure 5 noaa-16-amsub 96 Observation summary ob time 1 sound 85 global, 19 local synop 531 global, 0 local pilot 84 global, 18 local satem 78 global, 42 local geoamv 736 global, 103 local polaramv 0 global, 0 local airep 132 global, 18 local gpspw 183 global, 61 local gpsrf 0 global, 0 local metar 1043 global, 155 local ships 86 global, 13 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 61 global, 30 local buoy 216 global, 19 local bogus 0 global, 0 local pseudo 0 global, 0 local radar 0 global, 0 local radiance 44247 global, 834 local airs retrieval 0 global, 0 local sonde_sfc 85 global, 19 local mtgirs 0 global, 0 local tamdar 0 global, 0 local Set up background errors for regional application WRF-Var dry control variables are:psi, chi_u, t_u and psfc Humidity control variable is q/qsg Using the averaged regression coefficients for unbalanced part Vertical truncation for psi = 15( 99.00%) Vertical truncation for chi_u = 20( 99.00%) Vertical truncation for t_u = 29( 99.00%) Vertical truncation for rh = 22( 99.00%) Scaling: var, len, ds: 0.100000E+01 0.100000E+01 0.600000E+05 Scaling: var, len, ds: 0.100000E+01 0.100000E+01 0.600000E+05 Scaling: var, len, ds: 0.100000E+01 0.100000E+01 0.600000E+05 Scaling: var, len, ds: 0.100000E+01 0.100000E+01 0.600000E+05 Scaling: var, len, ds: 0.100000E+01 0.100000E+01 0.600000E+05 VARBC: Reading VARBC.in file VARBC: no matching for platform/satid/sensor 10 2 3 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 3 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 5 --> 0.20 VARBC: Cold-starting noaa-15-amsua 6 --> -1.10 VARBC: Cold-starting noaa-15-amsua 7 --> -0.20 VARBC: Cold-starting noaa-15-amsua 8 --> -0.40 VARBC: Cold-starting noaa-15-amsua 9 --> -0.90 VARBC: Cold-starting noaa-16-amsua 5 --> -0.10 VARBC: Cold-starting noaa-16-amsua 6 --> -0.50 VARBC: Cold-starting noaa-16-amsua 7 --> -0.90 VARBC: Cold-starting noaa-16-amsua 8 --> -0.80 VARBC: Cold-starting noaa-16-amsua 9 --> -0.70 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 noaa-18-amsua 9 --> 0.00 VARBC: Cold-starting noaa-15-amsub 3 --> 1.20 VARBC: Cold-starting noaa-15-amsub 4 --> 2.90 VARBC: Cold-starting noaa-15-amsub 5 --> 0.60 VARBC: Cold-starting noaa-16-amsub 3 --> 0.30 VARBC: Cold-starting noaa-16-amsub 4 --> -0.80 VARBC: Cold-starting noaa-16-amsub 5 --> -2.60 VARBC: Cold-starting noaa-17-amsub 3 --> -2.40 VARBC: Cold-starting noaa-17-amsub 4 --> -2.30 VARBC: Cold-starting noaa-17-amsub 5 --> -5.90 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 noaa-15-amsub VARBC: Applying bias correction for noaa-16-amsub VARBC: Applying bias correction for noaa-17-amsub VARBC: Estimate Hessian for preconditioning VARBC: 181 active observations for noaa-15-amsua channel 5 VARBC: 224 active observations for noaa-15-amsua channel 6 VARBC: 217 active observations for noaa-15-amsua channel 7 VARBC: 223 active observations for noaa-15-amsua channel 8 VARBC: 199 active observations for noaa-15-amsua channel 9 VARBC: 87 active observations for noaa-16-amsua channel 5 VARBC: 92 active observations for noaa-16-amsua channel 6 VARBC: 92 active observations for noaa-16-amsua channel 7 VARBC: 92 active observations for noaa-16-amsua channel 8 VARBC: 82 active observations for noaa-16-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 noaa-18-amsua channel 9 VARBC: 469 active observations for noaa-15-amsub channel 3 VARBC: 467 active observations for noaa-15-amsub channel 4 VARBC: 435 active observations for noaa-15-amsub channel 5 VARBC: 196 active observations for noaa-16-amsub channel 3 VARBC: 197 active observations for noaa-16-amsub channel 4 VARBC: 186 active observations for noaa-16-amsub channel 5 VARBC: 18 active observations for noaa-17-amsub channel 3 VARBC: 18 active observations for noaa-17-amsub channel 4 VARBC: 15 active observations for noaa-17-amsub channel 5 Minimize cost function using CG method Starting outer iteration : 1 Starting cost function: 2.30906223D+04, Gradient= 2.33547001D+02 For this outer iteration gradient target is: 2.33547001D+00 ---------------------------------------------------------- Iter Cost Function Gradient Step 1 2.09664575D+04 2.09873429D+02 7.78878333D-02 2 1.92411229D+04 1.90817683D+02 7.83408730D-02 3 1.79608188D+04 1.26005279D+02 7.03244057D-02 4 1.73665432D+04 8.96266468D+01 7.48583611D-02 5 1.70455918D+04 6.65528798D+01 7.99087697D-02 6 1.68505295D+04 5.34000940D+01 8.80781630D-02 7 1.67574102D+04 3.95370071D+01 6.53113728D-02 8 1.66806115D+04 3.31130659D+01 9.82587438D-02 9 1.66454627D+04 2.53684083D+01 6.41129586D-02 10 1.66207475D+04 1.88573130D+01 7.68067813D-02 11 1.66023933D+04 1.65580135D+01 1.03232828D-01 12 1.65901689D+04 1.42653113D+01 8.91742217D-02 13 1.65816374D+04 1.23395425D+01 8.38539582D-02 14 1.65737218D+04 1.13957167D+01 1.03971730D-01 15 1.65684794D+04 8.97792441D+00 8.07585155D-02 16 1.65640835D+04 8.29804597D+00 1.09053016D-01 17 1.65612184D+04 6.94985422D+00 8.32785642D-02 18 1.65588981D+04 6.28951972D+00 9.59448912D-02 19 1.65571265D+04 5.55890078D+00 8.97790793D-02 20 1.65557888D+04 4.88152139D+00 8.61620362D-02 21 1.65544483D+04 4.01269581D+00 1.13253086D-01 22 1.65534289D+04 3.65065079D+00 1.23663030D-01 23 1.65524728D+04 3.52267534D+00 1.56674707D-01 24 1.65512397D+04 3.75793036D+00 1.43819473D-01 25 1.65513915D+04 3.59450537D+00 1.17787381D-01 26 1.65485329D+04 3.80625062D+00 1.01880388D-01 27 1.65506689D+04 3.46268075D+00 7.27091839D-02 28 1.65484861D+04 2.99768249D+00 9.14178732D-02 29 1.65485819D+04 2.62355041D+00 1.18454146D-01 30 1.65478297D+04 2.31392756D+00 1.18278116D-01 ---------------------------------------------------------- Inner iteration stopped after 30 iterations Final: 30 iter, J= 1.73394853D+04, g= 2.31392756D+00 ----------------------------------------------------------