NOTICE: This version of the NSF Unidata web site (archive.unidata.ucar.edu) is no longer being updated.
Current content can be found at unidata.ucar.edu.
To learn about what's going on, see About the Archive Site.
NOTICE: This version of the NSF Unidata web site (archive.unidata.ucar.edu) is no longer being updated.
Current content can be found at unidata.ucar.edu.
To learn about what's going on, see About the Archive Site.
Hi All:I am working on creating a netcdf4 file with the structure we are trying to achieve. I have a file with Box10 groups and B0x2 subgroups, works great. Now I am trying to do ragged arrays. Below is the only fortran example I can find that does ragged arrays. It is only writing one,1-dimensional array, so I am unsure what I do when I am writing an N-dimensional ragged array. Do I just write each "column" in a loop - if so how does it know that it is writing to the same structure? Or if I can do multiple "columns" at once, how do I set that up. More realistic examples would be very helpful, including an example that writes an N-dimensional ragged array, and an example that adds one more column to an existing ragged array. Even just the fortran equivalent of the program tst_vl.c would be very helpful
It is clear how to do it in C (I set up an array of structures of the given form), but not in either version of FORTRAN.
Thanks, -Roy C This is part of the netCDF package.C Copyright 2008 University Corporation for Atmospheric Research/ Unidata.
C See COPYRIGHT file for conditions of use.
C This program tests netCDF-4 variable functions from fortran, even
C more, even more.
C $Id: ftst_vars4.F,v 1.9 2008/03/10 17:24:35 ed Exp $
program ftst_vars4
implicit none
include 'netcdf.inc'
C This is the name of the data file we will create.
character*(*) FILE_NAME
parameter (FILE_NAME='ftst_vars4.nc')
C NetCDF IDs.
integer ncid, vlen_typeid
integer max_types
parameter (max_types = 1)
C Need these to read type information.
integer num_types, typeids(max_types)
integer base_type, base_size, num_members, member_value
character*80 type_name, member_name
integer type_size, nfields, class
C Information for the vlen type we will define.
character*(*) vlen_type_name
parameter (vlen_type_name = 'vlen_type')
C Some data about and for the vlen.
integer vlen_len, vlen_len_in
parameter (vlen_len = 5)
integer data1(vlen_len), data1_in(vlen_len)
C These must be big enough to hold the stuct nc_vlen in netcdf.h.
integer vlen(10), vlen_in(10)
C Loop indexes, and error handling.
integer x, retval, index(1)
print *,'*** Testing VLEN types.'
do x = 1, vlen_len
data1(x) = x
end do
C Create the netCDF file.
retval = nf_create(FILE_NAME, NF_NETCDF4, ncid)
if (retval .ne. nf_noerr) call handle_err(retval)
C Create the vlen type.
retval = nf_def_vlen(ncid, vlen_type_name, nf_int, vlen_typeid)
if (retval .ne. nf_noerr) call handle_err(retval)
C Set up the vlen with this helper function, since F77 can't deal
C with pointers.
retval = nf_put_vlen_element(ncid, vlen_typeid, vlen,
& vlen_len, data1)
if (retval .ne. nf_noerr) call handle_err(retval)
C Write the vlen attribute.
retval = nf_put_att(ncid, NF_GLOBAL, 'att1', vlen_typeid, 1,
vlen)
if (retval .ne. nf_noerr) call handle_err(retval)
C Close the file.
retval = nf_close(ncid)
if (retval .ne. nf_noerr) call handle_err(retval)
**********************
"The contents of this message do not reflect any position of the U.S.
Government or NOAA."
********************** Roy Mendelssohn Supervisory Operations Research Analyst NOAA/NMFS Environmental Research Division Southwest Fisheries Science Center 1352 Lighthouse Avenue Pacific Grove, CA 93950-2097 e-mail: Roy.Mendelssohn@xxxxxxxx (Note new e-mail address) voice: (831)-648-9029 fax: (831)-648-8440 www: http://www.pfeg.noaa.gov/ "Old age and treachery will overcome youth and skill."
netcdfgroup archives: