[Thread Prev][Thread Next][Index]

[ferret_users] Correlation between one point and 2D spatial data



Dear Ferret users,

I am trying to compute a correlation between an observed current at one location (1D data: varies in time only) with SSH of a basin-wide region (2D data in x,y and time).

What I did is something like this:

-------------------
use daily-ssh.nc
use daily-current.nc
set mem/size=600

def axis/t=1-jan-2009:19-May-2010:1/unit=days tday
let ssh = sla[d=1, gt=tday]
let uvel = u[d=2, gt=tday]

let p = uvel
let q = ssh
go variance

shade correl
!
can data/a
can var/a
can mem
--------------------

But I got an error like:

-------
 *** NOTE: Ambiguous coordinates on X axis: P*Q-P*Q + 1
 *** NOTE: Ambiguous coordinates on X axis: Q*MASK_PQ
 *** NOTE: Ambiguous coordinates on X axis: COVAR / (P_VAR_MASK*Q_VAR_MASK)^.5
 *** NOTE: Ambiguous coordinates on X axis: P_DEV_MASK * Q_DEV_MASK
 **ERROR: dimensions improperly specified: must be a 2D region
--------

If I change between P and Q, I will get an error like:

--------
*** NOTE: Ambiguous coordinates on X axis: P*Q-P*Q + 1
 *** NOTE: Ambiguous coordinates on X axis: Q*MASK_PQ
 *** NOTE: Ambiguous coordinates on X axis: COVAR / (P_VAR_MASK*Q_VAR_MASK)^.5
 *** NOTE: Ambiguous coordinates on X axis: P_DEV_MASK * Q_DEV_MASK
 **ERROR: inconsistent sizes of data regions: X axis
          Q has 3 points (I=0:2)
          expression has 301 points (I=001:301)
---------

Do you have any suggestion?
Thank you for help.

Best regards,
Iskhaq



---------------------------------------------------------------- Dr. Iskhaq Iskandar Jurusan Fisika, Fakultas MIPA Universitas Sriwijaya Kampus Inderalaya, Ogan Ilir (OI) Sumatra Selatan, 30662 INDONESIA Phone: +62 (711) 580743 Fax: +62 (711) 580268 E-mail: iskhaq.iskandar@xxxxxxxxx http://iskhaq.wordpress.com/ ----------------------------------------------------------------


[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce / NOAA / OAR / PMEL / Ferret

Privacy Policy | Disclaimer | Accessibility Statement