[IDV #KVE-675145]: gfs - cdc grid

Hi Bill-

> I hadn't used sub before; I like it. And it works.

This came out of the GEMPAK grid diagnostic work.  I would recommend
using sub, mul, quo, add instead of -,*,/,+ since they call
GridMath routines which will handle the 2D/3D differences.

> BTW, the solution should be return sub(data,cgpm) not return sub(d2d,cgpm).
> Gotta have it right in the archive!

Thanks for pointing that out and getting the correction into the archives.

Don

> FW: [IDV #KVE-675145]: gfs - cdc grid
> 
> 
> 
> 
> -----Original Message-----
> From: Unidata IDV Support [mailto:address@hidden]
> Sent: Mon 2/9/2009 2:41 PM
> To: Fingerhut, William A @ LSC
> Cc: address@hidden
> Subject: [IDV #KVE-675145]: gfs - cdc grid
> 
> Hi Bill-
> 
> 
> > For some unknown reason, I decided to try again to solve the
> problem
> > described below. And I was successful
> > by writing the following jython:
> >
> > def zAnom250(data,clim):
> >    # slice at 250mb, every other point, every 4th time
> >    d2d=make2D(data)
> >    cgpm=newUnit(clim,"Geopotential_height","gpm")
> >    return d2d-cgpm
> >
> > I remember you once telling me to use newUnit, but that alone did
> not
> > fix my problem.
> > Success came when I used make2D, which was a new idea after I saw
> it in
> > the jython library.
> > I don't know if you created make2D after I asked for help and I
> didn't
> > receive the info, or
> > ....... ??? Anyway, I thought I would let you know that the
> problem has
> > been solved and
> > I am smiling.  :-)
> 
> You shouldn't need make2D if you change the return to be:
> 
>     return sub(d2d,cgpm)
> 
> sub actually will call what make2D does if it needs to.
> 
> The newUnit is needed if I recall because one is meters and the other is
> geopotential meters which are incompatible units for subtraction.
> 
> Don
> 
> > [IDV #KVE-675145]: gfs - cdc grid
> >
> >
> >
> > To: william.fingerhut@xxxxxxxxxxxxxxx
> > Subject: [IDV #KVE-675145]: gfs - cdc grid
> > From: "Unidata IDV Support"
> <support-idv@xxxxxxxxxxxxxxxx>
> > Date: Wed, 05 Mar 2008 07:58:33 -0700
> > Delivered-to: address@hidden by
> > laraine.unidata.ucar.edu (Postfix) with ESMTP id E5381CB199; Wed,
> 5 Mar
> > 2008 07:58:33 -0700 (MST) id C944AD5116; Wed, 5 Mar 2008 07:58:33
> -0700
> > (MST)
> >
> >
> >
> >
> >
> > Hi Bill-
> >
> > > Full Name: bill fingerhut
> > > Email Address: william.fingerhut@xxxxxxxxxxxxxxx
> > > Organization: lyndon state college
> > > Package Version: 2.4 build date:2007-12-20 16:34 UTC
> > > Operating System: Windows XP
> > > Hardware: Java: home: C:\Documents and
> Settings\waf01060\Program Files\IDV_2.4\jre version: 1.5.0_11 j3d:1.3.2
> fcs (build12)
> > > Description of problem: It did not work. :-(
> > > Is the problem with the resampling ?
> >
> > No, the problem is that one of the domains (GFS) is latitude,
> longitude and  the
> > other is (CDC) longitude, latitude.  I have no solution for that
> other than using
> > the Global GFS grids which are longitude, latitude.  A long time
> ago, we should
> > have had projections output longitude, latitude order instead of
> latitude, longitude
> > order, but that cat is long out of the bag.
> >
> > Don
> >
> > > Creating display: Color-Filled Contour Plan View
> > > org.python.core.PyException
> > >
> > > ******************
> > > Stack trace:
> > > Traceback (innermost last):
> > > File "<string>", line 1, in ?
> > > at visad.FlatField.resample(FlatField.java:4592)
> > >
> > > at visad.FlatField.binary(FlatField.java:2088)
> > >
> > > at visad.FieldImpl.binary(FieldImpl.java:835)
> > >
> > > at visad.DataImpl.binary(DataImpl.java:139)
> > >
> > > at visad.DataImpl.subtract(DataImpl.java:192)
> > >
> > > at visad.DataImpl.__sub__(DataImpl.java:1400)
> > >
> > > at sun.reflect.NativeMethodAccessorImpl.invoke0(Native
> Method)
> > >
> > > at
> sun.reflect.NativeMethodAccessorImpl.invoke(NativeMethodAccessorImpl.java:39)
> > >
> > > at
> sun.reflect.DelegatingMethodAccessorImpl.invoke(DelegatingMethodAccessorImpl.java:25)
> > >
> > > at java.lang.reflect.Method.invoke(Method.java:585)
> > >
> > >
> > > visad.VisADException: visad.VisADException: Sampling set
> is not compatible with domain
> > >
> > >
> > >
> >
> >
> > Ticket Details
> > ===================
> > Ticket ID: KVE-675145
> > Department: Support IDV
> > Priority: Normal
> > Status: Open
> >
> >
> > --
> > The Fingerhuts
> > Saint Johnsbury, Vt, 05819
> >
> > email: address@hidden
> >
> >
> >
> >
> 
> 
> Ticket Details
> ===================
> Ticket ID: KVE-675145
> Department: Support IDV
> Priority: Normal
> Status: Open
> 
> 
> 
> 
> 
> 
> 
> --
> The Fingerhuts
> Saint Johnsbury, Vt, 05819
> 
> email: address@hidden
> 
> 
> 
> 


Ticket Details
===================
Ticket ID: KVE-675145
Department: Support IDV
Priority: Normal
Status: Open