forrest@0: ;************************************************
forrest@0: ;    Read ascii, Write nc                         
forrest@0: ;************************************************
forrest@0: load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"  
forrest@0: load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"  
forrest@0: load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"   
forrest@0: ;************************************************
forrest@0: begin
forrest@0: ;---------------------------------------------------
forrest@0: ; final data
forrest@0:                                     
forrest@0:   c = addfile("amazon_mask_1x1.nc","c") 
forrest@0:        
forrest@0:   nlat  = 180
forrest@0:   nlon  = 360  
forrest@0: 
forrest@0: ;************************************************
forrest@0: ; read 1x1 mask data
forrest@0: ;************************************************
forrest@0: 
forrest@0:   diri  = "/fis/cgd/cseg/people/jeff/clamp_data/observed/biomass/"
forrest@0:   fili  = "mask_1x1.txt"
forrest@0:   t     = asciiread(diri+fili,(/180,360/),"integer")
forrest@0: 
forrest@0:   printVarSummary(t)
forrest@0: 
forrest@0:   print (min(t)+"/"+max(t))                      
forrest@0: 
forrest@0:   lon      = new((/nlon/),float)
forrest@0:   lat      = new((/nlat/),float)
forrest@0: 
forrest@0: ; data is (89.5N->89.5S) (-179.5->179.5)
forrest@0: 
forrest@0: ; sam result
forrest@0: ; lon = fspan(-179.5,179.5,nlon)
forrest@0: ; lat = fspan( -89.5, 89.5,nlat)
forrest@0:  
forrest@0:   do i = 0,nlon-1
forrest@0:      lon(i) = -179.5 + i*(360./nlon)
forrest@0:   end do
forrest@0:  
forrest@0:   do j = 0,nlat-1
forrest@0: ;    lat(j) = -89.5 + j*(180./nlat)
forrest@0:      lat(j) =  89.5 - j*(180./nlat)
forrest@0:   end do
forrest@0: 
forrest@0: ; print (lon)
forrest@0: ; print (lat)        
forrest@0: ;=============================
forrest@0:     ;  create lat and long coordinate variables
forrest@0:     ;============================                                         
forrest@0:      lon!0          = "lon"
forrest@0:      lon@long_name  = "lon"
forrest@0:      lon@units      = "degrees-east"
forrest@0:      lon&lon        = lon
forrest@0: 
forrest@0:      lat!0          = "lat"
forrest@0:      lat@long_name  = "lat"
forrest@0:      lat@units      = "degrees_north"
forrest@0:      lat&lat        = lat
forrest@0:     ;=============================
forrest@0:     ;  name dimensions of t and assign coordinate variables
forrest@0:     ;============================  
forrest@0:      t!0    = "lat"
forrest@0:      t!1    = "lon"
forrest@0:      t&lat  = lat
forrest@0:      t&lon  = lon
forrest@0:      t@long_name = "amazon mask"
forrest@0: ;    t@units     = ""
forrest@0:      t@_FillValue= -999         
forrest@0:      t@missing_value= -999 
forrest@0:         
forrest@0:   c->mask_amazon  = t
forrest@0:   c->lat  = lat
forrest@0:   c->lon  = lon
forrest@0: 
forrest@0: ;************************************************
forrest@0: ; create default plot
forrest@0: ;************************************************
forrest@0: 
forrest@0:   wks = gsn_open_wks("ps","xy")         ; open a ps file
forrest@0:   gsn_define_colormap(wks,"gui_default")     ; choose colormap
forrest@0: 
forrest@0:   res                     = True             ; Use plot options
forrest@0:   res@cnFillOn            = True             ; Turn on color fill
forrest@0:   res@gsnSpreadColors      = True            ; use full colormap
forrest@0: ; res@cnFillMode          = "RasterFill"     ; Turn on raster color
forrest@0: ; res@lbLabelAutoStride   = True
forrest@0:   res@cnLinesOn           = False            ; Turn off contourn lines
forrest@0:   res@mpFillOn            = False            ; Turn off map fill
forrest@0:   res@gsnAddCyclic        = False
forrest@0:   res@tiMainString        = "LC15_amazon_biomass"
forrest@0: 
forrest@0:   res@mpMinLatF            = -21.1      ; range to zoom in on
forrest@0:   res@mpMaxLatF            =  13.8
forrest@0:   res@mpMinLonF            =  277.28
forrest@0:   res@mpMaxLonF            =  326.43
forrest@0: 
forrest@0:   res@gsnSpreadColors     = True             ; use full colormap
forrest@0:   res@cnLevelSelectionMode = "ManualLevels"  ; Manual contour invtervals
forrest@0: ; res@cnMinLevelValF       = 1.              ; Min level
forrest@0: ; res@cnMaxLevelValF       = 11.             ; Max level
forrest@0: ; res@cnLevelSpacingF      = 1.              ; interval
forrest@0:   res@cnMinLevelValF       = 0.              ; Min level
forrest@0:   res@cnMaxLevelValF       = 1.             ; Max level
forrest@0:   res@cnLevelSpacingF      = 0.1              ; interval
forrest@0: ; pres                            = True     ; panel plot mods desired
forrest@0: ; pres@gsnMaximize                = True     ; fill the page
forrest@0: 
forrest@0:   plot = gsn_csm_contour_map_ce(wks,t,res)   ; for observed
forrest@0: 
forrest@0: end