Visual fortran , IMSL code

Visual fortran , IMSL code

2 posts / 0 new
Last post
For more complete information about compiler optimizations, see our Optimization Notice.

 Hello, can someone help me?!

I want to apply IMSL library in PGI visual fortran 2013, but my codes were been written in f90,

and after staring, the "This project is out of date" comes up, and also ask "Would you like to built it?".

and then after continuing again this error message is shown,"Visual studio cannot debug because the debug target 'C:\Users\wave\documents\visual studio 2013\projects\PVFProject5\PVFProject5\Win32\Debug\PVFProject5.exe'is missing". 

 

 USE MSIMSLMD
    EXTERNAL FCN
    Double precision v1,v2,v3,k12,k23,k13,MV1,MV2,MV3
    Double precision x,xguess,fnorm,errrel,v
    Dimension x(1),xguess(1),fnorm,errrel,v
    Common/const/k12,k23,k13,MV1,MV2,MV3,v3
    Data k12,k23,k13,MV1,MV2,MV3,v3
    Data k12,k23,k13,MV1,MV2,MV3/0.9,0.5,1.11,18.0,127.8,30532.0/
    ERRREL=1E-10
    ITMAX=300
    N=1
c
    open(unit=1,file='spin1.txt')
2 format(20x,16HVolume Fractions,//,12x,2Hv1,25x,2Hv2,22x,2Hv3)
    write(1,2)
    V(1)=0.4

C   Initial guess
10         write(*,*)'enter v3'
          read(*,*)v3
          xguess(1)=v(1)
          CALL DNEQNF(FCN,ERRREL,N,ITMAX,XGUESS,X,FNORM)
          v(1)=x(1)
          v2=1-v(1)-v3
          write(1,*)v(1),v2,V3
          GOTO 10
          clse(unit=1)
          end
       C  IMSL required subroutine FCN
             Subroutine fcn(x,f,n)
             Double precision v1,v2,v3,k12,k23,k13,MV1,MV2,MV3
             Double precision f,x,G22,G23,G33
             Dimension f(n),x(n)
             common/const/k12,k23,k13,MV1,MV2,MV3,v3
             v1=x(1)
             v2=1-v1-v3
             G22=1/V1+MV1/(MV2*V2)-2*K12
             G23=1/V1-(K12+K13)+MV1/MV2*K23
             G33=1V1+MV1/(MV3*V3)-2*K13
             f(1)=G22*G33-G23**2
             return
             end

Leave a Comment

Please sign in to add a comment. Not a member? Join today