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