| 123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302 | #include "f2c.h"#include "fio.h"#include "string.h"#ifndef NON_POSIX_STDIO#ifdef MSDOS#include "io.h"#else#include "unistd.h"	/* for access */#endif#endif#ifdef KR_headersextern char *malloc();#ifdef NON_ANSI_STDIOextern char *mktemp();#endifextern integer f_clos();#define Const /*nothing*/#else#define Const const#undef abs#undef min#undef max#include "stdlib.h"#ifdef __cplusplusextern "C" {#endifextern int f__canseek(FILE*);extern integer f_clos(cllist*);#endif#ifdef NON_ANSI_RW_MODESConst char *f__r_mode[2] = {"r", "r"};Const char *f__w_mode[4] = {"w", "w", "r+w", "r+w"};#elseConst char *f__r_mode[2] = {"rb", "r"};Const char *f__w_mode[4] = {"wb", "w", "r+b", "r+"};#endif static char f__buf0[400], *f__buf = f__buf0; int f__buflen = (int)sizeof(f__buf0); static void#ifdef KR_headersf__bufadj(n, c) int n, c;#elsef__bufadj(int n, int c)#endif{	unsigned int len;	char *nbuf, *s, *t, *te;	if (f__buf == f__buf0)		f__buflen = 1024;	while(f__buflen <= n)		f__buflen <<= 1;	len = (unsigned int)f__buflen;	if (len != f__buflen || !(nbuf = (char*)malloc(len)))		f__fatal(113, "malloc failure");	s = nbuf;	t = f__buf;	te = t + c;	while(t < te)		*s++ = *t++;	if (f__buf != f__buf0)		free(f__buf);	f__buf = nbuf;	} int#ifdef KR_headersf__putbuf(c) int c;#elsef__putbuf(int c)#endif{	char *s, *se;	int n;	if (f__hiwater > f__recpos)		f__recpos = f__hiwater;	n = f__recpos + 1;	if (n >= f__buflen)		f__bufadj(n, f__recpos);	s = f__buf;	se = s + f__recpos;	if (c)		*se++ = c;	*se = 0;	for(;;) {		fputs(s, f__cf);		s += strlen(s);		if (s >= se)			break;	/* normally happens the first time */		putc(*s++, f__cf);		}	return 0;	} void#ifdef KR_headersx_putc(c)#elsex_putc(int c)#endif{	if (f__recpos >= f__buflen)		f__bufadj(f__recpos, f__buflen);	f__buf[f__recpos++] = c;	}#define opnerr(f,m,s) {if(f) errno= m; else opn_err(m,s,a); return(m);} static void#ifdef KR_headersopn_err(m, s, a) int m; char *s; olist *a;#elseopn_err(int m, const char *s, olist *a)#endif{	if (a->ofnm) {		/* supply file name to error message */		if (a->ofnmlen >= f__buflen)			f__bufadj((int)a->ofnmlen, 0);		g_char(a->ofnm, a->ofnmlen, f__curunit->ufnm = f__buf);		}	f__fatal(m, s);	}#ifdef KR_headersinteger f_open(a) olist *a;#elseinteger f_open(olist *a)#endif{	unit *b;	integer rv;	char buf[256], *s;	cllist x;	int ufmt;	FILE *tf;#ifndef NON_UNIX_STDIO	int n;#endif	f__external = 1;	if(a->ounit>=MXUNIT || a->ounit<0)		err(a->oerr,101,"open")	if (!f__init)		f_init();	f__curunit = b = &f__units[a->ounit];	if(b->ufd) {		if(a->ofnm==0)		{		same:	if (a->oblnk)				b->ublnk = *a->oblnk == 'z' || *a->oblnk == 'Z';			return(0);		}#ifdef NON_UNIX_STDIO		if (b->ufnm		 && strlen(b->ufnm) == a->ofnmlen		 && !strncmp(b->ufnm, a->ofnm, (unsigned)a->ofnmlen))			goto same;#else		g_char(a->ofnm,a->ofnmlen,buf);		if (f__inode(buf,&n) == b->uinode && n == b->udev)			goto same;#endif		x.cunit=a->ounit;		x.csta=0;		x.cerr=a->oerr;		if ((rv = f_clos(&x)) != 0)			return rv;		}	b->url = (int)a->orl;	b->ublnk = a->oblnk && (*a->oblnk == 'z' || *a->oblnk == 'Z');	if(a->ofm==0)	{	if(b->url>0) b->ufmt=0;		else b->ufmt=1;	}	else if(*a->ofm=='f' || *a->ofm == 'F') b->ufmt=1;	else b->ufmt=0;	ufmt = b->ufmt;#ifdef url_Adjust	if (b->url && !ufmt)		url_Adjust(b->url);#endif	if (a->ofnm) {		g_char(a->ofnm,a->ofnmlen,buf);		if (!buf[0])			opnerr(a->oerr,107,"open")		}	else		sprintf(buf, "fort.%ld", (long)a->ounit);	b->uscrtch = 0;	b->uend=0;	b->uwrt = 0;	b->ufd = 0;	b->urw = 3;	switch(a->osta ? *a->osta : 'u')	{	case 'o':	case 'O':#ifdef NON_POSIX_STDIO		if (!(tf = FOPEN(buf,"r")))			opnerr(a->oerr,errno,"open")		fclose(tf);#else		if (access(buf,0))			opnerr(a->oerr,errno,"open")#endif		break;	 case 's':	 case 'S':		b->uscrtch=1;#ifdef NON_ANSI_STDIO		(void) strcpy(buf,"tmp.FXXXXXX");		(void) mktemp(buf);		goto replace;#else		if (!(b->ufd = tmpfile()))			opnerr(a->oerr,errno,"open")		b->ufnm = 0;#ifndef NON_UNIX_STDIO		b->uinode = b->udev = -1;#endif		b->useek = 1;		return 0;#endif	case 'n':	case 'N':#ifdef NON_POSIX_STDIO		if ((tf = FOPEN(buf,"r")) || (tf = FOPEN(buf,"a"))) {			fclose(tf);			opnerr(a->oerr,128,"open")			}#else		if (!access(buf,0))			opnerr(a->oerr,128,"open")#endif		/* no break */	case 'r':	/* Fortran 90 replace option */	case 'R':#ifdef NON_ANSI_STDIO replace:#endif		if (tf = FOPEN(buf,f__w_mode[0]))			fclose(tf);	}	b->ufnm=(char *) malloc((unsigned int)(strlen(buf)+1));	if(b->ufnm==NULL) opnerr(a->oerr,113,"no space");	(void) strcpy(b->ufnm,buf);	if ((s = a->oacc) && b->url)		ufmt = 0;	if(!(tf = FOPEN(buf, f__w_mode[ufmt|2]))) {		if (tf = FOPEN(buf, f__r_mode[ufmt]))			b->urw = 1;		else if (tf = FOPEN(buf, f__w_mode[ufmt])) {			b->uwrt = 1;			b->urw = 2;			}		else			err(a->oerr, errno, "open");		}	b->useek = f__canseek(b->ufd = tf);#ifndef NON_UNIX_STDIO	if((b->uinode = f__inode(buf,&b->udev)) == -1)		opnerr(a->oerr,108,"open")#endif	if(b->useek)		if (a->orl)			rewind(b->ufd);		else if ((s = a->oacc) && (*s == 'a' || *s == 'A')			&& FSEEK(b->ufd, 0L, SEEK_END))				opnerr(a->oerr,129,"open");	return(0);} int#ifdef KR_headersfk_open(seq,fmt,n) ftnint n;#elsefk_open(int seq, int fmt, ftnint n)#endif{	char nbuf[10];	olist a;	(void) sprintf(nbuf,"fort.%ld",(long)n);	a.oerr=1;	a.ounit=n;	a.ofnm=nbuf;	a.ofnmlen=strlen(nbuf);	a.osta=NULL;	a.oacc= (char*)(seq==SEQ?"s":"d");	a.ofm = (char*)(fmt==FMT?"f":"u");	a.orl = seq==DIR?1:0;	a.oblnk=NULL;	return(f_open(&a));}#ifdef __cplusplus}#endif
 |