Re: problem of WRAPIT

   Thanks Mary and Rick for your help very much!
   
   Now there are no problem on WRAPIT the rip_sub.f, but another problem coming. I got "Segmentation fault" when I calling the rip_sub.so, I have check my ncl script carefully but I can't find where s the error. I'm using the wrfout of WRFV2.2, and following is the information:
_________________________________________________________________________
[yangyh@gateway ncl]$ ncl cptp.ncl
 Copyright (C) 1995-2009 - All Rights Reserved
 University Corporation for Atmospheric Research
 NCAR Command Language Version 5.1.0
 The use of this software is governed by a License Agreement.
 See http://www.ncl.ucar.edu/ for more details.
(0) DX,DY = 9000 9000
(0) Working on time: 2010012505 3
(0) llz=34
(1) llz=264
(2) llz=264
(0) llz=34
(1) llz=264
(2) llz=264
(0) llz=34
(1) llz=264
(2) llz=264
(0) llz=34
(1) llz=264
(2) llz=264
(0) llt=264
(1) llt=264
(0) ll=34
(1) ll=264
(2) ll=264
(0) ll==34
(1) ll==264
(2) ll==264
(0) UUUUUU UUU
Segmentation fault
____________________________________________________________________________________________

   Could you please check the attachment cptp.ncl for me?
   
   By the way,following is the purpose of Why I using rip_sub.f :
   The "wrf_cape_2d" works well on my machine, but I'd like to get the temperature of Equilibrium level,which doesn't output by wrf_cape_2d, so I have to call the rip_sub.f in order to get extra information. Does the function wrf_cape_2d is coming from rip_cape.f?

  Thanks again!!

  Yuhua
   

  ----- Original Message -----
  From: Mary Haley
  To: Rick Brownrigg
  Cc: Yuhua Yang ; ncl-talk@ucar.edu
  Sent: Tuesday, January 26, 2010 1:10 AM
  Subject: Re: [ncl-talk] problem of WRAPIT


  Thanks for catching this Rick! I was not aware of this restriction.


  I've added this to the "special considerations" section of the WRAPIT documentation, which is a good document to read if you are having problems with WRAPIT:


     http://www.ncl.ucar.edu/Document/Tools/WRAPIT.shtml


  --Mary


  On Jan 25, 2010, at 9:43 AM, Rick Brownrigg wrote:


    Hi,


    The root of this error is that you declare the types of your two embedded functions as:


    FUNCTION TONPSADIABAT(THTE,PRS,...
            DOUBLE PRECISION TONPSADIABAT
    ...


    FUNCTION VIRTUAL(TEMP,RATMIX)
    DOUBLE PRECISION VIRTUAL
    ...


    While this is legitimate fortran code, its causing an error in the WRAPIT parsing code. Try declaring the types of the functions with one statement, as:


    DOUBLE PRECISION FUNCTION TONPSADIABAT(THTE,PRS,...
    ...
    DOUBLE PRECISION FUNCTION VIRTUAL(TEMP,RATMIX)




    Rick




    On Jan 25, 2010, at 3:12 AM, Yuhua Yang wrote:


      Hi,

        When I WRAPIT the attachment file(rip_cape.f) which coming from Mary, I got the following problem:
      ______________________________________________________
       WRAPIT rip_cape.f
       
      WRAPIT Version: 090115
      A syntax error occurred while parsing: TONPSADIABAT
      COMPILING rip_cape.f
      LINKING
      END WRAPIT
      ________________________________________________________

       The rip_cape.so file is created, but when I try to call the routine by my NCL script , I get the warning:
      Could not find Init() in external file ./rip_cape.so, file not loaded .

      the enviroment of my machine is linux,64-bit.

      Could you give me some hint to solve this problem? Thanks very much!

      Yuhua


        <rip_cape.f>
      _______________________________________________
      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

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk