pnf_test.F

John Tannahill tannahill1 at llnl.gov
Tue Sep 23 15:44:41 CDT 2003


Reiner,

I developed the original version of this test code, but it has been
edited significantly since then.  I took a look at pnf_test.F in
v0.9.0, and yes, it looks like your correction below needs to be
made.

I also noticed another place that I think needs to be fixed =>

        ierr = Nfmpi_Def_Dim
      &    (ncid, "level",     totsiz_3d(1), lev_id)
         ierr = Nfmpi_Def_Dim
      &    (ncid, "latitude",  totsiz_3d(2), lat_id)
         ierr = Nfmpi_Def_Dim
      &    (ncid, "longitude", totsiz_3d(3), lon_id)

I believe this should probably be =>

        ierr = Nfmpi_Def_Dim
      &    (ncid, "longitude", totsiz_3d(1), lon_id)
         ierr = Nfmpi_Def_Dim
      &    (ncid, "latitude",  totsiz_3d(2), lat_id)
        ierr = Nfmpi_Def_Dim
      &    (ncid, "level",     totsiz_3d(3), lev_id)

Also, there are a number of places where "d" exponent desciptors are
being used with numbers that don't need them.  These probably should
be changed to "e" for clarity, although I doubt that it affects much.

Regards,
John

Reiner Vogelsang wrote:
> Dear colleagues,
> 
> in pnf_test.F I found that filsiz was calculated via
> 
>       filsiz = (TOTSIZ_3D(1) * TOTSIZ_3D(2) * TOTSIZ_3D(3)) *
>      &         1.0d-6 *8.0d0
> 
> However, the code is writing/reading using NF_REAL, means 4-byte
> entities.
> The array tt is also specified as real*4.
> 
> So, I assume that 4 instead of 8 is correct.
> 
> Last but not least, are there any known restrictions to prevent me from
> running the code
> with 1024x1024x1024?
> 
> Best regards
>     Reiner
> 
> 
> 
> --
> --------------------------------------------------------------------------------
>                                              _
>                                            )/___                        _---_
>                                          =_/(___)_-__                  (     )
>                                         / /\\|/O[]/  \c             O   (   )
> Reiner Vogelsang                        \__/ ----'\__/       ..o o O .o  -_-
> Senior System Engineer
> 
> Silicon Graphics GmbH                   Home Office
> Am Hochacker 3
> D-85630 Grasbrunn                       52428 Juelich
> Germany
> 
> Phone   +49-89-46108-0                  +49-2461-939265
> Fax     +49-89-46108-222                +49-2461-939266
> Mobile  +49-171-3583208
> email    reiner at sgi.com
> 
> 
> 
> 


-- 
============================
John R. Tannahill
Lawrence Livermore Nat. Lab.
P.O. Box 808, M/S L-103
Livermore, CA  94551
925-423-3514
Fax:  925-423-4908
============================




More information about the parallel-netcdf mailing list