Re: help: when wrapit comes across allocatable array

Dennis,

So, as far as I know, the stub file just need to declare the parameters as
arguments of subroutine, and it doesn't have to declare all of the
variables in the subroutine. Right?

Regards,
Sha

> See: http://www.ncl.ucar.edu/Document/Tools/WRAPIT.shtml
>
> You would just need one f77 style stub file [ mystub.f].
> Stub files should be created only for subroutines
> invoked from within NCL explicitly.
>
> C NCLFORTSTART
> subroutine cstub(nlev,temp,pres,alt,pliml, plimu,plimlex,mis
> + ,numtp,tp,tt,ta)
> implicit none
> integer nlev
> real temp(nlev), pres(nlev), alt(nlev), pliml, plimu,plimlex,mis
> integer numtp
> real tp(nlev), tt(nlev), ta(nlev)
> C NCLEND
>
> %> WRAPIT mystub.f ShaFeng.f90
>
> external CSTUB "./mystub.so"
>
> CSTUB::cstub(...)
>
> On 9/29/11 3:32 PM, Sha Feng wrote:
>> Dennis,
>>
>> My code is .f90. So I don't have to use the same declaration with the
>> .f90 files in the .stub, do I?
>>
>> A snippet of my code is also follows.
>> !==================================
>> subroutine C(...)
>> implicit none
>> integer,intent(in) :: nlev
>> real,intent(in),dimension(nlev) :: temp
>> real,intent(in),dimension(nlev) :: pres
>> real,intent(in),dimension(nlev) :: alt
>> real, intent(in) :: pliml, plimu,plimlex
>> integer,intent(out) :: numtp
>> real,intent(in) ::mis
>> real,intent(out),dimension(nlev) :: tp,tt,ta
>>
>> integer :: i, j, validrow
>> real,dimension(nlev) :: t,p,a
>> real, dimension(:), allocatable :: t1,p1,a1
>> real :: trp
>> ....
>> call A(...)
>> ...
>> call B(...)
>> ...
>> end subroutine C
>>
>> subroutine A(...)
>> ...
>> end subroutine A
>>
>> subroutine B(....)
>> ...
>> end subroutine B
>> !==================================
>>
>>
>> BTW, this subroutine includes three subroutines (three subroutine in one
>> .f90 file). Should I make declaration for all of them in the .stub file.
>>
>>
>>
>>
>>
>> On Sep 29, 2011, at 5:10 PM, Dennis Shea wrote:
>>
>>> It would be good to see a snippet of your f90 code.
>>>
>>> As along as the 'allocateable' is not between the
>>> C NCLFORTSTART
>>> C NCLEND
>>>
>>> Also, is the f90 code .f or .f90?
>>>
>>> C ----------Sample for .f file -----------------------------------
>>> C NCLFORTSTART
>>> subroutine wgt_runave (x,npts,wgt,nwgt,kopt,xmsg,ier)
>>> integer npts, nwgt,kopt,ier
>>> real x(npts), wgt(nwgt), xmsg
>>> C NCLEND
>>> integer lwork
>>> real, allocatable, dimension(:) :: work
>>>
>>> nhalf = nwgt/2
>>> lwork = npts+2*nhalf
>>> allocate (work(lwork), stat=ier)
>>>
>>> call wrunavx77 (x,npts,wgt,nwgt,kopt,xmsg,work,lwork,ier)
>>>
>>> deallocate (work)
>>>
>>> return ! NCL return (wrunwgt[x])
>>> end
>>>
>>>
>>>
>>>
>>> On 9/29/11 12:19 PM, Sha Feng wrote:
>>>> Hi ncl-talker,
>>>>
>>>> How can I deal with the codes which are related with allocatable
>>>> arrays when I use WRAPIT to combine Fortran 90 subroutines with NCL
>>>> scripts?
>>>>
>>>> Any suggestion is very much appreciated.
>>>>
>>>> Sha
>>>> _______________________________________________
>>>> ncl-talk mailing list
>>>> List instructions, subscriber options, unsubscribe:
>>>> 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 Thu Sep 29 20:50:11 2011