[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
20030625: problem with GRDCOPY
- Subject: 20030625: problem with GRDCOPY
- Date: Wed, 25 Jun 2003 12:20:35 -0600
>From: Owen Cooper <address@hidden>
>Organization: Aeronomy Laboratory/NOAA
>Keywords: 200306251740.h5PHeMLd017765 McIDAS-X GRDCOPY
Hi Owen,
>I'm having trouble using MATH with the GRDCOPY command
>I copy a subsect of realtime AVN data to a work grid file
>
>TE AVNLEV "400 300 250
>GRDCOPY RTGRIDS/AVN MYDATA/GRIDS.1 PAR=U V Z T DAY=176 TIME=12 FHOUR=0
>GPRO=MERC LEV=#AVNLEV NUM=ALL SUBSECT=0 90 DEL=YES MAX=1065 TITLE='IPV
>work grid'
>
>This command is successful.
I replicated this sequence and was successful also.
>And I have no problem listing the copied grids:
>
>GRDLIST MYDATA/GRIDS.1 NUM=ALL
>Dataset position 1 Directory Title= IPV work grid
>PAR LEVEL DAY TIME SRC FHOUR FDAY FTIME GRID PR
> O
>---- --------- ------------ -------- ---- ----- ------------ -------- ----- --
> --
>T 400 MB 25 JUN 03176 12:00:00 AVN 0 25 JUN 03176 12:00:00 1 ME
> RC
>T 300 MB 25 JUN 03176 12:00:00 AVN 0 25 JUN 03176 12:00:00 2 ME
> RC
>T 250 MB 25 JUN 03176 12:00:00 AVN 0 25 JUN 03176 12:00:00 3 ME
> RC
>T 200 MB 25 JUN 03176 12:00:00 AVN 0 25 JUN 03176 12:00:00 4 ME
> RC
>U 400 MB 25 JUN 03176 12:00:00 AVN 0 25 JUN 03176 12:00:00 5 ME
> RC
>U 300 MB 25 JUN 03176 12:00:00 AVN 0 25 JUN 03176 12:00:00 6 ME
> RC
>U 250 MB 25 JUN 03176 12:00:00 AVN 0 25 JUN 03176 12:00:00 7 ME
> RC
>U 200 MB 25 JUN 03176 12:00:00 AVN 0 25 JUN 03176 12:00:00 8 ME
> RC
>V 400 MB 25 JUN 03176 12:00:00 AVN 0 25 JUN 03176 12:00:00 9 ME
> RC
>V 300 MB 25 JUN 03176 12:00:00 AVN 0 25 JUN 03176 12:00:00 10 ME
> RC
>V 250 MB 25 JUN 03176 12:00:00 AVN 0 25 JUN 03176 12:00:00 11 ME
> RC
>V 200 MB 25 JUN 03176 12:00:00 AVN 0 25 JUN 03176 12:00:00 12 ME
> RC
>Z 400 MB 25 JUN 03176 12:00:00 AVN 0 25 JUN 03176 12:00:00 13 ME
> RC
>Z 300 MB 25 JUN 03176 12:00:00 AVN 0 25 JUN 03176 12:00:00 14 ME
> RC
>Z 250 MB 25 JUN 03176 12:00:00 AVN 0 25 JUN 03176 12:00:00 15 ME
> RC
>Z 200 MB 25 JUN 03176 12:00:00 AVN 0 25 JUN 03176 12:00:00 16 ME
> RC
>DTDP 300 MB 25 JUN 03176 12:00:00 AVN 0 25 JUN 03176 12:00:00 17 ME
> RC
>SUMT 300 MB 25 JUN 03176 12:00:00 AVN 0 25 JUN 03176 12:00:00 18 ME
> RC
>Number of grids listed = 18
>
>GRDLIST - done
I got more or less the same list minus the DTDP and SUMT parameters and
minus the 200 mb level fields that you did not include in the AVNLEV
string you listed above. You probably added the 200 mb level to the
string later.
>I then check to see if there are any missing data points in the
>temperature fields.
>And I see that there are no missing data points.Please see below.....
>
>
>GRDINFO MYDATA/GRIDS.1 PAR=T LEV=400
>
>Statistics of
>Parameter: T
>Units: K
>Level of data: 400 MB
>Time of data: 120000 UTC
>Day of data: 2003176
>Data came from: AVN 0-hour forecast
>Dataset: MYDATA/GRIDS
>
>Minimum: 0.2336000E+03 occurred at [row,col]: [ 10, 1]
> [lat,lon]: [ 68.75, 100.00]
>Maximum: 0.2633100E+03 occurred at [row,col]: [ 46, 134]
> [lat,lon]: [ 23.75, -66.25]
>Mean: 0.2499987E+03
>SD : 0.7024199E+01
>Number of points analyzed: 9821
>Number of points missing: 0
>
>Row range: 1 to 61
>Col range: 1 to 161
>GRDINFO Done, Number of grids statistically analyzed=1
>GRDINFO - done
I got the same results.
>GRDINFO MYDATA/GRIDS.1 PAR=T LEV=300
>Statistics of
>Parameter: T
>Units: K
>Level of data: 300 MB
>Time of data: 120000 UTC
>Day of data: 2003176
>Data came from: AVN
>0-hour forecast
>
>Dataset: MYDATA/GRIDS
>Minimum: 0.2216500E+03 occurred at [row,col]: [ 15, 114]
> [lat,lon]: [ 62.50, -41.25]
>Maximum: 0.2489900E+03 occurred at [row,col]: [ 40, 158]
> [lat,lon]: [ 31.25, -96.25]
>Mean: 0.2354649E+03
>SD : 0.6825965E+01
>Number of points analyzed: 9821
>Number of points missing: 0
>
>Row range: 1 to 61
>Col range: 1 to 161
>GRDINFO Done, Number of grids statistically analyzed=1
>GRDINFO - done
Ditto on the same results.
>GRDINFO MYDATA/GRIDS.1 PAR=T LEV=250
>Statistics of
>Parameter: T
>Units: K
>Level of data: 250 MB
>Time of data: 120000 UTC
>Day of data: 2003176
>Data came from: AVN
>0-hour forecast
>
>Dataset: MYDATA/GRIDS
>Minimum: 0.2168500E+03 occurred at [row,col]: [ 15, 97]
> [lat,lon]: [ 62.50, -20.00]
>Maximum: 0.2397700E+03 occurred at [row,col]: [ 37, 156]
> [lat,lon]: [ 35.00, -93.75]
>Mean: 0.2280361E+03
>SD : 0.5235160E+01
>
>Number of points analyzed: 9821
>Number of points missing: 0
>
>Row range: 1 to 61
>Col range: 1 to 161
>GRDINFO Done, Number of grids statistically analyzed=1
>GRDINFO - done
Ditto on the same results.
>I then try to add the 400 and 300 hPa temperature fields.
>
>GRDCOPY MYDATA/GRIDS.1 MYDATA/GRIDS.1 G1='PARAM T ;LEV 400' G2='PARAM
>T ; LEV 300' MATH='G1+G2' NEWPAR=SUMT K 300 mb
>Total grids copied=1
>GRDCOPY - done
>
>No problem
This worked for me as well.
>But when I try to do it for the 400 and 250 T fields, it fails.....
>
>GRDCOPY MYDATA/GRIDS.1 MYDATA/GRIDS.1 G1='PARAM T ;LEV 400' G2='PARAM
>T ; LEV 250' MATH='G1+G2' NEWPAR=SUMT K 250 mb
>GRDCOPY: Set of grids to make grid is incomplete
>
>GRDCOPY - done
I got the same _bad_ results as you did, so I can verify that the
problem is not operator error :-)
>Seeing as there are no missing data points, I can't figure out why this
>is failing.
>Any ideas?
Since you have to get this working soon, I can suggest that you use the
grid number in the select clause for G2 in your second invocation:
GRDCOPY MYDATA/GRIDS.1 MYDATA/GRIDS.1 G1='PARAM T;LEV 400' G2='GRID 3'
MATH='G1+G2' NEWPAR=SUMT K 250 mb
Now, if it becomes a problem knowing exactly which grid the 250 mb T
field is written into in the GRDCOPY, you could change your initial
GRDCOPY to copy specific RTGRIDS/AVN grids to specific MYDATA/GRIDS.1
grids. This way, you would be sure of the grid number for every
parameter and could use the GRID # in the G1 or G2 specifications.
>I have to run this code as part of a script in support of a field
>mission that starts next week so if you can respond soon, I would really
>appreciate it!
Tracing down the error in GRDCOPY would take some time as I am heavily
loaded at the moment. Will the workaround be OK for now?
Tom