[ncl-talk] dividing by climatology
Jayant
jayantkp2979 at gmail.com
Mon Feb 10 13:03:43 MST 2020
Thank you Denis. I will discuss with my colleagues.
One more question:
in the example with nlat and nlon, "do loop" for nyrs or for each month on
top of lat-lon loop is not required?
best,
On Mon, Feb 10, 2020 at 12:53 PM Dennis Shea <shea at ucar.edu> wrote:
> re: "... not sure whether this means normalization?"
>
> You can "normalize by any non-zero quantity:
> eg: Max value, Mean, standard deviation,....
> You should consult with a colleague/advisor for which is appropriate for
> you.
> ---
> The following is a simple method of normalizing
>
> nmos=12
> nyrs=20
> ntim= nyrs*nmos
> monClim = random_uniform(10,100,nmos) ; [*]
> monData = random_uniform(10,100,ntim) ; [*]
> do nmo=0,nmos-1
> monData(nmo::nmos) = monData(nmo::nmos)/ monClim(nmo)
> end do
>
> nlat = 72
> mlon = 144
>
> monClim* :=* random_uniform(10,100,(/nmos,nlat,mlon/)) ;[*][*][*]
> monData *:=* random_uniform(10,100,(/ntim,nlat,mlon/))
> do nl=0,nlat-1
> do ml=0,mlon-1
> monData(nmo::nmos,nl,ml) = monData(nmo::nmos,nl,ml)/
> monClim(nmo,nl,ml)
> do
> end do
>
> On Mon, Feb 10, 2020 at 6:16 AM Jayant via ncl-talk <ncl-talk at ucar.edu>
> wrote:
>
>> Hi,
>> I have monthly fire data for 20 years. I have computed climatology. I
>> want to calculate ratio of each month and its respective climatology to
>> obtain weights. Is there a function that does this calculation? I am not
>> sure whether this means normalization?
>> Cheers,
>> _______________________________________________
>> 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/20200210/dc8f8540/attachment.html>
More information about the ncl-talk
mailing list