This archive contains answers to questions sent to Unidata support through mid-2025. Note that the archive is no longer being updated. We provide the archive for reference; many of the answers presented here remain technically correct, even if somewhat outdated. For the most up-to-date information on the use of NSF Unidata software and data services, please consult the Software Documentation first.
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 CHiswellFrom: =?ISO-8859-1?Q?Christian_Pag=E9?= <address@hidden> Organization: UCAR/Unidata Keywords: 200211151501.gAFF1KL26527Hi 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 gridsInterpolating HGHT 021115/0000F000 500 -1 PRESSo 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 HGTDParameters 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