[netcdfgroup] Reading in netcdf files using fortran 90

• To: netcdfgroup@xxxxxxxxxxxxxxxx
• Subject: [netcdfgroup] Reading in netcdf files using fortran 90
• From: SANJIV R <sanjiv@xxxxxxxx>
• Date: Wed, 6 May 2020 13:06:01 -0500

Hi,
   I am trying to get the dimension sizes from the gridfile for an ocean
model (ROMS). I am pasting below a short script that is trying to retrieve
the dimension id, dimension name and the dimension size for the vertical
coordinate. The script is returning junk values for all three variables. It
always prints out zero for the dimension id, a blank string for the
variable name and random numbers (e.g., -1067735120) for the dimension
size. Also, each time I run the code I get a different junk number for the
dimension size.

I am using netcdf-fortran version 4.4.4 on an intel machine. The libraries
have been compiled on my campus cluster.

I am compiling using:
ifort prog.f90 -I$NETCDF_HOME/include -L$NETCDF_HOME/lib -lnetcdff  -o
a.out,

where $NETCDF_HOME

Here is the script:
********************************************************************
program main

  use netcdf
  implicit none

  integer(kind=4) nz
  integer(kind=4) ncid,nzid,status

  character(len=200) gridfile, varname
!--------------------------------
  gridfile='grid.nc'

! Get the file id for the netcdf file being read
  status=  nf90_open(gridfile,nf90_nowrite,ncid)

! 's_rho' is the vertical coordinate
! Return the dimension id for s_rho in nzid
  status = nf90_inq_dimid(ncid, "s_rho", nzid)
! Using the dimension id to return the variable name (varname) no. of
vertical levels (nz)
  status = nf90_inquire_dimension(ncid, nzid, varname, nz)

! All three commands below are printing out junk!
  write(6,*) nzid
  write(6,*)' variable: ',varname
  write(6,*)'No. of vertical levels (should be 50)= ', nz

  status=nf90_close(ncid)

end program main
*************************************************************************

Does anybody know where am I making a mistake?

Thanks for your help.

Sanjiv,
Texas A & M University

• Follow-Ups:
  • Re: [netcdfgroup] [EXTERNAL] Reading in netcdf files using fortran 90
    • From: Warner, John C
  • Re: [netcdfgroup] Reading in netcdf files using fortran 90
    • From: Nick Papior