[Fortran] 纯文本查看 复制代码
PARAMETER (KCR =15)
PARAMETER (KCI =10)
PARAMETER (KNW =20)
PARAMETER (KDT =50)
PARAMETER (KNANG =60)
PARAMETER (KNSIZE=60)
PARAMETER (KKSZ=KNSIZE+KNW*2)
PARAMETER (KKDT=KDT*KNW+KNW)
PARAMETER (KNA1 =50)
PARAMETER (KNFI =50)
PARAMETER (KNLN =4)
PARAMETER (KNMEAS=KNANG*KNW+2*KNW,KNA1U=KNA1)
C
PARAMETER (PI=3.141592653589793,RAD=PI/180.0)
CHARACTER ERC*(*)
DIMENSION CONCA(KNLN),CONCM(KNLN)
DIMENSION WL(KNW),GA(KNW),SOLID(KNW),TAUO3(KNW)
& ,OPT(KNW),AUR(KNW,KDT),IAUR(KNW,KDT)
DIMENSION TM(KNW),WM(KNW)
DIMENSION VOL0(KNSIZE),VOL(KNSIZE)
DIMENSION VCR(KCR),VCI(KCI),AK(KNMEAS,KNSIZE,KCR,KCI)
& ,THETA(KNANG),PPM(KNANG),PRTR(KNMEAS),PRTR0(KNMEAS)
DIMENSION TH(KNW,KDT),FI(KNW,KDT),ISTH(KNW,KDT),ISFI(KNW,KDT)
DIMENSION AM1U(KNW,KNA1),FAI(KNW,KNFI),NA1U0(KNW),NFI0(KNW)
DIMENSION AURC(KNW,KDT),AURC0(KNW,KDT),OPTC(KNW),OPTC0(KNW)
& ,AURB(KNW,KDT),OPTB(KNW),AWGT(KNW,KDT),FWGT(KNW)
DIMENSION ERRV(KNW)
DIMENSION DLA(KKDT),DLAC(KKDT,KKSZ),DEV(KKSZ),DEV0(KKSZ),IEF(KKDT)
& ,H(KKSZ,KKSZ),DDV(KKSZ),SIG(KKSZ)
DIMENSION JCR(KNW),JCI(KNW),WCR(KNW),WCI(KNW),WCR0(KNW),WCI0(KNW)
& ,AKNEL(KNMEAS,KNSIZE)
C
DATA COEF /10./
DATA G0,VRATE /1., 0.01/
C
ERC=' '
DO IW=1,NW
WCR0(IW)=WCR(IW)
WCI0(IW)=WCI(IW)
ENDDO
CALL VLRTRN(IPLC,NLN,CONCA,CONCM,NW,WL,GA,SOLID
& ,TAUO3,TM,WM,TH0,NA,TH,FI,ISTH,ISFI,OPT,AUR,IAUR
& ,NA1U0,AM1U,NFI0,FAI,NSIZE,VOL0
& ,IPOL,NANG,THETA,AKNEL,PPM,OPTC0,AURC0,PRTR,ERC)
IF (ERC.NE.' ') GOTO 890
DO IW=1,NW
DO I=1,NA
AURB(IW,I)=AURC0(IW,I)
ENDDO
OPTB(IW)=OPTC0(IW)
ENDDO