如何将一个常量字符串从fortran函数传递给C

时间:2016-07-05 18:20:27

标签: c fortran filepath

我很难将字符串从Fortran传递给C.C代码将包含文件路径的字符串传递给Fortran函数。该函数读入文件并返回一个两个字符的字符串。

这是我的C代码

#include <stdlib.h>
#include <stdlib.h>
#include <stdio.h>
#include <string.h>
#include <math.h>
#include <sys/time.h>
#include <iostream>
#include <iomanip>
#include <fstream>
#include <sstream>
#include <time.h>
#include <ctype.h>
using namespace std;

extern "C" void read_id_type_(char *file_path, char *id_type, size_t *path_len_file);

int main(int argc, char **argv) {

char path[500],id_type[2];
char *dir_path, *file_path;
dir_path=path;
size_t path_len=strlen(dir_path);

file_path=strcat(path,"/rnspar_mpt1.dat");
size_t path_len_file=strlen(file_path);

read_id_type_(file_path,id_type,&path_len_file);
 printf("id_type is %s\n",id_type);

 return EXIT_SUCCESS ;
 }

这是我的Fortran功能

 integer function read_id_type(filename,id_type,path_len) 
 implicit none

 integer :: nrg, nrf, nrf_deform, nrgin,path_len
 character(400) :: filename
 character(2) ::  id_type,EQ_point

 filename=filename(1:path_len)
 write(*,*),"Reading file", filename

 open(unit = 1,file = trim(filename), status='old')
 read(1,'(4i5)') nrg
 read(1,'(4i5)') nrf
 read(1,'(2i5,2(3x,a2))') nrf_deform, nrgin, id_type, EQ_point
 close(1)
 print *,"Id_type=",id_type
 read_id_type = 0
 end function read_id_type

输出结果为:

Id_type=IR    (From the Fortran side)

id_type is IRV2?      (From the C side)

C代码中的输出应该只是前两个字符。 任何帮助将不胜感激

1 个答案:

答案 0 :(得分:0)

而不是

id_type[2]

至少

id_type[3]

并在调用Fortran例程后将3 rd 字符(在索引2处)设置为整数0,否则它将不是C理解的以零结尾的字符串。

由于当前代码尝试输出一个非零终止的字符串,因此它访问数组之外​​的内存并且具有正式的未定义行为。

UB什么都不是,或者可能发生什么。

你只是有一些奇怪的额外角色。

字符串在Fortran中正确显示,因为你已在那里指定了它的长度,

character(2) ::  id_type

C中没有相应的功能,但在C ++中,您可以使用长度为2的std::string并将其缓冲区传递给Fortran例程,

std::string id_type( 2, ' ' );

//...
read_id_type_( file_path, &id_type[0], &path_len_file );