[ncl-talk] Panel with Observed & Model Plots.

Alan Brammer abrammer at albany.edu
Tue Nov 17 10:59:17 MST 2015


You already know how to define a region in a plotted map. ? 

From your previous script 

 res at mpFillColor     = False 
 res at mpMinLatF       = 0.0  ; The data is to be Selected for Model Domain but not Happening
 res at mpMaxLatF       = 30.
 res at mpMinLonF       = 60.0
 res at mpMaxLonF       = 100.0






> On 17 Nov 2015, at 11:38, Geeta Geeta <geetag54 at yahoo.com> wrote:
> 
> Hi Alan. 
> thank you. 
> I defined the long variable and then defined rf_model variable once again. 
> 
> This is what I have done. 
>     latitude  = f->XLAT(0,:,0)
>     longitude = f->XLONG(0,0,:)
> 
>  latitude at units= "degree_north"
>  longitude at units= "degree_east"
>    print(latitude)
>    print("longitude =   "+longitude)
> 
> 
>    rf_model   = rainc28+rainnc28-rainc4-rainnc4
> ;   printVarSummary(rf_model)
> ;   print(dimsizes(rf_model))
> ;   print(getvardims(rf_model))
>    rf_model!0    = "latitude"
>    rf_model!1    = "longitude"
>    rf_model&latitude    = latitude
> 
> The print Summary of rf_model is Okay. 
> (137)    longitude =   79.8386
> (138)    longitude =   79.8836
> (139)    longitude =   79.9286
> (140)    longitude =   79.9735
> (141)    longitude =   80.0185
> (142)    longitude =   80.0635
> (143)    longitude =   80.1084
> (144)    longitude =   80.1534
> (145)    longitude =   80.1984
> (146)    longitude =   80.2434
> (147)    longitude =   80.2883
> (148)    longitude =   80.3333
> (149)    longitude =   80.3783
> (150)    longitude =   80.4233
> (151)    longitude =   80.4682
> (152)    longitude =   80.5132
> 
> Variable: rf_model
> Type: float
> Total Size: 80784 bytes
>             20196 values
> Number of Dimensions: 2
> Dimensions and sizes:    [latitude | 132] x [longitude | 153]
> Coordinates: 
>             latitude: [9.848549..15.59277]
>             longitude: [73.67728..80.51321]
> (0)    latitude
> (1)    longitude
> 
> Then I use gsn_csm_contour_map(wks,rf_model,res_model) I again get a map that is Global and has  data for my domain on that. 
> I also tried with res_model at gsnAddCyclic = False
> I am attaching both . Kindly Suggest.
>  
> Geeta.
> <trmm-model-v2.ncl><trmm_model-v2.ps>_______________________________________________
> 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/20151117/169028d7/attachment.html 


More information about the ncl-talk mailing list