Re: macos and linux ncl 6.1 time_step errors

From: Wei Huang <huangwei_at_nyahnyahspammersnyahnyah>
Date: Fri Dec 07 2012 - 15:10:38 MST

Barry,

Thanks for providing us the data to verify this issue.

It turns out to be an NCL attributes printing issue.
(We have fixed it in our current code, and it will be in the next release).

You may try to print this specific attribute by add 2 lines below to your test script:

 time = f->time

 print("time@time_step = " + time@time_step)

(its value should be 21600).

Thanks,

Wei

huangwei@ucar.edu
VETS/CISL
National Center for Atmospheric Research
P.O. Box 3000 (1850 Table Mesa Dr.)
Boulder, CO 80307-3000 USA
(303) 497-8924

On Dec 5, 2012, at 6:50 AM, Barry McInnes (NOAA Affiliate) wrote:

> For macos 10.8 the time_step is incorrect and changes
> For Linux its incorrect but constant.
> Mac
> Current ncl command linked to ncl.xq.fix
>
> [mac11:~/sjl] bmcinnes% ls -l /usr/local/bin/ncl
> lrwxr-xr-x 1 root bin 10 Nov 30 10:53 /usr/local/bin/ncl -> ncl.xq.fix
> [mac11:/usr/local/bin] bmcinnes% ldd ncl.xq.fix
> ncl.xq.fix:
> /usr/lib/libcurl.4.dylib (compatibility version 7.0.0, current version 7.0.0)
> /usr/lib/libiconv.2.dylib (compatibility version 7.0.0, current version 7.0.0)
> /usr/X11/lib/libXrender.1.dylib (compatibility version 5.0.0, current version 5.0.0)
> /usr/lib/libbz2.1.0.dylib (compatibility version 1.0.0, current version 1.0.5)
> /usr/X11/lib/libXpm.4.dylib (compatibility version 16.0.0, current version 16.0.0)
> /usr/X11/lib/libX11.6.dylib (compatibility version 10.0.0, current version 10.0.0)
> /usr/X11/lib/libXext.6.dylib (compatibility version 11.0.0, current version 11.0.0)
> /usr/lib/libSystem.B.dylib (compatibility version 1.0.0, current version 159.1.0)
> /usr/local/lib/libstdc++.6.dylib (compatibility version 7.0.0, current version 7.14.0)
> /usr/local/lib/libgcc_s.1.dylib (compatibility version 1.0.0, current version 1.0.0)
>
> Run 6.1 on a static file and search for time_step
>
> [mac11:~] bmcinnes% cd sjl
> [mac11:~/sjl] bmcinnes% ncl ncl.script | grep time_step
> time_step_setting : auto
> time_step : 115372808
> [mac11:~/sjl] bmcinnes% ncl ncl.script | grep time_step
> time_step_setting : auto
> time_step : 267854600
> [mac11:~/sjl] bmcinnes% ncl ncl.script | grep time_step
> time_step_setting : auto
> time_step : 117617416
> [mac11:~/sjl] bmcinnes% ncl ncl.script | grep time_step
> time_step_setting : auto
> time_step : 189281032
> [mac11:~/sjl] bmcinnes% ncl ncl.script | grep time_step
> time_step_setting : auto
> time_step : 120001288
> [mac11:~/sjl] bmcinnes% ncl ncl.script | grep time_step
> time_step_setting : auto
> time_step : 239285000
>
> Run the previous version now named ncl6.0
>
> [mac11:~/sjl] bmcinnes% ncl6.0 ncl.script | grep time_step
> time_step_setting : auto
> time_step : 21600
> [mac11:~/sjl] bmcinnes% ncl6.0 ncl.script | grep time_step
> time_step_setting : auto
> time_step : 21600
> [mac11:~/sjl] bmcinnes% ncl6.0 ncl.script | grep time_step
> time_step_setting : auto
> time_step : 21600
> [mac11:~/sjl] bmcinnes% ncl6.0 ncl.script | grep time_step
> time_step_setting : auto
> time_step : 21600
> [mac11:~/sjl] bmcinnes% ncl6.0 ncl.script | grep time_step
> time_step_setting : auto
> time_step : 21600
>
> [mac11:~/sjl] bmcinnes% cat ncl.script
> f=addfile("/Projects/CFSR/.pressure/air.201103.nc4","r")
> print(f)
>
> [mac11:~/sjl] bmcinnes% uname -a
> Darwin mac11.psd.esrl.noaa.gov 12.2.0 Darwin Kernel Version 12.2.0: Sat Aug 25 00:48:52 PDT 2012; root:xnu-2050.18.24~1/RELEASE_X86_64 x86_64
>
> Running direct ncl6.1 instead of the xq.fix version
>
> [mac11:~/sjl] bmcinnes% ncl6.1 ncl.script | grep time_step
> time_step_setting : auto
> time_step : 55198552
> [mac11:~/sjl] bmcinnes% ncl6.1 ncl.script | grep time_step
> time_step_setting : auto
> time_step : 193966936
> [mac11:~/sjl] bmcinnes% ncl6.1 ncl.script | grep time_step
> time_step_setting : auto
> time_step : 194491224
> [mac11:~/sjl] bmcinnes% ncl6.1 ncl.script | grep time_step
> time_step_setting : auto
> time_step : 171086680
> [mac11:~/sjl] bmcinnes% ncl6.1 ncl.script | grep time_step
> time_step_setting : auto
> time_step : 122868568
> [mac11:~/sjl] bmcinnes%
>
> ncdump -th /Projects/CFSR/.pressure/air.201103.nc4
> netcdf air.201103 {
> dimensions:
> latitude = 361 ;
> longitude = 720 ;
> time = UNLIMITED ; // (124 currently)
> level = 37 ;
> variables:
> double latitude(latitude) ;
> latitude:units = "degrees_north" ;
> latitude:long_name = "latitude" ;
> double longitude(longitude) ;
> longitude:units = "degrees_east" ;
> longitude:long_name = "longitude" ;
> double time(time) ;
> time:units = "seconds since 1970-01-01 00:00:00.0 0:00" ;
> time:long_name = "verification time generated by wgrib2 function verftime()" ;
> time:reference_time = 1298937600. ; // "2011-03-01"
> time:reference_date = "2011.03.01 00:00:00 UTC" ;
> time:time_step_setting = "auto" ;
> time:reference_time_type = 2 ; // "1970-01-01 00:00:2.000000"
> time:reference_time_description = "analyses, reference date is variable, min found reference date is given" ;
> time:time_step = 21600. ; // "1970-01-01 06"
> float level(level) ;
> level:lev_type = 100 ;
> level:units = "mb" ;
> level:long_name = "1 mb" ;
> level:_FillValue = 9.999e+20f ;
> float air(time, level, latitude, longitude) ;
> air:_FillValue = 9.999e+20f ;
> air:short_name = "air" ;
> air:long_name = "Temperature" ;
> air:level = "1 mb" ;
> air:units = "K" ;
>
> // global attributes:
> :Conventions = "COARDS" ;
> :History = "created by wgrib2" ;
> :GRIB2_grid_template = 0 ;
> }
>
> Linux
> [bmcinnes@batchb ~/sjl]$ uname -a
> Linux batchb.psd.esrl.noaa.gov 2.6.18-308.20.1.el5 #1 SMP Tue Nov 6 04:38:29 EST 2012 x86_64 GNU/Linux
>
> ncl is 6.0
> [bmcinnes@batchb ~/sjl]$ ncl ncl.script | grep time_step
> time_step_setting : auto
> time_step : 21600
> [bmcinnes@batchb ~/sjl]$ ncl ncl.script | grep time_step
> time_step_setting : auto
> time_step : 21600
> [bmcinnes@batchb ~/sjl]$ ncl ncl.script | grep time_step
> time_step_setting : auto
> time_step : 21600
> [bmcinnes@batchb ~/sjl]$ ncl ncl.script | grep time_step
> time_step_setting : auto
> time_step : 21600
>
>
> ncl6.1 is 6.1
> [bmcinnes@batchb ~/sjl]$ ncl6.1 ncl.script | grep time_step
> time_step_setting : auto
> time_step : 7256532
> [bmcinnes@batchb ~/sjl]$ ncl6.1 ncl.script | grep time_step
> time_step_setting : auto
> time_step : 7256532
> [bmcinnes@batchb ~/sjl]$ ncl6.1 ncl.script | grep time_step
> time_step_setting : auto
> time_step : 7256532
> [bmcinnes@batchb ~/sjl]$ ncl6.1 ncl.script | grep time_step
> time_step_setting : auto
> time_step : 7256532
> [bmcinnes@batchb ~/sjl]$ ncl6.1 ncl.script | grep time_step
> time_step_setting : auto
> time_step : 7256532
> [bmcinnes@batchb ~/sjl]$ which ncl6.1
> /usr/local/bin/ncl6.1
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sat Dec 8 07:10:56 2012

This archive was generated by hypermail 2.1.8 : Wed Dec 12 2012 - 19:07:01 MST