comp.lang.idl-pvwave archive: archive » Re: How to add 'd' to get the correct julian conversion ?

Home » Public Forums » archive » Re: How to add 'd' to get the correct julian conversion ?

Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend

Re: How to add 'd' to get the correct julian conversion ? [message #33725]

Thu, 23 January 2003 15:38

thompson
Messages: 584
Registered: August 1991

Senior Member

Kolbjorn Bekkelund <kolbjorn@arctic-linux.tnett.no> writes:

> Craig Markwardt wrote:
>> Kolbjorn Bekkelund <kolbjorn@arctic-linux.tnett.no> writes:
>>
>>
>>> How can I add the NEEDED d to get this:
>>>
>>> 2452662.305203d
>>>
>>> out of this:
>>> maxtime = jul2cal((data(0,maxgust_time)), /TO_STRING, /MDY)
>>>
>>> In my program (data(0,maxgust_time)) fetches 2452662.305203 out of the
>>> array, but if I don't add the d to the julian date it calculates the
>>> wrong time in the above statement.
>>
>>
>> You can use
>> double(data(0,maxgust_time)),
>> but the variable DATA should already be in double precision. At least
>> it should be if you expect 13 decimal digits of precision to be
>> maintained. When you type the number directly on the command line,
>> you probably do have to use the "D" to indicate double precision, but
>> you should not have to if the variable DATA is already double.
>>
>> Craig
>>

> I've checked my array a bit more and it seems as if there's something
> wrong with it. From the file I'm reading in with read-ascii I should
> have this:
> 2452662.499876 2.719500 6.216000 343.494000
> 955.793400 93.911600 -5.444307

> but the print, data in IDL shows:
> 2.45266e+06 2.71950 6.21600 343.494 955.793
> 93.9116 -5.44431

> If I replace the read-acsii with Reimar Bauers read_data_file I get:
> 2452662.5 2.7195000 6.2160000 343.49400
> 955.79340 93.911600 -5.4443070

> but as you see the julian date in the first element is wrong in both
> arrays. How can I do ensure that I get all digits inserted ?

I tried the following

IDL> a = 2452662.499876 ;Single precision
IDL> print,a
2.45266e+06
IDL> a = 2452662.499876d ;Double precision
IDL> print,a
2452662.5

It looks like read_data_file is reading the data correctly as double precision,
while read_ascii is apparently reading everything into as single precision. I
know that it looks like A is being rounded off in the second case, but that's
only because of the default format being used for printing. If you use an
explicit format, you can see more of the digits.

IDL> print,a,format='(F20.6)'
2452662.499876

If, on the other hand, the data was read in as single precision, it really will
be truncated.

IDL> a = 2452662.499876 ;Single precision
IDL> print,a,format='(F20.6)'
2452662.500000

Bill Thompson

Report message to a moderator

[Message index]

		Re: How to add 'd' to get the correct julian conversion ? By: thompson on Thu, 23 January 2003 15:38
		Re: How to add 'd' to get the correct julian conversion ? By: Paul Van Delst[1] on Thu, 23 January 2003 14:20
		Re: How to add 'd' to get the correct julian conversion ? By: David Fanning on Thu, 23 January 2003 14:15
		Re: How to add 'd' to get the correct julian conversion ? By: Kolbjorn Bekkelund on Thu, 23 January 2003 12:02
		Re: How to add 'd' to get the correct julian conversion ? By: Craig Markwardt on Thu, 23 January 2003 08:46
		Re: How to add 'd' to get the correct julian conversion ? By: David Fanning on Thu, 23 January 2003 08:19
		Re: How to add 'd' to get the correct julian conversion ? By: R.Bauer on Fri, 31 January 2003 10:54

Previous Topic:	How to eliminate trailing zeros in a string?
Next Topic:	PNGs with alpha channel?

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Fri Oct 10 15:00:48 PDT 2025

Total time taken to generate the page: 1.67933 seconds