equal
deleted
inserted
replaced
|
1 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" |
|
2 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" |
|
3 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" |
|
4 |
|
5 begin |
|
6 nlat = 3600 |
|
7 mlon = 7200 |
|
8 |
|
9 diri = "/fis/cgd/cseg/people/jeff/clamp_data/" |
|
10 fili = "MOD15A2_GEO_2000001.nc" |
|
11 a = addfile(diri+fili,"r") |
|
12 |
|
13 y = a->Lai_0_05deg ; byte |
|
14 y@_FillValue = inttobyte(253) |
|
15 y@_FillValue = inttobyte(255) |
|
16 |
|
17 x = byte2flt(y) |
|
18 |
|
19 x@_FillValue = 1.e+36 ; new fill |
|
20 |
|
21 lat = latGlobeFo(nlat, "lat", "latitude", "degrees_north") |
|
22 lat = (/ lat(::-1) /) ; make N->S |
|
23 lon = lonGlobeFo(mlon, "lon", "longitude", "degrees_east") |
|
24 lon = (/ lon - 180. /) ; subtract 180 from all values |
|
25 lon&lon = lon ; update coordinates |
|
26 |
|
27 x!0 = "lat" |
|
28 x!1 = "lon" |
|
29 x&lat= lat |
|
30 x&lon= lon |
|
31 |
|
32 y!0 = "lat" |
|
33 y!1 = "lon" |
|
34 y&lat= lat |
|
35 y&lon= lon |
|
36 |
|
37 iy = new (dimsizes(y), "integer", -999) |
|
38 iy = y |
|
39 |
|
40 iy!0 = "lat" |
|
41 iy!1 = "lon" |
|
42 iy&lat= lat |
|
43 iy&lon= lon |
|
44 |
|
45 print(y(:,{130})+" "+ iy(:,{130})+" "+x(:,{130})) |
|
46 |
|
47 |
|
48 |
|
49 ; diro = "/fis/cgd/cseg/people/jeff/clamp_data/" |
|
50 ; filo = "LAI_2000001_2.nc" |
|
51 ; c = addfile(diro+filo,"c") |
|
52 ; c->LAI = x |
|
53 end |
|
54 |