Fortran Coder

查看: 85|回复: 2

[求助] matlab代码转fortran代码

[复制链接]

2

帖子

1

主题

0

精华

新人

F 币
13 元
贡献
6 点
发表于 2019-9-27 19:34:41 | 显示全部楼层 |阅读模式
预定义一些变量的维度的matlab代码怎么写成fortran,matlab里面直接调用的zeros

         s_pf=zeros(layers,1); s_ps=zeros(layers,1);    s_sv=zeros(layers,1);     s_em=zeros(layers,1);  a_pf=zeros(layers,1);  a_ps=zeros(layers,1); a_sv=zeros(layers,1); a_em=zeros(layers,1);
               b_pf=zeros(layers,1); b_ps=zeros(layers,1);    b_sv=zeros(layers,1);     b_em=zeros(layers,1);
                K_w=zeros(layers,1);  L_w=zeros(layers,1); Sigma_w=zeros(layers,1);epsilon_w=zeros(layers,1); rou_w=zeros(layers,1); rou_t=zeros(layers,1);
               
            length_k=length(k_array);
                ur_0=zeros(1,length_k);  ut_0=zeros(1,length_k); uz_0=zeros(1,length_k);  
                Er_0=zeros(1,length_k);  Et_0=zeros(1,length_k); Ez_0=zeros(1,length_k);  
                Hr_0=zeros(1,length_k);  Ht_0=zeros(1,length_k); Hz_0=zeros(1,length_k);
               
                ur_1_1=zeros(1,length_k);     ur_1_2=zeros(1,length_k);     Er_1_1=zeros(1,length_k);    Er_1_2=zeros(1,length_k);   
                ut_1_1=zeros(1,length_k);     ut_1_2=zeros(1,length_k);     Et_1_1=zeros(1,length_k);    Et_1_2=zeros(1,length_k);  
                  uz_1=zeros(1,length_k);       Ez_1=zeros(1,length_k);                                          
                Hr_1_1=zeros(1,length_k);     Hr_1_2=zeros(1,length_k);   
                Ht_1_1=zeros(1,length_k);     Ht_1_2=zeros(1,length_k);   
                  Hz_1=zeros(1,length_k);     
               
                ur_1N_1=zeros(1,length_k);    ur_1N_2=zeros(1,length_k);    Er_1N_1=zeros(1,length_k);    Er_1N_2=zeros(1,length_k);  % 2*sin
                ut_1N_1=zeros(1,length_k);    ut_1N_2=zeros(1,length_k);    Et_1N_1=zeros(1,length_k);    Et_1N_2=zeros(1,length_k);  % 2*cos
                 uz_1N=zeros(1,length_k);        Ez_1N=zeros(1,length_k);                                        % 2*sin
                Hr_1N_1=zeros(1,length_k);    Hr_1N_2=zeros(1,length_k);  % 2*cos
                Ht_1N_1=zeros(1,length_k);    Ht_1N_2=zeros(1,length_k);  % 2*sin
                Hz_1N=zeros(1,length_k);        
                  
                  ur_2_1=zeros(1,length_k);    ur_2_2=zeros(1,length_k);    Er_2_1=zeros(1,length_k);    Er_2_2=zeros(1,length_k);  % 2*sin
                  ut_2_1=zeros(1,length_k);    ut_2_2=zeros(1,length_k);    Et_2_1=zeros(1,length_k);    Et_2_2=zeros(1,length_k);  % 2*cos
                   uz_2=zeros(1,length_k);       Ez_2=zeros(1,length_k);                                        % 2*sin
                  Hr_2_1=zeros(1,length_k);    Hr_2_2=zeros(1,length_k);  % 2*cos
                  Ht_2_1=zeros(1,length_k);    Ht_2_2=zeros(1,length_k);  % 2*sin
                  Hz_2=zeros(1,length_k);                                           % 2*cos         
                  
                           ur_2N_1=zeros(1,length_k);    ur_2N_2=zeros(1,length_k);   Er_2N_1=zeros(1,length_k);    Er_2N_2=zeros(1,length_k);  % 2*cos
                           ut_2N_1=zeros(1,length_k);    ut_2N_2=zeros(1,length_k);   Et_2N_1=zeros(1,length_k);    Et_2N_2=zeros(1,length_k);  % 2*sin
                             uz_2N=zeros(1,length_k);      Ez_2N=zeros(1,length_k);                                             % 2*cos
                           Hr_2N_1=zeros(1,length_k);     Hr_2N_2=zeros(1,length_k);  % 2*sin
                           Ht_2N_1=zeros(1,length_k);      Ht_2N_2=zeros(1,length_k);  % 2*cos
                             Hz_2N=zeros(1,length_k);   
                  
                  P_0=zeros(1,length_k);   P_1=zeros(1,length_k);  P_1N=zeros(1,length_k); P_2=zeros(1,length_k);  P_2N=zeros(1,length_k);
                  

回复

使用道具 举报

2

帖子

1

主题

0

精华

新人

F 币
13 元
贡献
6 点
 楼主| 发表于 2019-9-27 19:45:37 | 显示全部楼层
如何能够写一个循环或者类似的代码,将这种同类型的定义维度的内容进行编写

385

帖子

1

主题

0

精华

宗师

F 币
2030 元
贡献
1321 点
发表于 2019-9-28 16:20:32 | 显示全部楼层
matlab#       s_pf=zeros(layers,1)
fortran#       real::s_pf(layers,1)=0.0
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

QQ|捐赠本站|Archiver|关于我们 About Us|群聊|Fcode

GMT+8, 2019-10-21 03:59

Powered by Discuz! X3.2

© 2001-2017 Comsenz Inc.

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