<div dir="ltr"> RE: DECADAL SVD <br><div> Thanks Dennis. After getting rid from @_FillValue of filtered time series, the decadal svd works well.</div><div><br></div><div>Regards</div><div>Dickson</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Dec 12, 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: period mean removed (Debashis Paul)<br>
2. Re: How to average monthly wind direction and represent them<br>
as wind barbs (Kunal Dayal)<br>
3. GFS interpolation (vertical and horizontal) (Ehsan Taghizadeh)<br>
4. DECADAL SVD (dickson mbigi)<br>
5. help in remove period mean (Debashis Paul)<br>
6. Re: DECADAL SVD (Dennis Shea)<br>
7. Re: How to average monthly wind direction and represent them<br>
as wind barbs (Dennis Shea)<br>
8. Re: GFS interpolation (vertical and horizontal) (Dennis Shea)<br>
9. Re: help in remove period mean (Dennis Shea)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Wed, 11 Dec 2019 11:34:22 +0530<br>
From: Debashis Paul <<a href="mailto:debgeo668@gmail.com" target="_blank">debgeo668@gmail.com</a>><br>
To: Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>><br>
Cc: Ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Subject: Re: [ncl-talk] period mean removed<br>
Message-ID:<br>
<<a href="mailto:CAFV8iJ%2B9t8CYMDG4RWxPxhQBt3V8EQ%2BzjBLqBr2DDP0AedLaHQ@mail.gmail.com" target="_blank">CAFV8iJ+9t8CYMDG4RWxPxhQBt3V8EQ+zjBLqBr2DDP0AedLaHQ@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Thanks a lot, sir will work accordingly.<br>
<br>
With regards,<br>
DP<br>
**************************<br>
*Debashis Paul*<br>
<br>
*Research Scholar, NIT-RKL*<br>
*Ph:+916002713628*<br>
*E-mail: <a href="mailto:Debgeo668@gmail.com" target="_blank">Debgeo668@gmail.com</a> <<a href="mailto:Debgeo668@gmail.com" target="_blank">Debgeo668@gmail.com</a>>*<br>
<br>
<br>
On Wed, Dec 11, 2019 at 5:04 AM Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:<br>
<br>
> Please see examples:<br>
><br>
> *<a href="http://www.ncl.ucar.edu/Applications/mjoclivar.shtml*" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Applications/mjoclivar.shtml*</a><br>
> <<a href="http://www.ncl.ucar.edu/Applications/mjoclivar.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Applications/mjoclivar.shtml</a>><br>
> mjovlivar_2.ncl<br>
><br>
> *<a href="http://www.ncl.ucar.edu/Applications/climo.shtml*" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Applications/climo.shtml*</a><br>
> <<a href="http://www.ncl.ucar.edu/Applications/climo.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Applications/climo.shtml</a>><br>
> climo_5.ncl<br>
><br>
> ===<br>
> You must do the work. Use the above to guide your thought.<br>
><br>
> [1] compute climatology<br>
> [2] Cpmpute anomalies<br>
> [3] perform running means<br>
> [4] Compute areal averages<br>
><br>
><br>
> On Tue, Dec 10, 2019 at 6:40 AM Debashis Paul via ncl-talk <<br>
> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
><br>
>> Hello NCL community,<br>
>> For my work, I need to do a task in which I have to<br>
>> calculate a 5-day running mean of daily 200-hPa velocity potential<br>
>> anomalies and then average it 10degS-10degN with period mean removed. Is<br>
>> there any function or have anyone done removal of the period mean, please<br>
>> let me know.<br>
>> Your help will be appreciated.<br>
>><br>
>> With regards,<br>
>> DP<br>
>> *Debashis Paul*<br>
>><br>
>> *Research Scholar,NIT-RKL*<br>
>> *Ph:+916002713628*<br>
>> *E-mail: <a href="mailto:Debgeo668@gmail.com" target="_blank">Debgeo668@gmail.com</a> <<a href="mailto:Debgeo668@gmail.com" target="_blank">Debgeo668@gmail.com</a>>*<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/20191211/c0ebfdec/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191211/c0ebfdec/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Wed, 11 Dec 2019 23:24:07 +1300<br>
From: Kunal Dayal <<a href="mailto:kday202@aucklanduni.ac.nz" target="_blank">kday202@aucklanduni.ac.nz</a>><br>
To: Ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Subject: Re: [ncl-talk] How to average monthly wind direction and<br>
represent them as wind barbs<br>
Message-ID:<br>
<CAKca7rO376cPNZ4ZRxita2qD7_fEfj4mz+NX4S=<a href="mailto:E9qc9fiGg5g@mail.gmail.com" target="_blank">E9qc9fiGg5g@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Hi,<br>
<br>
I working on mapping wind resources from WRFOUT files.<br>
<br>
I can easily plot contour maps of average wind speed for a month but I am<br>
unable to represent average wind direction in the form of wind barbs.<br>
<br>
The code I am using for the average wind speed contour maps is as follows:<br>
"<br>
begin<br>
a =<br>
addfile("/scale_wlg_persistent/filesets/project/uoa02450/wrf-topo2sstpblmix/wrfout_d03_2017-06-28_00:00:<br>
<a href="http://00.nc" rel="noreferrer" target="_blank">00.nc</a>","r")<br>
type = "pdf"<br>
wks = gsn_open_wks(type,"Average Wind Speed - Jul 2017")<br>
opts = True<br>
opts@MainTitle = "WIND SPEED"<br>
pltres = True<br>
mpres = True<br>
<br>
times = wrf_user_getvar(a,"times",-1)<br>
ntimes = dimsizes(times)<br>
<br>
u = wrf_user_getvar(a,"ua",-1)<br>
v = wrf_user_getvar(a,"va",-1)<br>
spd = (u*u + v*v)^(0.5)<br>
spd@description = "Wind Speed"<br>
spd@units = "m s-1"<br>
<br>
spd_avg = dim_avg_n(spd(0:30,:,:,:),0)<br>
spd_avg@description = "Average Wind Speed"<br>
spd_avg@units = "m s-1"<br>
<br>
res = True<br>
opts = res<br>
opts@cnFillOn = True<br>
opts@cnLinesOn = True<br>
opts@ContourParameters = (/ 0., 12., 0.5/)<br>
opts@gsnSpreadColorEnd = -3<br>
contour_spd_avg = wrf_contour(a,wks,spd_avg(0,:,:),opts)<br>
delete(opts)<br>
<br>
plot = wrf_map_overlays(a,wks,(/contour_spd_avg/),pltres,mpres)<br>
<br>
end<br>
"<br>
I wish to add monthly average wind direction to the same maps in the form<br>
of wind barbs.<br>
<br>
Appreciate your assistance and advice.<br>
<br>
Regards<br>
Kunal<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191211/f6e1a4fc/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191211/f6e1a4fc/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Wed, 11 Dec 2019 11:42:49 +0000 (UTC)<br>
From: Ehsan Taghizadeh <<a href="mailto:ehsantaghizadeh@yahoo.com" target="_blank">ehsantaghizadeh@yahoo.com</a>><br>
To: Ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Subject: [ncl-talk] GFS interpolation (vertical and horizontal)<br>
Message-ID: <<a href="mailto:105133163.7794812.1576064569861@mail.yahoo.com" target="_blank">105133163.7794812.1576064569861@mail.yahoo.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Dears,I hope you're doing great.May I ask if anybody has a script to interpolate GFS (0.5 deg) on the desired station points and also the desired vertical level? I mean (both horizontal and vertical interpolation).Should I attach my script, however it faces many problems?I'll be thankful for any help.<br>
SincerelyEhsan<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191211/06da539d/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191211/06da539d/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Wed, 11 Dec 2019 20:44:35 +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] DECADAL SVD<br>
Message-ID:<br>
<<a href="mailto:CALUZNzN-vpgK_RvCBGbOn4WVaV%2BFCV%2Bdb-j_PqHUCOSe%2Bn06Jw@mail.gmail.com" target="_blank">CALUZNzN-vpgK_RvCBGbOn4WVaV+FCV+db-j_PqHUCOSe+n06Jw@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Dear concerned expert(s),<br>
<br>
I am running decadal SVD between rainfall and SST. The script (attached)<br>
works well without filtering(smoothing) but once Lanczos low pass filter is<br>
applied I get the following error. The problem occur in the place of<br>
dealing with missing values because the variable NMS1 has [1]dimensions<br>
and sizes . Without filtering NMS1 had [657]. As far I understand Lanczos<br>
filter by default replace missing values with fill value. Could be that be<br>
the reason and if that is the case how could I run the SVD without the<br>
aspect of missing values?<br>
<br>
Any help/guidance will highly be appreciated.<br>
<br>
<br>
<br>
[image: image.png]<br>
<br>
Regards,<br>
<br>
Dickson Mbigi,<br>
<br>
Institute of Atmospheric Physics,<br>
<br>
University of Chinese Academy of Sciences,<br>
<br>
Beijing, China.<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191211/37af9dc3/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191211/37af9dc3/attachment-0001.html</a>><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: image.png<br>
Type: image/png<br>
Size: 34054 bytes<br>
Desc: not available<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191211/37af9dc3/attachment-0001.png" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191211/37af9dc3/attachment-0001.png</a>><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: decadal.svd.ncl<br>
Type: application/octet-stream<br>
Size: 24022 bytes<br>
Desc: not available<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191211/37af9dc3/attachment-0001.obj" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191211/37af9dc3/attachment-0001.obj</a>><br>
<br>
------------------------------<br>
<br>
Message: 5<br>
Date: Wed, 11 Dec 2019 20:26:55 +0530<br>
From: Debashis Paul <<a href="mailto:debgeo668@gmail.com" target="_blank">debgeo668@gmail.com</a>><br>
To: Ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Subject: [ncl-talk] help in remove period mean<br>
Message-ID:<br>
<<a href="mailto:CAFV8iJKJioZstgK2_Q1mFY-C93a_Xtnq2T_NNQZ1OujtnKMtZg@mail.gmail.com" target="_blank">CAFV8iJKJioZstgK2_Q1mFY-C93a_Xtnq2T_NNQZ1OujtnKMtZg@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Hi NCL,<br>
This mail is in reference to the mail I have done yesterday<br>
regarding help in doing removal of the periodic mean. I did how I was<br>
advised to do by Dennis Sir, but am getting an error of coordinate<br>
mismatch. I have attached my ncl script, if anyone can say me where am I<br>
wrong or how to do remove periodic mean it will be much appreciated.<br>
Thanks in advance,<br>
With regards,<br>
*Debashis Paul*<br>
<br>
*Research Scholar,NIT-RKL*<br>
*Ph:+916002713628*<br>
*E-mail: <a href="mailto:Debgeo668@gmail.com" target="_blank">Debgeo668@gmail.com</a> <<a href="mailto:Debgeo668@gmail.com" target="_blank">Debgeo668@gmail.com</a>>*<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191211/46280d6d/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191211/46280d6d/attachment-0001.html</a>><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: anomaly.ncl<br>
Type: application/octet-stream<br>
Size: 5168 bytes<br>
Desc: not available<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191211/46280d6d/attachment-0001.obj" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191211/46280d6d/attachment-0001.obj</a>><br>
<br>
------------------------------<br>
<br>
Message: 6<br>
Date: Wed, 11 Dec 2019 10:09:16 -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] DECADAL SVD<br>
Message-ID:<br>
<CAOF1d_6Nv=ivZ=<a href="mailto:JwX-_kWKnq9yucV0n_27g2TSF7Lra35LVr7Q@mail.gmail.com" target="_blank">JwX-_kWKnq9yucV0n_27g2TSF7Lra35LVr7Q@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
The script would require a substantive time investment to understand it.<br>
<br>
---<br>
<br>
*filwgts_lanczos*<br>
<<a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/filwgts_lanczos.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Built-in/filwgts_lanczos.shtml</a>><br>
<br>
<br>
*---*<br>
<br>
*re: "As far I understand Lanczos filter by default replace missing values<br>
with fill value."*.<br>
<br>
NCL has a nice feature where upon input if a variable has a missing_value<br>
attribute and no _FillValue attribute, NCL will 'dynamically' add a<br>
_FillValue attribute. I believe that CRU datasets have both _FillValue and<br>
missing_value attributes so this language feature is not applied.<br>
<br>
The Lanczos filter only recognizes the _FillValue attribute.<br>
<br>
ncl 0> nwt = 9<br>
ncl 1> fca = 0.2<br>
ncl 2> ihp = 0<br>
ncl 3> nsigma = 1.<br>
ncl 4> wgt = *filwgts_lanczos* (nwt, ihp, fca, -999., nsigma)<br>
ncl 5> N = 120<br>
ncl 6> x = *random_normal*(0,5,N)<br>
ncl 7> xFilter = *wgt_runave* ( x, wgt, 0 )<br>
ncl 8> print(x+" "+xFilter)<br>
<br>
(0) 4.35164 *9.96921e+36*<br>
(1) 4.36462 9.96921e+36<br>
(2) 9.73397 9.96921e+36<br>
(3) 4.91265 9.96921e+36<br>
(4) 1.76982 1.85172<br>
(5) -0.634983 -0.730493<br>
(6) -4.29466 -0.420796<br>
<br>
[SNIP]<br>
<br>
(114) -0.780523 0.369537<br>
(115) -5.00307 -1.90461<br>
(116) 3.39464 *9.96921e+36*<br>
(117) -5.93985 9.96921e+36<br>
(118) -1.86471 9.96921e+36<br>
(119) -2.42371 9.96921e+36<br>
<br>
So the beginning and ending (nwgt/2=4) of the filtered series are padded<br>
with _FillValue.<br>
<br>
Since NCL's SVD does not allow _FillValue, you must input just the temporal<br>
segment that has no _FillValue.<br>
<br>
===========<br>
<br>
*PLEASE NOTE: * Cherry (1996) discusses Singular Value Decomposition (SVD)<br>
and Canonical Correlation Analysis (CCA). Cherry's summary comment is:<br>
"Both methods have a high potential to produce spurious spatial patterns.<br>
Caution is always called for in interpreting results from either method."<br>
Newman and Sardeshmukh (1995) came to a similar conclusion in their paper<br>
which focused on SVD: "These results suggest that any physical<br>
interpretation of SVD pairs may be unjustified."<br>
<br>
*REFERENCES:*<br>
<br>
Cherry, S. (1996): *Singular Value Decomposition Analysis and<br>
Canonical Correlation Analysis*<br>
<<a href="https://doi.org/10.1175/1520-0442(1996)009" rel="noreferrer" target="_blank">https://doi.org/10.1175/1520-0442(1996)009</a><2003:SVDAAC><a href="http://2.0.CO" rel="noreferrer" target="_blank">2.0.CO</a>;2>.<br>
J. Climate: pp 2003-2009.<br>
<br>
Newman, M. and Sardeshmukh, P.D. (1995): *A Caveat Concerning Singular<br>
Value Decomposition*<br>
<<a href="https://doi.org/10.1175/1520-0442(1995)008" rel="noreferrer" target="_blank">https://doi.org/10.1175/1520-0442(1995)008</a><0352:ACCSVD><a href="http://2.0.CO" rel="noreferrer" target="_blank">2.0.CO</a>;2><br>
J. Climate: pp 352-360.<br>
<br>
<br>
On Wed, Dec 11, 2019 at 5:44 AM dickson mbigi via ncl-talk <<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
<br>
> Dear concerned expert(s),<br>
><br>
> I am running decadal SVD between rainfall and SST. The script (attached)<br>
> works well without filtering(smoothing) but once Lanczos low pass filter is<br>
> applied I get the following error. The problem occur in the place of<br>
> dealing with missing values because the variable NMS1 has [1]dimensions<br>
> and sizes . Without filtering NMS1 had [657]. As far I understand Lanczos<br>
> filter by default replace missing values with fill value. Could be that be<br>
> the reason and if that is the case how could I run the SVD without the<br>
> aspect of missing values?<br>
><br>
> Any help/guidance will highly be appreciated.<br>
><br>
><br>
><br>
> [image: image.png]<br>
><br>
> Regards,<br>
><br>
> Dickson Mbigi,<br>
><br>
> Institute of Atmospheric Physics,<br>
><br>
> University of Chinese Academy of Sciences,<br>
><br>
> Beijing, China.<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/20191211/a80a1645/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191211/a80a1645/attachment-0001.html</a>><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: image.png<br>
Type: image/png<br>
Size: 34054 bytes<br>
Desc: not available<br>
URL: <<a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191211/a80a1645/attachment-0001.png" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191211/a80a1645/attachment-0001.png</a>><br>
<br>
------------------------------<br>
<br>
Message: 7<br>
Date: Wed, 11 Dec 2019 10:25:44 -0700<br>
From: Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>><br>
To: Kunal Dayal <<a href="mailto:kday202@aucklanduni.ac.nz" target="_blank">kday202@aucklanduni.ac.nz</a>><br>
Cc: Ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Subject: Re: [ncl-talk] How to average monthly wind direction and<br>
represent them as wind barbs<br>
Message-ID:<br>
<CAOF1d_6o2Zv9R1dr=<a href="mailto:YCorvhVKRXGNNAfQQ9e0qdWWYnTYXTXxA@mail.gmail.com" target="_blank">YCorvhVKRXGNNAfQQ9e0qdWWYnTYXTXxA@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
To compute average wind direction<br>
[1] compute average zonal [u] and meridional [v] components<br>
[2] use the u and v component averages as input to *wind_direction*<br>
<<a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/wind_direction.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Contributed/wind_direction.shtml</a>><br>
to commute the average direction(s)<br>
---<br>
NCL also has a function *wind_stats*<br>
<<a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/wind_stats.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Contributed/wind_stats.shtml</a>><br>
<br>
<br>
On Wed, Dec 11, 2019 at 3:24 AM Kunal Dayal via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
wrote:<br>
<br>
> Hi,<br>
><br>
> I working on mapping wind resources from WRFOUT files.<br>
><br>
> I can easily plot contour maps of average wind speed for a month but I am<br>
> unable to represent average wind direction in the form of wind barbs.<br>
><br>
> The code I am using for the average wind speed contour maps is as follows:<br>
> "<br>
> begin<br>
> a =<br>
> addfile("/scale_wlg_persistent/filesets/project/uoa02450/wrf-topo2sstpblmix/wrfout_d03_2017-06-28_00:00:<br>
> <a href="http://00.nc" rel="noreferrer" target="_blank">00.nc</a>","r")<br>
> type = "pdf"<br>
> wks = gsn_open_wks(type,"Average Wind Speed - Jul 2017")<br>
> opts = True<br>
> opts@MainTitle = "WIND SPEED"<br>
> pltres = True<br>
> mpres = True<br>
><br>
> times = wrf_user_getvar(a,"times",-1)<br>
> ntimes = dimsizes(times)<br>
><br>
> u = wrf_user_getvar(a,"ua",-1)<br>
> v = wrf_user_getvar(a,"va",-1)<br>
> spd = (u*u + v*v)^(0.5)<br>
> spd@description = "Wind Speed"<br>
> spd@units = "m s-1"<br>
><br>
> spd_avg = dim_avg_n(spd(0:30,:,:,:),0)<br>
> spd_avg@description = "Average Wind Speed"<br>
> spd_avg@units = "m s-1"<br>
><br>
> res = True<br>
> opts = res<br>
> opts@cnFillOn = True<br>
> opts@cnLinesOn = True<br>
> opts@ContourParameters = (/ 0., 12., 0.5/)<br>
> opts@gsnSpreadColorEnd = -3<br>
> contour_spd_avg = wrf_contour(a,wks,spd_avg(0,:,:),opts)<br>
> delete(opts)<br>
><br>
> plot = wrf_map_overlays(a,wks,(/contour_spd_avg/),pltres,mpres)<br>
><br>
> end<br>
> "<br>
> I wish to add monthly average wind direction to the same maps in the form<br>
> of wind barbs.<br>
><br>
> Appreciate your assistance and advice.<br>
><br>
> Regards<br>
> Kunal<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/20191211/dae31089/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191211/dae31089/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 8<br>
Date: Wed, 11 Dec 2019 10:43:03 -0700<br>
From: Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>><br>
To: Ehsan Taghizadeh <<a href="mailto:ehsantaghizadeh@yahoo.com" target="_blank">ehsantaghizadeh@yahoo.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] GFS interpolation (vertical and horizontal)<br>
Message-ID:<br>
<<a href="mailto:CAOF1d_7xy3XfZQQrbELas8nGtOkYfyURWjaUM064%2BqJHdeZW1w@mail.gmail.com" target="_blank">CAOF1d_7xy3XfZQQrbELas8nGtOkYfyURWjaUM064+qJHdeZW1w@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
ncl-talk does not want to see a script "with many problems."<br>
We do not have the time to sort/debug multiple issues.<br>
---<br>
GFS are ordered N->S<br>
[1] reorder S->N via NCL syntax: x(:,:,::-1,:)<br>
[2] Use *linint2_points_Wrap*<br>
<<a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/linint2_points_Wrap.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Contributed/linint2_points_Wrap.shtml</a>><br>
to interpolate to desired locations<br>
[3] Use *int2p_n_Wrap*<br>
<<a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/int2p_n_Wrap.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Contributed/int2p_n_Wrap.shtml</a>><br>
to interpolate to desired level(s)<br>
<br>
<br>
<br>
On Wed, Dec 11, 2019 at 4:42 AM Ehsan Taghizadeh via ncl-talk <<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
<br>
> Dears,<br>
> I hope you're doing great.<br>
> May I ask if anybody has a script to interpolate GFS (0.5 deg) on the<br>
> desired station points and also the desired vertical level? I mean (both<br>
> horizontal and vertical interpolation).<br>
> Should I attach my script, however it faces many problems?<br>
> I'll be thankful for any help.<br>
><br>
> Sincerely<br>
> Ehsan<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/20191211/e6dece39/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191211/e6dece39/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 9<br>
Date: Wed, 11 Dec 2019 11:07:44 -0700<br>
From: Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>><br>
To: Debashis Paul <<a href="mailto:debgeo668@gmail.com" target="_blank">debgeo668@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] help in remove period mean<br>
Message-ID:<br>
<CAOF1d_7mVa9WkyC00mJ2pMOZcPRWMT60PW=qeK8_-8a=<a href="mailto:2u5SmQ@mail.gmail.com" target="_blank">2u5SmQ@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
re: "...getting an error of coordinate mismatch"<br>
<br>
This does not provide any useful information.<br>
*You* should be using 'printVarSummary' on variables prior to the offending<br>
line.<br>
Look at the output. What coordinates could be a mismatch?<br>
<br>
<br>
<br>
On Wed, Dec 11, 2019 at 7:56 AM Debashis Paul via ncl-talk <<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br>
<br>
> Hi NCL,<br>
> This mail is in reference to the mail I have done yesterday<br>
> regarding help in doing removal of the periodic mean. I did how I was<br>
> advised to do by Dennis Sir, but am getting an error of coordinate<br>
> mismatch. I have attached my ncl script, if anyone can say me where am I<br>
> wrong or how to do remove periodic mean it will be much appreciated.<br>
> Thanks in advance,<br>
> With regards,<br>
> *Debashis Paul*<br>
><br>
> *Research Scholar,NIT-RKL*<br>
> *Ph:+916002713628*<br>
> *E-mail: <a href="mailto:Debgeo668@gmail.com" target="_blank">Debgeo668@gmail.com</a> <<a href="mailto:Debgeo668@gmail.com" target="_blank">Debgeo668@gmail.com</a>>*<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/20191211/fb0c3e2d/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20191211/fb0c3e2d/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 193, Issue 21<br>
*****************************************<br>
</blockquote></div>