爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4860|回复: 10

[求助] fortran

[复制链接]

新浪微博达人勋

发表于 2014-11-19 20:57:10 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
integer a(5,5)
integer i,j,k
do i=1,5
    a(i,j)=0,j=1,5
end do
a(2,3)=10
do i=1,5
write(*,*)(a(i,j),j=1,5)
end do
end

我想用隐式循环赋值,应该把
do i=1,5
    a(i,j)=0,j=1,5
end do
改成什么呢

密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 0
发表于 2014-11-19 21:14:19 | 显示全部楼层
你的这个运算就显式就行了。
do i=1,m
do j=1,n
a(i,j)=0
enddo
enddo

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-11-20 08:18:26 | 显示全部楼层
言深深 发表于 2014-11-19 21:14
你的这个运算就显式就行了。
do i=1,m
do j=1,n

直接a=0不就完了!!!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-11-20 09:03:16 | 显示全部楼层
我不明白你到底是要干嘛,你的意思好像是要把所有的值都设为0.除了a(2,3)=10
但是a(i,j)=0,j=1,5是有语法错误的
program ex
    implicit none
    integer :: a(5,5)
    integer :: i,j
   
    forall(i=1:5,j=1:5)
    a(i,j)=0
    end forall

a(2,3)=10

write(*,'(5(2x,i2),/)') ((a(i,j),j=1,5),i=1,5)
end

应该可以的
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-11-20 09:08:18 | 显示全部楼层
楼上正解,矩阵赋初值,直接矩阵名=初值即可
real a(n,m)
a=0
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 0
发表于 2014-11-20 09:23:13 | 显示全部楼层
lysx 发表于 2014-11-20 08:18
直接a=0不就完了!!!

是的呢···针对这个,确实可以这样,或者干脆不赋值也行。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-11-20 11:49:11 | 显示全部楼层
我不知道楼主的意思
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-3-9 21:49:50 | 显示全部楼层
言深深 发表于 2014-11-19 21:14
你的这个运算就显式就行了。
do i=1,m
do j=1,n

谢谢您的帮助,这是我自己想的题,就是把5*5的矩阵除a(2,3)以外的值都赋成0,a(2,3)赋成10。
就是想看看
   forall(i=1:5,j=1:5)
     a(i,j)=0
     end forall
用do的隐式循环能不能写?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-3-9 21:50:57 | 显示全部楼层
zcilly 发表于 2014-11-20 11:49
我不知道楼主的意思

我只是想看看
   forall(i=1:5,j=1:5)
     a(i,j)=0
     end forall
用do的隐式循环怎么写?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 0
发表于 2015-3-10 10:02:54 | 显示全部楼层
林明宇 发表于 2015-3-9 21:50
我只是想看看
   forall(i=1:5,j=1:5)
     a(i,j)=0

a=0;a(2,3)=10
就可以了,当然,非要循环的话,可以用显式循环,隐式循环报错。
do i=1,5
do j=1,5
if(i.eq.2.and.j.eq.3)then
a(i,j)=10
else
a(i,j)=0
endif
enddo
enddo
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表