C PROGRAM TRENDS JANUARY 1990 JRS PRIME FORTRAN 77 C DOUBLE PRECISION VARMAR, MARGM, WTXPX,BACKM,BLUES, DIST(28) DOUBLE PRECISION CMAX,CMIN,SUM,UMIN,DENOMIN INTEGER UNIT(50),IX,REPIN,NUMIN,NIN,inbird(1000) CHARACTER*2 HNTHOLD,HUNTST(50) CHARACTER*23 HEADER C C VARIABLES WITH THE SUFFIX 'NOW' REPRESENT THE CURRENT ROUTE. C THE OTHERS ARE FROM THE ROUTE MOST RECENTLY READ IN. C PARAMETER (MYRS=50,MCOLS=40) CHARACTER*5 OBS(MYRS), OBSNOW, OBSIN character*5 regi3,tmpfn*80,pid*12 character*4 caou CHARACTER*80 BLIST,OBSLIST,BBSUNIT, STALIST, STRLIST, COMMAND, & BBSOUT2,BBSZERO,SLOPFILE,STATLIST,STRALIST,REGILIST,tmplst, & INDEX,FNAME3,FNAMES,inlist2,indelist,inbirdf(1000),tmploes CHARACTER*80 TITLE,tmp*240,cmd*160 REAL STA, STANOW, STATE(50),AMAXIN,AMININ,AVEIN,ADENOM REAL YR,RT,STR,HRD,ANG,STATEE,STRAT REAL YRNOW,RTNOW,STRNOW,MIDYR,yr1 INTEGER*4 RNUM,ICNT INTEGER MFLAG, ALT, ALT2 INTEGER aou, IAOU, trim, inde integer*4 iseed,rand1 character*3 regi,regi2 character*100 fin,tmpindo character*25 name,regname character*26 sci fin='/home/httpd/htdocs/bbs/trend/scitab.95' tmpfn='/home/httpd/htdocs/bbs/trend/tmp/trendsel.' c print *,"htmtrend started" C c regi='ALA' c call getarg(1,caou) c print *,'tmp=',tmp read(caou,'(i4)') aou c call getarg(2,tmp) c print *,'tmp=',tmp c read(tmp,*) aou c call getarg(3,tmp) c print *,'tmp=',tmp c regi=tmp c print *,iyrf,aou,pid cx read(tmp,*) regi call getarg(2,tmp) c print *,'tmp=',tmp pid=tmp c print *, 'iyrf,aou,pid,regi=',iyrf,aou,pid,regi,'
' tmpfn(43:)=pid c print *,'
file:',tmpfn,'
' c print*, 'iyrf,aou,regi=',iyrf,aou,regi,'
' c write(20,*) 'aou is ',tmp, aou c c here, the regions are the state-stratum of interest cc call getarg(2,tmp) cc read(tmp,'(i4)')aou cc call getarg(3,regi) c c call getarg(4,tmp) c read(tmp,'(i4)') iyrl c call getarg(5,tmp) c read(tmp,'(i4)') linear c print *,aou,regi,iyrf,iyrl c c here is the read from the scitab file c c c open this if iyrf = 1 (an aou selection) c iyrf=1 if(iyrf.ne.1) go to 600 open(43,file=fin) mflag=0 10 read(43,12,end=400) iaou, name,sci 12 format(i4,7x,a24,a26) if(aou.ne.iaou) go to 10 mflag=1 400 continue close(43) 600 continue c c here is the regional option c if(iyrf.ne.2) go to 620 fin='/home/jrs/f/regnames.3' open(43,file=fin) mflag=0 610 read(43,612,end=404) regi2, name 612 format(a3,1x,a25) if(regi2.eq.'IA ') regi2=' IA' if(regi.ne.regi2) go to 610 mflag=1 404 continue close(43) 620 continue C open(22,file='/frank/jrs/f/t9999.trn') read(22,54,end=70) title(1:68) read(22,54,end=70) tmpindo(1:68) 54 format(4x,a68) c open(20,file=tmpfn) if(iyrf.ne.1) go to 630 if(mflag.eq.0) write(20,*)'No species name found for ',aouin if(mflag.eq.0) go to 299 write(20,40) name,sci 40 format('
') write(20,44) title(1:68) 44 format(27x,a68) write(20,46) tmpindo(1:68) 46 format(5x,'Region',13x,a68) 630 if(iyrf.ne.2) go to 640 write(20,434) name(1:23) 434 format('',a23,' Trend Results
') write(20,47) title(1:68) 47 format(25x,a68) write(20,48) tmpindo(1:68) 48 format(5x,'Species',13x,a68) 640 continue c c now read in the file c 50 read(22,55,end=70) inaou,regi2,name(1:23), ®name,title(1:69) 55 format(i4,a3,a23,a25,a69) c c now for species one iyrf=1 print out region c if(iyrf.ne.1 ) go to 66 if(inaou.ne.aou) go to 50 c if(inaou.lt.aou) go to 50 c if(inaou.gt.aou) go to 70 do 325 i=1,3 if(regi2(i:i).eq.'A') regi2(i:i)='a' if(regi2(i:i).eq.'B') regi2(i:i)='b' if(regi2(i:i).eq.'C') regi2(i:i)='c' if(regi2(i:i).eq.'D') regi2(i:i)='d' if(regi2(i:i).eq.'E') regi2(i:i)='e' if(regi2(i:i).eq.'F') regi2(i:i)='f' if(regi2(i:i).eq.'G') regi2(i:i)='g' if(regi2(i:i).eq.'H') regi2(i:i)='h' if(regi2(i:i).eq.'I') regi2(i:i)='i' if(regi2(i:i).eq.'J') regi2(i:i)='j' if(regi2(i:i).eq.'K') regi2(i:i)='k' if(regi2(i:i).eq.'L') regi2(i:i)='l' if(regi2(i:i).eq.'M') regi2(i:i)='m' if(regi2(i:i).eq.'N') regi2(i:i)='n' if(regi2(i:i).eq.'O') regi2(i:i)='o' if(regi2(i:i).eq.'P') regi2(i:i)='p' if(regi2(i:i).eq.'Q') regi2(i:i)='q' if(regi2(i:i).eq.'R') regi2(i:i)='r' if(regi2(i:i).eq.'S') regi2(i:i)='s' if(regi2(i:i).eq.'T') regi2(i:i)='t' if(regi2(i:i).eq.'U') regi2(i:i)='u' if(regi2(i:i).eq.'V') regi2(i:i)='v' if(regi2(i:i).eq.'W') regi2(i:i)='w' if(regi2(i:i).eq.'X') regi2(i:i)='x' if(regi2(i:i).eq.'Y') regi2(i:i)='y' if(regi2(i:i).eq.'Z') regi2(i:i)='z' 325 continue if(regi2(1:3).eq.'1ea') regi2(1:3)='ea ' if(regi2(1:3).eq.'1ce') regi2(1:3)='ce ' if(regi2(1:3).eq.'1we') regi2(1:3)='we ' if(regi2(1:3).eq.'1us') regi2(1:3)='us ' if(regi2(1:3).eq.'1ca') regi2(1:3)='ca ' if(regi2(1:3).eq.'1co') regi2(1:3)='sur' if(regi2(3:3).eq.' ') write(20,700) &inaou,regi2(1:2),regname,title(1:69) 700 format('' &,a23,'',1x,a69) if(regi2(3:3).ne.' ') write(20,701) &inaou,regi2,regname,title(1:69) 701 format('' &,a23,'',1x,a69) go to 50 66 if(iyrf.ne.2 ) go to 50 if(regi2.ne.regi) go to 50 write(20,710) name(1:23),title(1:69) 710 format(a23,1x,a69) go to 50 70 continue 299 continue close(20) 667 STOP END