2016-12-15 5 views
1

enter image description hereWie werden die Kanten des Plots in der NCAR-Befehlssprache (NCL) geglättet?

Wie die Ränder des Grundstücks durch gsn_csm_contour_map glatt stattdessen erstellt, um jene blauen Quadraten zu verwenden?

Ich habe verschiedene Füllmodus versucht, keiner von ihnen geholfen. Gibt es in ncl eine aktive Konturfunktion? Ist es möglich, eine Kontur zwischen Null und Nicht-Null-Wert zu erstellen?

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" 
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" 
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" 

;====================================================================== 
; The main code 
;====================================================================== 

begin 
;---Read desired data 
    sfile = addfile(f,"r") 
    var = sfile->var1 

;---For zooming in on map 
    minlat = foo 
    maxlat = bar 
    minlon = foo1 
    maxlon = bar1 

;---Get dimentions 
    dims = dimsizes(var) 
    nlev = dims(0) 
    time = var&time 
    date_str_i = getDate(time) 

;---Set some resources 
    res     = True 

    [email protected]   = True 
    ;[email protected]  = "RasterFill" 
    ;[email protected] =True 
    [email protected]   = False 
    [email protected] = True 

    [email protected] = "ManualLevels" 
    [email protected]  = -100 
    [email protected]  = 3000 
    [email protected]  = 200 ; 300 ; 50 ; 150 

    [email protected]   = minlat 
    [email protected]   = maxlat 
    [email protected]   = minlon 
    [email protected]   = maxlon 
    [email protected] = "HighRes" 

    [email protected]  = True 
    [email protected] = 0.005 
    [email protected] = 0.001 
    [email protected]  = (minlon+maxlon)*0.5 
    [email protected]  = (minlat+maxlat)*0.5 

    [email protected] = "Always" 

    [email protected] = 0.01 

    [email protected]  = False  ; this is regional data 

;---Loop across each level and plot to a different PNG file every time 
    do n=4,nlev-1 
    wks_type = "png" 
    wks_type[email protected] = 2000 
    [email protected] = 2000 
     wks = gsn_open_wks(wks_type,fname(0)+sprinti("%03i",n)) ; 
     [email protected] = "Time:" + date_str_i(n) 
     [email protected] = 0.010 
     plot = gsn_csm_contour_map_ce(wks,var(n,:,:),res) 
     delete(wks) 

    end do 
end 

Aber ich bezweifle, dass es Code-bezogen ist.

+0

Sie sollten den Code hinzufügen, die Sie verwendet haben, um dieses Ergebnis auf Ihre Frage zu erhalten. –

+0

danke, ich werde es hinzufügen. –

Antwort

Verwandte Themen