1.1 --- /dev/null Thu Jan 01 00:00:00 1970 +0000
1.2 +++ b/class/03.class_T31_write.ncl Mon Jan 26 22:08:20 2009 -0500
1.3 @@ -0,0 +1,61 @@
1.4 +;*************************************************
1.5 +; ce_1.ncl
1.6 +;************************************************
1.7 +load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
1.8 +load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
1.9 +load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
1.10 +;************************************************
1.11 +begin
1.12 +
1.13 +; final output
1.14 +
1.15 + diro = "/fis/cgd/cseg/people/jeff/clamp_data/model/"
1.16 + filo = "class_pft_T31.nc"
1.17 + c = addfile(diro+filo,"c")
1.18 +
1.19 +; read model data
1.20 +
1.21 + diri = "/fis/cgd/cseg/people/jeff/surface_data/"
1.22 + fili = "surfdata_48x96_c070501.nc"
1.23 + f = addfile(diri+fili,"r")
1.24 +
1.25 + y = f->PCT_PFT
1.26 +; printVarSummary(y)
1.27 +
1.28 +; read model grid data
1.29 +
1.30 + diri = "/fis/cgd/cseg/people/jeff/surface_data/"
1.31 + fili = "lnd_T31.nc"
1.32 + g = addfile(diri+fili,"r")
1.33 +
1.34 + landmask = g->landmask
1.35 + lat = g->lat
1.36 + lon = g->lon
1.37 +
1.38 + nlat = dimsizes(lat)
1.39 + nlon = dimsizes(lon)
1.40 +
1.41 + x = y(0,:,:)
1.42 +
1.43 + x!0 = "lat"
1.44 + x&lat = lat
1.45 + x!1 = "lon"
1.46 + x&lon = lon
1.47 + x@_FillValue = 1.e36
1.48 + x@long_name = "Model PFT Classes"
1.49 +
1.50 + do j= 0,nlat-1
1.51 + do i= 0,nlon-1
1.52 + x(j,i) = maxind(y(:,j,i))
1.53 + end do
1.54 + end do
1.55 +
1.56 +; print (x)
1.57 +
1.58 + x = where(landmask .lt. 1.,x@_FillValue,x)
1.59 +
1.60 + c->CLASS_PFT = x
1.61 + c->lat = lat
1.62 + c->lon = lon
1.63 +
1.64 +end
1.65 \ No newline at end of file