FORTRAN語言

--FORTRAN77結構化程序設計講師：徐紅波Email:第九章子程序

應用程序由一個主程序和若干個子程序組成。程序中可以只有主程序而沒有子程序，但是不能沒有主程序而只有子程序。以前程序舉例都只是由主程序組成的。程序在運行時總是從主程序開始執行，所以如果沒有主程序，程序就無法運行。 子程序有函數子程序、子例行程序和數據塊子程序三種。數據塊子程序僅用來給公用區中的變量賦初值，故將它放在后面與公用區的使用同時介紹。

子程序可以被主程序調用，也可以被其它子程序調用，但是子程序不能直接或間接地自己調用自己，不允許遞歸調用。把調用子程序的主程序和子程序統稱為調用程序，把被調用的子程序統稱為被調用程序。9.1函數子程序函數實例：f(x)=x4+x3+2x-7 +1 (x>0) Y=0 (x=0) -1 (x<0)函數子程序：

FUNCTIONY(X) IF(X.GT.0.0)THEN Y=1.0 ELSEIF(X.EQ.0.0)THEN Y=0.0 ELSE Y=-1.0 ENDIF END

三種函數：內部函數、語句函數、函數子程序(外部函數)。 函數子程序是一個單獨的程序單位，它不屬于調用程序所在的程序單位，因此它稱為外部函數。語句函數則屬于其所在的程序單位，只能被本程序單位中表達式引用。內部函數是FORTRAN提供的系統標準函數，可以為任何程序單位引用。[例9.1]編寫求實型數絕對值的函數 REALFUNCTIONMYABS(X) REALX IF(X.LT.0)THEN MYABS=-1*X RETURN ELSE MYABS=X ENDIF END PROGRAMMAIN REALMYABS,X,A WRITE(*,*)‘Enteranumber:’ READ(*,*)A X=MYABS(A) WRITE(*,*)‘A=’,A,‘X=’,X END9.1.1函數子程序的定義格式1：類型說明FUNCTION函數名(虛擬參數表)格式2：FUNCTION函數名(虛擬參數,…)

類型說明函數名格式3：類型說明FUNCTION函數名()虛參可以是變量名、數組名、子程序名。9.1.2函數子程序的調用

外部函數的調用形式與內部函數的調用形式完全相同，其形式如下：

函數名(實在參數,…)或者函數名()[9.2]編寫求1+2+…+n的函數子程序，調用此函數求以下表達式的值Y=((1+2+3)+(1+2+3+4)+(1+2+3+4+5))/((1+2+3+4+5+6)+

(1+2+3+4+5+6+7)) PROGRAMMAIN2 N=3 Y=(SUM(N)+SUM(N+1)+SUM(N+2))/(SUM(N+3)+SUM(N+4)) WRITE(*,*)‘Y=’,Y END FUNCTIONSUM(X) INTEGERX SUM=0 DO10I=1,X SUM=SUM+I10 CONTINUE END[例9.3]已知sinh(x)=(ex-e-x)/2，由程序輸入x的值，求出sinh(x)，其中ex=1+x+x2/2!+x3/3!+...+xn/n!

PROGRAMMAIN3 REALMYEXP READ(*,*)X Y=(MYEXP(X)-MYEXP(-X))/2.0 WRITE(*,*)‘SINH(‘,X,‘)=’,Y END FUNCTIONMYEXP(X) REALMYEXP N=0 T=1 S=1 DO10WHILE(ABS(T).GT.1E-5) N=N+1

T=T*X/N S=S+T10 CONTINUE MYEXP=S END9.2子例行程序

子例行程序和函數子程序都是子程序，它們的區別在于：函數子程序的名字代表一個值，在函數子程序中求出的函數值存放在函數名中。因此，函數名是函數值的體現者，對函數名應作類型說明。而子例行程序的名字只供調用，它不代表某個值，當然也不屬于某個類型。在子例行程序中求得的值不是由子程序名帶回調用單位，而是通過實參與虛參的聯系帶回調用單位。[例9.4]使用子例行程序求1+2+...+n SUBROUTINESUM(N,S) INTEGERN REALS S=0 DO10I=1,N S=S+I10 CONTINUE END PROGRAMMAIN4 CALLSUM(3,Y1) CALLSUM(4,Y2) CALLSUM(5,Y3) CALLSUM(6,Y4) CALLSUM(7,Y5) Y=(Y1+Y2+Y3)/(Y4+Y5) WRITE(*,*)‘Y=’,Y END11.2.1子例行程序的定義

子例行程序必須以SUBROUTINE語句開頭，以END語句結束。 形式如下：

SUBROUTINE子例行程序名(虛擬參數,…)

SUBROUTINE子例行程序名9.2.2子例程程序的調用

必須用一條獨立的CALL語句來調用子例行程序，CALL語句的形式如下：

CALL子例行程序名(實在參數,…)

當子例行程序沒有虛參時，則調用形式可以以下兩種形式之一：

CALL子例行程序名或CALL子例行程序名()子例行程序和函數子程序的主要差別： 1、名字的作用不同。子例行程序名只供調用時使用，不得在子程序體內賦值；函數子程序名除了供調用時使用外，還代表著函數值，因此必須在返回前在函數體內給函數名賦予函數值。

2、調用的方式不同。函數子程序的調用總是出現在表達式中，并將所得的函數值參與表達式的運算；子例行程序必須用一條獨立的CALL語句調用。

在程序設計中，如果只需得到一個計算值，一般采用函數子程序的形式。在設計函數子程序時通常不主張在函數體內改變形參的值，從而保證不使對應的實參發生變化。如果希望通過子程序的運算操作得到一批數組（例如矩陣相加的結果），則應該采用子例行程序，將得到的新值通過實參與虛參的結合傳送回調用單位，這是由在子例行程序體內給形參賦值實現的。另外，還通過使用子例行程序來完成一些特定的操作。[例9.5]編寫子例行程序用來在一行上打印20個星號SUBROUTINEASTERISKCHARACTER*20STARSTAR=‘********************’PRINT*,STARENDPROGRAMMAIN5CALLASTERISKCALLASTERISKEND[例9.6]編寫子例行程序分別求出5*5矩陣兩個對角線上元素值和 PROGRAMMAIN6 REALW(5,5) CALLREADIN(W) CALLOPP(W,X1,X2) END SUBROUTINEREADIN(A) REALA(5,5) WRITE(*,*)‘Enter5*5matrix:’ DO10I=1,5 READ(*,*)(A(I,J),J=1,5)10CONTINUE END SUBROUTINEOPP(A,S1,S2) REALA(5,5) S1=0 DO10I=1,5 S1=S1+A(I,I)10CONTINUE S2=0 DO20I=1,5 J=5-I+1 S2=S2+A(I,J)20 CONTINUE END9.3實參和虛參之間的數據傳送

函數子程序和子例行程序中的虛參可以是變量名、數組名、函數名、子例行程序名，在子例行程序中還可用星號作為虛參。 當函數子程序和子例行程序被調用之前，所有虛參都無定義，所有虛參既沒有具體的存儲單元也沒有具體的值，只是起到形式上的作用，所以人們也經常將虛擬參數稱為形式參數（或簡稱形參）。當執行的流程轉向子程序時，實參和虛參按地址結合，這時實參的存儲單元就是對應虛參的存儲單元。當執行的流程退出子程序時，該子程序的虛參又變成無定義的。9.3.1變量作為虛參

當虛參是變量時對應的實參可以是用一類型的常量、變量、數組元素和表達式。 如果實參是變量或數組元素，在調用子程序時，對應的虛參實際上將與之共用同一個存儲單元。因此，虛參的值就是實參的值，虛參的值改變時，對應實參的值也就同時改變。PROGRAMMAINSUBROUTINESUB(X,A)INTEGERA,C(3) INTEGERX,ADATAC/3*0/ A=2*XA=100 X=2*ACALLSUB(A,C(2)) ENDWRITE(*,*)‘A=’,A,‘C(2)=’,C(2)END

如果實參是常量或表達式，則子程序中對應的虛參變量的值不應當改變，也就是說在這種情況下虛參變量不應當在子程序中被賦值，否則會引起不可預料的結果。 如果虛參是字符變量，則它的長度定義應當遵循以下兩條規則之一：（1）虛參字符變量的長度定義必須小于等于對應實參變量的長度。（2）虛參字符變量的長度可用（*）來定義，表示長度不定。當調用子程序時，具有不定長度的虛參變量自動定義成為與對應實參具有同樣長度。例如以下程序中：PROGRAMMAIN SUBROUTINESUB(CH)CHARACTERSTR1*8,STR2*40 CHARACTER*(*)CH…. ….CALLSUB(STR1) ENDCALLSUB(STR2)….END9.3.2數組作為虛參

當虛參是數組名時，對應的實參可以是同一類型的數組名或數組元素。以下將區分各種情況加以說明。 （1）如果虛參數組是數值類型或邏輯類型，并且對應實參也是相同的類型，在調用子程序時將把實參數組第一個元素的地址傳送給子程序作為虛參數組的第一個元素的地址，從而導致了它們共用一個存儲單元。由于每個數組占一串連續的存儲單元，因此虛參數組的第二個元素必須按存放次序與實參數組的第二個元素共用存儲單元，其它依次類推。圖9.2給出了一維數組情況下虛參和實參數組結合的示意圖。

PROGRAMMAIN DIMENSIONA(8) CALLSUB(A) END

SUBROUTINESUB(B) DIMENSIONB(-1,5)

…. END圖9.3給出了不同維數情況下，虛參和實參數組結合的示意圖PROGRAMMAINDIMENSIONA(2,4)….CALLSUB(A)….ENDSUBROUTINESUB(B)DIMENSIONB(6)….END圖9.4給出了3*3的實參數組與2*2的虛參數組結合的情況PROGRAMMAINDIMENSIONA(3,3)……CALLSUB(A)……ENDSUBROUTINESUB(B)DIMENSIONB(2,2)……END

注意，在子程序中說明虛參數組時，它的元素個數必須小于等于對應實參數組中元素的個數。虛參數組總是按照內存排列次序與實參數組結合。 雖然虛實結合的數組允許維數不同，下標的上下界不同，但在這種情況下由于對應元素使用的下標完全不同使程序很難讀懂，也很容易造成一些隱藏的錯誤，因此應該盡量避免出現這種情況。（2）如果虛參數組是數值類型或邏輯類型；并且對應實參是一個相同類型的數組元素，在調用子程序時，將把該元素的地址傳送到子程序作為虛參數組第一個元素的地址，從而使它們共用一個存儲單元；并且虛參數組的其余元素將與該實參數組元素后的元素按排列順序一一對應結合。 在這種情況下，子程序中實際用到的虛參數組最后一個元素必須落在對應實參數組的范圍內，而且與數組維數無關。圖9.5給出了當虛參為數組名，對應實參為數組元素時，兩個數組結合的情況。 PROGRAMMAIN DIMENSIONA(8) …… CALLSUB(A(3)) …. END

SUBROUTINESUB(B) DIMENSIONB(4)

…… END圖9.6給出了當虛參為數組名，對應實參為不同維數的實參數組元素時兩個數組結合的情況。PROGRAMMAINDIMENSIONA(3,3)……CALLSUB(A(1,2))……ENDSUBROUTINESUB(B)DIMENSIONB(4)……END圖9.7中虛參數組的元素已超出對應實參數組的范圍。這將引起運行時出錯。PROGRAMMAINDIMENSIONA(6)……CALLSUB(A(3))……ENDSUBROUTINESUB(B)DIMENSIONB(6)……END（3）當數組是字符類型時，虛參和實參數組不是按數組元素的順序一一對應結合，而是按字符位置一一對應結合。虛參數組中允許的字符總數必須小于等于實參數組中允許的字符總數。在此條件下，虛、實數組的維數可以不同，下標的上、下界可以不同，數組元素的字符長度可以不同。圖9.8給出了字符類型數組結合的示意圖。 PROGRAMMAIN CHARACTER*4B(6) …… CALLSUB(B) …… END SUBROUTINESUB(C) CHARACTER*5C(4) …… END

除非特殊需要，虛參字符數組元素的長度應該與對應實參相同，這樣的程序不僅可讀性好，而且易于調試檢查。 與虛參字符數組對應的實參也可以是一個字符數組元素，虛參字符數組的第一個字符與該元素的第一個字符結合，其它按次序依次類推，只是虛參字符數組中最后一個字符必須落在對應實參數組的范圍之內。（4）在子程序中，允許虛參數組是可調數組。可調數組的使用大大提高了子程序的通用性和靈活性。讀者在了解數組虛實結合情況的基礎上應該充分利用可調數組這一強有力的工具來進行程序設計。

PROGRAMMAIN DIMENSIONA(10),B(4,3) M1=4 M2=3 CALLSUB(A,10,B,M1,M2) … END

SUBROUTINESUB(X,NX,Y,N1,N2)

DIMENSIONX(NX),Y(N1,N2) ... END 可調數組名必須是虛參。可調數組中每一維的上下界可以是整型虛參變量，其值通過對應的實參傳送過來；也可以是公用區中變量。為了使程序清晰易讀，建議采用虛參變量來說明可調數組的上、下界而不用公用區變量。另外可調數組也遵循前述1、2、3中的有關規定。對于初學者來說要特別注意，只能在子程序中使用可調數組，而且對于那些只是在子程序中局部使用的（并非通過虛實結合傳送的）數組不允許是可調的。（5）虛實數組中最后一維的上界可以用星號（*）表示。例如：

PROGRAMMAIN SUBROUTINESUB(X,Y) DIMENSIONA(10),B(4,3) DIMENSIONX(0,*),Y(2,*) ... ... CALLSUB(A,B) END ... END

用這種方式定義的數組稱為不定長數組。以上SUB子程序中的X、Y數組在調用之前其大小不固定，當調用該子程序時，它們的大小由對應的實參數組來確定。因此在上例中當執行了CALLSUB(A,B)調用語句后，SUB子程序中X、Y的實際說明為：

DIMENSIONX(0,9),Y(2,3)9.3.3子程序名作為虛參 在虛參表中允許出現函數名和子例行程序名。FORTRAN編譯程序完全根據某個虛參名字在子程序中出現時的上下文關系來確定它是函數名還是子例行程序名。函數名在必要時應該進行類型說明。 例如以下子例行程序中：

SUBROUTINESUB(A,N,X,F) DIMENSIONA(N) ... Y=F(I)+1 X=A(I)+1 ... END

虛參A由說明語句可知是一個數組名，子程序表達式中A(I)的形式是數組元素的引用。N和X分別是整型和實型變量名。虛參中，F是以F(I)的形式出現在表達式中，在FORTRAN中，F(I)的形式可能是數組元素的引用或者可能是函數的調用。由于在SUB子程序中并沒有說明F為數組名，因此FORTRAN編譯程序根據上下文關系確定F為實型函數名。又如以下子例行程序中：

SUBROUTINESUB(X,P) ... CALLP ... END

由于虛參P出現在CALL語句中，因此確定P為子例行程序名。 在虛參中的函數名、子例行程序名（如以上的F、P）只是一個虛設的名字，它并不代表程序中實際存在的某個函數名或子例行程序名，實際的函數名或子例行程序名要在調用時通過實參傳送給它。如果虛參是一個函數名，要求對應的實參代表一個在程序中實際存在的尾部函數名或內部函數名。如果虛參是一個子例行程序名，則要求對應的實參代表一個在程序中實際存在的子例行程序名。 當在調用的程序單位中，實在參數中出現外部函數名或子例行程序名時，必須在調用程序的說明部分用EXTERNAL語句說明這些名字。如果在實在參數表中出現內部函數名時，必須在調用程序的說明部分用INTRINSIC語句說明這些名字。以下我們通過例子來說明函數名和子例行程序名的虛實結合以及EXTERNAL和INTRINSIC語句的使用。[例9.7]編寫一個函數子程序，通過函數名的傳送，使之既能求正切值TAN(X)又能求余切值CTN(X)PROGRAMMAIN7INTRINSICSIN,COSX=3.14159/5.0Y1=TRIANGLE(SIN,COS,X)Y2=TRIANGLE(COS,SIN,X)WRITE(*,*)‘TAN(‘,X*180/3.14159,’)=‘,Y1WRITE(*,*)‘CTN(‘,X*180/3.14159,’)=’,Y2ENDFUNCTIONTRIANGLE(F1,F2,X)TRIANGLE=F1(X)/F2(X)END[例9.8]以下程序由你選擇對輸入的兩個數（X1,X2）進行加法（Add）運算還是減法（Subtract）。當你打入字母A（或a）時，則對隨后輸入的兩個數進行加法運算并輸出結果；當你打入字母S（或s）時，則對隨后輸入的兩個數進行減法運算并打印出結果。 PROGRAMMAIN8 EXTERNALADD,SUB CHARACTERCH WRITE(*,*)‘EnterA(Add)orS(Subtract):’ READ(*,‘(A)’)CH WRITE(*,*)‘EnterX1,X2:’ READ(*,*)X1,X2 IF(CH.EQ.‘A’.OR.CH.EQ.‘a’)THEN CALLTEST(ADD,X1,X2,Y)

WRITE(*,100)X1,‘+’,X2,‘=’,Y

ELSEIF(CH.EQ.‘S’.OR.CH.EQ.‘s’)THEN CALLTEST(SUB,X1,X2,Y)

WRITE(*,100)X1,‘-‘,X2,‘=’,Y ELSE WRITE(*,*)‘Selecterror!’ ENDIFFORMAT(1X,2(F6.1,A),F6.1)ENDSUBROUTINETEST(T,X2,X2,X)CALLT(X1,X2,X)ENDSUBROUTINEADD(A,B,C) C=A+BENDSUBROUTINESUB(A,B,C) C=A-B END 9.3.4星號作為虛參 當形參表中出現一個*號時，對應的實參應該是一個冠有*號的語句標號。例如：

PROGRAMMAIN SUBROUTINEEXAM(A,*，*) ... IF(...)THEN120X=X1+X2 ... ... RETURN1 CALLEXAM(X,*120,*140) ELSEIF(…)THEN … …140… RETURN2 END ENDIF … END

在CALLEXAM(X,*120,*140)語句中，與虛參第一個*號對應的語句標號為120，與虛參第二個*號對應的語句標號為140。在執行EXAM子例行程序時，如果遇到END語句，執行的流程將按正常情況返回到調用語句的后繼語句去繼續執行。當遇到RETURN1語句時，執行的流程返回主程序并跳到與第1個*號對應的語句標號120去繼續執行。當遇到RETURN2語句時，執行的流程返回主程序并跳到與第3個星號對應的語句標號140去繼續執行。 用星號作為形參將使子程序有一個入口而有多個出口，這種返回方式不符合結構化程序設計的要求，因此除非特殊需要，一般不主張采用。9.4利用子程序實現程序的模塊化設計

函數子程序和子例行程序都是獨立的程序單位，可以單獨進行編譯。除了函數子程序名和子例行程序名外，子程序內所用的語句標號、變量、數組、符號常量等都只局限于本程序單位；子程序中可以通過虛實結合使用可調數組，可以通過虛實結合傳送子程序名，這都為結構化程序設計提供了有力的工具。

在開發程序的時候，總是把一個大任務分解成若干個小子任務，并且更進一步把子任務分解成更小的子任務，一直細分到每一個子任務只完成比較簡單的操作為止。當把所有的子任務都由子程序來完成時，通常細分的程度以每個子程序不超過50-60語句行為宜。子程序可以多次調用，這使得一些相同的操作都可調用同一個子程序來完成而不必重復編程。用這種方法設計出現的程序不僅可讀性好而且容易調試、修改和維護。在調試程序時可以一個個模塊逐塊進行；當需要修改程序時，對一個模塊的修改往往不會影響到其它模塊。主程序子程序1子程序n……9.5在子程序中的SAVE語句和DATA語句

函數子程序或子例行程序中用到的所有變量，在被調用前通常都沒有確定的存儲單元，每當子程序被調用時才臨時分配存儲單元，而且在退出子程序時這些存儲單元又都被釋放并重新分配另作它用。所以這些變量的值不被保留。在下次進入子程序時，給這些變量分配的可能是另外的一些存儲單元，上次調用時這些變量所具有的值已不復存在。稱這些變量在子程序未調用時是無定義的。 在函數或子例行程序中可以使用SAVE說明語句來指定子程序中的某些變量的存儲單元不被釋放，它們的內容在退出子程序期間保持不變，在下次調用時仍可使用。SAVE語句的形式如下：

SAVE變量表或SAVE

變量表可以是本子程序中的變量名、數組名，若變量或數組名不止一個，它們之間用逗號隔開。在SAVE語句中不允許出現虛擬參數名。當SAVE語句中帶有變量表時，則為變量表中指定的變量和數組保留存儲單元，如果SAVE語句中沒有變量表，則表示讓子程序中所有的變量和數組在程序執行期間都占有確定的存儲單元，不因退出子程序而被釋放。[例9.9]求階乘的函數，每次調用此函數時，把最后一項乘數放在變量N0中，把階乘的值放在變量F中。當下次調用此函數時，如果N大于等于N0，則不必再從頭開始進行連乘運算，只需把上次的函數值再乘以N0+1、N0+2、…N即可。 PROGRAMMAIN9 WRITE(*,*) WRITE(*,*)‘ENTERN:’ READ(*,*)N DO10WHILE(N.GE.0) WRITE(*,*)’FUN=’,FUN(N) WRITE(*,*)‘ENTERN:’ READ(*,*)N10 CONTINUE END FUNCTIONFUN(N) SAVEN0,F DATAN0/0/,F/1/ IF(N.LT.0)THEN WRITE(*,*)‘DATAERROR!’ FUN=-1 RETURN ELSEIF(N.LE.N0)THEN F=1 N0=1 ENDIF DO10I=N0+1,N F=F*I10 CONTINUE FUN=F N0=N END

對于只需在第一次調用時賦初值，后續調用不必賦初值的情況，使用DATA語句是合適的。當用DATA語句給子程序中的變量或數組賦初值時，這些變量和數組元素只是在子程序第一次調用時才有初值。如果沒有SAVE語句說明，DATA語句中的變量和數組在第一次調用退出后就變成無定義的了。不要誤解為調用一次子程序DATA語句就賦一次初值，因為DATA是非執行語句。如果每調用一次就需要賦一次相同的初值，這應該用賦值語句來實現。FUNCTIONFAC(N)FUNCTIONFAC(N)

DATAF/1/ F=1IF(N.LT.0)THEN IF(N.LT.0)THEN WRITE(*,*)‘ERROR’WRITE(*,*)‘DATAERROR!’ FAC=-1 FAC=-1 RETURN RETURNELSE ELSE DO10I=1,N DO10I=1,N

F=F*I F=F*I10 CONTINUE 10CONTINUE ENDIF ENDIF FUN=F FUN=F END END[例9.10]給出一個產生偽隨機數的函數RAND()，函數沒有自變量，產生的隨機數在0到1之間，不包括0和1 PROGRAMMAIN10 DIMENSIONA(10) DO10I=1,10 A(I)=RAND()10 CONTINUE WRITE(*,*)A END

FUNCTIONRAND() SAVESEED INTEGERSEED,C1,C2,C3 PARAMETER(C1=29,C2=217,C3=1024) DATASEED/0/ SEED=MOD(SEED*C1+C2,C3) RAND=REAL(SEED)/C3 END從程序運行的結果可以看到，每次產生相同的隨機數序列。這是因為每次運行開始，隨機數的種子SEED的初值總是為0。如果每次運行時想從終端輸入種子SEED的初值，從而使每次運行時產生不同的隨機數序列，可把函數改寫如下： FUNCTIONRAND() SAVESEED INTEGERSEED,C1,C2,C3 PARAMETER(C1=29,C2=217,C3=1024) DATASEED/0/ IF(SEED.EQ.0)THEN WRITE(*,*)‘ENTERRANDNUMBERSEED:’ READ(*,*)SEED ENDIF SEED=MOD(SEED*C1+C2,C3) RAND=REAL(SEED)/C3 END9.6程序舉例

[例9.11]編寫函數求去掉尾部空格后字符串的長度 PROGRAMMAIN11 CHARACTER*80LINE READ(*,*)LINE L=LENGTH(LINE) WRITE(*,*)L END FUNCTIONLENGTH(LINE) CHARACTER*(*)LINE I=LEN(LINE) DO10WHILE(I.GT.0.AND.LINE(I:I).EQ.‘‘) I=I-110CONTINUELENGTH=IEND[例9.12]編寫子例行程序把字符串中的小寫字母轉換成大寫字母，其它字符不變

SUBROUTINETOUPCASE(STR) CHARACTER*(*)STR,CH*1 L=LEN(STR) DO10I=1,L,1 IF(STR(I:I).GE.‘a’.AND.STR(I:I).LE.‘z’)THEN CH=CHAR(ICHAR(STR(I:I))+ICHAR(‘A’)-ICHAR(‘a’)) STR(I:I)=CH ENDIF10 CONTINUE ENDPROGRAMMAIN12CHARA