[Thread Prev][Thread Next][Index]

Re: [ferret_users] Plot.



Hi,
First, I have a general suggestion (for all of us).  When you write a question, please use an informative subject. In the subject line, it is very helpful if you include a short description of the question.  For example with this question, instead of the subject that says "Plot", you might have said something like "error comparing datasets".


When comparing data from two variables that have different grids, you'll need to put the data from one dataset onto the grid of the other.  So, you will need to decide on the grid for your result, and then change your script, so that use substitute either

   POTTMP[d=2,g=temperature[d=1]]  instead of POTTMP[d=2]

or use

   TEMPERATURE[d=1,g=POTEMP[d=2]]  instead of TEMPERATURE[d=1]


Often you can plot two variables on different grids on the same plot. I think your commands that overlay a contour line on a fill plot will work, because the coordinate axes in the grids have the same units: calendar time, longitude, latitude, depth in meters.  So I can't test them without the data, but these commands are probably fine:

   yes? fill/nolabel TEMPERATURE[d=1]
   yes? contour/ov/nolabel/colour=7 TEMPERATURE[d=1]

but combining data from two different grids in one _expression_

   yes? let dT=TEMPERATURE[d=1,x=0.5]-POTTMP[d=2]

will not in general work correctly.

-Ansley

On 4/1/2015 2:12 AM, Satyam Srivastava wrote:
Hi Folks,

I am trying to compare two datasets and trying to plot their difference as well i am using the scrip as mentioned below but i am getting the error as:::

 """""NOTE: Ambiguous coordinates on Y axis: TEMPERATURE[D=1,X=0.5]-POTTMP[D=2]
 *** NOTE: Ambiguous coordinates on Z axis: TEMPERATURE[D=1,X=0.5]-POTTMP[D=2]
 *** NOTE: Ambiguous coordinates on T axis: TEMPERATURE[D=1,X=0.5]-POTTMP[D=2]
 **ERROR: illegal limits: POTTMP is not in the range T=13-MAY-3807 00:00:03-JUN-3807 00:00
          Axis extremes are T=16-DEC-2007 12:00:16-DEC-2008 00:00"""""
 Script is:::


yes? set mode metafile
yes? use pottm2008.nc
yes? use pottmp.2008.nc
yes? set reg/x=0.5E/y=35S/t=10-may-2008:31-may-2008/z=10:40
yes? let dT=TEMPERATURE[d=1]-POTTyes? use pottmp.2008.nc
yes? show grid pottmp
    GRID GCN2
 name       axis              # pts   start                end
 LON       LONGITUDE          360mr   0.5E                 0.5W
 LAT       LATITUDE           418 r   74.5S                64.499N
 LEVEL     DEPTH (m)           40 i-  5                    4478
 TIME      TIME                12 i   01-JAN-2008 00:00    01-DEC-2008 00:00
 normal    E
 normal    F
MP[d=2]
yes? let dT=TEMPERATURE[d=1,x=0.5]-POTTMP[d=2]
yes? set window/size=.9
yes? set viewport UL
yes? fill/nolabel TEMPERATURE[d=1]
yes? contour/ov/nolabel/colour=7 TEMPERATURE[d=1]
yes? set viewport UR
yes? fill/nolabel POTTMP[d=2]
yes? contour/ov/nolabel/colour=7 POTTMP[d=2]
yes? set viewport LL
yes? fill/nolabel/pal=mycol3 dT



the data pottm2008.nc grids are as::


yes? use pottm2008.nc
yes? show grid temperature
    GRID GOZ1
 name       axis              # pts   start                end
 LONN2551_N551 LONGITUDE     2001mr   130W(-130)           30E(30)
 LAT959_2051 LATITUDE        1093 i   39.976S              39.976N
 LEV       DEPTH (m)           10 i-  0                    200
 TIME1     TIME                32 r   08-MAY-2008 00:00    08-JUN-2008 00:00
 normal    E
 normal    F


And for the pottmp.2008.nc is as follows:::

yes? use pottmp.2008.nc
yes? show grid pottmp
    GRID GCN2
 name       axis              # pts   start                end
 LON       LONGITUDE          360mr   0.5E                 0.5W
 LAT       LATITUDE           418 r   74.5S                64.499N
 LEVEL     DEPTH (m)           40 i-  5                    4478
 TIME      TIME                12 i   01-JAN-2008 00:00    01-DEC-2008 00:00
 normal    E
 normal    F


could you guys tell me where i am going wrong or why i am getting this error.

Regards,
Satyam Srivastava.



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

Privacy Policy | Disclaimer | Accessibility Statement