[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[GEMPAK #YTR-838531]: Canadian RADAR data file



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