[ÇóÖú] ʦÐÖ¸øÁËÒ»¶Î´¦ÀíÆøÏóÊý¾ÝµÄ³ÌÐò ¿´²»Ì«¶® ¹òÇó´óÉñÖ¸µ¼½â»ó

ÕâÊÇ×Ó³ÌÐò²¿·Ö 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

»Ø¸´´ËÂ¥