EXAMPLE FORTRAN PROGRAM PRICLAN
DIMENSION ABSVAR(990),SVVAR(990),
*W(990) ,V(990) , (80) ,AKORRX (990) ,
* VEORRT (990) ,X (990) ,AEN(990) ,
*UZUN(990) , NÉ (990) ,NN(990)
* GRANIA (4000,40) ,SÜMGR(1000)
DO 44 KK=1,1
READ(5,20) TT
20 FORMAT(8O, A^)
PRINT 20,TT
READ(5,102)N
102 FORMAT(T 3) is
74 FORMAT(AR ,60(*-^))
DO 3875 I=1,N
READ 200,NN(I),NE(I),X(I),AEN(I),
*UZUN(I), (GRANIÄ(I,J),J=1,10)
BRINT 200,NN(I),NK(I),X(L),AEN(I),
*UZUN(I), SRANIACI,J) 424, 10))
3875 CONTINUE
200 FORMAT(13,12,F5.1,F5.1,10F5.1)
DO 623 I=1,N
SUMGR(I)=0.0
DO 574 J=1,10
SUMGR (I)-SÜMGR(I)4*GRANIA(I,J)
574 CONTINUE
623 CONTINUE
CALL FORMA(NN,NK,X,AEN,UZUN,
*SUMGR,N)
DO 66 IA=1,3
IF(IN.EQ.1960 TO 985
IF(IN:EQ:29G0 TO 986
IF(IN,EQ.39GO TO 987
986 DO 456 L2=1,N
X(L2)=AEN(L2)
456 CONTINUE
GO TO 985
987 DO 867 L2-4,N
X(L2)=UZUN(L2)
867 CONTINUE
985 SYM-O,
XMIN=X (1)
XMAX-XMIN
DO 77 I=1,N
Ip GE CD) Lb anre qr)
IF(X(1).G8T.XMAX OXMAX-X (I)
7? CONTINUE A
PRINT 88,XMIN,X à |
88 FORMAT(4H ,*"XMINz^,P40.2,*XMAXz^,
*F10.2)
DO 4 1=1,N
SYM=SYM+X (I)
4 CONTINUE
AN=N
SR=SYM/AN
WRITE (6,400)STM,SR | ;
400 FORMAT SYM *,P41.2, "HIDDIE^,
* *^ARTTHMETICAL^, "SRa*,F41.2)
AVAR=O °
AKV. AR=0,
AINFOR-O.
DO 23 I=1,N
DOLJÁ-X(1)/SYM
AINFOR=AINFOR+3. 3222* ALOG10(DOLJA)*
*DOLJA
ABSVAR(I)-X(I)-SR
AVAR=AVAR+ABSVAR(T)
SVVAR(I)-ABSVAR(I)**2
AKVARZAKVAR4SVVAR(I)
23 CONTINUE
AINFOR=ABS (AINFOR)
DO TK =SQRT(AKVAR/ (AN=1,0))
AMAXTN=3. 3222* ATOG1 O(AN)
OTINF=(AMAXIN-ATNFOR) /AMAXIN
HDO LI =AMAXIN-ATNFOR
HDELTASAMAXIN/AINFOR
WRITE(6,50)HDOLI,HDELTA
//P
50 FORMAT(1H ,"HDOLIz ^,P410.2,
** HDELTA= ",R40.2)
PRINT 75,AINFOR,AMAXIN,OTINF _ |
75 FORMAT(1H , “GENERAL ENTROPY HO=",
*F10.4/°MAXIMAL ENTROPY HMAX=",
*F10.4/ RELATIVE ENTROPY HOT=",
*F10.4)
CV=DOTK/SR
WRITE(6,24)AVAR,AKVAR, DOTK,CV
24 FORMAT(AH ,"AVAR-^ ,F40.2
* AKVARz" ,E40.2, ^DO Ts" ,FA0.2,
**'cyz* e.
TRITE(6,74
NM=
CALL KABKAZ(NM,X)
WRITE(6,74)
66 CONTINUE
WRITE(6,74)
CALL FORMA(X,AEN,UZUN,SUNGR,NM)
44. CONTINUE
STOP
END
SUBROUTINE FORMA(NN,NK,PLOSAD,
*SIRINA,DLINA, GRANIS ,N)
DIMENSION DLINA(1000),SIRINA(1000),
* GRANTS (1000) FIOSAD (1500) , E0077),
+5 (1000) ,NN (1000) , NE(1000
DO 3.1=1,N
GRANIS(19=GRANIS(1)*2.0
SIRINA(I)=SIRINA(I)* 2.0
DLINA(I)=DLINA(I)* 2,0
PLOSAD(1)=PLOSAD(1)* 4.0
3 CONTINUE
DO 4 I=1,N
S(I)-DLINA(I)/SIRINA(I)
SPK(1)=GRANIS(1)/3. 54* SORT(PLOSAD(I))
XS(1)24.24* PIOSAD(I)/DLINA(I)**2
4 CONTINUE
WRITE(6,5)(S(I),I=1 ,N)/
* (SPK(1),1=1,N)/(XS(1),1=1,N)
5 FORMAT(12840, 4/** 42840. 4/4240, 4)
SUBROUTINE KABKAZ(N,SQ)
*DIMENSION SQ(990) ,WWAC80) „QSQ(990)
*S (990) , AM(990)
AN=N
DO 3 I=1,N
3 WN=WN+SQ(I)
DO 33 1=1,N
S(1)=(SQ(I)/WN)* 100.0
AM(I)zS(1)/(100.0/AN)
33 CONTINUE
SUMzO,
N1=N-1
DO 7 I=1,M
I1=I41
DO 4 Jw I1,N
QSQ(I)zAM(I)*AM(J)
SUM=SUM+QSQ(I)
4 CONTINUE
7 CONTINUE
ODNOR=SUM/ (AN* (AN-1.0)/2.0)
RITE(6,5)WN,SUM,ODNOR | 2
5 FORMAT(SX,"WNz^,F40.5,0X, "SUMz^, |,
*F1105,4X, COEFFICIENT LANDSCAPIC “,
* ’UNHOMOGÉNIC -ODNOR= “,P10.4)
RETURN
END
EXAMPLE SYSTEMIC PROCEDUR
RICLAND JOB MSGLEVEL=(1,1),MSGCIASS=P
//GEOGRAF EXEC FORTGCLG,PARM.FORT=SOURCE
//FORT.SYSPRINT DD SYSOUT-P
//FORT.SYSIN DD UNIT=SYSDA,VOL=SER=GEOLAN,
/
DSN=PROGLAN(PRICLAND)
//G0.SYSIN DD UNIT-SYSDA,VOL-SER-GEOBAS,
//
175
DSN=LAN.DATA( LAND BASE)