| 查看: 914 | 回复: 2 | ||
[求助]
师兄给了一段处理气象数据的程序 看不太懂 跪求大神指导解惑
|
|
这是子程序部分 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(:,:,  :: varinclude '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 |
回复此楼» 猜你喜欢
青C资助名额大幅增加！
已经有15人回复
-
风把牡丹吹跑了
已经有8人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
评审感受-评审感受-评审感受
已经有18人回复
-
这年头没有找到涵评专家，还有中面上的可能吗
已经有10人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
教学课件你会给同学吗
已经有8人回复
-
26应届毕业生考博求助
已经有3人回复
-
重磅！青年科学基金项目（C类）资助增幅预计超过50%
已经有7人回复
2楼2016-03-31 14:19:45
yituoshiniao
至尊木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 23588.9
- 散金: 605
- 红花: 1
- 沙发: 1
- 帖子: 2808
- 在线: 96.6小时
- 虫号: 2150835
- 注册: 2012-11-27
- 性别: GG
- 专业: 石油、天然气地质学
3楼2016-08-02 21:11:22
精华III:
20