anchor.diff:  	README(IPP(I)) = "//$ELT_BC"
   	README(IPP(I))=  "//"//'"'//"THE KEYWORD ELT_BC IS USED TO DEFINE THE BOUNDARY CONDITION APPLIED ON THE ELEMENT GROUP."//'"'
   	README(IPP(I))="//NOTE THAT THE FINAL BC VALUE OF A NODE IS WHAT INPUT BY THE FIRST ELEMENT SHARING THE NODE."
  -	README(IPP(I))= "//{NELT_BC(I)} "
  -	README(IPP(I)) = "//{NO(I),GROUP(I),NDIM(I),DOF(I),STEPFUNC(I),VALUE(R)}  // NDIM=0,1,2,3 ֱʾ㡢ߡ桢Լ. Ϊ߽λƣˮͷȣߵԣNDIMȡֵԱ߽Ӱ."
  +	README(IPP(I))= "//{NELT_BC} "
  +	README(IPP(I)) = "//{NO,GROUPID,NDIM,DOF,STEPFUNC,VALUE}  // NDIM=0,1,2,3 ֱʾ㡢ߡ桢Լ. Ϊ߽λƣˮͷȣߵԣNDIMȡֵԱ߽Ӱ."
   	README(IPP(I))=  "//                                                      // DOF(I)=1,2,...,7,ֱʾԼX,Y,Z,H,MX,MY,MZ "
  -    README(IPP(I))=  "//{......}   //NELT_BC. "
  +    README(IPP(I))=  "//NOTE: FOR THE LAYERED MODEL, THE GROUP IDS OF EACH EXTRUDED LAYERED IS REFERED AS GROUP_BASEMESH*100+ILAYER.AS SHOWN BELOW." 
  +	README(IPP(I))=  "//   -----------------								"
  +	README(IPP(I))=  "//   |               |                                "
  +	README(IPP(I))=  "//   |    N-LAYER    | GROUPID=GROUP_BASEMESH*100+N   "
  +	README(IPP(I))=  "//   |               |                                "
  +	README(IPP(I))=  "//   -----------------                                "
  +	README(IPP(I))=  "//   |               |                                "
  +	README(IPP(I))=  "//   |    ..LAYER    |                                "
  +	README(IPP(I))=  "//   |               |                                "
  +	README(IPP(I))=  "//   -----------------                                "
  +	README(IPP(I))=  "//   |               |                                "
  +	README(IPP(I))=  "//   |    2-LAYER    | GROUPID=GROUP_BASEMESH*100+2   "
  +	README(IPP(I))=  "//   |               |                                "
  +	README(IPP(I))=  "//   -----------------                                "
  +	README(IPP(I))=  "//   |               |                                "
  +	README(IPP(I))=  "//   |    1-LAYER    | GROUPID=GROUP_BASEMESH         "
  +	README(IPP(I))=  "//   |               |                                "
  +	README(IPP(I))=  "//   -----------------                                "
  +	README(IPP(I))=  "//{......}   //NELT_BC. "
   	README(IPP(I)) = "//$ENDELT_BC"
   	
   	
  
positive.code: 	README(IPP(I)) = "//$ELT_LOAD"
  	README(IPP(I))=  "//"//'"'//"THE KEYWORD ELT_LOAD IS USED TO DEFINE THE LOAD APPLIED ON THE ELEMENT GROUP."//'"'
  	README(IPP(I))="//NOTE THAT, IF THE NDIM=1,2 OR 3, THE FINAL BC VALUE OF A NODE IS THE SUM OF ALL THE ELEMENT SHARING THE NODE."
  	README(IPP(I)) = "//{NELT_LOAD(I)} "
  	README(IPP(I)) = "//{NO(I),GROUP(I),NDIM(I),DOF(I),STEPFUNC(I),VALUE(R)}  // NDIM=0,1,2,3 ֱʾ㡢ߡ桢ĺ."
  	README(IPP(I)) =  "//{......}   //NELT_LOAD. "
  	README(IPP(I)) ="//Ŀǰɴ:NDIM=0(); "
  	README(IPP(I)) ="//			 NDIM=1Ŀǰɴ2ڵ㡢3ڵԼ5ڵĵԪߣҽھߺ;"
  	README(IPP(I)) ="//			 NDIM=2Ŀǰɴ3ڵ㣨Σ4ڵ㣨ıΣ6ڵ㣨Σ8ڵ㣨ıΣԼ15ڵ㣨ĴΣĵԪ棬ҽھ;"
  	README(IPP(I)) ="//			 NDIM=3Ŀǰɴ4ڵ㣨壩10ڵ㣨壩6ڵ(Prism)8ڵ()ҽھ;"
  	README(IPP(I)) = "//$ENDELT_LOAD"		
  	
  	README(IPP(I)) ="\N//******************************************************************************************************"C
  	README(IPP(I)) = "//$ELT_SPGFACE"
  	README(IPP(I))=  "//"//'"'//"THE KEYWORD ELT_SPGFACE IS USED TO DEFINE THE SEEPAGEFACE CONDITION APPLIED ON THE ELEMENT GROUP."//'"'
  	README(IPP(I)) = "//{NELT_SPGFACE(I)} "
  	README(IPP(I)) = "//{NO(I),GROUP(I),NDIM(I),DOF(I),STEPFUNC(I),VALUE(R)}  // NDIM=0,1,2,3 ɣVALUE=ֵڲΪYZ"
  	README(IPP(I)) =  "//{......}   //NELT_SPGFACE. "
  	README(IPP(I)) = "//$ENDELT_SPGFACE"	
  	
  
negative.code: SUBROUTINE GEN_LAYERED_ELEMENT()
      USE DS_Gmsh2Solver
      IMPLICIT NONE
      INTEGER::I,J,K,N1,N2,N3,N4,IGP1,NA1(100),ALLOC_ERR,NSUBLAYER1
      REAL(8)::T1
      INTEGER,ALLOCATABLE::SUMLAYER1(:)
      REAL(8),ALLOCATABLE::AR1(:)
  	character(32)::CH1
      type(element_type),allocatable::element1(:)
      type(node_type),allocatable::node1(:)
      
      ALLOCATE(SUMLAYER1(0:NLAYER))
      SUMLAYER1(0)=0
      DO I=1,NLAYER
          SUMLAYER1(I)=SUBLAYER(I)+SUMLAYER1(I-1)
      ENDDO
      NSUBLAYER1=SUMLAYER1(NLAYER)
      ALLOCATE(ELEMENT1(NSUBLAYER1*NEL),STAT=ALLOC_ERR)
   	
      DO CONCURRENT (I=1:NEL)
          N1=(I-1)*NSUBLAYER1
          N3=ELEMENT(I).NNODE
          N2=2*N3
          NA1(1:N3)=(ELEMENT(I).NODE-1)*(NSUBLAYER1+1)
          DO CONCURRENT (J=1:NLAYER)
              DO CONCURRENT (K=1:SUBLAYER(J))
              N4=N1+SUMLAYER1(J-1)+K
              ELEMENT1(N4).NNODE=N2
              ELEMENT1(N4).ILAYER=J
              ELEMENT1(N4).NTAG=ELEMENT(I).NTAG
  			SELECT CASE(ELEMENT(I).ET)
  				CASE(1)	
  					ELEMENT1(N4).ET=3					
  				CASE(2)
  					ELEMENT1(N4).ET=6
  				CASE(3)
  