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 #33838 is a reply to message #33729]

Fri, 31 January 2003 10:54

R.Bauer
Messages: 1424
Registered: November 1998

Senior Member

Kolbjorn Bekkelund wrote:

> 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 ?

Dear Kolbjorn

The problem I think you have is that's the default format for print is
defined for float numbers.
read_data_file uses as default double if you don't give a type.

So you should try something like

x=read_data_file('test.dat')
print, x.data[0],format='(F20.10)'
2452662.4998760000

best regards

Reimar

>
> Kolbjorn
>
>

--
Forschungszentrum Juelich
email: R.Bauer@fz-juelich.de
http://www.fz-juelich.de/icg/icg-i/
============================================================ ======
a IDL library at ForschungsZentrum Juelich
http://www.fz-juelich.de/icg/icg-i/idl_icglib/idl_lib_intro. html

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: Sun Oct 12 00:09:11 PDT 2025

Total time taken to generate the page: 0.72268 seconds