problem reading ncep grib file

From: James D Scott (James.D.Scott AT noaa.gov)
Date: Thu May 05 2005 - 12:15:57 MDT

  • Next message: Addy Fu: "Plot a round-shaped cloud on map"

    I'm reading a grib file containing ncep ensemble gfs forecast using the
    addfile function.
    There should be multiple ensemble members (12) in this file, but ncl
    indicates there is
    only one variable (time,lat, lon). Am I doing something wrong or is
    this a limitation in the
    addfile function?

    The data file I'm working with is on ftp if anyone wants to
    look at it: ftp://ftp.cdc.noaa.gov/Public/jds/t700.2005040400

    Thanks,
    -Jamie Scott

    Here is my ncl code :
    ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    begin
    name="t700.2005050400.grb"
    finfo=addfile(name,"r")
    print(finfo)
    z=finfo->TMP_2_ISBL_10
    printVarSummary(z)
    end
    ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    Here's the output from this code:

      Copyright (C) 1995-2004 - All Rights Reserved
      University Corporation for Atmospheric Research
      NCAR Command Language Version 4.2.0.a032
      The use of this software is governed by a License Agreement.
      See http://ngwww.ucar.edu/ncl/ for more details.
    Variable: finfo (file variable)

    filename: t700.2005050400
    path: /data/map/ensemble/mrfdata/t700.2005050400
        file global attributes:
        dimensions:
           forecast_time0 = 65
           lat_2 = 73
           lon_2 = 144
        variables:
           float TMP_2_ISBL_10 ( forecast_time0, lat_2, lon_2 )
              sub_center : NCEP Ensemble Products
              center : US National Weather Service - NCEP (WMC)
              long_name : Temperature
              units : K
              _FillValue : -999
              level_indicator : 100
              grid_number : 2
              parameter_number : 11
              model : Spectral Statistical Interpolation (SSI) analysis from
    "Aviation" run.
              level : 700
              initial_time : 05/04/2005 (00:00)

           float lat_2 ( lat_2 )
              grid_description : 10512-point (73x144) global
    longitude-latitude grid. (0,0) at 0E, 90N, latitude grid. (0,0) at
    0E, 90N, matrix layout. N.B.: prime meridian not duplicated.
              units : degrees_north
              long_name : latitude

           float lon_2 ( lon_2 )
              grid_description : 10512-point (73x144) global
    longitude-latitude grid. (0,0) at 0E, 90N, latitude grid. (0,0) at
    0E, 90N, matrix layout. N.B.: prime meridian not duplicated.
              units : degrees_east
              long_name : longitude

           integer forecast_time0 ( forecast_time0 )
              long_name : Forecast offset from initial time
              units : hours

    Variable: z
    Type: float
    Total Size: 2733120 bytes
                 683280 values
    Number of Dimensions: 3
    Dimensions and sizes: [forecast_time0 | 65] x [lat_2 | 73] x [lon_2 |
    144]
    Coordinates:
                 forecast_time0: [0..384]
                 lat_2: [90..-90]
                 lon_2: [ 0..357.5]
    Number Of Attributes: 11
       sub_center : NCEP Ensemble Products
       center : US National Weather Service - NCEP (WMC)
       long_name : Temperature
       units : K
       _FillValue : -999
       level_indicator : 100
       grid_number : 2
       parameter_number : 11
       model : Spectral Statistical Interpolation (SSI) analysis from
    "Aviation" run.
       level : 700
       initial_time : 05/04/2005 (00:00)

    _______________________________________________
    ncl-talk mailing list
    ncl-talk@ucar.edu
    http://mailman.ucar.edu/mailman/listinfo/ncl-talk

    This archive was generated by hypermail 2b29 : Thu May 05 2005 - 14:35:48 MDT