Re: Question for inverse_matrix

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Wed Apr 02 2014 - 10:10:37 MDT

> I have a problem about using inverse_matrix.
> The longitude and latitude of final output I
> got from my code was reversed.

NCL makes a direct call to the LAPACK routines DGETRF and DGETRI
which are written in fortran. NCL is row major, fortran is column major.

   NCL: x(lat,lon) <===> x(lon,lat): fortran

   Maybe you should reorder: eg

   xr = x(lon|:,lat|:) ; use cr

==
It would be better to include a simple example to illustrate the issue.

==

missing values

The LAPACK routines DGETRF and DGETRI do not contain code
to handle missing values.

http://www.netlib.no/netlib/lapack/double/dgetrf.f
http://www.netlib.no/netlib/lapack/double/dgetri.f

Note: this is not an NCL issue ... Matlab and other
software tools are unlikely to have an 'inverse_matrix'
function that handles missing values (NaN in Matlab).

You must fill in the values prior to calling the code.
Maybe setting them to 0.0 will work .. I don't know.

On 4/1/14, 12:11 PM, 김미애 wrote:
> Hello ncl team,
>
> I have a problem about using inverse_matrix.
> The longitude and latitude of final output I got from my code was reversed.
> I tried to use 'inverse_matrix' function. But, I got an error message as below.
> warning:inverse_matrix: info = 11; missing values not allowed
> What does this mean? How can I solve this problem?
>
> Thank you so much
>
> Sincerely,
>
>
> Miae Kim
>
>
>
> -------------------------------------------------
> Miae Kim
> School of Urban and Environment Engineering
> Ulsan National Institute of Science and Technology (UNIST)
> UNIST-gil 50, Ulsan 689-798, South Korea
> Mobile: +82-10-9421-8287
> E-mail: toa0710@naver.com
>
>
>
> _______________________________________________
> 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 Wed Apr 2 10:10:40 2014

This archive was generated by hypermail 2.1.8 : Thu Apr 03 2014 - 13:36:27 MDT