[ncl-talk] Error creating binary file

[Rick Brownrigg]

Hi Andrew,

That binary file looks correct to me, as far as I can tell. The begins with
5929560 values in the 100,000-ish range, followed by values in the -20/+20
ish range.

Is there possibly a mismatch between the array sizes on the NCL side, and
the array sizes on the Fortran side (which appear to be determined via a
parameter file?)?

Rick

On Fri, Feb 23, 2018 at 2:46 PM, Andrew Kren - NOAA Affiliate <
andrew.kren at noaa.gov> wrote:

> It says this in the fortran program. However, I am confused if ncl prints
> out bad wind values, so was thinking it was something on my end that I did
> in the ncl script since slp was fine. I attached the fortran code and the
> binary file and grads ctl file I made as google drive links.​
>  analysis_1.dat
> <https://drive.google.com/a/noaa.gov/file/d/1Y3JXqmxZ7ViBeGkGncFkkIS9wncmAVg4/view?usp=drive_web>
> ​​
>  create_binary.ncl
> <https://drive.google.com/a/noaa.gov/file/d/1fyDQeLSom9Za7r73xpeWvv6nMrVO0o7P/view?usp=drive_web>
> ​​
>  track.F
> <https://drive.google.com/a/noaa.gov/file/d/1Pn9ON5epsr6u_24jxi8eDXBjV8DgizqC/view?usp=drive_web>
> ​​
>  grads_ecm.ctl
> <https://drive.google.com/a/noaa.gov/file/d/1M9KUlm2HZwlaD_gw29Z7hn2qahYJ2Jhv/view?usp=drive_web>
> ​
>
> open(luin,file=filename,form='unformatted',status='unknown',
>
>      &       access='direct',recl=idim*jdim*4)
>
>
> On Fri, Feb 23, 2018 at 2:27 PM, Rick Brownrigg <brownrig at ucar.edu> wrote:
>
>> Hi Andrew,
>>
>> How are you opening the file from the fortran side? The docs for
>> fbindirwrite state that the fortran equivalent is a file opened with
>> "'form=unformatted, access=direct, recl= ...'.".
>>
>> If that doesn't help, I'd be happy to take a look at the file if you are
>> able to send it to me.
>>
>> Rick
>>
>> On Fri, Feb 23, 2018 at 2:03 PM, Andrew Kren - NOAA Affiliate <
>> andrew.kren at noaa.gov> wrote:
>>
>>> Dear ncl-talk,
>>>
>>> I wrote a program to convert grib data to binary format, as needed for a
>>> fortran program to track cyclones. I write slp, u at 500, and v at 500 to
>>> binary using fbindirwrite. I ran the fortran program, which gave erroneous
>>> wind values. So, I checked my binary program again and printed out the wind
>>> for a few points and they were erroneous. Is this because the wind is
>>> negative at times in direction? Here are some wind values I get when
>>> putting into binary:
>>>
>>> (0) -3.98
>>>
>>> (1) -3.98
>>>
>>> (2) -3.98
>>>
>>> (3) -3.99
>>>
>>> (4) -3.95
>>>
>>> (5) -3.96
>>>
>>> (6) -3.96
>>>
>>> (7) -3.96
>>>
>>> (8) -3.96
>>>
>>> (9) -3.96
>>>
>>> (10) -3.95
>>>
>>> And here is what I get when I read the binary file back in after created:
>>>
>>> (0) 102011.4
>>>
>>> (1) 102011.4
>>>
>>> (2) 102011.4
>>>
>>> (3) 102011.4
>>>
>>> (4) 102011.4
>>>
>>> (5) 102011.4
>>>
>>> (6) 102011.4
>>>
>>> (7) 102011.4
>>>
>>> (8) 102011.4
>>>
>>> (9) 102011.4
>>>
>>> (10) 102011.4
>>>
>>> Am I missing something. I also attached my code. Any help is much
>>> appreciated!
>>>
>>> Thanks,
>>>
>>> --
>>> Andrew Kren
>>> Assistant Scientist
>>> CIMAS & NOAA/AOML
>>> 325 Broadway, Boulder, CO 80305
>>> <https://maps.google.com/?q=325+Broadway,+Boulder,+CO+80305&entry=gmail&source=g>
>>> (303) 497-5418
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>>
>
>
> --
> Andrew Kren
> Assistant Scientist
> CIMAS & NOAA/AOML
> 325 Broadway, Boulder, CO 80305
> <https://maps.google.com/?q=325+Broadway,+Boulder,+CO+80305&entry=gmail&source=g>
> (303) 497-5418
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180223/4ae580be/attachment.html>