用ncl写nc文件，数据小数点后保留的位数

1368984969

大家好。我刚开始使用ncl，用ncl写nc文件，但是在把经纬度写到nc文件后，经纬度本来是小数后面2位数，lat:[-89.95, -89.85, ... , 89.85, 89.95]
lon:[-179.95,-179.85, ..., 179.85, 179.95]
不管我是用lat=fspan(-89.95, 89.95, 1800)
lon=fspan(-179.95, 179.95, 3600)
还是从原本的nc数据里面读取lat，lon重新写入新的nc文件中
或者从txt文件里面读取已经写好的lat，lon数据写入到nc文件中，
nc文件的经纬度数据从小数点后两位变成在小数点后多出很多位数

这个问题困扰我很多天，希望大家能帮我解决一下

1. lat_1 = asciiread("filepath/lat.txt", 1800, "float")
   
2. lon_1 = asciiread("filepath/lon.txt", 3600, "float")
   
3. time  = new(12,double,"No_FillValue")
   
4. YYYY = tointeger(year) ;将字符转为数值
   
5. do i = 0, ntim-1
   
6.     MM = i+1
   
7.     time(i) = cd_inv_calendar(YYYY, MM, 01, 00, 00, 00, "hours since 1985-01-01 00:00:00", 0)
   
8. end do
   
9. 
   
10. do i = 0, dimsizes(pollutant_names)-1
    
11.     emi_1 = emi_anthro_1(i,:,:,:)
    
12.     var = pollutant_names(i)
    
13.     fout_name = "PKU_" + pollutant_names(i) + "_" + year + "_anthro.nc"
    
14.     ;Remove the file if it already existed
    
15.     system("rm -f " + diro + fout_name)
    
16.    
    
17.     ;create the file:
    
18.     fout = addfile(diro + fout_name, "c")
    
19.    
    
20.     ;define some file properties:
    
21.     atts = True
    
22.     atts@title = "NetCDF file containing " + pollutant_names(i) + " anthropogenic emission dataset"
    
23.     fout@creation_date=systemfunc("date")
    
24.     fileattdef(fout, atts)
    
25. 
    
26.     ;define dimensions:
    
27.     dim_Names = (/"time","lat","lon"/)
    
28.     dim_Sizes = (/12, nlat, nlon/)
    
29.     dim_Unlim = (/True, False, False/)
    
30.     filedimdef(fout, dim_Names, dim_Sizes, dim_Unlim)
    
31. 
    
32.     ;define file variables:
    
33.     filevardef(fout, "time", "double", "time")
    
34.     filevardef(fout, "lat", "float", "lat")
    
35.     filevardef(fout, "lon", "float", "lon")
    
36.     filevardef(fout, var, "double", (/"time","lat","lon"/))
    
37.     ;define file variables attributes:
    
38.     ;;time
    
39.     tatts = True
    
40.     tatts@long_name = "Time"
    
41.     tatts@units = "hours since 1985-01-01 00:00:00"
    
42.     tatts@axis      = "T"
    
43.     filevarattdef(fout, "time", tatts)
    
44.     ;;latitude
    
45.     latatts = True
    
46.     latatts@long_name = "latitude"
    
47.     latatts@units     = "degree_north"
    
48.     latatts@axis      = "X"
    
49.     filevarattdef(fout, "lat", latatts)
    
50.     ;;longitude
    
51.     lonatts = True
    
52.     lonatts@long_name = "longitude"
    
53.     lonatts@units     = "degree_east"
    
54.     lonatts@axis      = "Y"
    
55.     filevarattdef(fout, "lon", lonatts)
    
56.     ;;emission data
    
57.     emiatts = True
    
58.     emiatts@long_name = "Gridded emission desities of " + pollutant_names(i)
    
59.     emiatts@units = "kg m-2 s-1"
    
60.     filevarattdef(fout, var, emiatts)
    
61. 
    
62.    
    
63.     ;write out  the data
    
64.     fout->lat   = (/lat_1/)
    
65.     fout->lon   = (/lon_1/)
    
66.     fout->time  = (/time/)
    
67.     fout->$var$ = (/emi_1/)
    
68. end do
    
69. print("Finish the Step 5")

复制代码