Re: tiegcm and timegcm .nc files

Don,

Well, here we are 2 months later ;-)  I really would like to
use your idv application if possible, especially since we went
to the trouble of adding a new vertical dimension coord to the 
CF conventions.

In my previous message, I left out part of the file name of the nc 
file I am trying to read w/ IDV. The file is:

http://download.hao.ucar.edu/pub/foster/murray/FOSTER.timegcm1.2.py2j001.nc

In our previous correspondence, you suggested it was having
problems with the units of time.  It is a double, with units and
attributes as in my previous message below. I thought these units
comply with the explanation in section 4.4 "Time coordinate" of
the CF conventions document. If you can suggest a modification 
of the units or any other part of the file, I can try it out..
Or maybe I should install a newer version..

BTW, I am more than willing to hike over to your office from CG1 
sometime if that would help...I could login to the linux workstation 
where I've installed IDV.  Thanks Don,

--Ben


>Date: Tue, 2 Jan 2007 10:50:18 -0700 (MST)
>From: Ben Foster <address@hidden>
>Subject: Re: tiegcm and timegcm .nc files
>To: address@hidden
>Cc: address@hidden
>MIME-Version: 1.0
>Content-MD5: jD9ojh0FgaixmOq7LoLtJQ==
>X-Keywords: 
>
>
>Hi Don,
>
>Here is the declaration for time:
>
>dimensions:
>        time = UNLIMITED ; // (6 currently)
>
>double time(time) ;
>                time:long_name = "time" ;
>                time:units = "minutes since    0-1-0 0:0:0" ;
>                time:initial_year = 2002 ;
>                time:initial_mtime = 80, 0, 0 ;
>                time:initial_day = 80 ;
>
>It would be easy to change the units if necessary -- I don't
>think any of our post-proc codes rely on that specific unit.
>
>I'm not sure why you could not find the files -- maybe the hao
>server was down related to the storms, or maybe I had a typo
>in the link I gave you. Anyway the files are still in
>
>http://download.hao.ucar.edu/pub/foster/murray
>
>This morning I added an ncdump of FOSTER.timegcm1.2.py2j001.nc.
>The above is from the .ncd file.
>
>Thanks for your help, and BTW, I can always walk over to your
>office (I'm at CG1) anytime if that would help,
>
>--Ben
>
>
>>Date: Tue, 26 Dec 2006 09:45:20 -0700
>>From: Don Murray <address@hidden>
>>User-Agent: Thunderbird 1.5.0.8 (Windows/20061025)
>>MIME-Version: 1.0
>>To: Ben Foster <address@hidden>
>>Subject: Re: tiegcm and timegcm .nc files
>>Content-Transfer-Encoding: 7bit
>>X-Virus-Scanned: amavisd-new at ucar.edu
>>X-Keywords: 
>>
>>Hi Ben-
>>
>>Sorry for the delay in responding.  I was out at AGU and then
>>the storm hit.
>>
>>My guess is that you have a bad unit specification for the time.
>>I tried accessing the file, but it said the link could not be
>>found.
>>
>>Don
>>
>>Ben Foster wrote:
>>> 
>>> Don,
>>> 
>>> I downloaded and installed it on my local linux workstation:
>>> 
>>> Linux oreo 2.6.18-1.2200.fc5smp #1 SMP Sat Oct 14 17:15:35 EDT 2006 i686 
>>> i686 i386 GNU/Linux
>>> 
>>> It installed fine (I think it is idv2.1RC3), and I set the Pressure to 
>>> Height
>>> option to vis5d. Then I added a new data source, namely a copy of
>>> http://download.hao.ucar.edu/pub/foster/murray/timegcm1.2.py2j001.nc.
>>> But it is "Unable to load field: TN", or any of the 3d fields. The error
>>> to stdout is as follows:
>>> 
>>> getGeoGridTimes() visad.VisADException: Error parsing 
>>> unit:visad.data.units.ParseException: Encountered " " at line 1, column 23.
>>> Was expecting one of:
>>>    <EOF>
>>>    <INTEGER> ...
>>>    <REAL> ...
>>>    <DIVIDE> ...
>>>    "." ...
>>> 
>>> I also went to View->Properties and set the vertical scale to
>>> -17 to +5 -- this is the new unitless vertical dimension (ln(p0/p))
>>> in the CF conventions)... I left the units in meters.  Maybe it doesnt
>>> like meters < 0?
>>> 
>>> Thanks,
>>> 
>>> --Ben
>>> 
>>> 
>>> Don Murray wrote:
>>>> Hi Ben-
>>>>
>>>> Ben Foster wrote:
>>>>> Don,
>>>>>
>>>>> How is this going??
>>>>
>>>> We're swamped getting ready for the workshop so this took
>>>> a back seat.  Here's where your stuff stands:
>>>>
>>>> Download the latest nightly release from:
>>>>
>>>> http://www.unidata.ucar.edu/software/idv/release/nightly
>>>>
>>>> and install it.
>>>>
>>>> Before loading in your dataset, go to the Edit-Preferences
>>>> menu and select the Formats & Data tab.  Set the
>>>> Pressure to Height option to be Vis5D instead of
>>>> Standard Atmosphere.
>>>>
>>>> Now, load in your datasets and create displays.  You'll
>>>> also need to set your vertical scale to be in the range of
>>>> your datasets.  Use the View->Properties dialog to do this.
>>>>
>>>> Let me know if you run into any problems.
>>>>
>>>> Unfortunately, our workshop is full so I can't guarantee
>>>> you a seat, but you could sit and audit it if you wanted
>>>> (if you had a laptop, you could work along with us on that).
>>>>
>>>> Don
>>>> *************************************************************
>>>> Don Murray                               UCAR Unidata Program
>>>> address@hidden                        P.O. Box 3000
>>>> (303) 497-8628                              Boulder, CO 80307
>>>> http://www.unidata.ucar.edu/staff/donm
>>>> *************************************************************
>>>>
>>>>
>>
>>-- 
>>*************************************************************
>>Don Murray                               UCAR Unidata Program
>>address@hidden                        P.O. Box 3000
>>(303) 497-8628                              Boulder, CO 80307
>>http://www.unidata.ucar.edu/staff/donm
>>*************************************************************
>>
>>
>
>-----------------------------------------------------------------------
>Ben Foster                     High Altitude Observatory (HAO)
>address@hidden                 phone: 303-497-1595  fax: 303-497-1589  
>Nat. Center for Atmos. Res.     P.O. Box 3000 Boulder CO 80307 USA
>-----------------------------------------------------------------------
>

-----------------------------------------------------------------------
Ben Foster                      High Altitude Observatory (HAO)
address@hidden                  phone: 303-497-1595  fax: 303-497-1589  
Nat. Center for Atmos. Res.     P.O. Box 3000 Boulder CO 80307 USA
-----------------------------------------------------------------------