怎么用wrf_map_overlays口令画出墨卡托投影方式的图？

veroniaca

想要出一张具有墨卡托投影方式的温度图，但是设置了 mpres@mpProjection ="Mercator"，是wrf_map_overlays这个函数的问题吗？如果要换成墨卡托投影，该怎么设置呢？恳请大神指导指导代码如下：
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
begin
       filelist=systemfunc("ls /mnt/d/test/wrfout*")     
        f=addfiles(filelist,"r")
        temdata=wrf_user_getvar(f,"T2",-1)
        ave=dim_avg_n(temdata,0)
        aave=ave-273.15
;-- define the workstation
                cmap                                  =read_colormap_file("matlab_hot")
                wks_type                                   ="pdf"                       
                wks_type@wkPaperWidthF           =15                               
                wks_type@wkPaperHeightF           =10                       
                wks                                     =gsn_open_wks(wks_type,"/mnt/f/dasixia/data/wrf/air/air_indian")

        res                                           =True
        res@gsnMaximize                        =True                       
        res@gsnPaperOrientation             ="portrait"               
        res@gsnPaperMargin                        =0.2                       
        res@cnRasterSmoothingOn                =True         
        res@cnLineLabelsOn                   =False                                       
        res@cnInfoLabelOn                    =False                                       
        res@pmTickMarkDisplayMode      ="Conditional"                       
        res@tmXTOn                             =False
        res@tmYROn                             =False

        gsn_define_colormap(wks,"matlab_hot")
        res@cnLevelSelectionMode        ="ExplicitLevels"
        res@cnLevels                    =(/20,20.5,21,21.5,22,22.5,23,23.5,24,24.5,25,25.5,26,26.5,27,27.5,28,28.5,29,29.5,30,30.5/)   
        res@cnFillColors                =(/60,59,54,53,50,48,46,44,42,40,38,36,34,32,30,28,26,24,22,20,18,16,12/)                     
        res@cnFillOn                    =True
        res@cnLinesOn                   =False
        res@InitTime                    =False
        res@ValidTime                   =False
        res@Footer                       =False
        res@FieldTitle                        ="WRF Temperature of the Indian Ocean in Oct. 2015"
        contour                            =wrf_contour(f[0],wks,aave,res)

        mpres = True
        mpres@mpNationalLineColor       ="Black"
        mpres@mpDataSetName               ="Earth..4"
        mpres@mpOutlineBoundarySets       ="National"
        mpres@mpProjection                ="Mercator"
        mpres@mpDataBaseVersion           ="HighRes"
        mpres@mpLeftCornerLatF            =-5
        mpres@mpRightCornerLatF           =25
        mpres@mpLeftCornerLonF            =65
        mpres@mpRightCornerLonF           =100
        mpres@mpGridAndLimbOn             =True
        mpres@mpGridLatSpacingF           =2.5
        mpres@mpGridLonSpacingF           =2.5
        mpres@mpGridLineThicknessF        =0.5
        mpres@mpDataSetName               ="Earth..4"               
       mpres@mpLimitMode                 ="Corners"               
        mpres@mpOutlineOn                 =True
        over_id=wrf_map_overlays(f[0],wks,(/contour/),True,mpres)
end