[ncl-talk] query

[Alan Brammer]

It seems you didn't copy over my last additional line.

You need to be aware of whether metadata is being transferred in the ncl
routines and check your array before passing them to a plotting routine.

copyVarCoords(data, land_only)
On 30 Aug 2014 00:01, "Kunal Bali" <kunal.bali9 at gmail.com> wrote:

> Dear Alan
>
> please find the attached file . Here you can see all scripts with out put
> result. The out put is not coming within the coordinates. How to fix it ??
>
> Kunal Bali
> Research Scholar
> Radio & Atmospheric Science Division
> CSIR - National Physical Laboratory
> New Delhi - 110012
>
>
>
>
>
>
>
> On Sat, Aug 30, 2014 at 3:40 AM, Alan Brammer <abrammer at albany.edu> wrote:
>
>> https://www.ncl.ucar.edu/Document/Functions/Shea_util/landsea_mask.shtml
>>
>> https://www.ncl.ucar.edu/Document/Functions/Built-in/mask.shtml
>>
>> This is based on a 1x1 land sea file, if you have higher res data then
>> you might want a more accurate mask.
>>
>> 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" load
>> "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl
>> <https://www.ncl.ucar.edu/Document/Functions/Shea_util/shea_util.shtml>"
>>
>> ncl 2> begin
>> ncl 3> in = addfile("/home/kunal/14-01. <http://14-01.nc/>nc
>> <http://14-01.nc/>","r")
>> ncl 4> data = in->TotCH4_A(:,:)
>> ncl 5>
>> printVarSummary(data) ; check coords are named lat lon
>>
>> a =*addfile*
>> <https://www.ncl.ucar.edu/Document/Functions/Built-in/addfile.shtml>
>> ("$NCARG_ROOT/lib/ncarg/data/cdf/landsea.nc","r") lsdata = a->LSMASK
>> lsm =*landsea_mask*(lsdata,data&lat,data&lon)
>> land_only =*mask*
>> <https://www.ncl.ucar.edu/Document/Functions/Built-in/mask.shtml>
>> (data,lsm.ge.1,True)
>> copyVarCoords(data, land_only) ; in case mask doesn't do this
>>
>> ncl 7>
>> ncl 8> wks = gsn_open_wks("x11","mask")
>> ncl 9> gsn_define_colormap(wks,"BlAqGrYeOrRe")
>> ncl 10> res = True
>> ncl 11> res at cnFillOn = True
>> ncl 12> res at cnLinesOn = False
>> ncl 13> res at gsnSpreadColors = True
>> ncl 14> res at gsnSpreadColorStart = 10
>> ncl 15> res at gsnSpreadColorEnd = 96
>> ncl 16> res at cnLevelSpacingF = 3
>> ncl 17> res at lbLabelString = 4
>> ncl 18> res at tiMainString = "Land Only"
>> ncl 19> plot = gsn_csm_contour_map_ce(wks,land_only,res)
>> ncl 20> end
>> On 29 Aug 2014 15:26, "Kunal Bali" <kunal.bali9 at gmail.com> wrote:
>>
>>> Could anyone provdies me the scripts for attached file.
>>>
>>> If i have variable TotCH4_A , Latitude, Longitude. So what could be the
>>> scripts that one. If you write here that would be better because i have
>>> tried the example from NCL tutorials but getting errors
>>>
>>> SO please provides me here if you can. Or i can send one file if anyone
>>> interested ??
>>>
>>> Kunal Bali
>>> Research Scholar
>>> Radio & Atmospheric Science Division
>>> CSIR - National Physical Laboratory
>>> New Delhi - 110012
>>>
>>>
>>>
>>>
>>>
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> 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/20140830/f1302b68/attachment.html