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

20010323: 20010322: 20010321: dcreanal update



Lodovica,

It really doesn't seem like that is the new dcreanal.c routine since 
I believe there should at least be compiler output warnings about
the <ulog.h> instead of "ulog.h" include in dcreanal.c (something
that changed between 5.4 and 5.6.

So, I created a completely new tar file for you with a test data set that
I downloaded from LDEO. The compiled program will be called dcreanal2.

In the gbuddy ftp account:
host: ftp.unidata.ucar.edu
login: gbuddy
password: XXXXXX

you will find:
nawips-5.4/contrib/dcreanal2.tar.Z

unpack this file from $NAWIPS with the command:
zcat dcreanal2.tar.Z | tar xvf -

Then build and test with:

cd $NAWIPS/unidata/ldmbridge/dcreanal2
make all
make install
make test


The "make test" should run: dcreanal2 datau.cdf YYYYMM_ldeo.gem
The datau.cdf file contains Jan 1-Jan2 1948 u component grids.

I added some extra print statements to dcreanal.c so you should see output like:

dcreanal2 datau.cdf YYYYMM_ldeo.gem
 Changing center table to cntrgrib1.tbl
 Changing vertical coord table to vcrdgrib1.tbl
 Changing WMO parameter table to ncarncep1.tbl
 Changing center parameter table to ncarncep1.tbl

********************
Vertical coordinate id P
Latitude coordinate Y  size 73
Longitude coordinate X  size 144
Time coordinate T
Variable parameter u
********************
opening 194801_ldeo.gem 480101/0000
 480101/0000         UREL        PRES                1000          -1
 480101/0000         UREL        PRES                 925          -1
 480101/0000         UREL        PRES                 850          -1
 480101/0000         UREL        PRES                 700          -1
 480101/0000         UREL        PRES                 600          -1
 480101/0000         UREL        PRES                 500          -1
 480101/0000         UREL        PRES                 400          -1
 480101/0000         UREL        PRES                 300          -1
 480101/0000         UREL        PRES                 250          -1
 480101/0000         UREL        PRES                 200          -1
 480101/0000         UREL        PRES                 150          -1
 480101/0000         UREL        PRES                 100          -1
 480101/0000         UREL        PRES                  70          -1
 480101/0000         UREL        PRES                  50          -1
 480101/0000         UREL        PRES                  30          -1
 480101/0000         UREL        PRES                  20          -1
 480101/0000         UREL        PRES                  10          -1
 480102/0000         UREL        PRES                1000          -1
 480102/0000         UREL        PRES                 925          -1
 480102/0000         UREL        PRES                 850          -1
 480102/0000         UREL        PRES                 700          -1
 480102/0000         UREL        PRES                 600          -1
 480102/0000         UREL        PRES                 500          -1
 480102/0000         UREL        PRES                 400          -1
 480102/0000         UREL        PRES                 300          -1
 480102/0000         UREL        PRES                 250          -1
 480102/0000         UREL        PRES                 200          -1

If the above works for you, then try one of your data sets. If you still 
continue to have trouble decoding your data, then I must be missing
something in your cdl output.

Sorry this is taking several iterations. Hope this solves the trouble.

Steve Chiswell
Unidata User Support



>From: Lodovica Illari <address@hidden>
>Organization: UCAR/Unidata
>Keywords: 200103231825.f2NIPXL05904

>Hi Steve,
>
>I did what you suggested.
>
>This the output of my compilation. Does it look right?
>
>Lodovica
>
>------------------------------------------------------------------------------
> ---
>make all
>ln -s /usr/local/nawips-5.4/gempak5.4/include/ulog.h .
>ln -s /usr/local/nawips-5.4/gempak5.4/include/gemprm.h .
>ln -s /usr/local/nawips-5.4/gempak5.4/include/geminc.h .
>ln -s /usr/local/nawips-5.4/gempak5.4/source/gemlib/na/nacmn.cmn .
>ln -s /usr/local/nawips-5.4/gempak5.4/source/gemlib/gb/gbcmn.h .
>cc -DUNDERSCORE -DOSF1 -O0   -c -I/usr/local/include -I. -DGEMPAK
>-D__STDC__ -DLONG_64_BITS=1 -DNO_REPLACE_SYSLOG -DLOGNAME_ISSOCK
>dcreanal.c
>f77 -O0  -c open_gemgrid.f
>f77 -O0  -c put_gemgrid.f
>/usr/bin/cc -DUNDERSCORE -DOSF1 -O0 -v -o dcreanal dcreanal.o
>open_gemgrid.o put_gemgrid.o /usr/local/nawips-5.4/lib/osf/gemlib.a
>/usr/local/nawips-5.4/lib/osf/cgemlib.a /usr/local/lib/libnetcdf.a
>/usr/local/nawips-5.4/lib/osf/bridge.a -lUfor -lfor -lFutil -lm -lots 
>/usr/lib/cmplrs/cc/ld -o dcreanal -g0 -O0 -call_shared
>/usr/lib/cmplrs/cc/crt0.o dcreanal.o open_gemgrid.o put_gemgrid.o
>/usr/local/nawips-5.4/lib/osf/gemlib.a
>/usr/local/nawips-5.4/lib/osf/cgemlib.a /usr/local/lib/libnetcdf.a
>/usr/local/nawips-5.4/lib/osf/bridge.a -lUfor -lfor -lFutil -lm -lots -lc 
>/usr/lib/cmplrs/cc/ld: 
>0.16u 0.10s 0:03 8% 0+29k 434+117io 23pf+0w 29stk+4560mem
>(null command)
>> make install
>updated /usr/local/nawips-5.4/bin/osf/dcreanal
>------------------------------------------------------------------------------
> ---
>
>On Thu, 22 Mar 2001, Unidata Support wrote:
>
>> 
>> Lodovica,
>> 
>> This is the same format of the files I tested under the CDAS directory
>> on LDEO.
>> 
>> When you downloaded the dcreanal.c file to your 
>> $NAWIPS/unidata/ldmbridge/dcreanal directory, did you rebuild dcreanal with:
>> 
>> cd $NAWIPS/unidata/ldmbridge/dcreanal
>> make clean
>> make all
>> make install
>> make clean
>> 
>> Did the compiliation issue any errors?  Is your $GEMEXE/dcreanal updated?
>> 
>> Steve Chiswell
>> Unidata User Support
>> 
>> 
>> 
>> netcdf u_i {
>> dimensions:
>>         Y = 73 ;
>>         X = 144 ;
>>         T = 9 ;
>>         P = 17 ;
>> variables:
>>         float Y(Y) ;
>>                 Y:long_name = "Latitude" ;
>>                 Y:gridtype = 0 ;
>>                 Y:units = "degree_north" ;
>>         float X(X) ;
>>                 X:long_name = "Longitude" ;
>>                 X:gridtype = 1 ;
>>                 X:units = "degree_east" ;
>>         int T(T) ;
>>                 T:long_name = "Time" ;
>>                 T:defaultvalue = 19435 ;
>>                 T:calendar = "standard" ;
>>                 T:gridtype = 0 ;
>>                 T:units = "days since 1948-01-01" ;
>>         int P(P) ;
>>                 P:long_name = "Pressure" ;
>>                 P:gridtype = 0 ;
>>                 P:units = "mb" ;
>>         float u(T, P, Y, X) ;
>>                 u:scale_max = 150.8f ;
>>                 u:center = "US Weather Service - National Met. Center" ;
>>                 u:subcenter = "NCEP Reanalysis Project" ;
>>                 u:gribleveltype = 100 ;
>>                 u:units = "m/s" ;
>>                 u:calendar = "standard" ;
>>                 u:grib_name = "UGRD" ;
>>                 u:gribparam = 33 ;
>>                 u:long_name = "zonal wind" ;
>>                 u:PTVersion = 2 ;
>>                 u:PDS_TimeRange = 113 ;
>>                 u:scale_min = -104.7f ;
>>                 u:GRIBgridcode = 2 ;
>>                 u:process = "62 wave triangular, 28 layer Spectral model fro
> m \"Medium Ran
>> ge Forecast\" run" ;
>>                 u:missing_value = 9.999e+20f ;
>>                 u:Compression.compressiontype = "null" ;
>> }
>> >From: Lodovica Illari <address@hidden>
>> >Organization: UCAR/Unidata
>> >Keywords: 200103222209.f2MM9cL27954
>> 
>> >Hi Steve,
>> >
>> >Thanks for your update version of dcreanal.
>> >
>> >We tried the Intrisic pressure files from LDEO, but still dcreanal did
>> >not work.
>> >
>> >Could you tell us more specifically where these "Intrisic" NetCDF files
>> >are located on the LDEO server?
>> >
>> >Thanks again
>> >
>> >
>> >
>> 
>> ****************************************************************************
>> Unidata User Support                                    UCAR Unidata Program
>> (303)497-8644                                                  P.O. Box 3000
>> address@hidden                                   Boulder, CO 80307
>> ----------------------------------------------------------------------------
>> Unidata WWW Service                        http://www.unidata.ucar.edu/     
>> ****************************************************************************
>> 
>