[ncl-talk] Custom panel plot

Appo derbetini appopson4 at gmail.com
Wed May 23 00:20:07 MDT 2018


Hi Mary,

See attached data and script we are using.
Regards,

Appo

2018-05-22 22:32 GMT+01:00 Mary Haley <haley at ucar.edu>:

> Hi James and Appo,
>
> Both of you have created some nice scripts that would be useful for the
> panel examples page. Would one or both of you be willing to contribute your
> script (and data file so I can run it) so we can add it to our page?
>
> We do ask for a clean script, which means one that has unsed commented
> code removed, somewhat consistent indentation, no references to personal
> directories, etc.  I'm happy to help clean it up as well, if there are
> things that you need help improving.
>
> If you are able to contribute the script and data, feel free to do this
> offline if you don't want to share the data publicly. We don't make the
> data file available, but I do like to have it so I can make sure the script
> runs with the latest version of NCL and that we create a nice PNG image
> from it.
>
> One small tip that might help with the subtitles not quite lining up.  If
> you use "degK" instead of just "K" for gsnRightString, then your titles
> will line up better because both titles will now have a character ('p' in
> the left title and 'g' in the right title) that goes below the text base
> line. If you don't want to repeat this title for every single plot, then I
> suggest setting this title just once at the top. See script panel_41.ncl
> for an example on how to do this:
>
> http://www.ncl.ucar.edu/Applications/panel.shtml#ex41
>
> Regards,
>
> --Mary
>
>
> On Tue, May 22, 2018 at 3:11 AM, Appo derbetini <appopson4 at gmail.com>
> wrote:
>
>> Hi James
>> I am on the right now as you can see on the attached plot.
>>
>> I am also sending script i used.
>>
>> Thank you
>>
>>
>>
>>
>>    fname = "tas_2006.nc"
>>     f     = addfile(fname,"r")
>>         t       = f->tas(:, {-15:15}, {5:35})
>>
>>     lon = t&lon
>>     lat = t&lat
>>
>>
>>
>>    pltDir = "./"                ; where output will be written
>>    pltName= "panel29"
>>    pltType= "png"               ; ps, eps, png, pdf, x11 .....
>>
>> ;***********************************************************
>> *****************
>> ; Miscellaneous
>> ;***********************************************************
>> *****************
>>   nrow = 3
>>   ncol = 4
>>
>>
>>
>> ;***********************************************************
>> *****************
>> ; Graphic resources
>> ;***********************************************************
>> *****************
>>   pltPath = pltDir+pltName
>>   wks     = gsn_open_wks(pltType ,pltPath)
>>
>>   res                       = True     ; plot mods desired
>>   res at gsnDraw               = False
>>   res at gsnFrame              = False
>>
>>   res at cnFillOn              = True     ; turn on color fill
>>   res at cnLinesOn             = False    ; turn of contour lines
>>   res at cnLineLabelsOn        = False    ; turn of contour line labels
>>   ;res at cnLevelSelectionMode  = "ManualLevels"    ; set manual contour
>> levels
>>
>>   res at lbLabelBarOn          = False
>>
>>
>>
>>  res at mpProjection                = "CylindricalEquidistant"       ;
>> projection
>>  res at mpLimitMode                 = "Corners"        ; method to zoom
>>  res at mpLeftCornerLatF            = min(lat)
>>  res at mpLeftCornerLonF            = min(lon)
>>  res at mpRightCornerLatF           = max(lat)
>>  res at mpRightCornerLonF           = max(lon)
>>  res at tfDoNDCOverlay              = True             ; do no
>>  res at gsnAddCyclic                = False            ; regional data,
>> don't add pt
>>
>>
>>   res at tmYLLabelsOn          = False    ; turn off lat labels
>>   res at tmXBLabelsOn          = False    ; turn off lon labels
>>
>>   res at mpDataBaseVersion     = "MediumRes"
>>   res at mpFillOn              = True
>>   res at mpOutlineOn           = True
>>   ;res at mpFillColors          = (/-1,0,-1,-1/)       ; mask out all
>> non-land areas
>>   res at mpFillOn                    = False        ; turn off gray fill
>>
>>   res at mpFillDrawOrder       = "PostDraw"
>>   res at mpProjection          = "CylindricalEquidistant"
>>
>>
>>   ;Ressource column 1
>>    res1 = True
>>    res1 = res
>>    res1 at cnFillPalette        = "WhiteYellowOrangeRed"
>>    res1 at lbOrientation      = "Horizontal"        ; orientation
>>    ;res1 at pmLabelBarWidthF   = 0.8
>>    ;res1 at pmLabelBarHeightF  = 0.03
>>
>>    res1 at lbLabelBarOn       = False
>>    res1 at lbBoxEndCapStyle      = "TriangleBothEnds"
>>
>>
>>    ;Ressource columns 2 & 3
>>      res2 = True
>>      res2 = res
>>      res2 at cnFillPalette        = "MPL_YlGn"
>>      res2 at lbOrientation      = "Horizontal"        ; orientation
>>      res2 at pmLabelBarWidthF   = 1.1
>>      ;res2 at pmLabelBarOrthogonalPosF  = -0.02
>>      res2 at pmLabelBarParallelPosF  = -0.02
>>
>>      res2 at lbBoxEndCapStyle      = "TriangleBothEnds"
>>
>>
>>     ;Ressource column 4
>>      res3 = True
>>      res3 = res
>>      res3 at cnFillPalette        = "MPL_RdPu"
>>      res3 at lbOrientation      = "Horizontal"        ; orientation
>>      ;res3 at pmLabelBarWidthF   = 0.8
>>
>>      res3 at lbBoxEndCapStyle      = "TriangleBothEnds"
>>
>>
>> ;***********************************************************
>> *****************
>> ; Generate panel entries
>> ;***********************************************************
>> *****************
>>                                ; only one time per file
>>   plot = new(nrow*ncol , "graphic")
>>
>>   ;res at cnFillPalette        = "WhiteYellowOrangeRed"
>>     ;res1 at tmYLLabelsOn          = True    ; turn off lat labels
>>    plot(0) = gsn_csm_contour_map(wks, t(0, :, :), res1)
>>
>>      res2 at tmYLLabelsOn          = False    ; turn off lat labels
>>    plot(1) = gsn_csm_contour_map(wks, t(1, :, :), res2)
>>    plot(2) = gsn_csm_contour_map(wks, t(2, :, :), res2)
>>    plot(3) = gsn_csm_contour_map(wks, t(3, :, :), res3)
>>
>>
>>     ;res1 at tmYLLabelsOn          = True
>>     plot(4) = gsn_csm_contour_map(wks, t(4, :, :), res1)
>>     res2 at tmYLLabelsOn          = False    ; turn off lat labels
>>     plot(5) = gsn_csm_contour_map(wks, t(5, :, :),  res2)
>>     plot(6) = gsn_csm_contour_map(wks, t(6, :, :), res2)
>>     plot(7) = gsn_csm_contour_map(wks, t(7, :, :),  res3)
>>
>>
>>     res1 at lbLabelBarOn          = True
>>     ;res2 at tmYLLabelsOn          = True
>>     ;res1 at tmXBLabelsOn          = True    ; turn off lon labels
>>     ;res2 at tmXBLabelsOn          = True    ; turn off lon labels
>>
>>     plot(8) = gsn_csm_contour_map(wks, t(8, :, :), res1)
>>     res2 at tmYLLabelsOn          = False    ; turn off lat labels
>>     plot(9) = gsn_csm_contour_map(wks, t(9, :, :),  res2)
>>
>>       res2 at lbLabelBarOn       = True
>>     plot(10) = gsn_csm_contour_map(wks, t(10, :, :),res2)
>>
>>     res3 at lbLabelBarOn       = True
>>    plot(11) = gsn_csm_contour_map(wks, t(11, :, :), res3)
>>
>>
>>
>>
>>
>> ;***********************************************************
>> *****************
>> ; Draw the panel and label bars
>> ;***********************************************************
>> *****************
>>   resP = True
>>   ;resP at gsnPanelMainString = "Land Only: Tmax: 2 Colorbars"
>>   resP at gsnMaximize = True
>>   resP at gsnFrame    = False
>>  ;resP at gsnPanelRight  = 0.98
>>   ;resP at gsnPanelBottom = 0.045
>>
>>    resP at gsnPanelFigureStrings= (/"1","1","1","1","2","2","2","2","3","3","3","3"/)
>> ; add strings to panel
>>     resP at amJust   = "BottomLeft"
>>
>>   gsn_panel(wks,plot,(/nrow,ncol/),resP)  ; plots
>>
>>
>>
>>   frame(wks)
>>
>>
>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180523/6241952e/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: tas_01.nc
Type: application/x-netcdf
Size: 996760 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180523/6241952e/attachment-0001.nc>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: panel_29.ncl
Type: application/octet-stream
Size: 6088 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180523/6241952e/attachment-0001.obj>


More information about the ncl-talk mailing list