[ncl-talk] Managing terrain height in the script wrf_CrossSection_add_smooth_terrain4.ncl

Wed Jun 1 23:46:00 MDT 2016

Hi Mary,
It have been using gsn_csm_xy to plot terrain .  The script is in the
attachment.

After your suggestion, I replaced gsn_csm_xy  by gsn_csm_contour but
the plot didnot show any terrain contour (see gsn_csm_contour.png in
the attachment).

I just  replaced   " contour_ter =
gsn_csm_xy(wks,X_plane,ter_plane,opts_ter) "   with

" contour_ter = gsn_csm_contour(wks,ter_plane,opts_ter)"     in the script.

Thanks.

Sang.

On Sat, May 28, 2016 at 1:00 AM, Mary Haley <haley at ucar.edu> wrote:
>
> Again, without a script, we can't tell what the problem might be.  Your tickmarks are pointing in, which is strange.  This seems to indicate that you are plotting your data with gsn_contour, and not gsn_csm_contour.
>
> One thing might be to try using gsn_csm_contour, and setting res at gsnYAxisIrregular2Linear = True
>
> http://www.ncl.ucar.edu/Document/Graphics/Resources/gsn.shtml#gsnYAxisIrregular2Linear
>
> If you continue to have issues with this, *please* include a script. We might also have to ask for the data.
>
> --Mary
>
>
> On Tue, May 24, 2016 at 11:39 PM, sangeeta maharjan <sangeetamaharjan at gmail.com> wrote:
>>
>> Hello Xiaoming and  Mary,
>>
>> Thanks for your reply.
>> I am aware that it is a west-east cross section plot. I tried but could not add the minor tick marks in vertical axes. The values in ter_plane used in the scripts yields (max value is red  in color (21th value)
>>
>> Variable: ter_plane
>> Type: float
>> Total Size: 272 bytes
>>             68 values
>> Number of Dimensions: 1
>> Dimensions and sizes:   [Horizontal | 68]
>> Coordinates:
>> Number Of Attributes: 3
>>   _FillValue :    9.96921e+36
>>   units :   m
>>   description :   Terrain Height
>> (0)   537.3745
>> (1)   556.2093
>> (2)   556.2093
>> (3)   552.6631
>> (4)   551.8929
>> (5)   588.4611
>> (6)   640.3762
>> (7)   713.1735
>> (8)   804.7409
>> (9)   854.6913
>> (10)  849.4119
>> (11)  805.8835
>> (12)  763.5001
>> (13)  735.0173
>> (14)  856.4299
>> (15)  1169.83
>> (16)  1492.322
>> (17)  1652.575
>> (18)  1694.905
>> (19)  1765.424
>> (20)  1889.178
>> (21)  1945.099
>> (22)  1869.276
>> (23)  1755.734
>> (24)  1730.924
>> (25)  1704.013
>> (26)  1546.452
>> (27)  1411.18
>> (28)  1363.609
>> (29)  1243.166
>> (30)  1051.444
>> (31)  936.7072
>> (32)  941.271
>> (33)  928.4125
>> (34)  861.8909
>> (35)  960.158
>> (36)  1142.136
>> (37)  1091.726
>> (38)  925.554
>> (39)  876.4838
>> (40)  894.1756
>> (41)  927.5086
>> (42)  1036.723
>> (43)  1188.228
>> (44)  1198.813
>> (45)  1186.619
>> (46)  1269.593
>> (47)  1163.258
>> (48)  913.6985
>> (49)  847.5559
>> (50)  987.7681
>> (51)  1175.981
>> (52)  1240.032
>> (53)  1103.708
>> (54)  881.4406
>> (55)  781.1543
>> (56)  893.9818
>> (57)  1047.922
>> (58)  1127.229
>> (59)  1131.302
>> (60)  1081.057
>> (61)  976.4305
>> (62)  924.3822
>> (63)  904.832
>> (64)  861.3491
>> (65)  830.9657
>> (66)  798.5671
>> (67)  748.0306
>>
>> With reference to plot in the  attachment, the temperature  contour  left and right  of maximum peak is close to 1500 m so i fell that the the terrain height is reduced by some hundreds.
>>
>> Regards
>> Sang.
>>
>> On Mon, May 23, 2016 at 8:47 PM, xiaoming Hu <yuanfangcan at hotmail.com> wrote:
>>>
>>> I don't see anything wrong with your plot.
>>>
>>> Are you questioning the interpolated terrain height is wrong?
>>> Since your angle is 90degree, basically west-east cross section.   It is very easy to check the terrain height yourself to see if the interpolated terrain height by wrf_CrossSection_add_smooth_terrain4.ncl   is reasonable.
>>>
>>> Cheers!
>>> Xiaoming
>>>
>>> ________________________________
>>> Date: Sun, 22 May 2016 09:38:27 +0545
>>> From: sangeetamaharjan at gmail.com
>>> To: ncl-talk at ucar.edu
>>> Subject: [ncl-talk] Managing terrain height in the script wrf_CrossSection_add_smooth_terrain4.ncl
>>>
>>> Dear All,
>>>
>>> I have used the wrf_CrossSection_add_smooth_terrain4.ncl to study distribution of  TKE along the cross section. The problem is  the height of terrain seems to lesser than the actual one. For example, the maximum terrain height with in the cross section has to be  about 1900 meters but it is lesser by 300 meters in the plot generated by the script.
>>>
>>>
>>> Please find the sample of plot showing the reduced terrain height  for your kind reference in the attachment.
>>>
>>> I would appreciate any insight to manage the terrain heights in the script.
>>>
>>> Thanks & Regards.
>>>
>>> Sang,
>>>
>>> _______________________________________________ ncl-talk mailing list ncl-talk at ucar.edu List instructions, subscriber options, unsubscribe: http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ncltalk_sang.ncl
Type: application/octet-stream
Size: 7578 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160602/3d1fe645/attachment.obj 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: gsn_csm_contour.png
Type: image/png
Size: 46657 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160602/3d1fe645/attachment.png 