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