24小时热门版块排行榜

返回列表

无知小白菜

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1.5
帖子: 2
在线: 2小时
虫号: 4525499
注册: 2016-03-21
专业: 水文、水资源

[求助] 师兄给了一段处理气象数据的程序看不太懂跪求大神指导解惑

这是子程序部分

subroutine readcmfd(ix,jx,iyear,imon,nmstep,dir,cmfd)  !fang

   implicit none

   character(len=256),intent(in) :: dir

   integer, intent(in) :: ix,jx
   integer, intent(in) :: iyear,imon
   integer, intent(in) :: nmstep
   real, intent(out),dimension(ix,jx,248,7) :: cmfd

   integer :: i,j,itm,ivar
   integer, parameter :: nvar = 7
   real, allocatable, dimension(:,:,

:: var

   include 'netcdf.inc'

   CHARACTER(len=256) :: ncfile
   CHARACTER(len=12)  :: varname(nvar)
   INTEGER STATUS
   INTEGER NCIDin
   INTEGER varID
   INTEGER START(3), COUNT(3)

   data (varname(ivar), ivar=1,nvar)          &
                  /'temp','shum','wind','pres',  &
                  'prec','srad','lrad'/

   allocate( var (ix,jx,nmstep) )

   START(1) = 1
   START(2) = 1
   START(3) = 1

   COUNT(1) = IX
   COUNT(2) = JX
   COUNT(3) = NMSTEP

   do ivar=1,7
   do itm=1,NMSTEP
   do J=1,JX
   do I=1,IX
         CMFD(i,j,itm,ivar)    = -999.
   end do
   end do
   end do
   end do

   write(*,*) '-------------------------------------------'

   do ivar = 1,  nvar
      write(*,*) varname(ivar)
      !write(*,'(I4.4,I2.2)') iyear,imon
      write(ncfile,100) iyear,varname(ivar),iyear,imon
      ncfile = TRIM(DIR)//ncfile

      write(*,*) 'READFORC: opening ',trim(ncfile)
      STATUS=NF_OPEN(ncfile,NF_NOWRITE,NCIDin)

      IF (STATUS .NE. NF_NOERR) then
         print *,'error when reading forcing data'
         print *,NF_STRERROR(status)
         stop
      ENDIF

      STATUS=NF_INQ_VARID (NCIDin, varname(ivar), varID)
      STATUS=NF_GET_VARA_real(NCIDin,varID,START,COUNT,var)
      !print *,status

      ! if (ivar .eq. 7) then
      ! !print *,var(223,209,1)
      ! !print *,var(223,209,1) * 1.e-6  + 0.025
      ! print *,var(223,209,1) * 2.5e-1 + 685.

      ! do j=1,400
      !    write(1215,"(700F10.3)"

(var(i,jx-j+1,1) * 2.5e-1 +685.,i=1,700)
      ! enddo

      !endif

      do I=1,IX
      do J=1,JX
      do itm=1,nmstep
         cmfd(i,jx-j+1,itm,ivar) = var(i,j,itm)
      end do
      end do
      end do

      STATUS=NF_CLOSE(NCIDin)
   end do

100  format('/',i4.4,'/',a4,                            &
         '_ITPCAS-CMFD_V0106_B-01_03hr_010deg_',i4.4,i2.2,'.nc')

   write(*,*) '-------------------------------------------'

  10 format(1x,5f15.5)

   return
   end subroutine readcmfd

回复此楼