[ncl-talk] WRAPIT, error in the fortran code
Mary Haley
haley at ucar.edu
Thu Jan 12 11:29:32 MST 2017
Ying,
Can you send me the location of the file on yellowstone so I can try
running it? You can send this offline if you prefer.
The error message you provided makes no sense for the script you sent,
because it's complaining about the X axis coordinates:
warning:ContourPlotInitialize: X coordinates out of data range: defaulting
warning:ContourPlotInitialize: Zero X coordinate span:
defaulting:[errno=1104]
but from your script, it doesn't look like you are setting any coordinates.
--Mary
On Thu, Jan 5, 2017 at 9:27 PM, Ying Song <ysong4 at slu.edu> wrote:
> Thanks for the reply. Sorry this is no longer a WRAPIT question, I will
> pay attention and write a new subject next time.
>
> Follow your suggestion, I did not set any trX/YMin/MaxF resources, and
> cut the "create plot" section into very short. However, the error is still
> the same.
> ;*********************
> ;create plot
> ;*********************
> res at cnFillOn = True ; turn on color
> res at gsnSpreadColors = True ; spread out color table
> res at gsnSpreadColorStart = 2
> res at gsnSpreadColorEnd = -3
>
> res at tiXAxisString = "Time"
> res at tiYAxisString = "vertical velocity (m/s)"
>
> plot = gsn_csm_contour(wks,w_bin,res)
> frame(wks)
> end
>
> ---------------
> Also, the plot includes 2 pages. But they are both blank. I only have one
> plot in NCL script. Don't understand why this happened. Please help. Thank
> you!
>
> On Thu, Jan 5, 2017 at 6:50 PM, Mary Haley <haley at ucar.edu> wrote:
>
>> Ying,
>>
>> Your w_bin array is 100 x 61. Since you didn't provide any X/Y
>> coordinates for this data, NCL will draw the plot with the Y axis ranging
>> from 0 to 99 and the X axis ranging from 0 to 60.
>>
>> You set:
>>
>> res at trXMaxF = 61
>>
>> which is out of the range of the actual X axis, so that's why you got the
>> "X coordinates out of data range" error. This may have also been the
>> cause of the second warning as well.
>>
>> You also set:
>>
>> res at trYMinF = 0.0 ; Limits for Y axis. percentage from 0
>> to 100.
>> res at trYMaxF = 1.0
>>
>> Your Y axis is going from 0 to 99, so if you set it to go from 0 to 1,
>> you are only seeing a very small subset of your data, and actually, this
>> may have caused the second warning.
>>
>> I suggest first plotting your data without setting *any* trX/YMin/MaxF
>> resources. Then, when you see what the X and Y axis looks like, you can
>> decide what values you want to set these resources to.
>>
>> --Mary
>>
>>
>> On Thu, Jan 5, 2017 at 3:31 PM, Ying Song <ysong4 at slu.edu> wrote:
>>
>>> Hello, Dennis,
>>>
>>> Thanks for helping me fix the WRAPIT error. It works fine now! But I
>>> have a new question here. My NCL script can not plot my output (w_bin). I
>>> use "gsn_csm_contour" to plot the velocity frequency changes with time and
>>> submit the job to yellowstone. The error message is:
>>> -------------------------------
>>> warning:ContourPlotInitialize: X coordinates out of data range:
>>> defaulting
>>> warning:ContourPlotInitialize: Zero X coordinate span:
>>> defaulting:[errno=1104]
>>> ----------------------------
>>> I think my output from WRAPIT is correct. And I use "gsn_csm_contour" to
>>> plot before. Never have this error before.
>>> I attached my scripts and .out file of the job.
>>> Please help me check the problem. Thank you very much!
>>>
>>> On Mon, Jan 2, 2017 at 10:09 AM, Dennis Shea <shea at ucar.edu> wrote:
>>>
>>>> Yes. WRAPIT does not need any 'local' fortran variables. They can be
>>>> placed after the 'C NCLEND' delimeter
>>>>
>>>> Good Luck
>>>>
>>>> On Fri, Dec 30, 2016 at 1:59 PM, Ying Song <ysong4 at slu.edu> wrote:
>>>>
>>>>> Thanks for the reply!
>>>>>
>>>>> You are right, Dennis. The error is caused by nbin-1, it is fixed now.
>>>>> I want to make sure I understand you correctly. You mean the local
>>>>> variables such as bins, bcnt, and npts need to be defined in .f script but
>>>>> not in stub file. Stub file only care about the variables pass into or out
>>>>> the fortran subroutine?
>>>>>
>>>>> Also, thanks for sending the binning function information also. Let me
>>>>> take a look!
>>>>>
>>>>> On Fri, Dec 30, 2016 at 2:02 PM, Dennis Shea <shea at ucar.edu> wrote:
>>>>>
>>>>>> Forgot: NCL has some binning functions
>>>>>>
>>>>>> http://www.ncl.ucar.edu/Applications/
>>>>>>
>>>>>> Under "Data Analysis" ... click "Binning"
>>>>>>
>>>>>> See the examples. They use the following functions:
>>>>>>
>>>>>> http://www.ncl.ucar.edu/Document/Functions/Built-in/bin_sum.shtml
>>>>>> http://www.ncl.ucar.edu/Document/Functions/Built-in/bin_avg.shtml
>>>>>>
>>>>>> On Fri, Dec 30, 2016 at 12:31 PM, Dennis Shea <shea at ucar.edu> wrote:
>>>>>>
>>>>>>> WRAPIT is not a complete fortran-77 parser. WRAPIT strictly looks
>>>>>>> at the arguments being passed between NCL (written in C) and the fortran
>>>>>>> subroutine. The types and sizes must match.
>>>>>>>
>>>>>>> WRAPIT does not understand fortran 'dynamic' array sizing. Also,
>>>>>>> WRAPIT does not allow things like redefining the array extent such as:
>>>>>>> foo(-N:N). A stub file or is one way of addressing these limitations. NOTE:
>>>>>>> the stub file could be a valid fortran subroutine containing the NCL
>>>>>>> delimeters [ C NCLFORTSTART and C NCLEND ] that handles assorted issues in
>>>>>>> the fortran world.
>>>>>>>
>>>>>>> That sadi:
>>>>>>>
>>>>>>> I speculate that the issue is the 'nbin-1' id the array
>>>>>>> declarations. Maybe, add and 'nbin1' argument where
>>>>>>> nbin1 = nbin-1
>>>>>>> and make the stub file be
>>>>>>>
>>>>>>> C NCLFORTSTART
>>>>>>> subroutine wacftd(nxx,nyx,nkx,ntx,nbin,nbin1, wa,bfra)
>>>>>>> integer nxx,nyx,nkx,ntx,nbin, nbin1
>>>>>>> real wa(nxx,nyx,nkx,ntx) ! input variable
>>>>>>>
>>>>>>> c c c real bins(ntx,nbin),bcnt(ntx,nbin-1), bfra(ntx,nbin-1)
>>>>>>> real bfra(ntx,nbin1) ! return
>>>>>>> variable
>>>>>>> C NCLEND
>>>>>>> C WRAPIT will not look at the following.
>>>>>>> integer npts
>>>>>>> real bins(ntx,nbin), bcnt(ntx,nbin-1) ! local ... fortran
>>>>>>> compiler handles this
>>>>>>>
>>>>>>> npts = nxx*nyx*nkx
>>>>>>>
>>>>>>> This would have to be added to the NCL script calling the subroutine
>>>>>>> as a shared object.
>>>>>>> Further, variable 'bins' is not being passed to-from the NCL script.
>>>>>>> I can an
>>>>>>>
>>>>>>> ---
>>>>>>> Re: PARAMETER statements
>>>>>>>
>>>>>>> WRAPIT strictly looks at the arguments being passed between NCL
>>>>>>> (written in C) and the fortran subroutine. The types and sizes must match.
>>>>>>>
>>>>>>> C NCLFORTSTART
>>>>>>> subroutine fooo(nxx,nyx,nkx,ntx,nbin, ........wa,.....)
>>>>>>> integer nxx,nyx,nkx,ntx,nbin, nbin1
>>>>>>> real wa(nxx,nyx,nkx,ntx) ! input variable
>>>>>>> C NCLEND
>>>>>>> integer ii, jj ! WRAPIT does not
>>>>>>> care about the following
>>>>>>> parameter (ii=30, jj=999) ! baecause they are not passed
>>>>>>> into or out the subroutine
>>>>>>> double precision work(ii,jj)
>>>>>>>
>>>>>>> ---
>>>>>>>
>>>>>>>
>>>>>>> You do not have to have a stub file. A minor rearrangement
>>>>>>>
>>>>>>>
>>>>>>> C NCLFORTSTART
>>>>>>> subroutine wacftd(nxx,nyx,nkx,ntx,nbin,nbin1, wa,bfra)
>>>>>>> integer nxx,nyx,nkx,ntx,nbin, nbin1
>>>>>>> real wa(nxx,nyx,nkx,ntx) ! input variable
>>>>>>> real bfra(ntx,nbin1) ! return
>>>>>>> variable
>>>>>>> C NCLEND
>>>>>>> real bins(ntx,nbin),bcnt(ntx,nbin-1)
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> Good luck
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> On Fri, Dec 30, 2016 at 11:33 AM, Ying Song <ysong4 at slu.edu> wrote:
>>>>>>>
>>>>>>>> Hello,
>>>>>>>>
>>>>>>>> I am using WRAPIT in NCL to calculate velocity frequency, following
>>>>>>>> the ex01. After running the command:
>>>>>>>> WRAPIT cftd.stub cftd.f
>>>>>>>> The error shows below:
>>>>>>>> ++++++++++++++++++++
>>>>>>>> WRAPIT Version: 120209
>>>>>>>> A syntax error occurred while parsing: 1
>>>>>>>> COMPILING cftd.f
>>>>>>>> LINKING
>>>>>>>> END WRAPIT
>>>>>>>> ++++++++++++++++++++
>>>>>>>> I think there are some problems in my cftd.f script but I can not
>>>>>>>> find it. So I write this email asking for help!
>>>>>>>>
>>>>>>>> cftd.f and cftd.stub are attached. Thank you very much!
>>>>>>>>
>>>>>>>> --
>>>>>>>> Ying Song, PhD, Post Doctoral Fellow
>>>>>>>> Department of Earth and Atmospheric Sciences
>>>>>>>> Saint Louis University
>>>>>>>>
>>>>>>>> ysong4 at slu.edu
>>>>>>>>
>>>>>>>>
>>>>>>>> _______________________________________________
>>>>>>>> ncl-talk mailing list
>>>>>>>> ncl-talk at ucar.edu
>>>>>>>> List instructions, subscriber options, unsubscribe:
>>>>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>>>>>
>>>>>>>>
>>>>>>>
>>>>>>
>>>>>
>>>>>
>>>>> --
>>>>> Ying Song, PhD, Post Doctoral Fellow
>>>>> Department of Earth and Atmospheric Sciences
>>>>> Saint Louis University
>>>>>
>>>>> ysong4 at slu.edu
>>>>>
>>>>>
>>>>
>>>
>>>
>>> --
>>> Ying Song, PhD, Post Doctoral Fellow
>>> Department of Earth and Atmospheric Sciences
>>> Saint Louis University
>>>
>>> ysong4 at slu.edu
>>>
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>>
>
>
> --
> Ying Song, PhD, Post Doctoral Fellow
> Department of Earth and Atmospheric Sciences
> Saint Louis University
>
> ysong4 at slu.edu
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170112/976fb15f/attachment-0001.html
More information about the ncl-talk
mailing list