[ungrib] data not read in when data not in sequence

Post here if you think you might have found any type of bug

[ungrib] data not read in when data not in sequence

Postby packard » Fri Oct 24, 2014 12:51 pm

I am using WPS3.6.1, and ECMWF YOTC surface analysis (CISL rda 629.1, grib1, subsetted for 2009/09/13 00UTC~2009/09/17 00UTC). When I run the ungrib program (for 13_00utc ~ 17_00utc), the sea ice (CI) and SST data for 13_00utc~14_00utc and 17_00utc is not read in.

Part of the ungrib.log is shown here:
Code: Select all
2014-10-24 22:13:24.601 --- Inventory for date = 2009-09-17 00:00:00
2014-10-24 22:13:24.601 ---      PRES     GEOPT       HGT        TT        UU        VV        RH     DEWPT   LANDSEA   SOILGEO   SOILHGT      PSFC      PMSL  SKINTEMP    SEAICE       SST   SNOW_EC      SNOW  ST000007  ST007028  ST028100  ST100255  SM000007  SM007028  SM028100  SM100255
2014-10-24 22:13:24.601 --- -------------------------------------------------
2014-10-24 22:13:24.601 ---    200100         O         O         X         O         O         O         O         O         O         O         O         O         O         O         O         O         O         O         O         O         O         O         O         O         O
2014-10-24 22:13:24.602 --- -------------------------------------------------
2014-10-24 22:13:24.765 --- Inventory for date = 2009-09-14 12:00:00
2014-10-24 22:13:24.765 ---      PRES     GEOPT       HGT        TT        UU        VV        RH     DEWPT   LANDSEA   SOILGEO   SOILHGT      PSFC      PMSL  SKINTEMP    SEAICE       SST   SNOW_EC      SNOW  ST000007  ST007028  ST028100  ST100255  SM000007  SM007028  SM028100  SM100255
2014-10-24 22:13:24.765 --- -------------------------------------------------
2014-10-24 22:13:24.765 ---    200100         O         O         O         O         O         O         O         O         O         O         O         O         O         X         O         O         O         O         O         O         O         O         O         O         O
2014-10-24 22:13:24.766 --- -------------------------------------------------
2014-10-24 22:13:24.803 --- Inventory for date = 2009-09-14 18:00:00
2014-10-24 22:13:24.803 ---      PRES     GEOPT       HGT        TT        UU        VV        RH     DEWPT   LANDSEA   SOILGEO   SOILHGT      PSFC      PMSL  SKINTEMP    SEAICE       SST   SNOW_EC      SNOW  ST000007  ST007028  ST028100  ST100255  SM000007  SM007028  SM028100  SM100255
2014-10-24 22:13:24.804 --- -------------------------------------------------
2014-10-24 22:13:24.804 ---    200100         O         O         O         O         O         O         O         O         O         O         O         O         O         X         O         O         O         O         O         O         O         O         O         O         O
2014-10-24 22:13:24.804 --- -------------------------------------------------
2014-10-24 22:13:24.842 --- Inventory for date = 2009-09-15 00:00:00
2014-10-24 22:13:24.842 ---      PRES     GEOPT       HGT        TT        UU        VV        RH     DEWPT   LANDSEA   SOILGEO   SOILHGT      PSFC      PMSL  SKINTEMP    SEAICE       SST   SNOW_EC      SNOW  ST000007  ST007028  ST028100  ST100255  SM000007  SM007028  SM028100  SM100255
2014-10-24 22:13:24.843 --- -------------------------------------------------
2014-10-24 22:13:24.843 ---    200100         O         O         O         O         O         O         O         O         O         O         O         O         O         X         O         O         O         O         O         O         O         O         O         O         O
2014-10-24 22:13:24.843 --- -------------------------------------------------
2014-10-24 22:13:24.880 --- Inventory for date = 2009-09-15 06:00:00
2014-10-24 22:13:24.880 ---      PRES     GEOPT       HGT        TT        UU        VV        RH     DEWPT   LANDSEA   SOILGEO   SOILHGT      PSFC      PMSL  SKINTEMP    SEAICE       SST   SNOW_EC      SNOW  ST000007  ST007028  ST028100  ST100255  SM000007  SM007028  SM028100  SM100255
2014-10-24 22:13:24.881 --- -------------------------------------------------
2014-10-24 22:13:24.881 ---    200100         O         O         O         O         O         O         O         O         O         O         O         O         O         X         O         O         O         O         O         O         O         O         O         O         O
2014-10-24 22:13:24.881 --- -------------------------------------------------
2014-10-24 22:13:24.919 --- Inventory for date = 2009-09-15 12:00:00
2014-10-24 22:13:24.919 ---      PRES     GEOPT       HGT        TT        UU        VV        RH     DEWPT   LANDSEA   SOILGEO   SOILHGT      PSFC      PMSL  SKINTEMP    SEAICE       SST   SNOW_EC      SNOW  ST000007  ST007028  ST028100  ST100255  SM000007  SM007028  SM028100  SM100255
2014-10-24 22:13:24.919 --- -------------------------------------------------
2014-10-24 22:13:24.919 ---    200100         O         O         O         O         O         O         O         O         O         O         O         O         O         X         O         O         O         O         O         O         O         O         O         O         O
2014-10-24 22:13:24.920 --- -------------------------------------------------
2014-10-24 22:13:24.957 --- Inventory for date = 2009-09-15 18:00:00
2014-10-24 22:13:24.957 ---      PRES     GEOPT       HGT        TT        UU        VV        RH     DEWPT   LANDSEA   SOILGEO   SOILHGT      PSFC      PMSL  SKINTEMP    SEAICE       SST   SNOW_EC      SNOW  ST000007  ST007028  ST028100  ST100255  SM000007  SM007028  SM028100  SM100255
2014-10-24 22:13:24.958 --- -------------------------------------------------
2014-10-24 22:13:24.958 ---    200100         O         O         O         O         O         O         O         O         O         O         O         O         O         X         O         O         O         O         O         O         O         O         O         O         O
2014-10-24 22:13:24.958 --- -------------------------------------------------
2014-10-24 22:13:24.996 --- Inventory for date = 2009-09-16 00:00:00
2014-10-24 22:13:24.996 ---      PRES     GEOPT       HGT        TT        UU        VV        RH     DEWPT   LANDSEA   SOILGEO   SOILHGT      PSFC      PMSL  SKINTEMP    SEAICE       SST   SNOW_EC      SNOW  ST000007  ST007028  ST028100  ST100255  SM000007  SM007028  SM028100  SM100255
2014-10-24 22:13:24.997 --- -------------------------------------------------
2014-10-24 22:13:24.997 ---    200100         O         O         O         O         O         O         O         O         O         O         O         O         O         X         O         O         O         O         O         O         O         O         O         O         O
2014-10-24 22:13:24.997 --- -------------------------------------------------
2014-10-24 22:13:25.034 --- Inventory for date = 2009-09-16 06:00:00
2014-10-24 22:13:25.034 ---      PRES     GEOPT       HGT        TT        UU        VV        RH     DEWPT   LANDSEA   SOILGEO   SOILHGT      PSFC      PMSL  SKINTEMP    SEAICE       SST   SNOW_EC      SNOW  ST000007  ST007028  ST028100  ST100255  SM000007  SM007028  SM028100  SM100255
2014-10-24 22:13:25.034 --- -------------------------------------------------
2014-10-24 22:13:25.035 ---    200100         O         O         O         O         O         O         O         O         O         O         O         O         O         X         O         O         O         O         O         O         O         O         O         O         O
2014-10-24 22:13:25.035 --- -------------------------------------------------
2014-10-24 22:13:25.073 --- Inventory for date = 2009-09-16 12:00:00
2014-10-24 22:13:25.073 ---      PRES     GEOPT       HGT        TT        UU        VV        RH     DEWPT   LANDSEA   SOILGEO   SOILHGT      PSFC      PMSL  SKINTEMP    SEAICE       SST   SNOW_EC      SNOW  ST000007  ST007028  ST028100  ST100255  SM000007  SM007028  SM028100  SM100255
2014-10-24 22:13:25.073 --- -------------------------------------------------
2014-10-24 22:13:25.073 ---    200100         O         O         O         O         O         O         O         O         O         O         O         O         O         X         O         O         O         O         O         O         O         O         O         O         O
2014-10-24 22:13:25.074 --- -------------------------------------------------
2014-10-24 22:13:25.111 --- Inventory for date = 2009-09-16 18:00:00
2014-10-24 22:13:25.111 ---      PRES     GEOPT       HGT        TT        UU        VV        RH     DEWPT   LANDSEA   SOILGEO   SOILHGT      PSFC      PMSL  SKINTEMP    SEAICE       SST   SNOW_EC      SNOW  ST000007  ST007028  ST028100  ST100255  SM000007  SM007028  SM028100  SM100255
2014-10-24 22:13:25.112 --- -------------------------------------------------
2014-10-24 22:13:25.112 ---    200100         O         O         O         O         O         O         O         O         O         O         O         O         O         X         O         O         O         O         O         O         O         O         O         O         O
2014-10-24 22:13:25.112 --- -------------------------------------------------
2014-10-24 22:13:25.277 --- Inventory for date = 2009-09-14 06:00:00
2014-10-24 22:13:25.278 ---      PRES     GEOPT       HGT        TT        UU        VV        RH     DEWPT   LANDSEA   SOILGEO   SOILHGT      PSFC      PMSL  SKINTEMP    SEAICE       SST   SNOW_EC      SNOW  ST000007  ST007028  ST028100  ST100255  SM000007  SM007028  SM028100  SM100255
2014-10-24 22:13:25.278 --- -------------------------------------------------
2014-10-24 22:13:25.278 ---    200100         O         O         O         O         O         O         O         O         O         O         O         O         O         X         O         O         O         O         O         O         O         O         O         O         O
2014-10-24 22:13:25.278 --- -------------------------------------------------
2014-10-24 22:13:25.399 --- Inventory for date = 2009-09-13 00:00:00
2014-10-24 22:13:25.399 ---      PRES     GEOPT       HGT        TT        UU        VV        RH     DEWPT   LANDSEA   SOILGEO   SOILHGT      PSFC      PMSL  SKINTEMP    SEAICE       SST   SNOW_EC      SNOW  ST000007  ST007028  ST028100  ST100255  SM000007  SM007028  SM028100  SM100255
2014-10-24 22:13:25.399 --- -------------------------------------------------
2014-10-24 22:13:25.399 ---    200100         O         O         O         O         O         O         O         X         O         O         O         O         O         O         O         O         O         O         O         O         O         O         O         O         O


Meanwhile, other sfc files seem to be correctly read in. Then I found that the sea-ice file is different from other sfc files that they are not in sequence according to util/g1print.exe
Code: Select all
[packard@jn ecmwf]$ ../WRF/WPS35/util/g1print.exe fc/yt.oper.an.sfc.regn400sc.ci_31.200909.chan93941
Copen: File = fc/yt.oper.an.sfc.regn400sc.ci_31.200909.chan93941                                                                     
Fortran Unit = 0
UNIX File descriptor: 3

----------------------------------------------------
 rec GRIB GRIB  Lvl  Lvl  Lvl         Time      Fcst
 Num Code name  Code one  two                   hour
----------------------------------------------------
   1  31 CI       1    0    0  2009-09-14_12:00 + 00
   2  31 CI       1    0    0  2009-09-14_18:00 + 00
   3  31 CI       1    0    0  2009-09-15_00:00 + 00
   4  31 CI       1    0    0  2009-09-15_06:00 + 00
   5  31 CI       1    0    0  2009-09-15_12:00 + 00
   6  31 CI       1    0    0  2009-09-15_18:00 + 00
   7  31 CI       1    0    0  2009-09-16_00:00 + 00
   8  31 CI       1    0    0  2009-09-16_06:00 + 00
   9  31 CI       1    0    0  2009-09-16_12:00 + 00
  10  31 CI       1    0    0  2009-09-16_18:00 + 00
  11  31 CI       1    0    0  2009-09-17_00:00 + 00
  12  31 CI       1    0    0  2009-09-13_00:00 + 00
  13  31 CI       1    0    0  2009-09-13_06:00 + 00
  14  31 CI       1    0    0  2009-09-13_12:00 + 00
  15  31 CI       1    0    0  2009-09-13_18:00 + 00
  16  31 CI       1    0    0  2009-09-14_00:00 + 00
  17  31 CI       1    0    0  2009-09-14_06:00 + 00

***** End-Of-File on C unit   3


So, I assume that when data are not recorded in time sequence in grib files, some data may be neglected by ungrib program, which I consider to be a bug. And I am forced to use the original monthly grib file (not subsetted), in which data are in sequence, and things are alright again.

Hope someone will check the source code and give some solution to this bug.

As a good practice, the namelist.wps is given here:
Code: Select all
&share
 wrf_core = 'ARW',
 max_dom = 4,
 start_date = '2009-09-13_00:00:00','2009-09-13_00:00:00','2009-09-13_00:00:00','2009-09-13_00:00:00',
 end_date   = '2009-09-17_00:00:00','2009-09-13_00:00:00','2009-09-13_00:00:00','2009-09-13_00:00:00',
 interval_seconds = 21600
 io_form_geogrid = 2,
/

&geogrid
 parent_id         =     1,    1,    2,    3,
 parent_grid_ratio =     1,    3,    3,    3,
 i_parent_start    =     1,   41,   33,   37,   70,
 j_parent_start    =     1,   29,   35,   37,   70,
 e_we              =   101,  109,  127,  163,  121,
 e_sn              =   101,  109,  127,  163,  121,
 geog_data_res     = '30s','30s','30s','30s','30s',
 dx = 67500,
 dy = 67500,
 map_proj = 'lambert',
 ref_lat   =  18.0,
 ref_lon   = 145.0,
 truelat1  =  10.0,
 truelat2  =  25.0,
 stand_lon = 145.0,
!geog_data_path = '/mnt/node71data1/NWP2013/geog'
 geog_data_path = '/lustre/jiangyx/geog'
/

&ungrib
 out_format = 'WPS',
 prefix = 'FILE',
/

&metgrid
 fg_name = 'FILE'
 io_form_metgrid = 2,
/
packard
 
Posts: 6
Joined: Fri Apr 25, 2014 10:33 pm

Return to Bugs

Who is online

Users browsing this forum: No registered users and 1 guest