program f90_NUM4seep !2次元飽和不飽和定常浸透流解析データrenumbering implicit none integer::i,j,k,l,ne character::strcom*256 !書出用コメント integer::NODT !節点総数 integer::NELT !要素総数 integer::MATEL !材料種類数 integer::KOH !全水頭指定節点数 integer::KOQ !流量指定節点数 integer::KOU !浸潤境界指定節点数 integer::nod=4 !1要素節点数 integer::nhen=1 !1節点自由度 integer::nt integer,allocatable::kakom(:,:) !要素構成節点番号 integer,allocatable::matno(:) !材料種別No real(8),allocatable::wm(:,:) !材料物性作業用 real(8),allocatable::x(:) !節点x座標 real(8),allocatable::z(:) !節点y座標 integer,allocatable::nokh(:) !全水頭指定節点番号 integer,allocatable::nokq(:) !流量指定節点番号 integer,allocatable::noku(:) !浸潤境界節点番号 real(8),allocatable::Hinp(:) !指定全水頭入力 real(8),allocatable::Qinp(:) !指定流量入力 real(8),allocatable::hvec(:) !全水頭入力 integer,allocatable::numr(:) !元の節店番号格納 integer,allocatable::jnd(:) !新しい節点番号格納 character :: linebuf*1000 !*************************** !data input !*************************** open(11,file='fnameR.csv',status='old') !!-------------------------------------------------------------------------------------- !!書出用コメント入力 !!-------------------------------------------------------------------------------------- read(11,'(a)') strcom !------------------------------------------------------------------------------------------ !節点数,要素数,材料種類数,全水頭指定節点数,流量指定節点数,浸潤境界節点数 !------------------------------------------------------------------------------------------ read(11,*) NODT,NELT,MATEL,KOH,KOQ,KOU !------------------------------------------------------------------------------------------ !配列寸法宣言 !------------------------------------------------------------------------------------------ nt=NODT*nhen allocate(kakom(1:NELT,1:nod)) allocate(matno(1:NELT)) allocate(wm(1:NELT,1:3)) allocate(x(1:NODT)) allocate(z(1:NODT)) if(0