[ncl-talk] plotting monthly climatology form hourly data

Dennis Shea shea at ucar.edu
Fri Feb 21 12:29:05 MST 2020 

You can *help yourself *by using print statements to isolate where the
messages are generated.
You must do a bettr job of examining the code.

At the top of your plot section you have


*plot  = new (nmos, graphic)  *

You never fill this array. It contains all missing

For example:

*print("BEFORE gsn_add_shapefile_polylines")*
  shape1 = gsn_add_shapefile_polylines(wks,*plot*,shapefile_name1,lnres)
  shape2 = gsn_add_shapefile_polylines(wks,*plot*,shapefile_name2,lnres)
*print("AFTER gsn_add_shapefile_polylines")*

Again, the *plot* argument contains missing. Likely, you want:

  shape1 = gsn_add_shapefile_polylines(wks,*plot*1,shapefile_name1,lnres)
  shape2 = gsn_add_shapefile_polylines(wks,*plot*,1shapefile_name2,lnres)

You can not just copy a section from some script and expect it to work.
You must examine the arguments to make sure they agree with you are using.

Good luck

On Fri, Feb 21, 2020 at 3:08 AM Setareh Rahimi <setareh.rahimi at gmail.com>
wrote:

> Dear Dennis,
>
> Thanks  so much for your help. However, I do not know why should I get the
> following warning (as appeared in previous errors):
>
> warning:_NclIRemovePrimitive: First parameter is a missing value,
> returning missing values
>
>
> How can I remove this, please?
>
> Best wishes,
>
>
>
> On Fri, Feb 21, 2020 at 3:51 AM Dennis Shea <shea at ucar.edu> wrote:
>
>> The error message is tellng you the issue. It says that variable *plot1
>> is undefined.*
>> *plot1 *and *plot2 *are arrays and must be explicitly declared as such
>> before use.
>>
>> Add the following::
>>
>> * plot1 = new( nmos, "graphic")*
>> * plot2= new( nmos, "graphic")*
>>
>>   do nmo=0,nmos-1                               ; loop over the months
>>      res at gsnCenterString   = months(nmo)
>>      plot1(nmo) = gsn_csm_contour_map(wks,xClm(nmo,:,:),res)  ; create
>> plot
>>      plot2(nmo) = gsn_csm_vector(wks,uClm(nmo,:,:),vClm(nmo,:,:),resv)
>>      overlay(plot1(nmo),plot2(nmo))
>>   end do
>>
>> ==================
>>
>> Also, you should be aware of the
>> *gsn_csm_vector_scalar_map function*
>> <http://ww.ncl.ucar.edu/Document/Graphics/Interfaces/gsn_csm_vector_scalar_map.shtml>
>>
>>
>> On Thu, Feb 20, 2020 at 5:17 AM Setareh Rahimi <setareh.rahimi at gmail.com>
>> wrote:
>>
>>> Dear all,
>>> I would like to plot a variable (dust) from model output [like the
>>> previous discussion] and also overlay wind vectors on the plot. I did some
>>> changes in the script but got some errors. Please have a look at the
>>> attached script and image showing the errors. I do not know why (plot1)
>>> is undefined! May I ask you please to guide me on how to address this issue?
>>> Many thanks in advance,
>>> Best wishes,
>>>
>>> On Fri, Nov 8, 2019 at 4:35 PM Setareh Rahimi <setareh.rahimi at gmail.com>
>>> wrote:
>>>
>>>> Dear Dennis,
>>>> Thank you so much for your explanations and guidance. However, I
>>>> noticed that the results for every year are similar to each other. Please
>>>> have look at the attached files. This is not reasonable.  I see no problems
>>>> in the script, so why this happened?
>>>> Best wishes,
>>>>
>>>> On Fri, Nov 8, 2019 at 4:36 AM Dennis Shea <shea at ucar.edu> wrote:
>>>>
>>>>> The issue is the variable "time". This a subtle issue.
>>>>>
>>>>>   f        = *addfile*("khuzetan_DUST04.2000010100.nc", "r")
>>>>>   time  = f->time
>>>>>   timeb = f->time_bnds
>>>>>
>>>>>   ymdh_tag = *cd_calendar*
>>>>> <http://www.ncl.ucar.edu/Document/Functions/Built-in/cd_calendar.shtml>(time
>>>>>      ,-3)
>>>>>   ymdh_Strt= *cd_calendar*
>>>>> <http://www.ncl.ucar.edu/Document/Functions/Built-in/cd_calendar.shtml>(timeb(:,0),-3)
>>>>> ; beginning time 'bound'
>>>>>   ymdh_Last= *cd_calendar*
>>>>> <http://www.ncl.ucar.edu/Document/Functions/Built-in/cd_calendar.shtml>(timeb(:,1),-3)
>>>>> ; ending time 'bound'
>>>>>
>>>>>   print(timeb(:,0)+"  "+timeb(:,1)+"  "+time+"  " \
>>>>>        +ymdh_Strt +"  "+ymdh_Last +"  "+ymdh_tag)
>>>>> =============
>>>>> One month of time values. Here 124 (31*4) January hours.
>>>>>
>>>>> (0)     timeb_0  timeb_1   time       ymdh_Strt    ymdh_Last
>>>>> ymdh_tag
>>>>>           439032  439038  439038  2000*010100*  2000*010106*  2000
>>>>> *010106*
>>>>> *                :            :             :
>>>>> :                     :                    :*
>>>>> (123) 439770  439776  439776  2000*013118*  2000*020100*  2000
>>>>> *020100*
>>>>> ===============
>>>>>
>>>>> *The last tag value [timeb(:,1)] for the month of JANUARY is FEBRUARY
>>>>> !!!  *
>>>>>
>>>>> *My personal view is that this is not the appropriate time tag. *
>>>>>
>>>>> *Rather, the timeb(:,0)]  should be used. *
>>>>> ================
>>>>>
>>>>> In a nutshell, the work-around is to reassign the time coordinate
>>>>> associated with each variable.
>>>>>
>>>>>   timebnd = a[:]->time_bnds
>>>>>   time0   = timebnd(:,0)                        ; extract the start
>>>>> bound
>>>>>   [SNIP]
>>>>>   x  = a[:]->emflx                     ; read emission flux from all
>>>>> files
>>>>>   x&time = time0                     ; reassign original 'tag' time
>>>>> [SNIP]
>>>>>   u = b[:]->uas(:,0,:,:)              ; read U from all files
>>>>>   u&time = time0                     ; reassign original 'tag' time
>>>>>
>>>>>   v = b[:]->vas(:,0,:,:)              ; read V from all files
>>>>>   v&time = time0                     ; reassign original 'tag' time
>>>>>
>>>>> =========
>>>>> Comments:
>>>>> [0] You MUST examine the time variable more carefully.
>>>>>       Using yyyymmdd start/stop dates is not appropriate for
>>>>> yyyymmddhh.
>>>>> [1] You do *not *have to compute the daily mean values
>>>>> [2] *calculate_monthly_values*
>>>>> <http://www.ncl.ucar.edu/Document/Functions/Contributed/calculate_monthly_values.shtml>
>>>>> will accept 6-hrly data
>>>>> [3] *http://www.ncl.ucar.edu/Applications/rcm.shtml*
>>>>> <http://www.ncl.ucar.edu/Applications/rcm.shtml>
>>>>> Examples 5,6,7: illustrate how to plot variables for the Regional
>>>>> Climate Model (RCM).
>>>>> =========
>>>>> See Attached
>>>>>
>>>>>
>>>>> On Thu, Nov 7, 2019 at 6:10 AM Setareh Rahimi via ncl-talk <
>>>>> ncl-talk at ucar.edu> wrote:
>>>>>
>>>>>> Dear all,
>>>>>>
>>>>>> I need to plot the monthly average from a model output, which gives
>>>>>> me 4times data per day( at 00, 06,12,18). So I first need to convert hourly
>>>>>> data tp daily data and then convert to monthly. And finally calculates the
>>>>>> climatology average of those monthly data from 2000-2014.
>>>>>> However, once I ran my script I faced some error:
>>>>>> 1. I did not have data for 2015, but when using print (yyyymmdd)   year
>>>>>> 2015 printed. ( please have a look at the attachment number 1)
>>>>>> 2. After reading variable x this warning comes up: Aggregated
>>>>>> dimension coordinate values are non-monotonic; check aggregated file
>>>>>> ordering (warning attached)
>>>>>> 3. Once using clmMonTLL function ( available on :
>>>>>> https://www.ncl.ucar.edu/Applications/Scripts/climo_3.ncl), NCL
>>>>>> tells that :  clmMonTLL: dimension must be a multiple of 12 . I just
>>>>>> wonder why this happened. I used 15 years with 12 months for each year!! (
>>>>>> attached)
>>>>>> I attached the script and data for the year 2000. Could you please
>>>>>> kindly guide me what are my mistakes?
>>>>>> Best wishes,
>>>>>>  str.zip
>>>>>> <https://drive.google.com/file/d/12-n6pQ4JF6nwIRBVW4iaTKaznfP1FGet/view?usp=drive_web>
>>>>>>
>>>>>> --
>>>>>> S.Rahimi
>>>>>>
>>>>>> _______________________________________________
>>>>>> ncl-talk mailing list
>>>>>> ncl-talk at ucar.edu
>>>>>> List instructions, subscriber options, unsubscribe:
>>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>>
>>>>>
>>>>
>>>> --
>>>> S.Rahimi
>>>>
>>>>
>>>
>>> --
>>> S.Rahimi
>>>
>>>
>
> --
> S.Rahimi
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200221/bcb45efa/attachment.html>

More information about the ncl-talk mailing list