equal
deleted
inserted
replaced
|
1 ; *********************************************** |
|
2 ; xy_4.ncl |
|
3 ; *********************************************** |
|
4 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" |
|
5 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" |
|
6 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" |
|
7 ;************************************************ |
|
8 begin |
|
9 ;************************************************ |
|
10 ; read in data: observed |
|
11 ;************************************************ |
|
12 diri = "/fis/cgd/cseg/people/jeff/clamp_data/co2/" |
|
13 fili = "co2_globalView_222.nc" |
|
14 g = addfile (diri+fili,"r") |
|
15 a = g->CO2_SEAS |
|
16 lon = g->LON |
|
17 lat = g->LAT |
|
18 |
|
19 ncase = dimsizes(lat) |
|
20 print (ncase) |
|
21 |
|
22 lat_tmp = lat |
|
23 lat_tmp@_FillValue = 1.e+36 |
|
24 |
|
25 do n = 0,ncase-1 |
|
26 if (.not. ismissing(lat_tmp(n))) then |
|
27 indexes = ind(lat(n) .eq. lat .and. lon(n) .eq. lon) |
|
28 if (dimsizes(indexes) .gt. 1) then |
|
29 lat_tmp(indexes(1:)) = lat_tmp@_FillValue |
|
30 end if |
|
31 delete (indexes) |
|
32 end if |
|
33 end do |
|
34 |
|
35 indexes = ind(.not. ismissing(lat_tmp)) |
|
36 print (indexes) |
|
37 |
|
38 lat_new = lat(indexes) |
|
39 lon_new = lon(indexes) |
|
40 print (lat_new +"/"+lon_new) |
|
41 |
|
42 exit |
|
43 |
|
44 end |