clamp: co2/11.test.ncl@0c6405ab2ff4

     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

    19  ncase = dimsizes(lat)

    20  print (ncase)

    22  lat_tmp = lat

    23  lat_tmp@_FillValue = 1.e+36

    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

    35  indexes = ind(.not. ismissing(lat_tmp))

    36  print (indexes)

    38  lat_new = lat(indexes)

    39  lon_new = lon(indexes)

    40  print (lat_new +"/"+lon_new)

    42 exit

    44 end

author	Forrest Hoffman <forrest@climatemodeling.org>
	Mon, 26 Jan 2009 22:08:20 -0500
changeset 0	0c6405ab2ff4
permissions	-rw-r--r--