Re: Problem calculating % of values in timeseries below a cutoff point - how to define a variable that you later calculate?

Hi - yes, I spotted the mis-named variable -so instead of copying the var
coords of precip_knt to precip_pct, I copied from the original precip file
(precip_mod1)... Thanks a lot!

-M.

On Wed, May 30, 2012 at 7:34 AM, Madeleine Patterson <
madeleine.patterson77@gmail.com> wrote:

> Hi -thanks for this; the script executes until the copy_varCoords part, I
> get an error message that the number of subscripts do not match number of
> dimensions of variable; 3 subscripts used, 2 subscripts expected...
>
> ? Can I plot a global map of these percentages without var coords?
>
> Thanks,
> M
>
>
> On Wed, May 30, 2012 at 6:50 AM, Dennis Shea <shea@ucar.edu> wrote:
>
>> [1] When reading an entire variable, it is more efficient
>> to not reference each dimension
>>
>> precip_mod1 = in1->rainfall ; more efficient
>>
>> precip_mod1 = in1->rainfall(:,:,:) ; less efficient
>>
>> If you want to indicate a variable's rank, then use a comment
>>
>> precip_mod1 = in1->rainfall ; (:,:,:) or 3D or (lat,lon,time)
>>
>> ===
>> [2]
>>
>> pcrt = 2.0
>> precip_knt = dim_num_n(precip_mod1.lt.pcrt, 2) ; (:,:)
>> precip_tot = dim_num_n(.not.ismissing(**precip_mod1), 2) ; (:,:)
>> precip_pcrt= (100.*precip_knt))/precip_tot ; %(:,:)
>>
>> copy_varCoords(precip_knt(;,:,**0),precip_pcrt) ; contributed.ncl
>> precip_pcrt@long_name = "..."
>> precip_pcrt@units = "%"
>> printVarSummary*precip_pcrt)
>>
>>
>> On 5/29/12 1:25 PM, Madeleine Patterson wrote:
>>
>>> Hi,
>>>
>>> I am trying to calculate % of values in a timeseries below a cutoff
>>> point, but keep getting error messages related to not having defined the
>>> variables that I am calculating (not reading in) properly. Just
>>> calculating the variable in the do loop doesn't seem to define them, and
>>> when I try to define them before I calculate them, I've tried a couple
>>> ways but haven't got it right so far. Can anyone help me out with this
>>> problem? My script is below...
>>>
>>> The error message is: undefined identifer: pct_drought_mod1 is
>>> undefined, can't continue.
>>>
>>> *********************************************************
>>> begin
>>>
>>> load "$NCARG_ROOT/lib/ncarg/**nclscripts/csm/gsn_code.ncl"
>>>
>>> load "$NCARG_ROOT/lib/ncarg/**nclscripts/csm/gsn_csm.ncl"
>>>
>>> load "$NCARG_ROOT/lib/ncarg/**nclscripts/csm/contributed.**ncl"
>>>
>>> load "$NCARG_ROOT/lib/ncarg/**nclscripts/csm/shea_util.ncl"
>>>
>>> ; use "$NCARG_ROOT/lib/ncarg/data/**cdf/landsea.nc <http://landsea.nc>"
>>>
>>>
>>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;**;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;**
>>> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;**;
>>>
>>> setfileoption("nc","**SuppressClose",False)
>>>
>>>
>>> ; open 1st pft map;
>>>
>>> in1 = addfile("rainfall.nc <http://rainfall.nc>","r")
>>>
>>>
>>> lat =in1->lat
>>> lon=in1->lon
>>>
>>> precip_mod1 = in1->rainfall(:,:,:) ;;; this is lat, lon time (91,144,
>>> 360)
>>>
>>> ; calculate the percentage of months in the timeseries where rainfall <
>>> 2 mm/d
>>>
>>>
>>> ;pct_drought_mod1(i,j)=0 ---> this sort of thing doesn't
>>> seem to work
>>> ;num_drought_mnth_mod1(i,j)=0
>>>
>>> ;num_drought_mnth_mod1 =0
>>> ;pct_drought_mod1 = 0
>>>
>>> count = 0
>>> do i=1,144
>>> do j=1,91
>>> do k=1,348
>>> if (k .lt. 2.0) then
>>> count = count +1
>>> end if
>>> end do
>>> print(count)
>>>
>>> num_drought_mnth_mod1(i,j) = count
>>> pct_drought_mod1(i,j) = count/348
>>> end do
>>> end do
>>>
>>>
>>> end
>>>
>>>
>>> ______________________________**_________________
>>> ncl-talk mailing list
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/**mailman/listinfo/ncl-talk<http://mailman.ucar.edu/mailman/listinfo/ncl-talk>
>>>
>>
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue May 29 15:58:37 2012