爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 2430|回复: 2

[求助] CVF中可以运行的程序在IVF中如何修改?

[复制链接]
发表于 2016-3-8 11:54:18 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
程序是关于读写FY-3C卫星资料上全球经纬度以及亮度温度。
我用的是IVF,在调试程序的过程中会触发一个断点,无法读取程序。
程序给我朋友帮忙在cvf上可以运行读写出来结果。
所以上来问问大神们,哪里需要改进。
  1.      program main
  2.      implicit none
  3. !-----------discription:---------------------------------------------------------------
  4. !----Observations-----------------------------
  5. !--------------------------------------------------------------------------------------

  6. !----1.0: defile local variables-------------------------------------------------------
  7.        integer           :: i,j,k,line,ilin,fov,siz,ch
  8.        real, allocatable ::  lat(:,:),lon(:,:)
  9.        integer*2, allocatable :: obs_T(:,:,:)
  10.        real, allocatable ::  obs(:,:,:)
  11.        character(len=100):: filelat,filelon,fileobs

  12. !--------------------------------------------------------------------------------------
  13.        OPEN(11,file='d:/ziliao/lat.list')
  14.        OPEN(12,file='d:/ziliao/lon.list')
  15.        OPEN(14,file='d:/ziliao/obs.list')
  16.       
  17.        DO WHILE(.not. eof(11))
  18.           READ(11,"(I6,1X,a)")siz,filelat
  19.           READ(12,"(a)")filelon
  20.           READ(14,"(a)")fileobs
  21.           line=siz/(90*4)
  22. print*,line
  23. !-----2: read out the MWTS data--------------------------------------------------------
  24.           allocate(lat(90,line))
  25.           allocate(lon(90,line))
  26.           allocate(obs_T(90,line,13))
  27.           allocate(obs(90,line,13))
  28.          

  29.          open(21,file=filelat,form='binary')
  30.          open(22,file=filelon,form='binary')
  31.          open(24,file=fileobs,form='binary')

  32.          read(21)((lat(fov,ilin),fov=1,90),ilin=1,line)   
  33.          read(22)((lon(fov,ilin),fov=1,90),ilin=1,line)      
  34.          !----1 land; 2 fresh water; 3 sea water; 5 coast
  35.          do ilin=1,line
  36.             do fov=1,90
  37.                do ch=1,13
  38.                   read(24)obs_T(fov,ilin,ch)
  39.                   if(obs_T(fov,ilin,ch)>0) then
  40.                      obs(fov,ilin,ch)=0.01*obs_T(fov,ilin,ch)
  41.                   else  
  42.                      obs(fov,ilin,ch)=0.01*(obs_T(fov,ilin,ch)+65536)
  43.                   endif
  44.                enddo
  45.             enddo
  46.          enddo
  47. !----------------------------------------------------------------------------
  48.             
  49.          close(21)
  50.          close(22)
  51.          close(24)

  52.          deallocate(lat)
  53.          deallocate(lon)
  54.          deallocate(obs_T)
  55.          deallocate(obs)
  56. !--------------------------------------------------------------------------------------
  57.       ENDDO

  58.      close(11)
  59.      close(12)
  60.      close(14)
  61. END
复制代码


会触发断点地方为:
         read(21)((lat(fov,ilin),fov=1,90),ilin=1,line)   
         read(22)((lon(fov,ilin),fov=1,90),ilin=1,line)   
C:\Users\宁\Desktop\QQ图片20160301112822.png

资料格式:
C:\Users\宁\Pictures\QQ截图20160308114829.png

C:\Users\宁\Pictures\QQ截图20160308115347.png
C:\Users\宁\Pictures\QQ截图20160308115308.png

p.s.我最近用ivf比cvf多,然而学校里系统学习的也是cvf的用法,编写程序的时候也是按照cvf的思路,故来求教大神们。
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2016-3-8 11:56:50 | 显示全部楼层
QQ截图20160308114829.png QQ截图20160308115308.png QQ截图20160308115347.png
嘎哦……图没发出来
密码修改失败请联系微信:mofangbao
发表于 2016-3-8 12:39:42 | 显示全部楼层
最好把用到的文件也放上来,能实际运行试试
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表