Re: escorc _FillValue£¿

The formula used is
    cor = SUM [(X(t)-Xave)*(Y(t)-Yave)}]/(Xstd*Ystd)

(a) either X or Y is all the same value then
either Xstd or Ystd will be zero.

(b) all X or all Y have _FillValue

In either case _FillValue is returned.

===
I am not sure what function you are using to print the results.
If you are using write_matrix, something like

   write_matrix (x, "20f9.3", False)

then, the *** indicate that the field width [here, 9] is not enough
to print the 1e20 value.

If, "r" is the correlation you computed then just set

     r@_FillValue = -9.999 ; any nunber which can fit the field width

http://www.ncl.ucar.edu/Document/Functions/Built-in/write_matrix.shtml

yanling_1220 wrote:
> Hello everyone,
>
> I use "escorc" to compute the "pattern correlation" between grid1 and
> grid2,which all with _FillValue = 1e+20,I used grid1@_FillValue
> <mailto:tas_gcm@_FillValue> = 1e+20
> grid2@_FillValue <mailto:tas_gcm@_FillValue> = 1e+20
> correl=new((/37,59/),float,1e+20)
> correl@_FillValue <mailto:correl@_FillValue> =1e+20
> and then wtite the reture array to a ascii file,the script run without
> any error,and produce a ascii file.
>
> But there is a doubt for me in the ascii file,the content is some like
> these:
> 0.974 0.974 0.974 0.975 0.975 0.976 0.976
> 0.974 0.978******************************
> **************************************************************************************************************
> 0.990
> 0.990 0.990 0.990 0.990 0.990 0.990
> 0.990 0.990 0.990 0.990 0.990 0.990
>
> I think "*****" means escorc was unable to compute maybe because of the
> _FillValue in one or the two grids,but why it not with the
> correl@_FillValue <mailto:correl@_FillValue> I set? Is there a mothod to
> make the reture array the filled with 1e+20 instead of "*****" in where
> correlation can't compute,because I want to deal these values in other
> software.
>
> Any suggestion will be appreciated.
>
> Regrads,
> sunflower
>
>
>
> ------------------------------------------------------------------------
> ûÓйã¸æµÄÖÕÉíÃâ·ÑÓÊÏä,www.yeah.net <http://www.yeah.net/?from=footer>
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> 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 Sun Sep 06 2009 - 07:41:18 MDT