error:conformable array is required ,请大侠帮忙

小芳林林妈

[Fortran] 纯文本查看 复制代码

查看代码

打印?

01	subroutine calcvis3(q,clw,rnw,ice,snow,t,prs,ista_num,vis3)

02	implicit none

03	integer :: ista_num

04	integer :: k

05	real,dimension(ista_num) :: t,q,clw,ice,rnw,snow,prs,vis3

06	real,dimension(ista_num) :: vovermd,celkel,tice,coeflc,coeflp,coeffc,coeffp,exponlc,                &

07	exponlp,exponfp,const1,rhoice,rhowat,tv,rhoair,beta,qprc,                &

08	qcld,qrain,qsnow,qclw,qclice,rgas,exponfc,concfc,concfp,                &

09	conclp,conclc

10	CHARACTER METH*1

11

12

13	----------------------------------------

14	celkel=273.15                         !K

15	rgas=287.04                          !J/K/kg

16	tice=celkel-10.

17	coeflc=144.7

18	coeflp=2.24

19	coeffc=327.8

20	coeffp=10.36

21	exponlc=0.8800

22	exponlp=0.7500

23	exponfc=1.0000

24	exponfp=0.7776

25	const1=-LOG(.02)

26	rhoice=917.

27	rhowat=1000.

28	DO k=1,ista_num

29

30	qprc(k)=rnw(k)+snow(k)

31	qcld(k)=clw(k)+ice(k)

32	qrain(k)=rnw(k)

33	qsnow(k)=snow(k)

34	qclw(k)=clw(k)

35	qclice(k)=ice(k)

36	rhoair(k)=prs(k)/(287.04*t(k)*(1.+0.608*q(k)))                        !air density 空气通量密度

37	vovermd(k)=(1.+q(k)/1000.)/rhoair(k)+(qclw(k)+qrain(k))&

38	/rhowat+(qclice(k)+qsnow(k))/rhoice

39	conclc(k)=qclw(k)/vovermd(k)*1000.

40	conclp(k)=qrain(k)/vovermd(k)*1000.

41	concfc(k)=qclice(k)/vovermd(k)*1000.

42	concfp(k)=qsnow(k)/vovermd(k)*1000.

43

44	beta=coeffc*concfc(k)**exponfc+coeffp*concfp(k)**exponfp                        &

45	+coeflc*conclc(k)**exponlc+coeflp*conclp(k)**exponlp                        &

46

+1.E-10

47	vis3(k)=min(25.,const1/beta)

48

ENDDO

49

50

RETURN

51

END

fcode

这个应该解决了，是数组不对应造成的。

const1,rhoice,rhowat,beta 等不应该定义为数组。