读出来的数非常小负的230次方，为什么？
SUBROUTINE USDFLD(FIELD,STATEV,PNEWDT,DIRECT,T,CELENT,TIME,DTIME,
1 CMNAME,ORNAME,NFIELD,NSTATV,NOEL,NPT,LAYER,KSPT,KSTEP,KINC,
2 NDI,NSHR,COORD,JMAC,JMTYP,MATLAYO,LACCFLG)
C
INCLUDE 'ABA_PARAM.INC'
C
CHARACTER*80 CMNAME,ORNAME
CHARACTER*3 FLGRAY(15)
DIMENSION FIELD(NFIELD),STATEV(NSTATV),DIRECT(3,3),T(3,3),TIME(2),
* COORD(*),JMAC(*),JMTYP(*)
DIMENSION ARRAY(15),JARRAY(15)
C
C GET TEMPERATURES FROM PREVIOUS INCREMENT
CALL GETVRM('S',ARRAY,JARRAY,FLGRAY,JRCD,
$ JMAC, JMTYP, MATLAYO, LACCFLG)
STATEV(1) = ARRAY(1)
WRITE(6,*)'the S11 is:', STATEV(1)
RETURN
END
C=======================================================================
SUBROUTINE FRIC(LM,TAU,DDTDDG,DDTDDP,DSLIP,SED,SPD,
1 DDTDDT,PNEWDT,STATEV,DGAM,TAULM,PRESS,DPRESS,DDPDDH,
2 SLIP,KSTEP,KINC,TIME,DTIME,NOEL,CINAME,SLNAME,
3 MSNAME,NPT,NODE,NPATCH,COORDS,RCOORD,DROT,TEMP,
4 PREDEF,NFDIR,MCRD,NPRED,NSTATV,CHRLNGTH,PROPS,NPROPS)
C
INCLUDE 'ABA_PARAM.INC'
C
CHARACTER*80 CINAME,SLNAME,MSNAME
DIMENSION TAU(NFDIR),DDTDDG(NFDIR,NFDIR),DDTDDP(NFDIR),
1 DSLIP(NFDIR),DDTDDT(NFDIR,2),STATEV(*),
2 DGAM(NFDIR),TAULM(NFDIR),SLIP(NFDIR),TIME(2),
3 COORDS(MCRD),RCOORD(MCRD),DROT(2,2),TEMP(2),
4 PREDEF(2,*),PROPS(NPROPS)
C
CALL SUBROUTINE USDFLD(FIELD,STATEV,PNEWDT,DIRECT,T,CELENT,TIME,DTIME,
1 CMNAME,ORNAME,NFIELD,NSTATV,NOEL,NPT,LAYER,KSPT,KSTEP,KINC,
2 NDI,NSHR,COORD,JMAC,JMTYP,MATLAYO,LACCFLG)
RETURN
END |
发表于 2006-3-7 18:11
| 
