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

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



Hi again,

This was the error. It works very well now. Thanks!!

On Vendredi, nove 15, 2002, at 13:14 America/Montreal, Unidata Support wrote:


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


*********************************************************************** ***** < Unidata User Support UCAR Unidata Program < (303)497-8643 P.O. Box 3000 < address@hidden Boulder, CO 80307 < ----------------------------------------------------------------------- ----- < Unidata WWW Service http://www.unidata.ucar.edu/ < *********************************************************************** ***** <


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