Re: [gembud] viewing reanalysis data in nmap2. . .

Thanks for the response, Kevin.  Probably part of the problem is that
I'm having a hard time finding what all of the fields represent in the
various table files.

Here is my entry for datatype.tbl:
NNRPPGB      $GEMDATA/model/nnrppgb    YYYYMMDD_nnrppgb.gem
CAT_GRD  SCAT_ANL  124  44640    360      OFF/0/0      4

I've tried changing the file name to just a non-date defined name (ie,
nnrppgb.gem), I've tried changing the subcategory to SCAT_FCT, I've
tried changing the other columns - nothing has worked so far.  I assume
#FRM is the number of frames in the file (or maximum # of frames?).  In
this case it would be 124 (31 days times 4x per day).  I'm guessing the
RANGE is the time span that the file covers in minutes (31 days times
1440 minutes per day).  I assume INTRVL is the data interval in minutes
(6 hours times 60 minutes per hour).  I don't know what BIN HRS refers
to, but they all look the same, and it looks like TIME MATCH shouldn't
matter in this case.

Here's an entry for mod_res.tbl for plotting PMSL:
NNRPPGB_PMSL                     $NMAP_RESTORE/IBR/nnrppgb/nnrppgb_pmsl
A0.nnrppgb       nnrppgb

That looks pretty straightforward - I think.

And here is the current restore file that I've been experimenting with
for PMSL (I've tried tweaking this as well):
!---------------------------------------------------------------------
!
! Restore File : PMSL
!
GDATTIM all
GLEVEL  0
GVCORD  none
PANEL   0        
SKIP    0
SCALE   0
GDPFUN  pmsl
TYPE    c
CONTUR  1        
CINT    1
LINE    1
FINT    
FLINE   
HILO    
HLSYM        
CLRBAR  1        
WIND    
REFVEC      
TITLE   5/-1/~ @ PMSL OA!0
TEXT    1/21//hw        
CLEAR   yes        
STNPLT        
SATFIL        
RADFIL        
STREAM        
POSN    4        
COLORS  2        
MARKER  2        
GRDLBL  5        
LUTFIL  none
FILTER  yes
PROJ    mer
GAREA   or
!---------------------------------------------------------------------

Here's what happens when I try to load the reanalysis data into nmap2.
I select GRID in the top level scroll list in the "Data Source" window.
Then I select nnrppgb in the next scroll list (titled GRID).  But then
in the next scroll list (titled nnrppgb) I get the listing of all the
times which are contained within my single gem file - not just the first
time as I would like to see.

I can select any of the times, and the same thing happens - I'm
presented with my GROUP (A0.nnrppgb).  I select that and get the name
for the field I've defined (NNRPPGB_PMSL).  I select that and click on
"Accept" - and nothing.  The pane that should show the available times
that I can select to plot is empty.

It really seems like it wants to treat each analysis time point
separately - even though they are all in a single file.  If I could
trick it into thinking that the first time was an analysis, and the rest
were forecasts then I think it might work.

Kevin, are you sure that the data in your file is all labeled with time
tags indicating that they are analysis data (like mine below)?  If so,
I'm not sure why mine isn't working.  Does anybody know how I might be
able to trick gempak into thinking my analysis data is a forecast?  Or
another way to make this work?

Thanks,
Mike



-----Original Message-----
From: Kevin R. Tyle [mailto:ktyle@xxxxxxxxxxxxxxxx] 
Sent: Monday, October 19, 2009 6:41 PM
To: Zulauf, Michael
Cc: gembud@xxxxxxxxxxxxxxxx
Subject: Re: [gembud] viewing reanalysis data in nmap2. . .

Hi Michael,

What is your table entry for this data set in the datatypes.tbl file,
and 
what did you name the GEMPAK file?

You may need to change the default #frames, default range, and default 
interval columns for this data type.  For example, this is what we have 
for a loop of GFS F000 grids that are cover a period of a month but live

in a single file:

FNL          /data1/gfsanl             gfsanl.gem        CAT_GRD
SCAT_FCT   30  44640    720      OFF/0/0      4

--Kevin 
______________________________________________________________________
Kevin Tyle, Systems Administrator               **********************
Dept. of Atmospheric & Environmental Sciences   ktyle@xxxxxxxxxxxxxxxx
University at Albany, ES-235                    518-442-4578 (voice)
1400 Washington Avenue                          518-442-5825 (fax)
Albany, NY 12222                                **********************
______________________________________________________________________

On Mon, 19 Oct 2009, Zulauf, Michael wrote:

> Hi all. . .
>
> I'm having a bit of a problem successfully viewing global reanalysis
> data (NNRP) using nmap2.  I believe the problem is that I'm not
> correctly handling the fact that the data is a series of analyses
> (f000), rather than forecasts.
>
> Here's what I'm doing.  First, I'm subsetting the original pgb grib
> files using wgrib, so I'm creating grib files with just the variables
I
> want to convert.  Then I use dcgrib2 to convert these into monthly
> gempak files with the fields I want.  It appears to me that the time
> tags are getting entered correctly - at least when I use something
like
> gdinfo.  But nothing I do is letting me visualize the data in nmap2.
>
> Here's an example of the output from gdinfo for PMSL, from a file
> containing data for January 2009:
>  NUM       TIME1              TIME2           LEVL1 LEVL2  VCORD PARM
>   28     090101/0000F000                          0         NONE PMSL
>
>   60     090101/0600F000                          0         NONE PMSL
>
>   92     090101/1200F000                          0         NONE PMSL
>
>  124     090101/1800F000                          0         NONE PMSL
>
>  156     090102/0000F000                          0         NONE PMSL
>
>  188     090102/0600F000                          0         NONE PMSL
>
>  220     090102/1200F000                          0         NONE PMSL
>
>  252     090102/1800F000                          0         NONE PMSL
>
>  284     090103/0000F000                          0         NONE PMSL
>
>  316     090103/0600F000                          0         NONE PMSL
>
>  348     090103/1200F000                          0         NONE PMSL
>
>  380     090103/1800F000                          0         NONE PMSL
>
>  412     090104/0000F000                          0         NONE PMSL
>
>  444     090104/0600F000                          0         NONE PMSL
>
>  476     090104/1200F000                          0         NONE PMSL
>
>  508     090104/1800F000                          0         NONE PMSL
>
> (etc etc)
>
> We're editing the table files for nmap2, and what is currently
happening
> is that _all_ of the times show up in the "Data Source" dialog window,
> as if they were separate "forecasts" in separate gem files.  But if I
> select any of them, I can't get any of them to display - and in
reality
> there is only a single gem file that contains all those times.  And I
> can successfully plot the data using gdcntr/gdplot/etc.
>
> I'm thinking maybe what I should do is somehow convince nmap2 that
each
> month is a single "forecast" (even though they are really all
analyses),
> with the first time being labeled as f000, 6 hours later as f006, 12
> hours later as f012, etc, etc.
> 
> But I don't know how to go about doing that - or even if that is the
> best approach.  Is there any utility that will let me change the time
> tags in a gem file?
>
> Any thoughts?
>
> Thanks,
> Mike
>
> -- 
> Mike Zulauf
> Meteorologist
> Wind Asset Management
> Iberdrola Renewables
> 1125 NW Couch, Suite 700
> Portland, OR 97209
> Office: 503-478-6304  Cell: 503-913-0403
>
>
>
>
> This message is intended for the exclusive attention of the
address(es) indicated.  Any information contained herein is strictly
confidential and privileged, especially as regards person data,
> which must not be disclosed.  If you are the intended recipient and
have received it by mistake or learn about it in any other way, please
notify us by return e-mail and delete this message from
> your computer system. Any unauthorized use, reproduction, alteration,
filing or sending of this message and/or any attached files to third
parties may lead to legal proceedings being taken. Any
> opinion expressed herein is solely that of the author(s) and does not
necessarily represent the opinion of Iberdrola. The sender does not
guarantee the integrity, speed or safety of this
> message, not accept responsibility for any possible damage arising
from the interception, incorporation of virus or any other manipulation
carried out by third parties.
>
> _______________________________________________
> gembud mailing list
> gembud@xxxxxxxxxxxxxxxx
> For list information or to unsubscribe,  visit:
http://www.unidata.ucar.edu/mailing_lists/
>