NSF Unidata

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

20010323: 20010322: 20010321: dcreanal update

This archive contains answers to questions sent to Unidata support through mid-2025. Note that the archive is no longer being updated. We provide the archive for reference; many of the answers presented here remain technically correct, even if somewhat outdated. For the most up-to-date information on the use of NSF Unidata software and data services, please consult the Software Documentation first.



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/     
>> ****************************************************************************
>> 
>