On Mon, Aug 2, 2010 at 1:59 PM, Matt Brewer <span dir="ltr"><<a href="mailto:mcbrewer83@yahoo.com">mcbrewer83@yahoo.com</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div><div style="font-family:arial,helvetica,sans-serif;font-size:10pt;color:#000000">Folks,<br>I am using lats4d to convert NCEP's NARR data (grib) to netCDF. It only kind of worked. This is what I type into the command line: <br>
<br>lats4d -i narr-a_221_19790101_0000_000.sub -o <a href="http://narr2009.nc" target="_blank">narr2009.nc</a> -ftype ctl -lon -160 -90 -lat 25 65 -time 0z1jan2009 21z31dec2009 -mxtimes 2920 -levs 925mb -xvars ugrdprs vgrdprs ugrd10m vgrd10m msletmsl prmslmsl<br>
<br>This works fine except for one issue. It seems to be totally messing up the axis dimensions, causing data to be in incorrect lat/lon locations.Notice how I asked for data from 160W to 90W and 25N to 65N, so with the longitude, that is a span of 70 degrees, and for latitude, it is 40 degrees. So I should get back a geopot. height matrix of 121x211 (there are data points every 1/3rd degree). It seems, in matlab at least, that I get
an array of 122x212. And so what happens is the first point on the second array actually gets saved in the 212th place in the first array, so the second line starts out with what is supose to be the second point in the array. Make sense? When I put the nc file into grads, it has similar issues as well....so something wierd is going on when I run the lats4d script....I just don't know what. So if anyone has any ideas of what is goin on, please help!!!<br>
<br></div></div></blockquote><div><br></div><div>The first step is to make sure GrADS is reading your GRIB file correctly. Go ahead, start GrADS, open your GRIB file, plot it and see if it is doing it correctly:</div><div>
<br></div><div>ga-> open <span class="Apple-style-span" style="font-family: arial, helvetica, sans-serif; font-size: 14px; ">narr-a_221_19790101_0000_000.sub</span></div><div><span class="Apple-style-span" style="font-family: arial, helvetica, sans-serif; font-size: 14px; ">ga-> d something</span></div>
<div><span class="Apple-style-span" style="font-family: arial, helvetica, sans-serif; font-size: 14px; "><br></span></div><div><span class="Apple-style-span" style="font-family: arial, helvetica, sans-serif; font-size: 14px; ">If this is pre-projected data, make sure your PDEF is correct. Once you sort this out, the netcdf version should produce the same results. (Our regression tests do just that: reads a GRIB file, converts to netcdf, then reads the netcdf file back and checks if it gets exactly the same values.)</span></div>
<div><span class="Apple-style-span" style="font-family: arial, helvetica, sans-serif; font-size: 14px; "><br></span></div><div><span class="Apple-style-span" style="font-family: arial, helvetica, sans-serif; font-size: 14px; ">One important detail: when asking for help, you need to let us know the particular version of grads, lats4d, and your operating system; basically, the otput of "q config".</span></div>
<div><span class="Apple-style-span" style="font-family: arial, helvetica, sans-serif; font-size: 14px; "><br></span></div><div><span class="Apple-style-span" style="font-family: arial, helvetica, sans-serif; font-size: 14px; "> Arlindo</span></div>
<div><br></div><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div><div style="font-family:arial,helvetica,sans-serif;font-size:10pt;color:#000000">
Matt<br></div><br>
</div><br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br>Arlindo da Silva<br><a href="mailto:dasilva@alum.mit.edu">dasilva@alum.mit.edu</a><br>