<div dir="ltr"><div dir="ltr">Hi Dennis,<div>Thank you for the reply and clarification as well.</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Nov 6, 2019 at 3:00 AM <<a href="mailto:ncl-talk-request@ucar.edu">ncl-talk-request@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Send ncl-talk mailing list submissions to<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:ncl-talk-request@ucar.edu" target="_blank">ncl-talk-request@ucar.edu</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:ncl-talk-owner@ucar.edu" target="_blank">ncl-talk-owner@ucar.edu</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of ncl-talk digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. Re: climo_3.ncl + loop + conventional subscripts (zoe jacobs)<br>
2. Re: Assignment type mismatch, right hand side can't be<br>
coerced to type of left hand side (zoe jacobs)<br>
3. convert string to fortran character (Zhifeng Yang)<br>
4. Remove Global mean SST from the SST data (dickson mbigi)<br>
5. Re: convert string to fortran character (Rick Brownrigg)<br>
6. Re: convert string to fortran character (Dennis Shea)<br>
7. Re: Remove Global mean SST from the SST data (Dennis Shea)<br>
8. Need a data file of NCL examples (Feiyu Wang)<br>
9. Re: Amplitude and Phase of the seasonal cycle<br>
(Sri.durgesh Nandini-Weiss)<br>
10. Re: convert string to fortran character (Zhifeng Yang)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Tue, 5 Nov 2019 00:59:53 +0330<br>
From: zoe jacobs <<a href="mailto:zoejacobs1990@gmail.com" target="_blank">zoejacobs1990@gmail.com</a>><br>
To: Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>>, <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
Subject: Re: [ncl-talk] climo_3.ncl + loop + conventional subscripts<br>
Message-ID:<br>
<CAL-eqBk1biqZmXhjo80FgsYWZ6BrBUEUZMq=<a href="mailto:DeipmdYqjukZ0g@mail.gmail.com" target="_blank">DeipmdYqjukZ0g@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Dear Dennis,<br>
It is ok now. Many thanks for your help and advice.<br>
Best wishes,<br>
<br>
On Mon, Nov 4, 2019 at 6:13 PM Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:<br>
<br>
> A quick overview of the script. It will 'work '<br>
><br>
> However, I would suggest you use<br>
> res@cnLevelSelectionMode = "ExplicitLevels"<br>
> and<br>
> res@cnFillMode = "RasterFill"<br>
> ===<br>
> Also, your script will need to be changed to accommodate the precipitation<br>
> variable plot values<br>
> There are many examples of how to do this.<br>
><br>
> eg: *<a href="http://www.ncl.ucar.edu/Applications/HiResPrc.shtml*" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Applications/HiResPrc.shtml*</a><br>
> <<a href="http://www.ncl.ucar.edu/Applications/HiResPrc.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Applications/HiResPrc.shtml</a>><br>
><br>
> On Sat, Nov 2, 2019 at 11:26 PM zoe jacobs <<a href="mailto:zoejacobs1990@gmail.com" target="_blank">zoejacobs1990@gmail.com</a>><br>
> wrote:<br>
><br>
>> Dear Dennis,<br>
>> Many thanks for your help. I used the attached script. I just want to<br>
>> check with you to make sure that the script is correct. Would you please<br>
>> check it?<br>
>> Best wishes,<br>
>><br>
>> On Mon, Oct 28, 2019 at 12:31 AM Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:<br>
>><br>
>>> The following is untested. However, it gives you the idea.<br>
>>> It is your responsibility to look up documentation and look at the<br>
>>> printed output.<br>
>>><br>
>>><br>
>>><br>
>>> fils = systemfunc ("ls precip.*.nc") ; file paths<br>
>>> a = addfiles (fils, "r")<br>
>>> printVarSummary (a)<br>
>>><br>
>>> ;************************************************<br>
>>> ; Read the file<br>
>>> ;************************************************<br>
>>><br>
>>> ymdStrt = 19790101 ; climatology start<br>
>>> year/month/day<br>
>>> ymdLast = 20151231 ; " last<br>
>>> year/month/day<br>
>>> time = a[:]->time<br>
>>> yyyymmdd= cd_calendar(time, -2) ; note -2<br>
>>> print(yyyymmdd)<br>
>>> ntStrt = ind(yyyymmdd.eq.ymdStrt) ; start time index<br>
>>> ntLast = ind(yyyymmdd.eq.ymdLast) ; last time index<br>
>>> print("ntStrt="+ntStrt+" ntLast="+ntLast)<br>
>>> prc = a[:]->precip(ntStrt:ntLast,:,:) ; read precip<br>
>>> from all files<br>
>>><br>
>>> printVarSummary (prc)<br>
>>> printMinMax (prc,0)<br>
>>> print("=====")<br>
>>> ; you decide what you need<br>
>>> prcMonth = *calculate_monthly_values(*<br>
>>> <<a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/calculate_monthly_values.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Contributed/calculate_monthly_values.shtml</a>>prc,<br>
>>> "avg", 0, False) ; monthly mean<br>
>>> ;;prcMonth = *calculate_monthly_values(*<br>
>>> <<a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/calculate_monthly_values.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Contributed/calculate_monthly_values.shtml</a>>prc,<br>
>>> "sum", 0, False) ; monthly total<br>
>>><br>
>>> printVarSummary (prcMonth)<br>
>>> printMinMax (prcMonth,0)<br>
>>> print("=====")<br>
>>> opt = True ; examine data distribute<br>
>>> [information only]<br>
>>> opt@PrintStat = True<br>
>>> stat_prc = stat_dispersion(prcMonth, opt )<br>
>>><br>
>>> prcClm = *clmMonTLL*<br>
>>> <<a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/clmMonTLL.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Contributed/clmMonTLL.shtml</a>>(prcMonth)<br>
>>> ; Compute monthly climatology<br>
>>> printVarSummary (prcClm)<br>
>>> printMinMax (prcClm,0)<br>
>>><br>
>>> ====<br>
>>> Add the plot code.<br>
>>> Look carefully at: <a href="http://www.ncl.ucar.edu/Applications/HiResPrc.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Applications/HiResPrc.shtml</a><br>
>>><br>
>>> EG: The following may be inappropriate for your dataset but this is the<br>
>>> idea.<br>
>>><br>
>>> res@cnLevelSelectionMode = "ExplicitLevels"<br>
>>> res@cnLevels = (/ 0.01, 0.02, 0.04, 0.08, 0.16, \<br>
>>> 0.32, 0.64, 0.96/)<br>
>>> res@cnFillColors = (/"white","cyan", "green","yellow",\<br>
>>> "darkorange","red","magenta","purple",\<br>
>>> "black"/)<br>
>>><br>
>>><br>
>>> On Sun, Oct 27, 2019 at 2:48 PM Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:<br>
>>><br>
>>>><br>
>>>> [0]<br>
>>>> As noted on ncl-talk many time the first rule of data processing is<br>
>>>> *Look at your data!"<br>
>>>> The temperature data were *monthly means* [K].<br>
>>>> The precipitation data are* daily* totals [mm/day]. How do I know? I<br>
>>>> looked at the README file at the URL. You should do that.<br>
>>>><br>
>>>> [1]<br>
>>>> a =* addfiles*<br>
>>>> <<a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/addfiles.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Built-in/addfiles.shtml</a>><br>
>>>> (fils, "r")<br>
>>>><br>
>>>> You have the following which is appropriate for an *addfile*<br>
>>>> <<a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/addfile.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Built-in/addfile.shtml</a>><br>
>>>> reference<br>
>>>> time = a->time<br>
>>>><br>
>>>> However, you are using *addfiles*<br>
>>>> <<a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/addfiles.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Built-in/addfiles.shtml</a>><br>
>>>> It should be:<br>
>>>> time = a*[:]*->time ; *[:]* syntax is appropriate for an<br>
>>>> *addfiles* reference<br>
>>>><br>
>>>> [2]<br>
>>>> Did you read the documentation for *cd_calendar *<br>
>>>> <<a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/cd_calendar.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Built-in/cd_calendar.shtml</a>>as<br>
>>>> I suggested? The whole point of using *cd_calendar* is to return<br>
>>>> 'human-readable' start and end times from which the actual start/stop time<br>
>>>> index values could be determined.. Your script uses the actual time values<br>
>>>> [ 692496 , 1016808 ] mixed with human-readable yyyymm. Further, you used<br>
>>>><br>
>>>> yyyymm = cd_calendar(time, *-1*) ; -1 returns YYYYMM<br>
>>>><br>
>>>> The precipitation files contain [as noted] daily values. You must use a<br>
>>>> human readable time like YYYYMMDD [eg: 19790101]. The us<br>
>>>><br>
>>>> yyyymmdd = cd_calendar(time, *-2*) ; -2 returns YYYYMMDD<br>
>>>><br>
>>>> This converts the actual time values [eg: 692496] to<br>
>>>> 'human-readable' yyyymmdd [year=>yyyy, month=>mm, day=dd>*<br>
>>>> 692496=>197901*01] values.<br>
>>>><br>
>>>><br>
>>>> ===<br>
>>>> Also, your script will need to be changed to accommodate the<br>
>>>> precipitation variable plot values<br>
>>>> There are many examples of how to do this.<br>
>>>><br>
>>>> eg: *<a href="http://www.ncl.ucar.edu/Applications/HiResPrc.shtml*" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Applications/HiResPrc.shtml*</a><br>
>>>> <<a href="http://www.ncl.ucar.edu/Applications/HiResPrc.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Applications/HiResPrc.shtml</a>><br>
>>>> ===<br>
>>>> In a previous email, you mentioned that you converted degrees Kelvin to<br>
>>>> C.<br>
>>>> It is best to self-document the variable with its new units.<br>
>>>><br>
>>>> tmp = tmp-273.15<br>
>>>> tmp@units = "degC" ; this will be automatically plotted<br>
>>>> ===<br>
>>>> Good Luck<br>
>>>><br>
>>>><br>
>>>><br>
>>>><br>
>>>><br>
>>>><br>
>>>> On Sat, Oct 26, 2019 at 8:19 AM zoe jacobs <<a href="mailto:zoejacobs1990@gmail.com" target="_blank">zoejacobs1990@gmail.com</a>><br>
>>>> wrote:<br>
>>>><br>
>>>>> Dear all,<br>
>>>>> I would like to have climatology precipitation like temperature<br>
>>>>> climatology discussed through previous emails. However, there are<br>
>>>>> precipitation files for each year from 1979- 2015 (<br>
>>>>> <a href="ftp://ftp.cdc.noaa.gov/Datasets/cpc_global_precip/" rel="noreferrer" target="_blank">ftp://ftp.cdc.noaa.gov/Datasets/cpc_global_precip/</a>). In the first<br>
>>>>> step, I tried to merge all those files as one, and then use Dennis script<br>
>>>>> to plot a map. I faced some errors and need your advice please. the script<br>
>>>>> I am running has been attached.Many thanks in advance,<br>
>>>>> Best regardsn Sat, Oct 26, 2019 at 3:20 PM zoe jacobs <<br>
>>>>> <a href="mailto:zoejacobs1990@gmail.com" target="_blank">zoejacobs1990@gmail.com</a>> wrote:<br>
>>>>><br>
>>>>>> I just noticed that I could use res@gsnRightString = " deg C" .<br>
>>>>>> Once again, many thanks for your guide.<br>
>>>>>> Best wishes,<br>
>>>>>><br>
>>>>>><br>
>>>>>> On Fri, Oct 25, 2019 at 11:07 PM Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:<br>
>>>>>><br>
>>>>>>> I forgot to mention that plots can be done in 'portrait' [default]<br>
>>>>>>> or 'landscape' mode.<br>
>>>>>>><br>
>>>>>>> You should experiment to see which is appropriate for your needs.<br>
>>>>>>><br>
>>>>>>> resP@gsnPaperOrientation = "landscape" ; "portrait" is<br>
>>>>>>> default<br>
>>>>>>><br>
>>>>>>> do nmo=0,11 ; loop over the<br>
>>>>>>> months<br>
>>>>>>> res@gsnLeftString = months(nmo)<br>
>>>>>>> plot(nmo) = gsn_csm_contour_map(wks,tmpClm(nmo,:,:), res) ;<br>
>>>>>>> create plot<br>
>>>>>>> end do<br>
>>>>>>><br>
>>>>>>> if (resP@gsnPaperOrientation .eq."landscape") then<br>
>>>>>>><br>
>>>>>>> gsn_panel(wks,plot,(/3,4/),resP) ; 3 rows<br>
>>>>>>> x 4 columns<br>
>>>>>>> else<br>
>>>>>>> gsn_panel(wks,plot,(/4,3/),resP) ; portrait<br>
>>>>>>> ; 4 rows x 3 columns<br>
>>>>>>> end if<br>
>>>>>>><br>
>>>>>>> On Fri, Oct 25, 2019 at 12:15 PM Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:<br>
>>>>>>><br>
>>>>>>>> Hello,<br>
>>>>>>>><br>
>>>>>>>> We do not usually do this but this is offline from ncl-talk.<br>
>>>>>>>><br>
>>>>>>>> [1]<br>
>>>>>>>> Conventional subscripts start at 1 for Fortran and Matlab. Hence,<br>
>>>>>>>> for 12 months [nmos=12]<br>
>>>>>>>> Fortran: do nmo=1,nmos<br>
>>>>>>>> Matlab: for 1:nmos or 1:1:nmos [I am not a Matlab user so<br>
>>>>>>>> this is a guess.]<br>
>>>>>>>> NCL 'conventional' subscripts start at 0 [like C/C++/IDL/Python]<br>
>>>>>>>> NCL: do nmo=0,ntim-1<br>
>>>>>>>><br>
>>>>>>>> [2] Dave illustrated 'coordinate subscripting' which uses the *{*<br>
>>>>>>>> ...*}* syntax.<br>
>>>>>>>><br>
>>>>>>>> The script I am attaching uses conventional subscripting. Just like<br>
>>>>>>>> Matlab/Fortran except it uses a range fro 0 to 11 [12 elements].<br>
>>>>>>>><br>
>>>>>>>> [3] Please carefully read the documentation for *clmMonTLL*<br>
>>>>>>>> <<a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/clmMonTLL.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Contributed/clmMonTLL.shtml</a>>,<br>
>>>>>>>> *cd_calendar*<br>
>>>>>>>> <<a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/cd_calendar.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Built-in/cd_calendar.shtml</a>><br>
>>>>>>>> ,* ind*<br>
>>>>>>>> <<a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/ind.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Built-in/ind.shtml</a>><br>
>>>>>>>><br>
>>>>>>>> [4] Look at the output from '*printVarSummary*<br>
>>>>>>>> <<a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/printVarSummary.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Built-in/printVarSummary.shtml</a>>'.<br>
>>>>>>>> USE *printVarSummary *frequently.<br>
>>>>>>>> Note the dimensions and added attributes:<br>
>>>>>>>> time_op_ncl : Climatology: 29 years<br>
>>>>>>>> info : function clmMonTLL: contributed.ncl<br>
>>>>>>>><br>
>>>>>>>> [5] NCL offer a large number of color tables [palettes]:<br>
>>>>>>>><br>
>>>>>>>> *<a href="http://www.ncl.ucar.edu/Document/Graphics/color_table_gallery.shtml*" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Document/Graphics/color_table_gallery.shtml*</a><br>
>>>>>>>> <<a href="http://www.ncl.ucar.edu/Document/Graphics/color_table_gallery.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Document/Graphics/color_table_gallery.shtml</a>><br>
>>>>>>>> Using a color map that emphasizes features can be very useful.<br>
>>>>>>>><br>
>>>>>>>> [6] You must invest the time to learn any new language.<br>
>>>>>>>> Karin Meier-Fliesher and Michael Bottinger [DKRZ] wrote a wonderful<br>
>>>>>>>> tutorial.<br>
>>>>>>>><br>
>>>>>>>> *<a href="http://www.ncl.ucar.edu/Document/Manuals/NCL_User_Guide/*" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Document/Manuals/NCL_User_Guide/*</a><br>
>>>>>>>> <<a href="http://www.ncl.ucar.edu/Document/Manuals/NCL_User_Guide/" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Document/Manuals/NCL_User_Guide/</a>><br>
>>>>>>>><br>
>>>>>>>> I suggest you read it.<br>
>>>>>>>><br>
>>>>>>>> Good luck<br>
>>>>>>>> ==============<br>
>>>>>>>> *%>* ncl zoe_jacobs.hgcn_cams.ncl<br>
>>>>>>>><br>
>>>>>>>> will produce a png file.<br>
>>>>>>>><br>
>>>>>>>><br>
>>>>>>>><br>
>>>>>>>> On Fri, Oct 25, 2019 at 11:16 AM Dave Allured - NOAA Affiliate via<br>
>>>>>>>> ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
>>>>>>>><br>
>>>>>>>>> 1. In your example, please notice that the month loop has a<br>
>>>>>>>>> stride of 3. Therefore it is making plots for only four months: January,<br>
>>>>>>>>> April, July, and October.<br>
>>>>>>>>><br>
>>>>>>>>> If you want all 12 months, remove the stride. Then make sure the<br>
>>>>>>>>> array "plot" is dimensioned 12. I think you can interchangeably use either<br>
>>>>>>>>> a 1-D or 2-D plot array (i.e. 3*4), your choice.<br>
>>>>>>>>><br>
>>>>>>>>> The secondary index "i" is used to write the four plots into<br>
>>>>>>>>> positions 0, 1, 2, 3 in the "plot" graphics array. If you are making 12<br>
>>>>>>>>> plots, you do not really need to use a secondary index.<br>
>>>>>>>>><br>
>>>>>>>>> Panel plotting uses a graphics array containing multiple plots.<br>
>>>>>>>>> Please study basic examples and documentation for making panel plots.<br>
>>>>>>>>><br>
>>>>>>>>> 2. That data file has full coordinates. Therefore you can use<br>
>>>>>>>>> either conventional or coordinate subscripting, or mixed, your choice.<br>
>>>>>>>>><br>
>>>>>>>>> I suggest using one of the date functions with coordinate<br>
>>>>>>>>> subscripting, to index the time subset that you want. Something like this:<br>
>>>>>>>>><br>
>>>>>>>>> time_units = f->air&time@units<br>
>>>>>>>>> time1 = cd_inv_calendar (year1, 1, 1, 0, 0, 0, time_units, 0)<br>
>>>>>>>>> time2 = cd_inv_calendar (year2, 12, 31, 23, 0, 0, time_units,<br>
>>>>>>>>> 0)<br>
>>>>>>>>> air_subset = f->air({time1:time2},:,:)<br>
>>>>>>>>><br>
>>>>>>>>> Use printVarSummary to ensure that the subset has the dimension<br>
>>>>>>>>> sizes that you expect. Then you can proceed to compute the climatology.<br>
>>>>>>>>><br>
>>>>>>>>><br>
>>>>>>>>> On Fri, Oct 25, 2019 at 10:14 AM zoe jacobs via ncl-talk <<br>
>>>>>>>>> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
>>>>>>>>><br>
>>>>>>>>>> Hi all,<br>
>>>>>>>>>> regarding climo_3.ncl (<br>
>>>>>>>>>> <a href="https://www.ncl.ucar.edu/Applications/Scripts/climo_3.ncl" rel="noreferrer" target="_blank">https://www.ncl.ucar.edu/Applications/Scripts/climo_3.ncl</a>) ,<br>
>>>>>>>>>> I have 2 questions:<br>
>>>>>>>>>> 1. I need to plot all months on one panel (say 3*4), and cannot<br>
>>>>>>>>>> understand the logic behind the below loop, which used in the climo_3.ncl<br>
>>>>>>>>>> script :<br>
>>>>>>>>>><br>
>>>>>>>>>> i = -1 ; Climatologies<br>
>>>>>>>>>> do nmo=0,11,3 ; loop over the months<br>
>>>>>>>>>> i = i+1<br>
>>>>>>>>>> res@gsnCenterString = months(nmo)+":"+time(0)/100 +"-"+ time(ntim-1)/100<br>
>>>>>>>>>> plot(i) = gsn_csm_contour_map(wks,prcClm(nmo,:,:), res) ; create plot<br>
>>>>>>>>>> end do<br>
>>>>>>>>>><br>
>>>>>>>>>> How does it work??/<br>
>>>>>>>>>><br>
>>>>>>>>>> 2. I would like to show climotology temperature from 1987- 2015. Data which I am using is <a href="http://air.mon.mean.nc" rel="noreferrer" target="_blank">air.mon.mean.nc</a> ( <a href="https://www.esrl.noaa.gov/psd/data/gridded/data.ghcncams.html" rel="noreferrer" target="_blank">https://www.esrl.noaa.gov/psd/data/gridded/data.ghcncams.html</a>) and it is using conventional subscripts. I am not familiar with conventional subscripts. So how can I convert coordinate subscripting to conventional subscripts?!!<br>
>>>>>>>>>><br>
>>>>>>>>>> Please kindly advice me .<br>
>>>>>>>>>><br>
>>>>>>>>>> Many thanks in advance,<br>
>>>>>>>>>><br>
>>>>>>>>>> Best regards,<br>
>>>>>>>>>><br>
>>>>>>>>>> _______________________________________________<br>
>>>>>>>>> ncl-talk mailing list<br>
>>>>>>>>> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
>>>>>>>>> List instructions, subscriber options, unsubscribe:<br>
>>>>>>>>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
>>>>>>>><br>
>>>>>>>><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191105/1b9f6820/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191105/1b9f6820/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Tue, 5 Nov 2019 01:44:32 +0330<br>
From: zoe jacobs <<a href="mailto:zoejacobs1990@gmail.com" target="_blank">zoejacobs1990@gmail.com</a>><br>
To: Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>>, <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
Subject: Re: [ncl-talk] Assignment type mismatch, right hand side<br>
can't be coerced to type of left hand side<br>
Message-ID:<br>
<<a href="mailto:CAL-eqBkqvmkZRkLp76M0asTFBjnudRT9-J%2Be5HvoXCk0atQZrg@mail.gmail.com" target="_blank">CAL-eqBkqvmkZRkLp76M0asTFBjnudRT9-J+e5HvoXCk0atQZrg@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Dear Rick,<br>
Many thanks for your help. it is working now. However, I noticed when I<br>
plot the data for two different dates ( 1979-2015 & 2000-2015), I get the<br>
same result. What is my mistake? Would you please advise me in this regard?<br>
Best wishes,<br>
<br>
On Mon, Nov 4, 2019 at 5:55 PM Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:<br>
<br>
> Hi,<br>
><br>
> Docs for clmMonTLL says that the input variable can be of type float or<br>
> "numeric", but it really must be either float or double. Notice that your<br>
> variable "u" is of type short, and that further there are "add_offset" and<br>
> "scale_factor" attributes that must be applied to the data to get the<br>
> actual values. "short2flt()" is exactly the function you want; I think this<br>
> will do it:<br>
><br>
> uClm = clmMonTLL (short2flt(u))<br>
> vClm = clmMonTLL (short2flt(v))<br>
><br>
> <a href="http://ncl.ucar.edu/Document/Functions/Contributed/short2flt.shtml" rel="noreferrer" target="_blank">http://ncl.ucar.edu/Document/Functions/Contributed/short2flt.shtml</a><br>
><br>
> There are numerous examples of how to write NetCDF file at the following<br>
> link. I recommend the "inefficient" way. Its much easier and not that<br>
> inefficient unless you are dealing with lots of really big variables or are<br>
> writing files on a production basis.<br>
><br>
> <a href="http://ncl.ucar.edu/Applications/write_netcdf.shtml" rel="noreferrer" target="_blank">http://ncl.ucar.edu/Applications/write_netcdf.shtml</a><br>
><br>
> Hope that helps...<br>
> Rick<br>
><br>
><br>
> On Mon, Nov 4, 2019 at 5:42 AM zoe jacobs via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
> wrote:<br>
><br>
>> Dear All,<br>
>> I attempted to plot climatological wind speed during 1979-2015. But I<br>
>> faced an error and have no idea from where this problem arose. ( I checked<br>
>> the NCL website but was not enough:<br>
>> <a href="https://www.ncl.ucar.edu/Document/Language/error_messages.shtml#AssignMismatch" rel="noreferrer" target="_blank">https://www.ncl.ucar.edu/Document/Language/error_messages.shtml#AssignMismatch</a><br>
>> ).Please have a look at the attached script, error message, and data file.<br>
>> I also have another question, after having climatological wind, how can I<br>
>> have the result as nc file?<br>
>> I am waiting for your helpful suggestions.<br>
>> Many thanks in advance.<br>
>> adaptor.mars.internal-1572848815.075499-2543-15...<br>
>> <<a href="https://drive.google.com/file/d/1DAUj5O8_HRB80p6f2ieiubdjqGlUokTL/view?usp=drive_web" rel="noreferrer" target="_blank">https://drive.google.com/file/d/1DAUj5O8_HRB80p6f2ieiubdjqGlUokTL/view?usp=drive_web</a>><br>
>> ,<br>
>> _______________________________________________<br>
>> ncl-talk mailing list<br>
>> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
>> List instructions, subscriber options, unsubscribe:<br>
>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
><br>
><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191105/ccd3c492/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191105/ccd3c492/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Mon, 4 Nov 2019 21:40:53 -0500<br>
From: Zhifeng Yang <<a href="mailto:yangzf01@gmail.com" target="_blank">yangzf01@gmail.com</a>><br>
To: Talk NCL <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Subject: [ncl-talk] convert string to fortran character<br>
Message-ID:<br>
<CAF92WDMVbWr+S60Qn=5W1BCqvGpGNyBoBfCy_re=8HqO=<a href="mailto:QhGuQ@mail.gmail.com" target="_blank">QhGuQ@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Hi<br>
I am reading wrf output which has variable Times as character. When I read<br>
it using NCL, the Times variable is still character, not string. I want to<br>
change the Times variable using NCL. For example,<br>
<br>
Variable: times<br>
<br>
Type: character<br>
<br>
Total Size: 19 bytes<br>
<br>
19 values<br>
<br>
Number of Dimensions: 2<br>
<br>
Dimensions and sizes: [Time | 1] x [DateStrLen | 19]<br>
<br>
Coordinates:<br>
<br>
Number Of Attributes: 0<br>
<br>
(0,0) 2<br>
<br>
(0,1) 0<br>
<br>
(0,2) 1<br>
<br>
(0,3) 4<br>
<br>
(0,4) -<br>
<br>
(0,5) 0<br>
<br>
(0,6) 5<br>
<br>
(0,7) -<br>
<br>
(0,8) 2<br>
<br>
(0,9) 5<br>
<br>
(0,10) _<br>
<br>
(0,11) 0<br>
<br>
(0,12) 0<br>
<br>
(0,13) :<br>
<br>
(0,14) 0<br>
<br>
(0,15) 0<br>
<br>
(0,16) :<br>
<br>
(0,17) 0<br>
<br>
(0,18) 0<br>
<br>
I want to change 4 to 5. I tried<br>
<br>
times(0, 3) = "5"<br>
"5" is a string, not a character. I got the following error.<br>
<br>
fatal:Assignment type mismatch, right hand side can't be coerced to type of<br>
left hand side<br>
<br>
Do you have any suggestions? Thank you<br>
<br>
Zhifeng<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191104/fbfaace7/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191104/fbfaace7/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Tue, 5 Nov 2019 11:00:42 +0800<br>
From: dickson mbigi <<a href="mailto:dickson.mbigi@gmail.com" target="_blank">dickson.mbigi@gmail.com</a>><br>
To: <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
Subject: [ncl-talk] Remove Global mean SST from the SST data<br>
Message-ID:<br>
<CALUZNzO1rVr9izVxcjB=oeK8m7eLAqfp89fhpP=<a href="mailto:3EN2qYM7PLg@mail.gmail.com" target="_blank">3EN2qYM7PLg@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Dear concerned expert,<br>
I am preparing a script to perfom decadal MAM EOF over Indian oceans using<br>
ERSSTV5 data set.<br>
Steps prior to do EOF:<br>
1) Select the three months of March-May and do the seasonal mean;<br>
2) Do a 9-year running average (Lanczos filter)<br>
3) Do the linear regression on the global mean SST and keep the residual<br>
for latter EOF analysis (i.e. remove the global mean component at each grid<br>
point).<br>
<br>
I cannot figure out the third step on how to remove the global mean<br>
component at each grid point.<br>
<br>
Kindly I am writing to ask for any assistance or directive regarding the<br>
said concern.<br>
<br>
Regards,<br>
Dickson.<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191105/ece8bdc9/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191105/ece8bdc9/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 5<br>
Date: Mon, 4 Nov 2019 20:05:11 -0700<br>
From: Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>><br>
To: Zhifeng Yang <<a href="mailto:yangzf01@gmail.com" target="_blank">yangzf01@gmail.com</a>><br>
Cc: Talk NCL <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Subject: Re: [ncl-talk] convert string to fortran character<br>
Message-ID:<br>
<<a href="mailto:CAGKRhbEoazaqKupXVsgRCtL61rrLpcbPY9JA7cvxZok0ozJxUA@mail.gmail.com" target="_blank">CAGKRhbEoazaqKupXVsgRCtL61rrLpcbPY9JA7cvxZok0ozJxUA@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Hi,<br>
<br>
You can convert characters to strings via the tostring() function:<br>
<br>
<a href="http://ncl.ucar.edu/Document/Functions/Built-in/tostring.shtml" rel="noreferrer" target="_blank">http://ncl.ucar.edu/Document/Functions/Built-in/tostring.shtml</a><br>
<br>
With regard to "...change 4 to 5." I'm not sure what you are trying to do<br>
there, but NCL is strongly typed, in both data type and array dimensions.<br>
Perhaps use tostring() in conjunction with a new array name (?)<br>
<br>
Rick<br>
<br>
On Mon, Nov 4, 2019 at 7:41 PM Zhifeng Yang via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
wrote:<br>
<br>
> Hi<br>
> I am reading wrf output which has variable Times as character. When I read<br>
> it using NCL, the Times variable is still character, not string. I want to<br>
> change the Times variable using NCL. For example,<br>
><br>
> Variable: times<br>
><br>
> Type: character<br>
><br>
> Total Size: 19 bytes<br>
><br>
> 19 values<br>
><br>
> Number of Dimensions: 2<br>
><br>
> Dimensions and sizes: [Time | 1] x [DateStrLen | 19]<br>
><br>
> Coordinates:<br>
><br>
> Number Of Attributes: 0<br>
><br>
> (0,0) 2<br>
><br>
> (0,1) 0<br>
><br>
> (0,2) 1<br>
><br>
> (0,3) 4<br>
><br>
> (0,4) -<br>
><br>
> (0,5) 0<br>
><br>
> (0,6) 5<br>
><br>
> (0,7) -<br>
><br>
> (0,8) 2<br>
><br>
> (0,9) 5<br>
><br>
> (0,10) _<br>
><br>
> (0,11) 0<br>
><br>
> (0,12) 0<br>
><br>
> (0,13) :<br>
><br>
> (0,14) 0<br>
><br>
> (0,15) 0<br>
><br>
> (0,16) :<br>
><br>
> (0,17) 0<br>
><br>
> (0,18) 0<br>
><br>
> I want to change 4 to 5. I tried<br>
><br>
> times(0, 3) = "5"<br>
> "5" is a string, not a character. I got the following error.<br>
><br>
> fatal:Assignment type mismatch, right hand side can't be coerced to type<br>
> of left hand side<br>
><br>
> Do you have any suggestions? Thank you<br>
><br>
> Zhifeng<br>
><br>
><br>
> _______________________________________________<br>
> ncl-talk mailing list<br>
> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
> List instructions, subscriber options, unsubscribe:<br>
> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191104/0e3f70f8/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191104/0e3f70f8/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 6<br>
Date: Mon, 4 Nov 2019 21:25:03 -0700<br>
From: Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>><br>
To: Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>><br>
Cc: Talk NCL <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Subject: Re: [ncl-talk] convert string to fortran character<br>
Message-ID:<br>
<<a href="mailto:CAOF1d_4uMJ_VMpgAn6o_cXQwHjukpSnhrPJNmVFRth-tJiVz4w@mail.gmail.com" target="_blank">CAOF1d_4uMJ_VMpgAn6o_cXQwHjukpSnhrPJNmVFRth-tJiVz4w@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
The question is not clear. Possibly:<br>
<br>
times = f->Times<br>
times(0,3) = tochar(tostring(5))<br>
print(times)<br>
<br>
On Mon, Nov 4, 2019 at 8:10 PM Rick Brownrigg via ncl-talk <<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
<br>
> Hi,<br>
><br>
> You can convert characters to strings via the tostring() function:<br>
><br>
> <a href="http://ncl.ucar.edu/Document/Functions/Built-in/tostring.shtml" rel="noreferrer" target="_blank">http://ncl.ucar.edu/Document/Functions/Built-in/tostring.shtml</a><br>
><br>
> With regard to "...change 4 to 5." I'm not sure what you are trying to do<br>
> there, but NCL is strongly typed, in both data type and array dimensions.<br>
> Perhaps use tostring() in conjunction with a new array name (?)<br>
><br>
> Rick<br>
><br>
> On Mon, Nov 4, 2019 at 7:41 PM Zhifeng Yang via ncl-talk <<br>
> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
><br>
>> Hi<br>
>> I am reading wrf output which has variable Times as character. When I<br>
>> read it using NCL, the Times variable is still character, not string. I<br>
>> want to change the Times variable using NCL. For example,<br>
>><br>
>> Variable: times<br>
>><br>
>> Type: character<br>
>><br>
>> Total Size: 19 bytes<br>
>><br>
>> 19 values<br>
>><br>
>> Number of Dimensions: 2<br>
>><br>
>> Dimensions and sizes: [Time | 1] x [DateStrLen | 19]<br>
>><br>
>> Coordinates:<br>
>><br>
>> Number Of Attributes: 0<br>
>><br>
>> (0,0) 2<br>
>><br>
>> (0,1) 0<br>
>><br>
>> (0,2) 1<br>
>><br>
>> (0,3) 4<br>
>><br>
>> (0,4) -<br>
>><br>
>> (0,5) 0<br>
>><br>
>> (0,6) 5<br>
>><br>
>> (0,7) -<br>
>><br>
>> (0,8) 2<br>
>><br>
>> (0,9) 5<br>
>><br>
>> (0,10) _<br>
>><br>
>> (0,11) 0<br>
>><br>
>> (0,12) 0<br>
>><br>
>> (0,13) :<br>
>><br>
>> (0,14) 0<br>
>><br>
>> (0,15) 0<br>
>><br>
>> (0,16) :<br>
>><br>
>> (0,17) 0<br>
>><br>
>> (0,18) 0<br>
>><br>
>> I want to change 4 to 5. I tried<br>
>><br>
>> times(0, 3) = "5"<br>
>> "5" is a string, not a character. I got the following error.<br>
>><br>
>> fatal:Assignment type mismatch, right hand side can't be coerced to type<br>
>> of left hand side<br>
>><br>
>> Do you have any suggestions? Thank you<br>
>><br>
>> Zhifeng<br>
>><br>
>><br>
>> _______________________________________________<br>
>> ncl-talk mailing list<br>
>> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
>> List instructions, subscriber options, unsubscribe:<br>
>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
><br>
> _______________________________________________<br>
> ncl-talk mailing list<br>
> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
> List instructions, subscriber options, unsubscribe:<br>
> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191104/f9677955/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191104/f9677955/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 7<br>
Date: Mon, 4 Nov 2019 21:30:25 -0700<br>
From: Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>><br>
To: dickson mbigi <<a href="mailto:dickson.mbigi@gmail.com" target="_blank">dickson.mbigi@gmail.com</a>><br>
Cc: Ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Subject: Re: [ncl-talk] Remove Global mean SST from the SST data<br>
Message-ID:<br>
<CAOF1d_4WD+gTnpF5sC4r2xK=<a href="mailto:7nTGMukWZmi1TKtD-UULqznzgA@mail.gmail.com" target="_blank">7nTGMukWZmi1TKtD-UULqznzgA@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
One approach:<br>
<br>
<br>
;--- Global: only 60S-60N<br>
<br>
sstGlb = f->SST(iStrt:iLast,{-60:60},:) ; read Global<br>
data (limit 60S to 60N)<br>
printVarSummary(sstGlb)<br>
<br>
sstGlbClm = clmMonTLL(sstGlb) ; climatology<br>
for clm period<br>
printVarSummary(sstGlbClm)<br>
<br>
sstGlbAnom= calcMonAnomTLL ( sstGlb , sstGlbClm ) ; anomalies at<br>
each grid point<br>
printVarSummary(sstGlbAnom) ; from base climatology<br>
<br>
<br>
On Mon, Nov 4, 2019 at 8:00 PM dickson mbigi via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
wrote:<br>
<br>
> Dear concerned expert,<br>
> I am preparing a script to perfom decadal MAM EOF over Indian oceans<br>
> using ERSSTV5 data set.<br>
> Steps prior to do EOF:<br>
> 1) Select the three months of March-May and do the seasonal mean;<br>
> 2) Do a 9-year running average (Lanczos filter)<br>
> 3) Do the linear regression on the global mean SST and keep the residual<br>
> for latter EOF analysis (i.e. remove the global mean component at each grid<br>
> point).<br>
><br>
> I cannot figure out the third step on how to remove the global mean<br>
> component at each grid point.<br>
><br>
> Kindly I am writing to ask for any assistance or directive regarding the<br>
> said concern.<br>
><br>
> Regards,<br>
> Dickson.<br>
><br>
> _______________________________________________<br>
> ncl-talk mailing list<br>
> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
> List instructions, subscriber options, unsubscribe:<br>
> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191104/0229480c/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191104/0229480c/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 8<br>
Date: Tue, 5 Nov 2019 15:37:24 +0800 (CST)<br>
From: "Feiyu Wang" <<a href="mailto:leightonw@126.com" target="_blank">leightonw@126.com</a>><br>
To: <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
Subject: [ncl-talk] Need a data file of NCL examples<br>
Message-ID: <<a href="mailto:3e952b3c.5190.16e3a7f309f.Coremail.leightonw@126.com" target="_blank">3e952b3c.5190.16e3a7f309f.Coremail.leightonw@126.com</a>><br>
Content-Type: text/plain; charset="gbk"<br>
<br>
<br>
<br>
Dear Sir or Madam, <br>
When using the ESMF_regrid.ncl, it dosen't work without a weight file. Would you please provide the wgt_file such as CAMSE_2_FV.nc and HOMME_to_FV.nc for me?<br>
Much thanks and best wishes.<br>
<br>
<br>
Feiyu Wang<br>
TEL:18681849808<br>
E-mail: <a href="mailto:leightonw@126.com" target="_blank">leightonw@126.com</a><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191105/6ce78c71/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191105/6ce78c71/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 9<br>
Date: Tue, 5 Nov 2019 12:06:01 +0100<br>
From: "Sri.durgesh Nandini-Weiss"<br>
<<a href="mailto:sri.durgesh.nandini-weiss@uni-hamburg.de" target="_blank">sri.durgesh.nandini-weiss@uni-hamburg.de</a>><br>
To: <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
Subject: Re: [ncl-talk] Amplitude and Phase of the seasonal cycle<br>
Message-ID: <<a href="mailto:4d29b832-6b26-67c3-4c3d-e063b83bf1bf@uni-hamburg.de" target="_blank">4d29b832-6b26-67c3-4c3d-e063b83bf1bf@uni-hamburg.de</a>><br>
Content-Type: text/plain; charset=utf-8; format=flowed<br>
<br>
<br>
On 11/4/19 3:27 PM, Sri.durgesh Nandini-Weiss via ncl-talk wrote:<br>
> Hello everyone,<br>
><br>
> I am learning, and trying to compute and plot the phase and amplitude <br>
> of the seasonal cycle, e.g. for Sea surface height over a period.<br>
><br>
> Can someone please let me know if the script fanal_4.ncl using Fourier <br>
> ( transformation) analysis is the correct way to do this?<br>
><br>
> Thanx in advance<br>
><br>
> Sri<br>
><br>
><br>
><br>
-- <br>
Dr. Sri Nandini-Wei?<br>
Research Scientist<br>
<br>
Universit?t Hamburg<br>
Center for Earth System Research and Sustainability (CEN)<br>
Cluster of Excellence 'Climate, Climatic Change, and Society' (CLICCS)<br>
<br>
Bundesstrasse 53, 20146 Hamburg<br>
Tel: +49 (0) 40 42838 7472<br>
<br>
<br>
<br>
------------------------------<br>
<br>
Message: 10<br>
Date: Tue, 5 Nov 2019 07:55:47 -0500<br>
From: Zhifeng Yang <<a href="mailto:yangzf01@gmail.com" target="_blank">yangzf01@gmail.com</a>><br>
To: Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>><br>
Cc: Talk NCL <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Subject: Re: [ncl-talk] convert string to fortran character<br>
Message-ID:<br>
<CAF92WDNjULUR-ZRjfkM+EusKa0yvnxnPqLdx07dQ9P=<a href="mailto:VNmz6Cw@mail.gmail.com" target="_blank">VNmz6Cw@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Hi Dennis and Rick,<br>
I found tochar will give me an array with 2 dimensions. Now I can use the<br>
following to get what I want<br>
times = f->Times<br>
times_tmp = tochar(tostring(5))<br>
times(0,3) = times_tmp(0)<br>
<br>
then it works. Thank you for your help.<br>
<br>
Best,<br>
Zhifeng<br>
<br>
On Mon, Nov 4, 2019 at 11:25 PM Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:<br>
<br>
> The question is not clear. Possibly:<br>
><br>
> times = f->Times<br>
> times(0,3) = tochar(tostring(5))<br>
> print(times)<br>
><br>
> On Mon, Nov 4, 2019 at 8:10 PM Rick Brownrigg via ncl-talk <<br>
> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
><br>
>> Hi,<br>
>><br>
>> You can convert characters to strings via the tostring() function:<br>
>><br>
>> <a href="http://ncl.ucar.edu/Document/Functions/Built-in/tostring.shtml" rel="noreferrer" target="_blank">http://ncl.ucar.edu/Document/Functions/Built-in/tostring.shtml</a><br>
>><br>
>> With regard to "...change 4 to 5." I'm not sure what you are trying to<br>
>> do there, but NCL is strongly typed, in both data type and array<br>
>> dimensions. Perhaps use tostring() in conjunction with a new array name (?)<br>
>><br>
>> Rick<br>
>><br>
>> On Mon, Nov 4, 2019 at 7:41 PM Zhifeng Yang via ncl-talk <<br>
>> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
>><br>
>>> Hi<br>
>>> I am reading wrf output which has variable Times as character. When I<br>
>>> read it using NCL, the Times variable is still character, not string. I<br>
>>> want to change the Times variable using NCL. For example,<br>
>>><br>
>>> Variable: times<br>
>>><br>
>>> Type: character<br>
>>><br>
>>> Total Size: 19 bytes<br>
>>><br>
>>> 19 values<br>
>>><br>
>>> Number of Dimensions: 2<br>
>>><br>
>>> Dimensions and sizes: [Time | 1] x [DateStrLen | 19]<br>
>>><br>
>>> Coordinates:<br>
>>><br>
>>> Number Of Attributes: 0<br>
>>><br>
>>> (0,0) 2<br>
>>><br>
>>> (0,1) 0<br>
>>><br>
>>> (0,2) 1<br>
>>><br>
>>> (0,3) 4<br>
>>><br>
>>> (0,4) -<br>
>>><br>
>>> (0,5) 0<br>
>>><br>
>>> (0,6) 5<br>
>>><br>
>>> (0,7) -<br>
>>><br>
>>> (0,8) 2<br>
>>><br>
>>> (0,9) 5<br>
>>><br>
>>> (0,10) _<br>
>>><br>
>>> (0,11) 0<br>
>>><br>
>>> (0,12) 0<br>
>>><br>
>>> (0,13) :<br>
>>><br>
>>> (0,14) 0<br>
>>><br>
>>> (0,15) 0<br>
>>><br>
>>> (0,16) :<br>
>>><br>
>>> (0,17) 0<br>
>>><br>
>>> (0,18) 0<br>
>>><br>
>>> I want to change 4 to 5. I tried<br>
>>><br>
>>> times(0, 3) = "5"<br>
>>> "5" is a string, not a character. I got the following error.<br>
>>><br>
>>> fatal:Assignment type mismatch, right hand side can't be coerced to type<br>
>>> of left hand side<br>
>>><br>
>>> Do you have any suggestions? Thank you<br>
>>><br>
>>> Zhifeng<br>
>>><br>
>>><br>
>>> _______________________________________________<br>
>>> ncl-talk mailing list<br>
>>> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
>>> List instructions, subscriber options, unsubscribe:<br>
>>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
>><br>
>> _______________________________________________<br>
>> ncl-talk mailing list<br>
>> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
>> List instructions, subscriber options, unsubscribe:<br>
>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
><br>
><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191105/cf1921e9/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191105/cf1921e9/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Subject: Digest Footer<br>
<br>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br>
------------------------------<br>
<br>
End of ncl-talk Digest, Vol 192, Issue 5<br>
****************************************<br>
</blockquote></div></div>