program wwcdf
CCC---------------------------------------------------------------------
CCC
CC declarations of variables
CC =========================
C
      IMPLICIT NONE
CC
CC    jpi = longitude
CC    jpj = latitude
CC    jpk = depth
CC    jit = number of time step
CC
      INTEGER jpi,jpj,jpk,jit
      PARAMETER(jpi=20,jpj=20,jpk=1,jit=1)
      INTEGER ji,jj,jk,it
CC
      REAL wflx(jpi,jpj,jpk,jit)
CC
      REAL lat(jpj), lon(jpi), depth(jpk), time(jit)
      REAL fix
CC
      INCLUDE 'netcdf.inc'
c
      INTEGER ncid, iret, varid
      INTEGER varlat, varlon, vardep, vartim
      INTEGER nvdims
      PARAMETER(nvdims = 4)
      INTEGER vardims(nvdims),start(nvdims),count(nvdims),stride(nvdims)
      INTEGER status
c
      DATA start/1,1,1,1/
      DATA count/jpi,jpj,jpk,jit/
      DATA stride/1,1,1,1/
      DATA fix/1.e35/
CC
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c  open the netCDF file
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c
C NOTE: the variable status should always return a 0 value. When this
C is not 0, then the operation failed (check the name of the file and
C variables, watch for capital letters)
C
       status = nf_open("fileou.nc",nf_noclobber,ncid)
       write(*,*) "File open", status
C      if (status.ne.0) write(*,*) nf_strerror(status)
c
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c  inquire attributes for the grid and variables
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c
       status = nf_inq_varid(ncid,'LONGITUDE',varlon)
       if (status.ne.0) write(*,*) "Problem in varlon",status
C
       status = nf_inq_varid(ncid,'LATITUDE',varlat)
       if (status.ne.0) write(*,*) "Problem in varlat",status
C
       status = nf_inq_varid(ncid,'DEPTH',vardep)
       if (status.ne.0) write(*,*) "Problem in vardep",status
C
       status = nf_inq_varid(ncid,'TIME',vartim)
       if (status.ne.0) write(*,*) "Problem in vartim",status
C
       status = nf_inq_varid(ncid,'co2',varid)
       if (status.ne.0) write(*,*) "Problem in varid",status
C
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c  read the data
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c
       status = nf_get_var_real(ncid,varlat,lat)
       write(*,*) "this is latitude",lat
C
       status = nf_get_var_real(ncid,varlon,lon)
       write(*,*) "this is longitude",lon
C
       status = nf_get_var_real(ncid,vardep,depth)
       write(*,*) "this is depth",depth
C
       status = nf_get_var_real(ncid,vartim,time)
       write(*,*) "this is time",time                    
C
       status = nf_get_var_real(ncid,varid,wflx)              
       write(*,*) "this is co2_flx",wflx                    
c   
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c  close the netCDF file
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c
      status = nf_close(ncid)
      write(*,*) 'all data read'
C
      END