anchor.diff:  module DS_Gmsh2Solver
       INTEGER::NLAYER=1
  +    INTEGER,ALLOCATABLE::SUBLAYER(:)
       LOGICAL::ISGENLAYER=.FALSE.
   	type node_type
   		integer::inode=0
  
positive.code: 	
  	
  	type physicalGroup_type
  		integer::ndim
  		logical::ismodel=.FALSE. !is a part of the model
  		integer::ET_GMSH=0 
  
negative.code: 	README(IPP(I)) ="\N//******************************************************************************************************"C
  	README(IPP(I)) = "//$MODELDIMENSION"
  	README(IPP(I))=  "//"//'"'//"THE KEYWORD MODELDIMENSION IS USED TO DEFINE THE MODEL DIMENSION."//'"'
  	README(IPP(I)) = "//{MODELDIMENSION}   // MODELDIMENSION=2, 3"  
  	README(IPP(I)) = "//$ENDMODELDIMENSION"	
  	
  	README(IPP(I)) ="\N//******************************************************************************************************"C
  	README(IPP(I)) = "//$GROUPPARAMETER"
  	README(IPP(I))=  "//"//'"'//"THE KEYWORD GROUPPARAMETER IS USED TO DEFINE THE MATERIAL AND ELEMENT TYPE IN THE MODELGROUP ELEMENTS."//'"'
  	README(IPP(I)) = "//{NGROUPPARAMETER}"  
  	README(IPP(I)) = "//{GROUDID,MATID(NLAYER),ET)}  //ET=CPE3,CPE6,TET4,TET10,... "
  	README(IPP(I)) = "//{......}  // NGROUPPARAMETER "     
  	README(IPP(I)) = "//$ENDGROUPPARAMETER"	
      
  	README(IPP(I)) ="\N//******************************************************************************************************"C
  	README(IPP(I)) = "//$WSP"
  	README(IPP(I))=  "//"//'"'//"THE KEYWORD WSP IS USED TO DEFINE THE BOUNDRAY LINE WHERE A WATER SURFACE PROFILE IS CALCULATED."//'"'
  	README(IPP(I)) = "//{NWSP}"  
  	README(IPP(I)) = "//{GROUPID,CHARACTERPOINTGROUP,STARTPOINTGROUP}*{NWSP}  // ÿ3NWSPУSTARTPOINTGROUPΪ(GROUPֻһ)㰴˳ϸ㡣"
  	README(IPP(I)) = "//עCHARACTERPOINTGROUPΪε㡣Щ㽫ֳɲͬʵ"
  	README(IPP(I)) = "//$ENDWSP"	
  	
  	README(IPP(I)) ="\N//******************************************************************************************************"C
  	README(IPP(I)) = "//$DATAPOINT"
  	README(IPP(I))=  "//"//'"'//"THE DataPoint IS USED TO OUTPUT THE DATA POINTS ALONG A POLYLINE YOU WANT TO LIST THE COMPUTATED RESULT."//'"'
  	README(IPP(I)) = "//{NDataPoint)}"  
  	README(IPP(I)) = "//{GROUPID [,ORDER=0,STARTPOINTGROUP]}*{NDataPoint(I)}  // ÿ3NDataPoint.STARTPOINTGROUPΪ(GROUPֻһ)㰴˳ϸ㡣"
  	README(IPP(I)) = "//ORDER=0,ʾɢʽ.=1ʾSTARTPOINTGROUPΪ㰴˳"
  	README(IPP(I)) = "//$ENDDATAPOINT"		
  	
  	README(IPP(I)) ="\N//******************************************************************************************************"C
  	README(IPP(I)) = "//$INCLUDEFILE"
  	README(IPP(I))=  "//"//'"'//"THE INCLUDEFILE IS USED TO READ IN DATA IN ANOTHER FILE."//'"'
  	README(IPP(I)) = "//{NINCLUDEFILE}"  
  	README(IPP(I)) = "//{ABSOLUTED PATH OF THE FILES}*{NINCLUDEFILE(I)}  // NINCLUDEFILE."
  	README(IPP(I)) = "//$ENDINCLUDEFILE"    
  	
  	README(IPP(I)) ="\N//******************************************************************************************************"C
  	README(IPP(I)) = "//$COPYFILE"
  	README(IPP(I))=  "//"//'"'//"THE COPYFILE IS USED TO COPY DATA IN THE FILES TO THE OUTPUT FILE."//'"'
  	README(IPP(I)) = "//{NCOPYFILE}"  
  	README(IPP(I)) = "//{ABSOLUTED PATH OF THE FILES}*{NCOPYFILE(I)}  // NCOPYFILE."
  	README(IPP(I)) = "//$ENDCOPYFILE" 
   
  