[求助][已解决]自定义经纬度数据的陆地填充问题

xuebiz

数据格式：
variables:
    double LON4(LON4 ;
        LON4:long_name = "Longitude" ;
        LON4:units = "degrees_east" ;
        LON4:modulo = 360. ;
        LON4:point_spacing = "even" ;
        LON4:axis = "X" ;
    double LAT4(LAT4) ;
        LAT4:long_name = "Latitude" ;
        LAT48:units = "degrees_north" ;
        LAT4:point_spacing = "even" ;
        LAT4:axis = "Y" ;
    double TIME(TIME) ;
        TIME:long_name = "Center time of the day" ;
        TIME:units = "days since 1978-01-01 00:00:00" ;
        TIME:time_origin = "01-JAN-1978 00:00:00" ;
        TIME:axis = "T" ;
    float SST(TIME, ZLEV, LAT4, LON4) ;
        SST:_FillValue = -9.99f ;
        SST:long_name = "Daily sea surface temperature" ;
        SST:units = "degrees C" ;
data:
LON4 = 117.125, 117.375, 。。。。。。
LAT4 = 32.125, 32.375,。。。。。。
---------------------------------------------------------------------------------
NCL程序：
  f = addfile("../data/SST_YS.nc", "r")
  wks = gsn_open_wks("png","sst")
   lon = f->LON4
   lat = f->LAT4
   sst = f->SST(0,0,:,:)
   res = True
  res@cnFillOn    = True                        ; turn on color fill
  res@cnLinesOn   = False                       ; turn off the contour lines
  res@cnLineLabelsOn = False                    ; turn the line labels off
  res@cnInfoLabelOn = False
  plot = gsn_csm_contour(wks,sst,res)
---------------------------------------------------------------------------------
出图如下：


求助：怎么增加XY轴的经纬度标示个数
          怎么添加陆地填充

海盗船长

plot = gsn_csm_contour_map(wks,sst,res)

xuebiz

海盗船长 发表于 2016-4-11 14:57
plot = gsn_csm_contour_map(wks,sst,res)

蠢死了，自己  一开始用的plot = gsn_csm_contour，后来换用了plot = gsn_csm_contour_map，这个测试里居然忘了改

还应该再加上下面这些，搞定！

  res@gsnAddCyclic=False              ;regional data

  res@mpMinLonF = 117.125
  res@mpMaxLonF = 127.875
  res@mpMinLatF = 32.125
  res@mpMaxLatF = 41.875

xuebiz

经度只有一个，看着很别扭，用什么参数可以增加标注，
这个数据的经纬度不是我最常见的那种，唉。。。

还有个比它更麻烦的经纬度的，等我搞定这个再来研究那个。。。