[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