Deleted attribute reappears

I've got a file with 32 years of daily data that I'm breaking into files of
individual years. The time variable in the original file has an incorrect
"description" attribute that I am deleting. But when I write the new yearly
file, the attribute reappears! I can figure out a work around, but I'd like
to understand what I'm doing wrong. I'm running 5.2.0. Here is a code
snippet followed by the screen output:

CODE:

    times = fin->day
    delete_VarAtts(times,"description")
    printVarSummary(times)
    t = ut_calendar(times,-5)
    strs = str_split(fin,"_")
    vname = strs(0)
    do y = 1979,2010
        ii = ind(t(:,0).eq.y)
        fname = strs(0)+"_"+y+".nc"
        fout = addfile(fname,"c")
        t2 = times(ii)
        printVarSummary(t2)
        dat = fin->data(ii,:,:)
        dat!0 = "time"
        dat!1 = "lat"
        dat!2 = "lon"
        dat&time = t2
        fout->$vname$ = dat
        delete(ii)
        delete(dat)
        delete(t2)
        system("ncdump -h "+fname)
    end do

OUTPUT:

Variable: times
Type: float
Total Size: 46752 bytes
            11688 values
Number of Dimensions: 1
Dimensions and sizes: [day | 11688]
Coordinates:
            day: [722816..734503]
Number Of Attributes: 2
  units : days since 0001-01-01
  calendar : gregorian

Variable: t2
Type: float
Total Size: 1460 bytes
            365 values
Number of Dimensions: 1
Dimensions and sizes: [day | 365]
Coordinates:
            day: [722816..723180]
Number Of Attributes: 2
  calendar : gregorian
  units : days since 0001-01-01

netcdf pr_1979 {
dimensions:
    time = 365 ;
    lat = 200 ;
    lon = 464 ;
variables:
    float pr(time, lat, lon) ;
        pr:units = "mm" ;
        pr:description = "Daily Accumulated Precipitation" ;
        pr:_FillValue = -9999.f ;
    float time(time) ;
        time:units = "days since 0001-01-01" ;
        time:calendar = "gregorian" ;
        *time:description = "days since 1900-01-01" ;*
    float lat(lat) ;
        lat:units = "degrees_north" ;
        lat:description = "Latitude of the center of the grid cell" ;
    float lon(lon) ;
}

Any insight?

-- 
Bridget Thrasher, PhD
Independent Contractor, Research Scientist

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Jun 30 16:24:22 2011