NOTICE: This version of the NSF Unidata web site (archive.unidata.ucar.edu) is no longer being updated.
Current content can be found at unidata.ucar.edu.
To learn about what's going on, see About the Archive Site.
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 >