[ncl-talk] WRAPIT, error in the fortran code
Ying Song
ysong4 at slu.edu
Thu Jan 5 21:27:51 MST 2017
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/20170105/0ad7e672/attachment.html
-------------- next part --------------
A non-text attachment was scrubbed...
Name: wa_evo_test591855.out
Type: application/octet-stream
Size: 5706 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170105/0ad7e672/attachment.obj
-------------- next part --------------
A non-text attachment was scrubbed...
Name: plt_wa_time_evo_Aug27_test.ps
Type: application/postscript
Size: 138045 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170105/0ad7e672/attachment.ps
More information about the ncl-talk
mailing list