SUBROUTINE XAJ(   FILELEN    ,&
                NODE    ,&		! 单元出流数组大小								//输入变量
                M		,&		! 无因次单位数组大小							//输入变量
                WUM		,&		! 上层张力水容量wum							//输入变量
                W1		,&		! 下层张力水容量wl							//输入变量
                WDM		,&		! 深层张力水容量wdm							//输入变量
                KC	    ,&		!蒸发能力折算系数KC							//输入变量
                C	    ,&		!深层蒸发系数c							    //输入变量
                B	    ,&		!张力水蓄水容量系数b						    //输入变量
                IMP1    ,&		!不透水面积比率imp1							//输入变量
                SM	    ,&		!自由水蓄水容量sm								//输入变量
                EX	    ,&		!自由水蓄水容量指数ex							//输入变量
                KG	    ,&		!地下水出流系数kg								//输入变量
                KSS	    ,&		!壤中流出流系数kss								//输入变量
                KKG	    ,&		!地下水出流系数kkg								//输入变量
                KKSS    ,&		!壤中流出流系数kkss								//输入变量
                NAREA	,&		! 单元面积个数                                  //输入变量
                AREA	,&		! 单元面积										//输入变量
                UH		,&		! 无因次单位线									//输入变量
                DT		,&		! 时段步长										//输入变量
                P		,&		! 降雨系列										//输入变量
                EP		,&		! 蒸发皿蒸发能力								//输入变量
                W		,&		! 土壤含水层	1.上层	2.下层	3.深层			//输入变量
                FR		,&		! 初始产流面积									//输入变量
                S		,&		! 初始自由水深									//输入变量
                QRSS0	,&		! 初始壤中流流量								//输入变量
                QRG0	,&		! 初始地下水径流量								//输入变量
                Q20     ,&      ! 出流流量初始值              浮点数            //输入     
                X       ,&      ! 流量比重因子          浮点数                  //输入参数
                K       ,&      ! 流量模数              浮点数                  //输入参数
                DETAT   ,&      ! 时间步长      （小时）    浮点数              //输入参数
                C0      ,&      ! 系数                      浮点数              //输出参数 
                C1      ,&      !                           浮点数              //输出参数 
                C2      ,&      !                          浮点数               //输出参数 
                QOUT      )       ! 出流流量                 浮点数               //输出 

    IMPLICIT NONE

    INTEGER::NODE	


    REAL::X
    REAL::K
    REAL::DETAT
    
    REAL::C0
    REAL::C1
    REAL::C2    
  
    INTEGER:: FILELEN
	INTEGER::M					!
	REAL::PAR(13)				! 1.上层张力水容量wum	2.下层张力水容量wl	3.深层张力水容量wdm
								! 4.蒸发能力折算系数KC.深层蒸发系数c		6.张力水蓄水容量系数b
								! 7.不透水面积比率imp1	8.自由水蓄水容量sm	9.自由水蓄水容量指数ex
								!10.地下水出流系数kg 11.壤中流出流系数kss	12.地下水出流系数kkg
								!13.壤中流出流系数kkss
    
    REAL::WUM                    !上层张力水容量wum	
 	REAL::W1                     !下层张力水容量wl	
 	REAL::WDM                   ! 深层张力水容量wdm
	REAL::KC					! 蒸发能力折算系数
	REAL::C						! 深层蒸发系数	
	REAL::B						! 张力水蓄水容量系数
 	REAL::IMP1					! 不透水面积比率		
	REAL::SM					! 自由水蓄水容量
	REAL::EX					! 自由水蓄水容量指数
	REAL::KG					! 地下水出流系数
	REAL::KSS					! 壤中流出流系数
	REAL::KKG					! 地下水出流系数	
	REAL::KKSS					! 壤中流出流系数		
    
    integer::NAREA			! 单元面积个数                                  //输入变量
    REAL::AREA(NAREA)					! 单元面积
    REAL::Q1(NAREA,NODE)
    REAL::UH(M)					! 无因次单位线
	REAL::DT					! 时段步长
	REAL::P(NAREA,NODE)					! 降雨系列
	REAL::EP(NAREA,NODE)					! 蒸发皿蒸发能力
	REAL::QR(NAREA,NODE)					! 单元出流
	REAL::W(3)					! 土壤含水层	1.上层	2.下层	3.深层
	REAL::FR					! 初始产流面积
	REAL::S						! 初始自由水深
	REAL::QRSS0					! 初始壤中流流量
	REAL::QRG0					! 初始地下水径流量
    REAL::Q2(NAREA,NODE)   
    INTEGER::I,J
    REAL::QOUT(NODE)
    REAL::Q20(NAREA)    
    
    DO I = 1,NODE
        QOUT(I) = 0.0
    END DO
    
    DO I =1,NAREA
    
        CALL XAJMX(   FILELEN    ,&
                        NODE		,&		! 单元出流数组大小								//输入变量
					    M		,&		! 无因次单位数组大小							//输入变量
					    WUM		,&		! 上层张力水容量wum							//输入变量
                        W1		,&		! 下层张力水容量wl							//输入变量
                        WDM		,&		! 深层张力水容量wdm							//输入变量
                        KC	    ,&		!蒸发能力折算系数KC							//输入变量
                        C	    ,&		!深层蒸发系数c							    //输入变量
                        B	    ,&		!张力水蓄水容量系数b						    //输入变量
                        IMP1    ,&		!不透水面积比率imp1							//输入变量
                        SM	    ,&		!自由水蓄水容量sm								//输入变量
                        EX	    ,&		!自由水蓄水容量指数ex							//输入变量
                        KG	    ,&		!地下水出流系数kg								//输入变量
                        KSS	    ,&		!壤中流出流系数kss								//输入变量
                        KKG	    ,&		!地下水出流系数kkg								//输入变量
                        KKSS    ,&		!壤中流出流系数kkss								//输入变量
					    AREA(I)	,&		! 单元面积										//输入变量
					    UH		,&		! 无因次单位线									//输入变量
					    DT		,&		! 时段步长										//输入变量
					    P(I,:)		,&		! 降雨系列										//输入变量
					    EP(I,:)		,&		! 蒸发皿蒸发能力								//输入变量
					    W		,&		! 土壤含水层	1.上层	2.下层	3.深层			//输入变量
					    FR		,&		! 初始产流面积									//输入变量
					    S		,&		! 初始自由水深									//输入变量
					    QRSS0	,&		! 初始壤中流流量								//输入变量
					    QRG0	,&		! 初始地下水径流量								//输入变量
					    QR(I,:)		)		! 单元出流										//输出变量
    
    DO J = 1,NAREA
        
        Q1(J,:)=QR(J,:)
        
    END DO
    
    CALL Mc_method( NODE           ,&          ! 入流流量数据点个数     整数         //输入参数
                    Q1(I,:)          ,&          ! 入流流量              浮点数            //输入参数
                    Q20(I)         ,&        ! 出流流量初始值              浮点数         //输入     
                    X           ,&          ! 流量比重因子          浮点数           //输入参数
                    K           ,&          ! 流量模数              浮点数             //输入参数
                    DETAT       ,&          ! 时间步长      （小时）    浮点数             //输入参数
                    C0          ,&          ! 系数                      浮点数          //输出参数  
                    C1          ,&          !                           浮点数         //输出参数  
                    C2          ,&           !                          浮点数           //输出参数  
                    Q2(I,:)          )           ! 出流流量                 浮点数                 //输出

    END DO  
    
    DO I = 1,NAREA
        DO J = 1,NODE
            QOUT(J)= Q2(I,J) + QOUT(J)
        END DO
    END DO
    
END SUBROUTINE XAJ