求各位大神有没有能看懂这个程序的

jeffggg

lqouc 发表于 2013-6-2 17:32

                               
登录/注册后可看大图

你贴出来的程序最后说了会调用这个的。下载以后你看看就知道了。

第二次贴出来的是就是get_scat_daily_v4.f
我没想明白 这两个分开独立的  怎么调用 而且读取数据的路径加在哪里他也没有说

lqouc

第二个程序你根本就不用管，你只需要把第一个的路径写好就好了，第二个这个是个子程序，在主程序里面第二个call的就是它，在call之前你已经把数据文件打开了。
这个子程序的调用和普通的完全一样，就是call就可以了，至于为什么分开的，之前告诉你了，在主程序的最后有一个include，这句话就是说有这么个子程序或者函数，我不写在主程序后面，但是我要调用，原因可能是太长，而且你用不着改这部分内容，所以没必要让你看到里面的代码，就分开放了。这样的话，你把这个用来读取数据的子程序放在和主程序相同的路径就好了。

jeffggg

lqouc 发表于 2013-6-2 18:51

                               
登录/注册后可看大图

第二个程序你根本就不用管，你只需要把第一个的路径写好就好了，第二个这个是个子程序，在主程序里面第二个 ...

嗯 我把两个放到一个目录下面了  我再把第一个程序发一下 里面一些地方大神再帮我看下吧 运行的时候提示错误
Program  Use_Scat_Daily_Example        
        character(120) filename

        real(4) mingmt(1440,720,2)
        real(4) windspd(1440,720,2)
        real(4) winddir(1440,720,2)
        real(4) radrain(1440,720,2)
        integer(4) scatflag(1440,720,2)

        integer lyear,jday,imonth,idaymonth,ierr,iasc,ilat,ilon
        
        DO lyear =2000,2000        !qscat verification 2000,2000        !seawinds verification  2003,2003        !full qscat data range  1999,2004 (or current year)
          DO jday=11,11         !qscat verification        11,11                !seawinds verification  115,115                !full year  1,366
                IF(lyear.eq.1999.and.jday.lt.200)                         cycle        !before qscat in orbit 如果等于1999中的某一天 cycle之后的语句不再被执行，而是直接开始新一轮的循环
                IF(lyear.ne.4*INT(lyear/4).and.jday.eq.366)  cycle        !not a leap year 如果是闰年的某一天，cycle之后的语句不再被执行，而是直接开始新一轮的循环。
               
            CALL find_month_day(lyear,jday,imonth,idaymonth) ! 调用子函数

                    ! Construct file name  建立文件名称
                  ! ADJUST FOR YOUR DIRECTORY STRUCTURE and INSTRUMENT 调整你的目录结构!
                ! This subroutine assumes you have unzipped files with no extention 这个子程序假定你解压后的文件没有扩展

                WRITE(filename,9001) lyear,imonth,lyear,imonth,idaymonth
9001                FORMAT('H:\Qscat\y',i4.4,'\m',i2.2,'\qscat_',i4.4,i2.2,i2.2,'v4')
                WRITE(*,*) filename

                                                                                                               
                CALL GET_SCAT_DAILY_V4(filename, mingmt,windspd,winddir,scatflag,radrain, IERR)
                IF (IERR.ne.0) THEN
                        WRITE(*,*) 'no scat data: ', filename
                        STOP
                ENDIF


                !!  write out scat verification file data for comparison 写出SCAT验证文件数据的比较
                DO iasc=1,1                !range !1,2(asc,dsc)    !seawinds and qscat verification !1,1
                  DO ilat=274,278           !1,720                                                   !274,278
                    DO ilon=170,175           !1,1440                                                   !170,175

                !                XLAT=0.25*ILAT-90.125
                !                XLON=0.25*ILON-0.125

                                write(*,'(2i5,f6.0,2f8.2,2x,i4,f11.2)') ilon,ilat,mingmt(ilon,ilat,iasc),windspd(ilon,ilat,iasc), winddir(ilon,ilat,iasc),scatflag(ilon,ilat,iasc), radrain(ilon,ilat,iasc)

                        ENDDO  !ilon
                  ENDDO  !ilat
                ENDDO  !iasc

          ENDDO        !jday
        ENDDO        !lyear

        END

          !! write the location of the qscat data subroutine 写QSCAT数据子程序的位置
            !! include 'YourDrive:\YourDirectory\GET_SCAT_DAILY_v4.F'  驱动器 目录
        include 'GET_SCAT_DAILY_V4.F'   !调入函数或自定义的内容




      SUBROUTINE FIND_MONTH_DAY(LYEAR,IDAYJL, IMON,IDAY)

        INTEGER(4) IDAYFX(12,0:1)
      DATA IDAYFX/1,32,60,91,121,152,182,213,244,274,305,335,1,32,61,92,122,153,183,214,245,275,306,336/           
                                                  
      ILEAP=0
        IF(LYEAR.EQ.4*INT(LYEAR/4)) ILEAP=1

      DO 10 JMON=2,12
        IF(IDAYFX(JMON,ILEAP).GT.IDAYJL) THEN
      IMON=JMON-1
        GO TO 20
        ENDIF
   10 CONTINUE
      IMON=12
   20 CONTINUE

      IDAY=1+IDAYJL-IDAYFX(IMON,ILEAP)
      RETURN
        END

jeffggg

lqouc 发表于 2013-6-2 18:51

                               
登录/注册后可看大图

第二个程序你根本就不用管，你只需要把第一个的路径写好就好了，第二个这个是个子程序，在主程序里面第二个 ...

编译连接都没有错误  但是运行后  提示我错误 错误是H：\Qscat\y2000\m01\qscat_20000111v4

jeffggg

lqouc 发表于 2013-6-2 18:51

                               
登录/注册后可看大图

第二个程序你根本就不用管，你只需要把第一个的路径写好就好了，第二个这个是个子程序，在主程序里面第二个 ...

错误提示的这个路径下面我按要求放进去了文件啊 为什么还提时错误呢？

lqouc

完整错误信息呢？这个不全啊。

jeffggg

lqouc 发表于 2013-6-2 19:43

                               
登录/注册后可看大图

完整错误信息呢？这个不全啊。

你看嘛  第二个程序我根本没敢动  第一个程序 我只是加进去了一个路径  其余都没动 。。错误提示就只有这一点  图片上有的 你的意思是哪里错误信息不全呢？

jeffggg

lqouc 发表于 2013-6-2 19:43

                               
登录/注册后可看大图

完整错误信息呢？这个不全啊。

然后我选择调试程序  就是图片提示处发生错了 FALSE

lqouc

好吧，之前没看见你的报错信息图片，那就是你的fortran里面堆栈设置太小了，去调一下，在link里面。

言深深

关于此贴，说两点：

1、你发的找资料的帖子有人回复了的，怎么可以说是没有人回复呢？？再者我可以很负责任的告诉你，基本上你求的资料，论坛都有人发过，找找别人发的帖子看看；

2、随便这么一贴程序就发问是对自己，也是对别人的不负责任。光看你的帖子内容完全不知道在看什么，看了回复才了解到，原来还有前续内容···

【以上是表达一下情绪，下面说点儿有用的吧】

7#的说法是对的，这个里面有好些个程序，相互调用的，如果你对fortran不熟悉，或者不懂的话，不建议你跑这么长的程序，子程序“ GET_SCAT_DAILY_V4(filename,mingmt,windspd,winddir,scatflag,radrain, IERR)”就是缺失的。

这个程序里面看起来像是根据日序挑选具体日期数据的，因为涉及到具体的平年闰年判断子程序“SUBROUTINE FIND_MONTH_DAY(LYEAR,IDAYJL, IMON,IDAY)”实现该功能。

其余的地方就是简单的输入输出，需要修改的地方主要在：

DO lyear =2000,2000
    DO jday=11,11  

两句