C Create the netCDF file. retval = nf_create(file_name, NF_NETCDF4, ncid) if (retval .ne. nf_noerr) call handle_err(retval) C Create a group and a subgroup. retval = nf_def_grp(ncid, group_name, grpid) if (retval .ne. nf_noerr) call handle_err(retval) retval = nf_def_grp(grpid, sub_group_name, sub_grpid) if (retval .ne. nf_noerr) call handle_err(retval) C Create a two dims and two vars. retval = nf_def_dim(sub_grpid, dim1_name, 0, dimids(1)) if (retval .ne. nf_noerr) call handle_err(retval) retval = nf_def_dim(sub_grpid, dim2_name, 0, dimids(2)) if (retval .ne. nf_noerr) call handle_err(retval) retval = nf_def_var(sub_grpid, var1_name, NF_UINT64, 2, dimids, & varids(1)) if (retval .ne. nf_noerr) call handle_err(retval) retval = nf_def_var(sub_grpid, var2_name, NF_UINT64, 2, dimids, & varids(2)) 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)