[ncl-talk] Error on relative units conversion, string = day as %Y%m%d.%f

[李嘉鹏]

Hi, I don't think cd_calendar supports units attribute "day as %Y%m%d.%f", see its documentation. Based on your data description, you can use

    ymd = round(time, 3)

---
Regards

于 2014年10月31日 格林尼治标准时间+0800下午11时30分20秒, Thomas Tobian <thomastobian at gmail.com> 写到:
>Hallo,
>
>Anyone know how to solve the problem with my historical model data,
>1850 to
>1875:
>'
>      institution : Max-Planck-Institute for Meteorology
>      CDO : Climate Data Operators version 1.5.9 (
>http://code.zmaw.de/projects/cdo)
>   dimensions:
>      lon = 192
>      lat = 96
>      lev = 1
>      time = 9496  // unlimited
>   variables:
>      double lon ( lon )
>         standard_name :        longitude
>         long_name :    longitude
>         units :        degrees_east
>         axis : X
>
>      double lat ( lat )
>         standard_name :        latitude
>         long_name :    latitude
>         units :        degrees_north
>         axis : Y
>
>      double lev ( lev )
>         standard_name :        air_pressure
>         long_name :    pressure
>         units :        Pa
>         positive :     down
>         axis : Z
>
>      double time ( time )
>         standard_name :        time
>         units :        day as %Y%m%d.%f
>         calendar :     proleptic_gregorian
>
>      float var131 ( time, lev, lat, lon )
>         table :        128
>         grid_type :    gaussian
>
>I got errors once I tried to subset the data using NCL as follows:
>
>f1 = addfile ("day_data.1850.to.1875.131.25000.ano.nc","r")
>ymdStrt  = 18500101
>ymdLast = 18751231
>time    = f2->time
>
>
>*ymd     = cd_calendar(time,-2)CDMS error: Error on relative units
>conversion, string = day as %Y%m%d.%f*
>
>
>*CDMS error: Error on time conversion: invalid month = 3664*
>;======
>
>ncl 9> print(time)
>
>
>Variable: time
>Type: double
>Total Size: 75968 bytes
>            9496 values
>Number of Dimensions: 1
>Dimensions and sizes:   [time | 9496]
>Coordinates:
>            time: [18500101.75..18751231.75]
>Number Of Attributes: 3
>  standard_name :       time
>  units :       day as %Y%m%d.%f
>  calendar :    proleptic_gregorian
>(0)     18500101.75
>(1)     18500102.75
>(2)     18500103.75
>(3)     18500104.75
>(4)     18500105.75
>(5)     18500106.75
>(6)     18500107.75
>(7)     18500108.75
>(8)     18500109.75
>(9)     18500110.75
>(10)    18500111.75
>(11)    18500112.75
>(12)    18500113.75
>
>
>Thank you
>-- Thomas
>
>
>------------------------------------------------------------------------
>
>_______________________________________________
>ncl-talk mailing list
>List instructions, subscriber options, unsubscribe:
>http://mailman.ucar.edu/mailman/listinfo/ncl-talk