[交流] 【求助】求助 程序我问题 解决送80金币 已有3人参与

Originally posted by bohdshen at 2010-05-10 14:53:23:
您好
这个PORI是初始孔隙度
最后POR应该在0~1之间
我做的这段深度内孔隙度POR实际值已经知道在0.09到0.25之间
我开是没用wasat直接写 问题也一样是做了1米就不做了
无加深处是FoRWARD测井软件的固定部分
它 ...

Originally posted by snoopyzhao at 2010-05-10 19:59:23:
你把 WASAT 需要的参数给出来，你现在给出的程序没有办法运行啊……或者你给一个可以运行的，可以反映你的问题的程序……

      program main
      open(10,file='zzz',status='old')
      do 5 i=1,509
        read(10,*) ac, den
        pori = 0.029
        AS=92414*DEN/(AC**2)
        CALL WASAT(AS,PORI,POR, ind)
        if(ind .ne. 1) write(*,*) por
5     continue
      END
      SUBROUTINE WASAT(AS,PORI,POR, ind)
      ind = 0
9     AA=(1.0-PORI)**(3.0/(1.0-PORI))
      AM=10.17*((1.0-AA)**2.0)
      AN=AS-47.317*AA
      AT=(AM/AN)-PORI
      IF(ABS(AT) .LE. 0.001) THEN
        POR=PORI*100
      ELSE
        pori=pori+0.000001
        if (pori .ge. .3) then
          write(*,*) 'i can not get the correct value!'
          ind = 1
          return
        endif
        GOTO 9
      ENDIF                                                        
      RETURN
      END