scatter plot and regression line for 2D vs 2D array

Good morning!

I am plotting scatter plots of two 2D datasets and finding their regression
line. The 2D arrays are:

V_obs(90,144)
V_obs@_FillValues = -999

V_model(90,144)

What I'm doing firstly is to change these two 2D arrays to 1D, and then use
the following method to create plot:

 plot = gsn_csm_xy (wks,data_obs,data,res)

The 1D array data_obs will be the X-axis, 1D data_model will be the Y-axis.
 Where,

rc = regline(data_obs,data_model)
data(0,:) = data_model
data(1,:) = rc*(data_obs - rc@xave) + rc@yave

To change the datasets from 2D to 1D, I used:

  i_data = - 1
  do isat = 0,90-1
   do jsat = 0,144-1
     i_data = i_data + 1
     data_model(i_data) = V_model(isat,jsat)
     data_obs(i_data) = V_obs(isat,jsat)
   end do
  end do

The error message I got is:

fatal:Number of dimensions on right hand side do not match number of
dimension in left hand side
fatal:["Execute.c":7556]:Execute: Error occurred at or near line 5263 in
file $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl

fatal:["Execute.c":7556]:Execute: Error occurred at or near line 177 in
file scatter_sat_model.ncl

The line 177 is:
    gsn_text(wks,plot,text_s,25,23,res_text)

Could you please assist me to find the problems in my code? I'm also not
sure if the places with _FillValues in the observation are automatically
treated as missing value in model data too.

Thank you!

Hongyan

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Sep 17 08:34:18 2012