[Met_help] Re: ncmod_gallus script - question (fwd)
John Halley Gotway
johnhg at rap.ucar.edu
Tue Jul 29 15:59:57 MDT 2008
Bill,
OK, try the version of ncmod_gallus in /d1/johnhg/ncmod_gallus.
I'm just looking at the first dimension in the NetCDF file. If it's 169 I'm using the original grid definition. If it's 99, I'm using the new one.
Let me know how it goes.
And please keep in mind that MET requires the forecast and observation files to be on the same grid for the Grid-Stat and MODE tools.
I'll be out of town for 10 days starting Monday, so if you need help setting anything up, please let me know by Friday. But Randy should be here to answer any questions that come up.
Good luck.
John
wgallus at iastate.edu wrote:
> My grad student who produced the data gives some info below. He didn't
> know what r_km would be. It seems to me it might be the radius of the
> earth, or something?
>
> Bill
>
> ------------------------------------------------
> | William A. Gallus, Jr. |
> | Professor of Meteorology |
> | Professor-in-charge, Meteorology Program |
> | 3025 Agronomy |
> | Iowa State University |
> | 515-294-2270 |
> ------------------------------------------------
>
> ---------- Forwarded message ----------
> Date: Tue, 29 Jul 2008 16:07:20 -0500 (CDT)
> From: Adam Clark <clar0614 at iastate.edu>
> To: wgallus at iastate.edu
> Subject: Re: ncmod_gallus script - question (fwd)
>
> The stuff below:
>
> // Set up the Lambert Conformal data for the grid to be used
> LambertData lc_data = {
> "gallus", // name
> 42.0, // secant lat 1
> 42.0, // secant lat 2
> 34.821, // lower-left lat
> 104.924, // lower-left lon
> 95.0, // lcen
> 10.379, // d_km
> 6367.47, // r_km
> 169, // nx
> 149 // ny
> };
>
> should be changed to:
>
> 40.5, // secant lat 1
> 40.5, // secant lat 2
> 31.9, // lower-left lat
> 106.50, // lower-left lon
> 97.0, // lcen
> 20.000, // d_km
> ??????, // r_km
> 99, // nx
> 92 // ny
> };
>
>
> The only thing I don't know is the "r_km" variable. Maybe John would
> know what that would be based on the other variables??
>
> I got most of the info above from the GrADS .ctl file I have to read the
> 4 and 20-km data. The line in the .ctl file defining the grid is:
>
> pdef 99 92 lcc 31.900 -106.50 1.000 1.000 40.500 40.500 -97.000
> 20000.000 20000.000
>
> Hope that helps!
>
> Adam
>
>
>
>
>
> On Tue, 29 Jul 2008, wgallus at iastate.edu wrote:
>
>> Adam,
>>
>> I tried to run the code to convert things to the type of nc file that
>> MODE could use this morning, and didn't realize the script had been
>> hardwired for the 15 km runs last year. See the email below. Can you
>> provide all the info that would be valid for your 20 km ensemble
>> datasets?
>>
>> Bill
>>
>> ------------------------------------------------
>> | William A. Gallus, Jr. |
>> | Professor of Meteorology |
>> | Professor-in-charge, Meteorology Program |
>> | 3025 Agronomy |
>> | Iowa State University |
>> | 515-294-2270 |
>> ------------------------------------------------
>>
>> ---------- Forwarded message ----------
>> Date: Tue, 29 Jul 2008 14:02:46 -0600
>> From: John Halley Gotway <johnhg at rap.ucar.edu>
>> To: wgallus at iastate.edu
>> Subject: Re: ncmod_gallus script - question
>>
>> Bill,
>>
>> I was able to find that file. I had been looking in /d1 and forgot
>> all about /d2.
>>
>> Anyway, yes it is "hardwired" in there. The previous data we were
>> using resided on 169x149 Lambert Conformal grid defined below:
>>
>> // Set up the Lambert Conformal data for the grid to be used
>> LambertData lc_data = {
>> "gallus", // name
>> 42.0, // secant lat 1
>> 42.0, // secant lat 2
>> 34.821, // lower-left lat
>> 104.924, // lower-left lon
>> 95.0, // lcen
>> 10.379, // d_km
>> 6367.47, // r_km
>> 169, // nx
>> 149 // ny
>> };
>>
>> The data in this new file resides on a grid that's 99x92. We'll need
>> to figure out the grid definition for this data and add it into the
>> tool. Any ideas how this grid should be defined?
>>
>> John
>>
>>
>> wgallus at iastate.edu wrote:
>>> It is on ebony in /d2/gallus/4km.
>>>
>>> Let me know if you can't find it there and I can email it.
>>>
>>> Bill
>>>
>>> ------------------------------------------------
>>> | William A. Gallus, Jr. |
>>> | Professor of Meteorology |
>>> | Professor-in-charge, Meteorology Program |
>>> | 3025 Agronomy |
>>> | Iowa State University |
>>> | 515-294-2270 |
>>> ------------------------------------------------
>>>
>>> On Tue, 29 Jul 2008, John Halley Gotway wrote:
>>>
>>> > Bill,
>>> > > Could you please send me that file: 20070502_p106_27.dat.nc
>>> > > Or let me know if it's sitting on one of the DTC machines.
>>> > > And I'll try running the tool to see what needs to be adjusted.
>>> > > John
>>> > > wgallus at iastate.edu wrote:
>>> > > John,
>>> > > > > I'm getting ready for my DTC visit the next 2 weeks, and
>>> tried to run > > the ncmod_gallus script you created for me, but
>>> this time, for some new > > output from 4 km CAPS runs, that is on a
>>> 20 km grid. I believe I have > > the output in the same .nc format
>>> as the output I used last year, > > except the dimensions of the
>>> grid are different. In any event, when I > > tried to run
>>> ncmod_gallus on the data, I got this error...
>>> > > > > Reading input netCDF file: ../4km/20070502_p106_27.dat.nc
>>> > > ncvarget: ncid 3; varid 2: Start+count exceeds dimension bound
>>> > > > > Is your ncmod_gallus script hardwired for the dimensions of
>>> the grid I > > had been using previously? If so, I guess I'll need
>>> an adjusted > > version that I can use now for these other datasets.
>>> > > > > Bill
>>> > > > > ------------------------------------------------
>>> > > | William A. Gallus, Jr. |
>>> > > | Professor of Meteorology |
>>> > > | Professor-in-charge, Meteorology Program |
>>> > > | 3025 Agronomy |
>>> > > | Iowa State University |
>>> > > | 515-294-2270 |
>>> > > ------------------------------------------------
>>> >
>>
>
More information about the Met_help
mailing list