load  "/home/server/student/homes/naraligi/contributed/gsn_code.ncl"
load  "/home/server/student/homes/naraligi/contributed/gsn_csm.ncl"
load  "/home/server/student/homes/naraligi/contributed/shea_util.ncl"
load  "/home/server/student/homes/naraligi/contributed/contributed.ncl"



begin

   rnccf=asciiread("ltrncc.dat",(/49,1/),"float")
   lat1 = (/-8.0,-5.0,-2.0,0,2.0,5.0,8.0/)
   lon1 = (/137.0,147.0,156.0,165.0,180.0,190.0,205.0/)
  
   print(dimsizes(rnccf))
   ccf = new((/7,7/),"float")
   gg=0
   do lt=0,6
	  do ln=0,6
		  print(gg) 
		  ccf(lt,ln)=rnccf(gg,0)
		  gg=gg+1
      end do
   end do


   ccf!0="lat"
   ccf!1="lon"

   ccf&lat = (/-8.0,-5.0,-2.0,0,2.0,5.0,8.0/)
   ccf&lon = (/137.0,147.0,156.0,165.0,180.0,190,205/)

    ccf@_Fillvalue = -999.0
	wks = gsn_open_wks("ps","proj")
	gsn_define_colormap(wks,"gui_default")


	res  = True
    res@mpProjection  = "Mercator"
	res@mpGridAndLimbon = False
	res@mpPerimOn       = True
	res@mpFillOn        = True
	res@mpCenterLonF    = 180

	res@cnFillOn        = True
	res@cnLineLabelsOn  = False
	res@cnLinesOn       = False
	res@lbLabelBarOn    = True

	res@mpLimitMode     = "LatLon"
	res@mpMinLatF       = -8
	res@mpMaxLatF       = 8
	res@mpMinLonF       = 130
	res@mpMaxLonF       = 205

	res@vpWidthF  = 0.8
	res@vpHeightF = 1.0

	res@pmTickMarkDisplayMode = "Always"
	res@pmTickMarkZone        = 0.5

    res@gsnAddCyclic          = False

	plot = gsn_csm_contour_map(wks,ccf,res)
	draw(plot)


end