Little_r format bogus data

Problems with data in 3DVAR.

Little_r format bogus data

Postby bhj9793 » Fri Feb 02, 2018 1:46 am

Hello, I want to use TC bogus data in Little_r format file.

So, I read Little_r format help in WRFDA online tutorial.
(http://www2.mmm.ucar.edu/wrf/users/wrfda/OnlineTutorial/Help/littler.html)

According to Little_r format help, Mandatory variables of Bogus data (FM-135) are pressure[Pa] and height[m] data.

Thus, I only input pressure and height data and executed obsproc.exe. However, bogus data wasn't writed in ob.ascii file.

Meanwhile, I included arbitrary temperature data in data record of Bogus data (FM-135) in Little_r file.

Then, bogus data was writed in ob.ascii file.

Why wasn't bogus data writed in ob.ascii file when I only input pressure and height data?

What was the problem?


- Below is my records of bogus data in Little_r format file -

Code: Select all
   
            23.80000           127.60000             Bogus observation            BOGUS DATA FROM ??????????? SOURCE    FM-135 BOGUS                                                                                 0.00000         1         0         0         0         0         F         T         F   -888888   -888888      20161003000000  92500.00000      0-888888.00000      0-888888.00000      0-888888.00000      0-888888.00000      0-888888.00000      0-888888.00000      0-888888.00000      0-888888.00000      0-888888.00000      0-888888.00000      0-888888.00000      0-888888.00000      0-888888.00000      0
  92500.00000      0      0.00000      0-888888.00000      0-888888.00000      0-888888.00000      0-888888.00000      0-888888.00000      0-888888.00000      0-888888.00000      0-888888.00000      0
-777777.00000      0-777777.00000      0-888888.00000      0-888888.00000      0-888888.00000      0-888888.00000      0-888888.00000      0-888888.00000      0-888888.00000      0-888888.00000      0
      1      0      0



- The results of obsproc.exe is as follows -

Code: Select all

 READ NAMELIST FILE: namelist.obsproc
 ------------------
 &RECORD1
 OBS_GTS_FILENAME = OBS_2016100300                                                                  ,
 OBS_ERR_FILENAME = obserr.txt                                                                      ,
 FG_FORMAT = WRF,
 GTS_FROM_NCAR_ARCHIVE =  F
 /
 &RECORD2
 TIME_WINDOW_MIN = 2016-10-02_23:00:00,
 TIME_ANALYSIS = 2016-10-03_00:00:00,
 TIME_WINDOW_MAX = 2016-10-03_01:00:00
 /
 &RECORD3
 MAX_NUMBER_OF_OBS =       400000,
 FATAL_IF_EXCEED_MAX_OBS =  T
 /
 &RECORD4
 QC_TEST_VERT_CONSISTENCY =  T,
 QC_TEST_CONVECTIVE_ADJ =  T,
 QC_TEST_ABOVE_LID =  T,
 REMOVE_ABOVE_LID =  T,
 DOMAIN_CHECK_H =  T,
 THINING_SATOB =  F,
 THINING_SSMI =  F,
 THINING_QSCAT =  F
 /
 &RECORD5
 PRINT_GTS_READ =  T,
 PRINT_GPSPW_READ =  T,
 PRINT_RECOVERP =  T,
 PRINT_DUPLICATE_LOC =  T,
 PRINT_DUPLICATE_TIME =  T,
 PRINT_RECOVERH =  T,
 PRINT_QC_VERT =  T,
 PRINT_QC_CONV =  T,
 PRINT_QC_LID =  T,
 PRINT_UNCOMPLETE =  T
 /
 &RECORD6
 PTOP =    1000.000    ,
 PS0 =    100000.0    ,
 TS0 =    290.0000    ,
 TLP =    50.00000    ,
 PIS0 =    20000.00    ,
 TIS0 =    215.0000    ,
 BASE_PRES =    100000.0    ,
 BASE_TEMP =    290.0000    ,
 BASE_LAPSE =    50.00000    ,
 BASE_STRAT_TEMP =    215.0000    ,
 BASE_TROPO_PRES =    20000.00
 /
 &RECORD7
 IPROJ =            3,
 PHIC =    24.80000    ,
 XLONC =    127.2000    ,
 TRUELAT1 =    30.00000    ,
 TRUELAT2 =    60.00000    ,
 MOAD_CEN_LAT =    24.80000    ,
 STANDARD_LON =    127.2000
 /
 &RECORD8
 IDD =            1,
 MAXNES =            1,
 NESTIX =          100,
                   200,
                   136,
                   181,
                   211,
                     0,
                     0,
                     0,
                     0,
                     0,
 NESTJX =          100,
                   200,
                   181,
                   196,
                   211,
                     0,
                     0,
                     0,
                     0,
                     0,
 DIS =    10.00000    ,
          10.00000    ,
          3.300000    ,
          1.100000    ,
          1.100000    ,
          0.000000    ,
          0.000000    ,
          0.000000    ,
          0.000000    ,
          0.000000    ,
 NUMC =            1,
                   1,
                   2,
                   3,
                   4,
                   0,
                   0,
                   0,
                   0,
                   0,
 NESTI =            1,
                   40,
                   28,
                   35,
                   45,
                   45,
                    0,
                    0,
                    0,
                    0,
                    0,
 NESTJ =            1,
                   60,
                   25,
                   65,
                   55,
                    0,
                    0,
                    0,
                    0,
                    0
 /
 &RECORD9
 PREPBUFR_OUTPUT_FILENAME = obs_3dvar                                                                       ,
 PREPBUFR_TABLE_FILENAME = obs_3dvar_table                                                                 ,
 OUTPUT_OB_FORMAT =            2,
 USE_FOR = 3DVAR,
 NUM_SLOTS_PAST =            3,
 NUM_SLOTS_AHEAD =            3,
 WRITE_SYNOP =  T,
 WRITE_SHIP =  T,
 WRITE_METAR =  T,
 WRITE_BUOY =  T,
 WRITE_PILOT =  T,
 WRITE_SOUND =  T,
 WRITE_AMDAR =  T,
 WRITE_SATEM =  T,
 WRITE_SATOB =  T,
 WRITE_AIREP =  T,
 WRITE_GPSPW =  T,
 WRITE_GPSZTD =  T,
 WRITE_GPSREF =  T,
 WRITE_GPSEPH =  T,
 WRITE_SSMT1 =  T,
 WRITE_SSMT2 =  T,
 WRITE_SSMI =  T,
 WRITE_TOVS =  T,
 WRITE_QSCAT =  T,
 WRITE_PROFL =  T,
 WRITE_BOGUS =  T,
 WRITE_AIRS =  T,
 WRITE_TAMDAR =  T
 /
 &RECORD10
 WIND_SD =  F,
 WIND_SD_SYNOP =  F,
 WIND_SD_SHIPS =  F,
 WIND_SD_METAR =  F,
 WIND_SD_BUOY =  F,
 WIND_SD_SOUND =  F,
 WIND_SD_QSCAT =  F,
 WIND_SD_PILOT =  F,
 WIND_SD_AIREP =  F,
 WIND_SD_GEOAMV =  F,
 WIND_SD_TAMDAR =  F,
 WIND_SD_PROFILER =  F
 /

               === 3DVAR_OBSPROC is used for WRF 3DVAR ===
          {number of slots = 1,  length of the full-slot = 7200 sec.}

Setting up MERCATOR map...
 map_info:            1            3
 code,lat1,lon1,dx,dlat,dlon,stdlon,truelat1,truelat2,hemi,cone,polei,polej,rsw,rebydx,knowni,knownj,init
            1    24.80000        127.2000        10000.00        0.000000
   1.8127168E-03    127.2000        30.00000        60.00000
    1.000000       -999.9000       -999.9000       -999.9000
    246.6061        637.0000        50.00000        50.00000      T

MERCATOR PROJECTION IS ONLY SUPPORTED AT 0 DEGREE TRUE LATITUDE.
TRUELAT1 IS RESET TO 0

 COARSE GRID CENTER IS AT X =      0.0 KM AND Y =   2847.6 KM.   DIS(1)=    10.000

         LLXY:lat, lon, (dot) xxi, yyj:    24.80000   127.20000    50.50000    50.50000
latlon_to_ij:lat, lon, (cross)xxi, yyj:    24.80000   127.20000    50.00000    50.00000

------------------------------------------------------------------------------
READ GTS OBSERVATIONS IN FILE OBS_2016100300

Diagnostics in file obs_gts_read.diag

 ==> obs_num=    1607 INVALID WMO CODE: platform=             ID=
Troubles with first line
Have reached end of observations file.

------------------------------------------------------------------------------
GTS OBSERVATIONS READ:

 SYNOP reports: 63963
 SHIPS reports:  4668
 BUOYS reports: 40315
 BOGUS reports:     1
 METAR reports:136692
 PILOT reports:   108
 SOUND reports:  1399
 AMDAR reports:     0
 SATEM reports:     0
 SATOB reports:391530
 GPSPW reports:   462
 GPSZD reports:     0
 GPSRF reports:     0
 GPSEP reports:     0
 AIREP reports:405349
TAMDAR reports:     0
 SSMT1 reports:     0
 SSMT2 reports:     0
 SSMI  reports:     0
 TOVS  reports:     0
 QSCAT reports:240034
 PROFL reports:  2093
 AIRST reports:     0
 OTHER reports:    56
 Total reports:******

Number of observations read:           1286670
Number of empty observations:               15
Number of observations out of domain:  1285049                                  Number of observations for ingestion:     1606
 
 
AIRCRAFT DATA: Total=     3  Above cut_height=     3
------------------------------------------------------------------------------
PRESSURE RECOVERED FROM HEIGHTS:
 
Diagnostics in file obs_recover_pressure.diag
 
------------------------------------------------------------------------------
CHECK PRESENCE OF HEIGHT OR/AND PRESSURE
 
Diagnostics in file obs_check_pressure.diag
 
------------------------------------------------------------------------------
REMOVE DUPLICATE STATIONS BY LOCATION:
 
Diagnostics in file obs_duplicate_loc.diag
 
Found   524 location duplicate stations that have been merged.
 
------------------------------------------------------------------------------
REMOVE DUPLICATE STATIONS BY TIME:
 
Diagnostics in file obs_duplicate_time.diag_2016-10-03_00:00:00
 
Found   323 time duplicate stations that have been removed.
 
------------------------------------------------------------------------------
DIAGNOSTICS: U, V, RH, QV etc...:
 
------------------------------------------------------------------------------
HEIGHT RECOVERED FROM P, T, Q,..:
 
Diagnostics in file obs_recover_height.diag
 
------------------------------------------------------------------------------
CHECK PRESENCE OF HEIGHT OR/AND PRESSURE
 
Diagnostics in file obs_check_height.diag
 
------------------------------------------------------------------------------
 
<AFWA> OBSERVATIONAL ERRORS PER TYPE:
 
HEIGHT SENSOR ERRORS:  RAOBS, PIBALS, AIREPS, TAMDARS, SURFACE SHIP, BUOY, DMSP - A RETRIEVAL,
 DMSP - A RETRIEVAL,
 NOAA - CLEAR PATH RETRIEVAL,
 NOAA - A RETRIEVAL, PROFL, AIRSRET,
 
UV SENSOR ERRORS:  RAOBS, PIBALS, AIREPS, TAMDARS, US LOW LEVEL WINDS, PROFL, BUOY,

SPD SENSOR ERRORS:  RAOBS, PIBALS, AIREPS, TAMDARS, US LOW LEVEL WINDS, PROFL, BUOY,
   
DIR SENSOR ERRORS:  RAOBS, PIBALS, TAMDARS, AIREPS, US LOW LEVEL WINDS, PROFL, BUOY,

TEMP SENSOR ERRORS:  RAOBS, PIBALS, AIREPS, TAMDARS, SURFACE LAND, SURFACE SHIP, BUOY, DMSP - A RETRIEVAL, DMSP - A RETRIEVAL, NOAA - A RETRIEVAL, AIRSRET,

RH SENSOR ERRORS ( % ):  RAOBS, AIREPS, TAMDARS, PIBALS, DMSP - A RETRIEVAL, NOAA - A RETRIEVAL, SURFACE LAND, AIRSRET,   

PRESSURE SENSOR ERRORS:  RAOBS, PIBALS, AIREPS, TAMDARS, SURFACE LAND, SURFACE SHIP, DMSP - A RETRIEVAL, NOAA - A RETRIEVAL, AIRSRET,

Number of processed stations:              758 =   1080 levels.
Number of processed surface stations:      731 =    731 surface levels.
Number of processed upper-air stations:     27 =    349 upper-air levels.

------------------------------------------------------------------------------
VERTICAL CONSISTENCY TEST QC:

Diagnostics in file obs_qc1.diag

------------------------------------------------------------------------------
CONVECTIVE ADJUSTEMENT TEST QC:

Diagnostics in file obs_qc1.diag
 
------------------------------------------------------------------------------
CHECK OBS HEIGHT...:
 
Diagnostics in file obs_qc2.diag
 
------------------------------------------------------------------------------
LOOK FOR UNCOMPLETE DATA:
 
Diagnostics in file obs_uncomplete.diag
Remove      0 surface stations.
Remove     47 upper-air levels.

------------------------------------------------------------------------------
INGESTED OBSERVATION AFTER CHECKS:

                 READ   EMPTY OUTSIDE LOCDUPL TIMDUPL UNCOMPL INGESTD
 SYNOP reports: 63963       0   63870      31      40       0      22
 SHIPS reports:  4668       0    4662       2       0       0       4
 BUOYS reports: 40315      15   40300       0       0       0       0
 BOGUS reports:     1       0       0       0       0       0       1
 METAR reports:136692       0  136563      44      71       0      14
 PILOT reports:   108       0     108       0       0       0       0
 SOUND reports:  1399       0    1393       3       0       0       3
 AMDAR reports:     0       0       0       0       0       0       0
 SATEM reports:     0       0       0       0       0       0       0
 SATOB reports:391530       0  390410     441     212       0     467
 AIREP reports:405349       0  405340       3       0       0       6
 GPSPW reports:   462       0     462       0       0       0       0
 GPSRF reports:     0       0       0       0       0       0       0
 GPSEP reports:     0       0       0       0       0       0       0
 SSMT1 reports:     0       0       0       0       0       0       0
 SSMT2 reports:     0       0       0       0       0       0       0
 SSMI  reports:     0       0       0       0       0       0       0
 TOVS  reports:     0       0       0       0       0       0       0
 QSCAT reports:240034       0  239792       0       0       0     242
 RPOFL reports:  2093       0    2093       0       0       0       0
 AIRS  reports:     0       0       0       0       0       0       0
TAMDAR reports:     0       0       0       0       0       0       0
 OTHER reports:    56       0      56       0       0       0       0
 Total reports:******      15  ******     524     323       0     759

------------------------------------------------------------------------------
REDUCE QC FROM 7 TO 2 DIGITS:



slot= 1  time_min, time_fg, time_max:  2016-10-02_23:00:00  2016-10-03_00:00:00  2016-10-03_00:59:59
------------------------------------------------------------------------------
OBSERVATIONS FOR OUTPUT SLOT 01

 SYNOP reports:    22
 SHIPS reports:     4
 BUOY  reports:     0
 BUGUS reports:     0
 METAR reports:    14
 PILOT reports:     0
 SOUND reports:     3
 AMDAR reports:     0
 SATEM reports:     0
 SATOB reports:   467
 AIREP reports:     6
TAMDAR reports:     0
 GPSPW reports:     0
 GPSZD reports:     0
 GPSRF reports:     0
 GPSEP reports:     0
 SSMT1 reports:     0
 SSMT2 reports:     0
 SSMI  reports:     0
 TOVS  reports:     0
 QSCAT reports:     0
 PROFL reports:     0
 AIRST reports:     0
 OTHER reports:     0
 Total reports:   516 =    489 single +     27 multi-level reports.

------------------------------------------------------------------------------
Write 3DVAR 3DVAR GTS observations in file obs_gts_2016-10-03_00:00:00.3DVAR (wrfvar version 3.0)

Wrote    1307 lines of data in file: obs_gts_2016-10-03_00:00:00.3DVAR

No SSMI observations available.

Warning: ieee_divide_by_zero is signaling
Warning: ieee_inexact is signaling
99999
bhj9793
 
Posts: 1
Joined: Sun Aug 20, 2017 1:34 pm

Return to Data Issues

Who is online

Users browsing this forum: No registered users and 4 guests