>
>
> The vertical integration routine "vibeta" seems to be causing a floating point
> error on our COMPAQ Alpha machines running Red Hat Linux. The same code works
> fine on our Solaris machines.
>
> How can I further debug this?
>
> Are there any other routines that do vertical integrals in pressure
> coordinates?
>
> Is it possible to get the (FORTRAN? C?) code to vibeta so that I can wrap it
> myself and see what's going on?
>
> Joe Barsugli :: jjb@cdc.noaa.gov :: 303-497-6042
> Climate Diagnostics Center :: University of Colorado
Hi Joe,
This problem, as you discovered, turned out to be something particular
to Compaq systems. It was choking on a missing value test. . I fixed
the problem, and put a new NCL binary on our download site:
http://ngwww.ucar.edu/ncl/download.html
Once you fill out the form and click on "Submit", a table with a list
of binaries will be presented. Choose the "Compaq Alpha running OSF
5.1" binaries, and follow the instructions for the installation, and
you should be set.
--Mary
>
> -- Example
>
> Here is the screen output:
>
> [jjb@cr01 ~/output]$ ncl < daily_psfc.ncl
> Copyright (C) 1995-2002 - All Rights Reserved
> University Corporation for Atmospheric Research
> NCAR Command Language Version 4.2.0.a023
> The use of this software is governed by a License Agreement.
> See http://ngwww.ucar.edu/ncl/ for more details.
> Floating exception (core dumped)
>
> Here is a test program bombs:
> ;
> load "/usr/local/ncarg/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "/usr/local/ncarg/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "/usr/local/ncarg/lib/ncarg/nclscripts/csm/shea_util.ncl"
> ;
> begin
> ccm3 = addfile("daily_debug.nc","r" )
> pm = ccm3->pm
> hr = ccm3->hr
> ps = ccm3->ps
> interpolation = ccm3->interp
> Top = 150.
> vhr = vibeta(pm(time|:,lat|:,lon|:,lev|:),hr(time|:,lat|:,lon|:,lev|:),interpolation,ps,110000.,Top*100.)
> end
>
> I'll put the input file at http://www.cdc.noaa.gov/~jjb/ncl/daily_debug.nc if
> anyone wants to test this.
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk@ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
_______________________________________________
ncl-talk mailing list
ncl-talk@ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
This archive was generated by hypermail 2b29 : Thu Aug 01 2002 - 11:14:08 MDT