Re: wrapit trouble

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Fri Aug 05 2011 - 10:20:47 MDT

[1] I did not see any
C NCLFORTSTART
...
C NCLEND
      in your code. I assume you sent me the original code.

[2] NCL's builtin f77 parser which is used by WRAPIT
      is not as sophisticated as a fortran compiler. It
      can not evaluate expressions like the following:

dimension xw((-(m-1)/2+1):(n+(m-1)/2)),ck(-(m-1)/2:(m-1)/2)

[3] A simple work around is to create a 'stub' or interface subouttine

C NCLFORTSTART
        subroutine gfiltncl(n,m,x,y)
        implicit none
        integer n,m
        real x(n),y(n)
C NCLEND
        call gfilt(n,m,x,y)
        return
        end

        subroutine gfilt(n,m,x,y)
        dimension x(n),y(n)
        dimension xw((-(m-1)/2+1):(n+(m-1)/2)),ck(-(m-1)/2:(m-1)/2)
[snip]

[4] %> WRAPIT Gfil.f

WRAPIT Version: 110125
COMPILING Gfil.f
LINKING
END WRAPIT

[5] NCL code

external GFILT Gfil.so

GFILT::gfiltncl(n,m,x,y)

On 08/05/2011 09:49 AM, Coly SAR wrote:
> please find the file below:
>
> Coly
>
> 2011/8/5 Dennis Shea <shea@ucar.edu <mailto:shea@ucar.edu>>
>
> You should provide a sample of thye fortran code you are wrapping.
>
>
> On 08/05/2011 09:40 AM, Coly SAR wrote:
>
> Hi all,
>
> I am learning wrapit with a fortran subroutine but I am having : A
> syntax error occurred while parsing:
>
[SNIP]

>
> Any hint, on this useful topic I am learning.
>
> Best Regards,
>
> Coly
>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Aug 5 10:20:52 2011

This archive was generated by hypermail 2.1.8 : Fri Aug 05 2011 - 14:53:54 MDT