[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[GEMPAK #YTR-838531]: Canadian RADAR data file
- Subject: [GEMPAK #YTR-838531]: Canadian RADAR data file
- Date: Fri, 10 Nov 2006 15:35:53 -0700
Stephen,
You must not have built the tar file I posted this morning when I sent you the
new binary
since the radar_grid.c routine uses rint() rather than lround(), so that is
the one fundemental problem. You should make sure that $GEMLIB/gdradrc.a is
removed
from your previous build attempt (which "make clean" will do) to ensure you are
getting
the new library code. The binary I sent doesn't use the prodmmap that the
previous
version did, so that should have alleviated the read errors under linux.
The code uses gtrans to convert the lat and lon of each data point to
a location within the KX and KY space. In my example then, the largest value
would be 1700. Each data point is read, converted, and tested for
insertion into the grid. The number of data vlues in the radar data values
will therefore be inconsequential.
While the program ran much faster on my slow SGI than the Linux PC,
it still was under 5 seconds.
% time gdradrc<<EOF
? r
?
? e
? EOF
GRDAREA Area covered by grid 44;-128;47;-48
PROJ Map projection/angles/margins|dr LCC/50.0;-110.0;50.0
KXKY Number of grid points in x;y 1700;950
GDPFUN Scalar grid or vector grid funct n1r
GDFILE Grid file 20061104_radr.gem
RADTIM Radar composite current/dattim 061104/1530
RADDUR Radar time window (minutes prior 30
RADFRQ Update Frequency
CPYFIL Grid file whose navigation is to
STNFIL Station information file cradar.tbl
MAXGRD Maximum number of grids 1000
NDVAL Replacement value for "ND" (defa -9999
Parameters requested: GRDAREA,PROJ,KXKY,GDPFUN,GDFILE,RADTIM,RADDUR,RADFRQ,
CPYFIL,STNFIL,MAXGRD,NDVAL.
GEMPAK-GDRADRC> [GDRADRC 0] Using: $RAD/cradar/WBI/N1R/N1R_20061104_1521
[GDRADRC 0] Using: $RAD/cradar/WGJ/N1R/N1R_20061104_1522
[GDRADRC 0] Using: $RAD/cradar/WHK/N1R/N1R_20061104_1521
[GDRADRC 0] Using: $RAD/cradar/WHN/N1R/N1R_20061104_1522
[GDRADRC 0] Using: $RAD/cradar/WKR/N1R/N1R_20061104_1522
[GDRADRC 0] Using: $RAD/cradar/WMB/N1R/N1R_20061104_1520
[GDRADRC 0] Using: $RAD/cradar/WSO/N1R/N1R_20061104_1521
[GDRADRC 0] Using: $RAD/cradar/WTP/N1R/N1R_20061104_1530
[GDRADRC 0] Using: $RAD/cradar/WUJ/N1R/N1R_20061104_1521
[GDRADRC 0] Using: $RAD/cradar/WVY/N1R/N1R_20061104_1520
[GDRADRC 0] Using: $RAD/cradar/WWW/N1R/N1R_20061104_1522
[GDRADRC 0] Using: $RAD/cradar/XRA/N1R/N1R_20061104_1530
[GDRADRC 0] Using: $RAD/cradar/XBE/N1R/N1R_20061104_1521
[GDRADRC 0] Using: $RAD/cradar/XFT/N1R/N1R_20061104_1521
[GDRADRC 0] Using: $RAD/cradar/XWL/N1R/N1R_20061104_1521
[GDRADRC 0] Using: $RAD/cradar/XME/N1R/N1R_20061104_1521
[GDRADRC 0] Using: $RAD/cradar/XGO/N1R/N1R_20061104_1520
[GDRADRC 0] Using: $RAD/cradar/XLA/N1R/N1R_20061104_1521
[GDRADRC 0] Using: $RAD/cradar/XAM/N1R/N1R_20061104_1521
[GDRADRC 0] Using: $RAD/cradar/XSM/N1R/N1R_20061104_1530
[GDRADRC 0] Using: $RAD/cradar/XSI/N1R/N1R_20061104_1522
[GDRADRC 0] Using: $RAD/cradar/XBU/N1R/N1R_20061104_1530
[GDRADRC 0] Using: $RAD/cradar/XFW/N1R/N1R_20061104_1521
[GDRADRC 0] Using: $RAD/cradar/XMB/N1R/N1R_20061104_1521
[GDRADRC 0] Using: $RAD/cradar/XDR/N1R/N1R_20061104_1522
[GDRADRC 0] Using: $RAD/cradar/XNI/N1R/N1R_20061104_1522
[GDRADRC 0] Using: $RAD/cradar/XSS/N1R/N1R_20061104_1521
[GDRADRC 0] Using: $RAD/cradar/XPG/N1R/N1R_20061104_1521
[GDRADRC 0] Using: $RAD/cradar/XNC/N1R/N1R_20061104_1530
[GDRADRC 0] Using: $RAD/cradar/XTI/N1R/N1R_20061104_1521
[GDRADRC 4] Write grid 061104/1530
Parameters requested: GRDAREA,PROJ,KXKY,GDPFUN,GDFILE,RADTIM,RADDUR,RADFRQ,
CPYFIL,STNFIL,MAXGRD,NDVAL.
GEMPAK-GDRADRC> GEMPAK-GDRADRC>
1.812u 0.112s 0:02.54 75.5% 0+0k 0+0io 0pf+0w
This is on an old 32 bit FC4 machine, which is what I make the linux binaries
on.
% uname -a
Linux rip.unidata.ucar.edu 2.6.17-1.2142_FC4smp #1 SMP Tue Jul 11 22:57:02 EDT
2006 i686 athlon i386 GNU/Linux
Steve Chiswell
Unidata User Support
> Thanks Chiz.
>
> I have tried using the binary for i386 that you attached and get the same
> results, all values greater than 10**31.
>
> I have installed the binary release and libraries for GEMPAK5.9.4 and copied
> the $NAWIPS/config from my source installation and tried a make of gdradrc.
> Same results. There must be something fundamentally different about my
> hardware/os.
>
> I have tried changing the KXKY values to reduce the size of my grid,...
>
> I did notice that it took over ten minutes to process 10 radar files, which
> does not match with what you saw (0.5 seconds for the sample I sent).
>
> There was a warning for the lroundf function in radar_grid.c in all cases:
> radar_grid.c:121: warning: incompatible implicit declaration of built-in
> function 'lroundf'
>
> If there are more than 64K points in the radar file, how do process that
> number with type int?
>
> I am learning ddd to try and walk through the logic.
>
> I have included the output from the latest run in case there might be some
> clue that you might see ...
>
> Thanks again for all your assistance.
>
> [ssinnis@gempak3 canadian_radar]$ ../gdradrc
> GRDAREA Area covered by grid 44;-128;47;-48
> PROJ Map projection/angles/margins|dr lcc/50;-110;50
> KXKY Number of grid points in x;y 1700;950
> GDPFUN Scalar grid or vector grid funct n1r
> GDFILE Grid file 20061109_cradr.gem
> RADTIM Radar composite current/dattim 061109/1715
> RADDUR Radar time window (minutes prior 25
> RADFRQ Update Frequency
> CPYFIL Grid file whose navigation is to
> STNFIL Station information file cradar.tbl
> MAXGRD Maximum number of grids 2
> NDVAL Replacement value for "ND" (defa -1
> Parameters requested: GRDAREA,PROJ,KXKY,GDPFUN,GDFILE,RADTIM,RADDUR,RADFRQ,
> CPYFIL,STNFIL,MAXGRD,NDVAL.
> GEMPAK-GDRADRC>kxky=500;200
> GEMPAK-GDRADRC>r
> [GDRADRC 4] Write grid 061109/1715
> Parameters requested: GRDAREA,PROJ,KXKY,GDPFUN,GDFILE,RADTIM,RADDUR,RADFRQ,
> CPYFIL,STNFIL,MAXGRD,NDVAL.
> GEMPAK-GDRADRC>e
> [ssinnis@gempak3 canadian_radar]$ ll
> total 652
> drwxr-xr-x 2 ssinnis twnwrk 4096 Nov 10 15:34 .
> drwxrwxrwx 3 ssinnis twnwrk 4096 Nov 10 13:06 ..
> -rw-rw-rw- 1 ssinnis twnwrk 204800 Nov 10 15:34 20061109_cradr.gem
> -rw-r--r-- 1 ssinnis twnwrk 2541 Nov 10 15:31 cradar.tbl
> -rw-r--r-- 1 ssinnis twnwrk 138 Nov 10 11:26 datatype.tbl
> -rwxr-xr-x 1 ssinnis twnwrk 234 Nov 6 12:28 dounzip.csh
> -rw-r--r-- 1 ssinnis twnwrk 1718 Nov 2 15:42 gdcanl.f
> -rw-r--r-- 1 ssinnis twnwrk 3406 Nov 2 15:42 gdcnav.f
> -rw-r--r-- 1 ssinnis twnwrk 4370 Nov 2 15:42 gdctbl.f
> -rw-r--r-- 1 ssinnis twnwrk 5480 Nov 6 16:00 gdplot2.nts
> -rwxrwxrwx 1 ssinnis twnwrk 305487 Nov 10 15:23 gdradrc
> -rw-r--r-- 1 ssinnis twnwrk 522 Nov 6 15:53 gdradrc.err
> -rw-r--r-- 1 ssinnis twnwrk 9706 Nov 6 15:22 gdradrc.f
> -rw-r--r-- 1 ssinnis twnwrk 4499 Nov 6 15:53 gdradrc.hlp
> -rw-r--r-- 1 ssinnis twnwrk 1644 Nov 10 11:29 gdradrc.nts
> -rw-r--r-- 1 ssinnis twnwrk 906 Nov 6 15:46 gdradrc.pdf
> -rw-r--r-- 1 ssinnis twnwrk 80 Nov 6 15:46 gdradrc.prm
> -rw-r--r-- 1 ssinnis twnwrk 3272 Nov 6 15:55 gdradrc.txt
> -rw-r--r-- 1 ssinnis twnwrk 1781 Nov 6 14:50 gdradr.nts
> -rw-r--r-- 1 ssinnis twnwrk 23153 Nov 10 15:34 gemglb.nts
> -rw-r--r-- 1 ssinnis twnwrk 1454 Nov 6 15:51 gpinp.f
> -rw-r--r-- 1 ssinnis twnwrk 1644 Nov 10 15:34 last.nts
> -rw-r--r-- 1 ssinnis twnwrk 1002 Nov 10 09:08 Makefile
> -rw-r--r-- 1 ssinnis twnwrk 3365 Nov 10 15:23 radar_grid.c
> -rw-rw-rw- 1 ssinnis twnwrk 10080 Nov 10 15:23 radar_grid.o
> -rw-r--r-- 1 ssinnis twnwrk 620 Nov 2 15:42 waittime.c
> [ssinnis@gempak3 canadian_radar]$ gdinfo
> GDFILE Grid file 20061109_cradr.gem
> LSTALL Full list flag yes
> OUTPUT Output device/filename t
> GDATTIM Grid date/time 061109/1715
> GLEVEL Grid level 0
> GVCORD Grid vertical coordinate none
> GFUNC Scalar grid n1r
> Parameters requested: GDFILE,LSTALL,OUTPUT,GDATTIM,GLEVEL,GVCORD,GFUNC.
> GEMPAK-GDINFO>r
>
> GRID FILE: 20061109_cradr.gem
>
>
> GRID NAVIGATION:
> PROJECTION: LCC
> ANGLES: 50.0 -110.0 50.0
> GRID SIZE: 500 200
> LL CORNER: 44.00 -128.00
> UR CORNER: 47.00 -48.00
>
> GRID ANALYSIS BLOCK:
> ANALYSIS TYPE: BARNES
> DELTAN: 0.298
> DELTAX: -9999.000
> DELTAY: -9999.000
> GRID AREA: 33.00 -138.00 64.00 -47.00
> EXTEND AREA: 32.00 -138.00 65.00 -47.00
> DATA AREA: 32.00 -138.00 65.00 -47.00
>
> Number of grids in file: 1
>
> Maximum number of grids in file: 2
>
> NUM TIME1 TIME2 LEVL1 LEVL2 VCORD PARM
> 1 061109/1715 0 NONE N1R
>
> Parameters requested: GDFILE,LSTALL,OUTPUT,GDATTIM,GLEVEL,GVCORD,GFUNC.
> GEMPAK-GDINFO>e
> [ssinnis@gempak3 canadian_radar]$ gdlist
> GDATTIM Grid date/time 061109/1715
> GLEVEL Grid level 0
> GVCORD Grid vertical coordinate none
> GFUNC Scalar grid n1r
> GDFILE Grid file 20061109_cradr.gem
> GAREA Graphics area #52.52;-107.44;1;1
> PROJ Map projection/angles/margins|dr lcc/50;-110;50
> SCALE Scalar scale / vector scale 0
> OUTPUT Output device/filename t
> Parameters requested: GDATTIM,GLEVEL,GVCORD,GFUNC,GDFILE,GAREA,PROJ,SCALE,
> OUTPUT.
> GEMPAK-GDLIST>r
>
>
> GDLIST PARAMETERS:
>
> Grid file: 20061109_cradr.gem
>
> GRID IDENTIFIER:
> TIME1 TIME2 LEVL1 LEVL2 VCORD PARM
> 061109/1715 0 NONE N1R
>
> GAREA: #52.52;-107.44;1;1
> SCALE FACTOR : 10** 0
> OUTPUT: TERM/
>
> MINIMUM AND MAXIMUM VALUES******************
> Enter <cr> to accept parameters or type EXIT:e
> Parameters requested: GDATTIM,GLEVEL,GVCORD,GFUNC,GDFILE,GAREA,PROJ,SCALE,
> OUTPUT.
> GEMPAK-GDLIST>
>
>
Ticket Details
===================
Ticket ID: YTR-838531
Department: Support GEMPAK
Priority: Urgent
Status: Closed