1 subroutine dia_sizref() 2c 3c $Id$ 4c 5 implicit none 6c 7#include "dia_common.fh" 8#include "mafdecls.fh" 9#include "msgids.fh" 10c 11 character*1 cdum 12 integer i,nsm 13c 14 integer ibroad(10),nhist 15 logical lforces 16c 17 if(me.eq.0) then 18 open(unit=lfnref,file=filref(1:index(filref,' ')-1), 19 + status='old',err=9999) 20c 21 read(lfnref,1000) cdum 22 read(lfnref,1000) cdum 23 read(lfnref,1000) cdum 24 1000 format(a) 25 read(lfnref,1001) nhist,lforces 26 1001 format(32x,i5,4x,l1) 27 if(nhist.gt.0) then 28 do 1 i=1,nhist 29 read(lfnref,1000) cdum 30 1 continue 31 endif 32c 33 read(lfnref,1000) cdum 34c 35 read(lfnref,1000) cdum 36 read(lfnref,1000) cdum 37 read(lfnref,1000) cdum 38c 39 read(lfnref,1000) cdum 40 read(lfnref,1000) cdum 41c 42 read(lfnref,1002) nwm,nwa,nsm,nsai 43 1002 format(4i10) 44c 45 close(unit=lfnref) 46 endif 47c 48 msai=nsai 49 nwmi=nwm 50 mwmi=nwm 51 nwai=nwa 52 mwai=nwa 53c 54c broadcast options to all nodes 55c 56 if(np.gt.1) then 57 ibroad(1)=nsai 58 ibroad(2)=msai 59 ibroad(3)=nwmi 60 ibroad(4)=mwmi 61 ibroad(5)=nwai 62 ibroad(6)=mwai 63 call ga_brdcst(mag_d01,ibroad,6*ma_sizeof(mt_int,1,mt_byte),0) 64 nsai=ibroad(1) 65 msai=ibroad(2) 66 nwmi=ibroad(3) 67 mwmi=ibroad(4) 68 nwai=ibroad(5) 69 mwai=ibroad(6) 70 endif 71c 72 return 73c 74 9999 continue 75 call md_abort('Unable to open reference file',0) 76 return 77 end 78 79