calculate the number of the cells

! ! calculate the number of the cells
! neall=0
! itotal=0
! do j=1,NJ-1
! if(j>1)then
! itotal=itotal+NI(j-1)
! endif
! if(NI(j)==NI(j+1))then ! NI(j) = NI(j+1)
! do i=1,NI(j)-1
! neall=neall+1
!
! p1(neall)=itotal+i
! p2(neall)=p1(neall)+1
! p3(neall)=p2(neall)+NI(j)
! p4(neall)=p3(neall)-1
! ! write(11,"(4(X,I5))")p1(neall),p2(neall),p3(neall),p4(neall)
! enddo
! elseif(NI(j) < NI(j+1))then ! NI(j) < NI(j+1)
! do i=1,NI(j)-1
! neall=neall+1
!
! p1(neall)=itotal+i
! p2(neall)=p1(neall)+1
! p3(neall)=p2(neall)+NI(j)
! p4(neall)=p3(neall)-1
! ! write(11,"(4(X,I5))")p1(neall),p2(neall),p3(neall),p4(neall)
! enddo
! do i=NI(j),NI(j+1)-1
! neall=neall+1
!
! p1(neall)=itotal+NI(j)
! p2(neall)=p1(neall)
! p3(neall)=p2(neall)+i+1
! p4(neall)=p3(neall)-1
! ! write(11,"(4(X,I5))")p1(neall),p2(neall),p3(neall),p4(neall)
! enddo
! else ! NI(j) > NI(j+1)
! do i=1,NI(j+1)-1
! neall=neall+1
!
! p1(neall)=itotal+i
! p2(neall)=p1(neall)+1
! p3(neall)=p2(neall)+NI(j)
! p4(neall)=p3(neall)-1
! ! write(11,"(4(X,I5))")p1(neall),p2(neall),p3(neall),p4(neall)
! enddo
! do i=NI(j+1),NI(j)-1
! neall=neall+1
!
! p1(neall)=itotal+i
! p2(neall)=p1(neall)+1
! p3(neall)=itotal+NI(j)+NI(j+1)
! p4(neall)=p3(neall)
! ! write(11,"(4(X,I5))")p1(neall),p2(neall),p3(neall),p4(neall)
! enddo
! endif
! enddo
! ! output total
! open(11,file="output/"//TRIM(fname)//"/"//"contour_1.dat")
! write(11,"(A)") 'TITLE = " 3D mesh " '
! write(11,"(A)") 'VARIABLES ="x(mm)", "y(mm)", "z(mm)", "beta", "gamma",'//&
! '"Ps_5T_Pa", "Pt_5T_Pa", "rho_5T", "V_5T", "Cps_pt_5T", "omega_pt_5T", "Unor_pt_5T","Ux_flow","Uy_flow","Uz_flow" '
! write(11,"(A,I3,A,I3,A)")' ZONE T="'//TRIM(fname)//'" n=',ntotal,',E=',neall,',ZONETYPE=FEQuadrilateral, DATAPACKING=POINT'
! do j=1,NJ
! do i=1,NI(j)
! write(11,"(15(E19.10e3))") &
! dx,&
! RALL(i,j,PC_REAL_L_mm), RALL(i,j,PC_REAL_T_mm), &
! RALL(i,j,beta_5T), RALL(i,j,gamma_5T), &
! RALL(i,j,Ps_5T_Pa),RALL(i,j,Pt_5T_Pa), &
! RALL(i,j,rho_5T), RALL(i,j,V_5T), &
! RALL(i,j,Cps_pt_5T),RALL(i,j,omega_pt_5T),RALL(i,j,Unor_pt_5T),&
! RALL(i,j,Ux_flow),RALL(i,j,Uy_flow),RALL(i,j,Uz_flow)
! enddo
! enddo
! do ieall=1,neall
! write(11,"(4(X,I5))")p1(ieall),p2(ieall),p3(ieall),p4(ieall)
! enddo
! close(11)

fortran读取文件的行书

integer::i,N,stat

!get data from file

stat=0 !这个判断读取是否成功，

N=1 ！这个是用来计数行数的。

open (99,file=\'myfile.dat\') !打开你想要的文件，

do while(stat==0) ！当读取成功的话，继续读取

read(99,*,iostat=stat)g ！这里我假定你一行数据只有一个实数

N=N+1 ！增加行数

end do

rewind(99) ！回到这个文件的开头

ln -s TARGET LINK_NAME

TARGET LINK_NAME

原文件 链接文件

需要完整地址