[ncl-talk] Fw: scatsat-1 data plotting

[Mary Haley]

Please, this is the second time I'm asking you to post your questions to
ncl-talk and not directly to me. By emailing me, you are expecting me to
individually answer your questions. I have my own deadlines I'm trying to
meet next week.

In the script I sent you yesterday, I showed you how to zoom in on an area
of interest.  I also included a PDF showing the zoomed in area.

You said you wanted to plot both ascending and descending winds in one plot
as a contour plot, and then that you wanted to plot over India. But, the
data for descending wind speed appears to not be over India (see the
attached asc_des_wspd_maps.000002.png file).

I have provided two more scripts, but I am not going to help any more
beyond this.  Again, I ask you to look these scripts over carefully. You
must try to make an effort to understand what the script is doing and
modify it for your purpose.

If you have further questions, start a new email thread and send the email
to ncl-talk at ucar.edu. You may need to resubscribe to ncl-talk, because as I
said yesterday, you do not appear in the membership list.

--Mary

On Tue, Jan 30, 2018 at 11:40 PM, jayanti pal <jiban_samudra18 at yahoo.com>
wrote:

> Thanks you for your reply. I don't know why my mail is not posted in
> ncl_talk though I have subscribe it.
>
> Sorry for wrong attachments. I have modified code that you have send me
> previously. This code has been working properly and it plots one variable
> (Ascending wind speed)over global region.
>
> But I want to plot over Indian region. Will it be possible to plot two
> variables (Ascending and descending wind speed) as contour in a single plot?
>
> Sincerely
>
>
> Dr. Jayanti Pal
> Post-doc Fellow
> Department of Atmospheric Sciences
> Institute of Atmospheric and Environmental Sciences
> University of Calcutta
> Contact: jiban_samudra18 at yahoo.com, jayantibright at gmail.com
> Ph. 8013906997 <(801)%20390-6997>
>
>
> On Tuesday 30 January 2018, 10:35:58 PM IST, Mary Haley <haley at ucar.edu>
> wrote:
>
>
> Please subscribe to ncl-talk if you want help with NCL, and send your
> questions there. I did not see your "jiban_samudra18 at yahoo.com" in the
> ncl-talk membership list.
>
> Dave, Rick, and I all sent you scripts a few weeks ago, but it looks like
> you are not using the fixed scripts we sent you. The script you sent has
> several things wrong with it, that we had fixed for you previously:
>
> [1] You have several resources misspelled:
>
>  res at cnfillon=True
>  rescnLineson=False
>  res at lbLabAutoStride=True
>
> should be:
>
>  res at cnFillOn=True
>  rescnLinesOn=False
>  res at lbLabelAutoStride=True
> ​     ; This is the default, so you don't need to set it.​
>
>
> ​[2] You didn't use the contour levels that we created for you before.  It
> looks like you have unequally spaced contour levels, so you must provide
> these for NCL.​
>
> res at cnLevelSelectionMode = "ExplicitLevels"
> res at cnLevels             = (/1,1.5,2,2.25,2.5,2.75,3,3.5,
> 4,4.5,5.2,5.6,6,6.4,6.8,7.2,7.6,8.5,11.8,12.4,14.0,14.5,15,30/)
>
>
> [3] You created the lat/lon coordinates, but you didn't attach them to the
> variable:
>
> lon = (360./numlon)*(ispan(0,numlon-1,1)+0.5)
> lat = (180./numlat)*(ispan(0,numlat-1,1)+0.5)-90
>
> lat at units="degrees_north"
> lon at units="degrees_east"
> data_actual!0 = "lat"
> data_actual!1 = "lon"
> data_actual&lat = lat
> data_actual&lon = lon
>
> Note that you don't need a do loop to calculate lat/lon.  This is another
> thing we tried to show you in a previous script.
>
> [4] If you want the value 0.0 to be treated as missing, use the "where"
> function:
>
> data_actual at _FillValue = -999.
> data_actual = where(data_actual.eq.0,data_actual at _FillValue,data_actual)
>
> ​I've attached the modified script and the resultant PNG. I didn't have
> time to clean it up, so I'll leave that to you. Please look over this
> script carefully and try to understand it this time.
>
> As for adding vectors, you will need to add code like this at the end of
> the contouring code, and you need to make sure u and v have lat/lon
> attached to them the same way I did for asc_wind_speed in the attached
> script.
>
> vcres = True
> ;;vcres at vcMinDistanceF  = 0.005    ; thins the vectors
> ;;vcres at vcRefLengthF    = 0.01      ; changes the NDC length of the
> reference vector
>
> vector_plot=gsn_csm_vector(xwks,u,v,vcres)
>
> overlay(plot,vector_plot)
> draw(plot)
> frame(xwks)
>
> To see more about customizing vectors, see the examples at:
>
> http://www.ncl.ucar.edu/Applications/vector.shtml
>
>
> I don't know what u and v are on the file, so this is something you will
> figure out yourself.
>
> If you have more questions, please subscribe to ncl-talk and post a new
> message.
>
> Thanks for your consideration, and good luck.
>
>
> --Mary
>
>
>
> On Tue, Jan 30, 2018 at 4:04 AM, jayanti pal <jiban_samudra18 at yahoo.com>
> wrote:
>
>
>
> Dr. Jayanti Pal
> Post-doc Fellow
> Department of Atmospheric Sciences
> Institute of Atmospheric and Environmental Sciences
> University of Calcutta
> Contact: jiban_samudra18 at yahoo.com, jayantibright at gmail.com
> Ph. 8013906997
>
>
> ----- Forwarded message -----
> *From:* jayanti pal <jiban_samudra18 at yahoo.com>
> *To:* ncl-talk at ucar.edu <ncl-talk at ucar.edu>
> *Sent:* Tuesday 30 January 2018, 4:24:31 PM IST
> *Subject:* Fw: scatsat-1 data plotting
>
> Sir
>
> I am new to ncl. I tried to plot scatsat-1 data with ncl. But figure is
> not correctly plotted.
> I have attached the code I have used, Hdf file, actual plot (in color) and
> ncl plot (in black).
> In this hdf file there are four variables: Asending wind speed, asending
> wind direction , descending wind speed, desending wind direction.
> I want to plot this four variable in a single plot where speeds are in
> contour and directions are in vector. Also I want to plot this over Indian
> region (Lat: 0-40 and Lon : 60-100).
>
> Please help me. I have to plot those data urgently as I have to submit
> project review.
>
> Sincerely
>
>
>
> Dr. Jayanti Pal
> Post-doc Fellow
> Department of Atmospheric Sciences
> Institute of Atmospheric and Environmental Sciences
> University of Calcutta
> Contact: jiban_samudra18 at yahoo.com, jayantibright at gmail.com
> Ph. 8013906997
>
>
> ----- Forwarded message -----
> *From:* jayanti pal <jiban_samudra18 at yahoo.com>
> *To:* ncl-talk at ucar.edu <ncl-talk at ucar.edu>
> *Sent:* Tuesday 2 January 2018, 12:47:33 PM IST
> *Subject:* Fw: scatsat-1 data plotting
>
> Sir
>
> I am new to ncl. I tried to plot scatsat-1 data with ncl. But figure is
> not correctly plotted.
> I have attached the code I have used, Hdf file, actual plot (in color) and
> ncl plot (in black).
>
> Please help me. I have to plot those data urgently.
>
>
>
> Dr. Jayanti Pal
> Post-doc Fellow
> Department of Atmospheric Sciences
> Institute of Atmospheric and Environmental Sciences
> University of Calcutta
> Contact: jiban_samudra18 at yahoo.com, jayantibright at gmail.com
> Ph. 8013906997
>
>
> ----- Forwarded message -----
> *From:* jayanti pal <jiban_samudra18 at yahoo.com>
> *To:* eoshelp at hdfgroup.org <eoshelp at hdfgroup.org>
> *Sent:* Friday 22 December 2017, 12:28:59 PM IST
> *Subject:* scatsat-1 data plotting
>
> Sir
>
> I am new to ncl. I tried to plot scatsat-1 data with ncl. But figure is
> not correctly plotted.
> I have attached the code I have used, Hdf file, actual plot (in color) and
> ncl plot (in black).
>
> Please help me. I have to plot those data urgently.
>
> Thanking you,
>
> Dr. Jayanti Pal
> Post-doc Fellow
> Department of Atmospheric Sciences
> Institute of Atmospheric and Environmental Sciences
> University of Calcutta
> Contact: jiban_samudra18 at yahoo.com, jayantibright at gmail.com
> Ph. 8013906997
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180131/5b35ce83/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: asc_des_wspd_overlay.png
Type: image/png
Size: 169495 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180131/5b35ce83/attachment-0003.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: asc_des_wspd_maps.000001.png
Type: image/png
Size: 161811 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180131/5b35ce83/attachment-0004.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: asc_des_wspd_maps.000002.png
Type: image/png
Size: 155335 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180131/5b35ce83/attachment-0005.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: plot_ncl_mod_asc_des_maps.ncl
Type: application/octet-stream
Size: 3862 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180131/5b35ce83/attachment-0002.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: plot_ncl_mod_asc_des_overlay.ncl
Type: application/octet-stream
Size: 4558 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180131/5b35ce83/attachment-0003.obj>