[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

20021115: oabsnd, first guess and missing grid again...

Christian,

Try taking SHOW out of the STNDEX list of parameters to
grid in OABSND. I thinking that the error message is
related to not finding a guess grid for that. If that
isn't the problem, I'll have to set up a test.

Steve CHiswell


>From: =?ISO-8859-1?Q?Christian_Pag=E9?= <address@hidden>
>Organization: UCAR/Unidata
>Keywords: 200211151501.gAFF1KL26527

>Hi again,
>
>I am trying to also use a first guess for 500 hPa geopotential height 
>barnes analysis, but oabsnd never finds the grid, while gdinfo and 
>gdlist do list the grid I ask for.
>
>Here is a summary:
>
>I run oagrid to create a grid file for 500 hPa hgtd barnes analysis:
>$GEMEXE/oagrid  << EOF
>GDFILE  = /io/httpd/gempak/grids/nam_up_t.grid
>DELTAN  = 0
>DELTAX  = 0
>DELTAY  = 0
>GRDAREA = #45.30;-73.35;11;14
>EXTEND  = 2;2;2;2
>DTAAREA = #45.30;-73.35;11;14
>SOURCE  = SN
>SNFILE  = $GEMUPA/20021115_upa.gem
>SFFILE  =
>SNPARM  = hgtd;hght
>SFPARM  =
>DATTIM  = 021115/00
>LEVELS  = 500
>MAXGRD  = 1000
>VCOORD  = PRES
>run
>
>exit
>EOF
>
>                     OBJECTIVE ANALYSIS PARAMETERS
>
>  GRIDAREA           34.30    -87.35     55.30    -60.05  (       56.30  
>-59.35)
>  DELTAN              4.200
>  DELTAX              2.100
>  DELTAY              2.100
>
>  GRID POINTS     X :  14  ( 2, 2 )        Y :  11  ( 2, 2 )
>  EXTENDAREA           30.10    -91.55     59.50    -55.85
>  DATAAREA             34.30    -87.35     56.30    -59.35
>
>
>  CALCULATIONS:
>    FILE      $GEMUPA/20021115_upa.gem
>  SOURCE    SN               34.30    -87.35     56.30    -59.35       
>NSTN   23
>    PARM      HGTD        LEVEL    0.0        DATTIM       021115/0000
>
>      COMPUTED DATA SPACING           2.912
>      UNIFORM  DATA SPACING           5.486
>      RECOMMENDED   DELTAN =  4.199   DELTAX =  2.985   DELTAY =  2.099
>
>Then I run gdbiint to get the proper model field:
>  GDFILE    Grid file                        
>$GEMMOD/2002111500_cmcreg.gem
>  GDOUTF    Output grid file                  
>/io/httpd/gempak/grids/namumodt.grid
>  GFUNC     Scalar grid                       hght
>  GLEVEL    Grid level                        500
>  GVCORD    Grid vertical coordinate          pres
>  GDATTIM   Grid date/time                    F000
>  GDNUM     Grid numbers
>
>and I get:
>Building interpolation array for    14 columns    11 rows
>  Input grid file contains          565 grids
>  Interpolating HGHT    021115/0000F000                    500    -1 PRES
>
>So far so good.
>
>Then gddiag to change hght to hgtd:
>GDFILE = /io/httpd/gempak/grids/namumodt.grid
>GDOUTF = /io/httpd/gempak/grids/namumodt.grid
>GFUNC = QUO(hght,10)
>GLEVEL = 500
>GVCORD = pres
>GDATTIM = F000
>GRDNAM = hgtd
>
>with the output:
>     TIME1             TIME2         LEVL1 LEVL2   VCORD PARM
>021115/0000F000                       500          PRES HGTD
>
>Then the oabsnd run:
>  OABSND PARAMETERS:
>
>  Output grid file:        /io/httpd/gempak/grids/nam_up_t.grid
>  Number of x, y points:    14    11
>  Grid bnds:     34.30   -87.35    55.30   -60.05 Grid xtnd:    2   2   
>2   2
>  Station spacing:           4.20
>  Guess grid file:  /io/httpd/gempak/grids/namumodt.grid
>  Guess time:  F000
>  Sounding file:           $GEMUPA/20021115_upa.gem
>  Data area:               34.30;-87.35;55.30;-60.05
>  Level parameters:        HGTD
>  QC threshold:             0.0
>  Levels:                   500
>  Vertical coordinate:     PRES
>  Station parameters:      SHOW
>  Time:                   021115/0000
>  GAMMA: 0.30   SEARCH:  20.00  NPASS:  2
>  Data extrapolation may occur.
>
>  Processing data at time 021115/0000
>  [DG -7]  Input grid ... cannot be found.
>  [OA -5]  Guess grid could not be found.
>  SNFILE    Sounding data file                $GEMUPA/20021115_upa.gem
>  GDFILE    Grid file                         
>/io/httpd/gempak/grids/nam_up_t.grid
>  SNPARM    Sounding parameter list           hgtd
>  STNDEX    Stability indices                 SHOW
>  LEVELS    Vertical levels                   500
>  VCOORD    Vertical coordinate type          PRES
>  DATTIM    Date/time                         021115/00
>  DTAAREA   Data area for OA                  dset
>  GUESS     Guess file*time                   
>/io/httpd/gempak/grids/namumodt.grid*F000
>  GAMMA     Convergence parameter             0.3
>  SEARCH    Search radius/Extrapolation       20/EX
>  NPASS     Number of passes                  2
>  QCNTL     Quality control threshold
>
>
>I tried with the full guess time with no difference.
>
>Here is the gdinfo output:
>  GDFILE   = /io/httpd/gempak/grids/namumodt.grid
>  LSTALL   = YES
>  OUTPUT   = t
>  GDATTIM  = F000
>  GLEVEL   = 500
>  GVCORD   = pres
>  GFUNC    = all
>  GEMPAK-GDINFO>r
>
>  GRID FILE: /io/httpd/gempak/grids/namumodt.grid
>
>  GRID NAVIGATION:
>      PROJECTION:          CED
>      GRID SIZE:           14  11
>      LL CORNER:              34.30    -87.35
>      UR CORNER:              55.30    -60.05
>
>  GRID ANALYSIS BLOCK:
>      ANALYSIS TYPE:        BARNES
>      DELTAN:               4.200
>      DELTAX:               2.100
>      DELTAY:               2.100
>      GRID AREA:            34.30  -87.35   55.30  -60.05
>      EXTEND AREA:          30.10  -91.55   59.50  -55.85
>      DATA AREA:            34.30  -87.35   56.30  -59.35
>
>  Number of grids in file:     2
>
>  Maximum number of grids in file:   1000
>
>   NUM       TIME1              TIME2           LEVL1 LEVL2  VCORD PARM
>     1     021115/0000F000                        500         PRES HGHT
>     2     021115/0000F000                        500         PRES HGTD
>  Parameters requested: GDFILE,LSTALL,OUTPUT,GDATTIM,GLEVEL,GVCORD,GFUNC.
>
>So the guess grid is there, but oabsnd doesn't find it. Is there a way 
>to have more verbose output from oabsnd?
>Is there something I miss to have oabsnd fail?
>
>Christian Pagé
>address@hidden
>http://meteocentre.com/    http://meteoalerte.com/
>
>Etudiant au Doctorat en Sciences de l'environnement UQAM
>+1 514 987 3000 ext. 2376
>