小编Jer*_*amu的帖子

我有九个矩阵,其维数为(N乘N) A1(i,j),A2(i,j),A3(i,j),A4(i,j),A5(i,j),A6(i,j),A7(i,j),A8(i,j),A9(i,j)

然后我想构建一个更大的矩阵(3N乘3N),包括这九个矩阵:

A = [A1 A2 A3
     A4 A5 A6
     A7 A8 A9]

在fortran中,我可以使用命令行吗？

do i=1,FN
   do j=1,FML
      A(i,j) = [A1(i,j),A2(i,j),A3(i,j);A4(i,j),A5(i,j),A6(i,j);A7(i,j),A8(i,j),A9(i,j)]
   end do
end do

我有一个关于如何在颜色条上方添加单位的问题.我的代码如下:

hc=colorbar;
xlabel(hc,'psi');

但是,它显示单位不在彩条之上.

Integer :: NBE,ierr,RN,i,j
Real(kind=8), allocatable :: AA1(:,:),AA2(:,:)

NBE=40 
RN=3*NBE-2
Allocate(AA1(3*NBE,3*NBE),AA2(3*NBE,RN),stat=ierr)

If (ierr .ne. 0) Then
    print *, 'allocate steps failed 1'
    pause
End If

Do i=1,3*NBE
  Do j=1,3*NBE
     AA1(i,j)=1
  End Do
End Do

我想从矩阵中删除列97和113 AA1,然后这个矩阵变为AA2.我只想知道Fortran的任何命令是否可以实现此操作？

当我编译下面的程序时,我在call Coor_Trans命令行中出现错误和警告

警告:行被截断

错误:参数列表中的语法错误

我多次编译程序,但它不起作用.也许我的呼叫命令有问题.

program 3D

      implicit none
      integer :: i,j,k
      integer, parameter :: FN=2,FML=5,FMH=5
      integer, parameter :: NBE=FN*FML*FMH
      real, parameter :: pi = 4*atan(1.0)
      real(kind=4), dimension(1:FN,1:FML+1,1:FMH+1) :: BEXL,BEYL,BEZL
      real(kind=4), dimension(1:FN,1:FML,1:FMH) :: BEXC,BEYC,BEZC,BE2A,BE2B,ANGLE
      real(kind=4), dimension(1:NBE,1:1,1:1) :: BEXC1,BEYC1,BEZC1,BE2A1,BE2B1,ANGLE1
      real(kind=4), dimension(1:NBE,1:NBE) :: LOC_PTS1,LOC_PTS2,LOC_PTS3
      real :: LOC_1,LOC_2,LOC_3

      do i=1,FN
        do j=1,FML
           do k=1,FMH

            BEXC(i,j,k) = 0.5*(BEXL(i,j,k) + BEXL(i,j+1,k))
            BEYC(i,j,k) = 0.5*(BEYL(i,j,k) + BEYL(i,j+1,k))
            BEZC(i,j,k) = 0.5*(BEZL(i,j,k) + BEZL(i,j,k+1))
            BE2A(i,j,k) = FL(i)/FML + j*0 + k*0
            BE2B(i,j,k) = FH(i)/FMH + j*0 …

当我编译下面的程序时，我有一个错误

函数'j17'的返回类型不匹配（未知/真实（4））”

program Jfun

    implicit none
    real :: R = 1.5, X1 = 2.5, X3 = 4.8, EPS1 = 0.5
    write(*,*) J17(R,X1,X3,EPS1)

end program Jfun

real function J17(R,X1,X3,EPS1)

    implicit none
    real, intent(in) :: R,X1,X3,EPS1
    J17 = -((R+X1-EPS1)*(X3**2-R**2)+(X3**2)*R)/(R**3)/((R+X1-EPS1)**2)

end function J17

我该如何解决这个问题？

我知道Python中的sorted()函数可以在数字列表或字典上运行,并返回已排序的数字或键的列表.例如:

w = [1,4,6,8,2,0,3,2]
sorted(w)

然后我们可以获得重新排列的排序数字列表为[0,1,2,2,3,4,6,8]

我的问题是:sorted()函数有没有任何其他功能？我希望你们能给我一些简单的例子.:)