| Re: Contour data with two different lon/lat grids onto the same map projection [message #71300 is a reply to message #71239] |
Tue, 08 June 2010 18:59   |
tessp
Messages: 8
Registered: June 2010
|
Junior Member |
|
|
> Well, I take it by "script" you mean an IDL main
> program (that is, not a procedure or function).
> There is plenty in just this small code sample
> to cause me to well believe such a script might
> cause you trouble.
>
> In fact, in spite of Mike Galloy proclaiming here the
> other day that he often uses main programs, I rarely
> do. And for exactly this reason. They can get you
> in a world of strange hurt. You would probably find
> out in a minute or two why you can't run the program
> more than once if you made a procedure out of it.
>
> I'm mostly confused by statements like this, though:
>
> "if I load the variables from netcdf to idl variables..."
>
> When you read the netCDF variables out of the file they ARE
> IDL variables. No need to load anything. Any why is the
> graphic being drawn over and over again in the same window
> in a loop!?
>
> I'm just not sure what to make of all this. Especially when
> I see the limits on the MAP_SET command, which have no logical
> connection to the latitude and longitude coordinates discussed
> in your previous communications.
>
> Perhaps this is just too small a code fragment to visualize the
> rest of it. But, I would say, if it works at least once, I'd be
> tempted to let it go as a good show. :-)
>
> Cheers,
>
> David
>
> --
> David Fanning, Ph.D.
David - Yes, script is my word for main program. I write a heap of
different programs to do things, and haven't ever got around to making
them into proceudres; although I agree with Mike that I wouldn't
change the code much if at all, were I to do so....This is mainly due
to time pressures, and the fact that once I have done something and
got the plots I need it's generally on to the next thing on the list.
The code is generally not pretty, but it works - usually. :o)
When I talk about loading the netcdf variables to IDL variables that
is just my clumsy terminology for reading the data out of the netcdf
file.
The map is being drawn over and over again to a postscript file in a
loop because I am plotting each timestep of the data, usually every 6
hours, to get a sequence of maps showing the development of the
atmospheric variables I'm looking at. Stepping through the .ps file
then shows how weather systems develop over the time period of
interest.
The limits on the map_set command are to map just the Australian
region. In my full code I have options to do this, or the whole globe,
or just the southern hemisphere.
Andy from Reading has my full code and data, and is taking a look at
this for me....so he may have some insight into the issues I'm having
when running it more than once.
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
