[绘图界面库] 读取nc文件进行计算

program readsilt
      implicit none
      external read_silt_nc,handle_err


      !integer,parameter::lon_2=43200,lat_2=16800,depth_2=4
      real SILT_2
      character*299 silt_file
      integer status
      real E
      call read_silt_nc(SILT_2,silt_file)
      call handle_err(status)
      E=1.01*SILT_2**(3.0/5.0)
      write(*,*) E
      end program

      !number 1
      subroutine read_silt_nc(SILT_2,silt_file)
      implicit none

      integer ierr,ncid,varid,len_file
      integer,parameter::lon=43200,lat=16800,depth=4
      real lon_2(lon)
      real lat_2(lat)
      real depth_2(depth)
      real SILT_2(depth,lat,lon)
      !real E
      character*299 silt_file
      include 'netcdf.inc'

      silt_file="/home/iga_qtong/fortran_xcao/gae_c/ap-42/SILT1.nc"
      len_file=len_trim(silt_file)
      ierr=nf_open(trim(silt_file),nf_write,ncid)

      ierr=nf_inq_varid(ncid,'lon',varid)
      ierr=nf_get_var_real(ncid,varid,lon_2)

      ierr=nf_inq_varid(ncid,'lat',varid)
      ierr=nf_get_var_real(ncid,varid,lat_2)


      ierr=nf_inq_varid(ncid,'depth',varid)
      ierr=nf_get_var_real(ncid,varid,depth_2)

      ierr=nf_inq_varid(ncid,'SILT',varid)
      ierr=nf_get_var_real(ncid,varid,SILT_2)
      write(*,*) 'ierr'
      write(*,*) SILT_2(1,31.56247,-179.9958)


      end subroutine read_silt_nc

      subroutine handle_err(status)
      integer status
      if (status .NE. nf_noerr) then
           write(*,*) nf_strerror(status)
           stop 'stopped'
      end if
      end subroutine handle_err