NCL绘图

zys_cdut

大家好，我有这样一个问题想请教。
我是做深时气候模拟的，现在我有输出数据是180 Ma年前的180ma.cam.h0.0001-01.nc文件，现在我想做关于T、TS、U等变量的图，但是我绘制的图他不显示古地理的陆地轮廓，而显示现代的地理轮廓，我该怎么解决呢？

zys_cdut

begin
;************************************************
; read in netCDF file
;************************************************
  a = addfile("180ma.cam.h0.0001-01.nc","r")
;************************************************
; read in zonal winds
;************************************************
  t = a->TS(0,:,:)                                ; read July zonal winds
;************************************************
; reshape the data to 2D if needed
;************************************************
  ; Assuming TS is a 4D array (time, lev, lat, lon), you can reshape it to 2D
  ; by selecting a specific time step and flattening the lat and lon dimensions.
  ; If TS is already 2D, you can skip this step.
  ; t = reshape(t, dims(1), dims(2))  ; reshape t to be 2D

;************************************************
; create plot
;************************************************
  wks = gsn_open_wks("png","Surface_temperature")          ; send graphics to PNG file

  res                   = True
  res@mpProjection      = "Mollweide"       ; choose projection
  res@mpGridAndLimbOn   = True              ; turn on lat/lon lines
  res@mpPerimOn         = False             ; turn off box around plot
  res@mpGridLatSpacingF = 30.               ; spacing for lat lines
  res@mpGridLonSpacingF = 30.               ; spacing for lon lines
  res@mpFillOn          = False

  res@cnFillOn          = True              ; color plot desired
  res@cnFillPalette     = "gui_default"     ; set color map
  res@cnLineLabelsOn    = False             ; turn off contour lines
  res@txFontHeightF     = 0.015

  res@vpXF            = 0.1                 ; make plot bigger
  res@vpYF            = 0.9         
  res@vpWidthF        = 0.8
  res@vpHeightF       = 0.8

  res@lbLabelFontHeightF  = 0.015           ; label bar font height

  res@tiMainString       = "Example of Mollweide Projection"  ; add a title
  res@tiMainFontHeightF  = .018                               ; font height

  ; Turn off the boundary outlines of the modern map
  res@mpOutlineOn         = False          ; turn off the boundary outlines
  res@mpLandFillColor     = "white"        ; set land fill color to white or transparent
  res@mpOceanFillColor    = "lightblue"    ; set ocean fill color
  
  contour = gsn_csm_contour_map(wks,t,res)  ; create the plot
end

AuuuTi

用gsn_add_shapefile_polylines这个类似函数去画特定的shpfile（比如你要的古地理的陆地轮廓，把这个.shp文件读进去）
https://www.ncl.ucar.edu/Applications/shapefiles.shtml
https://www.ncl.ucar.edu/Documen ... ile_polylines.shtml

AuuuTi

用gsn_add_shapefile_polylines这个类似函数去画特定的shpfile（比如你要的古地理的陆地轮廓，把这个.shp文件读进去） https://www.ncl.ucar.edu/Applications/shapefiles.shtml https://www.ncl.ucar.edu/Document/Graphics/Interfaces/gsn_add_shapefile_polylines.shtml

zys_cdut

但是我没有这样的一个古地理轮廓shpfile，可以利用已有的古地理topo文件.nc吗