options memory=3000; options limerr=50000; options limwarn=50000; options crt; name KON; freq none; smpl 1 7192; pertes=0; read(file='estdatareg8JUmars2012ut.asc',format=binary) regno regnony contr categ mainarea ratef depth conlength leadtime build year month loep riggser9 oilfl0 oilfl1-oilfl143 gassprl0 gassprl1-gassprl143 dumposnull; contes=conlength; deftes=; ldepth=log(depth); year1986=0; year1987=0; year1988=0; year1989=0; year1990=0; year1991=0; year1992=0; year1993=0; year1994=0; year1995=0; year1996=0; year1997=0; year1998=0; year1999=0; year2000=0; year2001=0; year2002=0; year2003=0; year2004=0; year2005=0; year2006=0; year2007=0; year2008=0; year2009=0; smplif year=1986; year1986=1; smpl 1 7192; smplif year=1987; year1987=1; smpl 1 7192; smplif year=1988; year1988=1; smpl 1 7192; smplif year=1989; year1989=1; smpl 1 7192; smplif year=1990; year1990=1; smpl 1 7192; smplif year=1991; year1991=1; smpl 1 7192; smplif year=1992; year1992=1; smpl 1 7192; smplif year=1993; year1993=1; smpl 1 7192; smplif year=1994; year1994=1; smpl 1 7192; smplif year=1995; year1995=1; smpl 1 7192; smplif year=1996; year1996=1; smpl 1 7192; smplif year=1997; year1997=1; smpl 1 7192; smplif year=1998; year1998=1; smpl 1 7192; smplif year=1999; year1999=1; smpl 1 7192; smplif year=2000; year2000=1; smpl 1 7192; smplif year=2001; year2001=1; smpl 1 7192; smplif year=2002; year2002=1; smpl 1 7192; smplif year=2003; year2003=1; smpl 1 7192; smplif year=2004; year2004=1; smpl 1 7192; smplif year=2005; year2005=1; smpl 1 7192; smplif year=2006; year2006=1; smpl 1 7192; smplif year=2007; year2007=1; smpl 1 7192; smplif year=2008; year2008=1; smpl 1 7192; smplif year=2009; year2009=1; smpl 1 7192; read(file='../ser9div.asc',NCOL=2,NROW=299) rigg9; do k=14 to 298 by 1; smplif loep=k; riggser9=rigg9[k,1]; riggser9l1=rigg9[k,2]; smpl 1 7192; enddo; riggser9l1tr=log(1-riggser9l1); highutil=0; smplif riggser9l1>=0.98; riggser9l1tr=0; highutil=1; smpl 1 7192; msd loep; dregnony=regnony-regnony(-1); dcateg=categ-categ(-1); smplif dregnony=0 & dcateg^=0; smpl 1 7192; smplif regnony=22; print categ year month loep; smpl 1 7192; msd loep; ?print regnony ratef year month loep; conlength1=conlength; conlength=conlength/100; msd categ leadtime; smpl 1 7192; set test1=1; print test1; print @nob; dumcat12=0; dumcat13=0; dumcat14=0; dumcat15=0; smplif categ=12; msd categ; print @nob; dumcat12=1; smpl 1 7192; smplif categ=13; msd categ; print @nob; dumcat13=1; smpl 1 7192; smplif categ=14; msd categ; print @nob; dumcat14=1; smpl 1 7192; smplif categ=15; msd categ; print @nob; dumcat15=1; smpl 1 7192; trend t1; lratef=log(ratef); loilf0=log(oilfl0); loilf1=log(oilfl1); loilf2=log(oilfl2); loilf3=log(oilfl3); loilf4=log(oilfl4); loilf5=log(oilfl5); loilf6=log(oilfl6); loilf7=log(oilfl7); loilf8=log(oilfl8); loilf9=log(oilfl9); loilf10=log(oilfl10); loilf11=log(oilfl11); loilf12=log(oilfl12); loilf13=log(oilfl13); loilf14=log(oilfl14); loilf15=log(oilfl15); loilf16=log(oilfl16); loilf17=log(oilfl17); loilf18=log(oilfl18); loilf19=log(oilfl19); loilf20=log(oilfl20); loilf21=log(oilfl21); loilf22=log(oilfl22); loilf23=log(oilfl23); loilf24=log(oilfl24); loilf25=log(oilfl25); loilf26=log(oilfl26); loilf27=log(oilfl27); loilf28=log(oilfl28); loilf29=log(oilfl29); loilf30=log(oilfl30); loilf31=log(oilfl31); loilf32=log(oilfl32); loilf33=log(oilfl33); loilf34=log(oilfl34); loilf35=log(oilfl35); loilf36=log(oilfl36); loilf37=log(oilfl37); loilf38=log(oilfl38); loilf39=log(oilfl39); loilf40=log(oilfl40); loilf41=log(oilfl41); loilf42=log(oilfl42); loilf43=log(oilfl43); loilf44=log(oilfl44); loilf45=log(oilfl45); loilf46=log(oilfl46); loilf47=log(oilfl47); loilf48=log(oilfl48); loilf49=log(oilfl49); loilf50=log(oilfl50); loilf51=log(oilfl51); loilf52=log(oilfl52); loilf53=log(oilfl53); loilf54=log(oilfl54); loilf55=log(oilfl55); loilf56=log(oilfl56); loilf57=log(oilfl57); loilf58=log(oilfl58); loilf59=log(oilfl59); loilf60=log(oilfl60); loilf61=log(oilfl61); loilf62=log(oilfl62); loilf63=log(oilfl63); loilf64=log(oilfl64); loilf65=log(oilfl65); loilf66=log(oilfl66); loilf67=log(oilfl67); loilf68=log(oilfl68); loilf69=log(oilfl69); loilf70=log(oilfl70); loilf71=log(oilfl71); loilf72=log(oilfl72); loilf73=log(oilfl73); loilf74=log(oilfl74); loilf75=log(oilfl75); loilf76=log(oilfl76); loilf77=log(oilfl77); loilf78=log(oilfl78); loilf79=log(oilfl79); loilf80=log(oilfl80); loilf81=log(oilfl81); loilf82=log(oilfl82); loilf83=log(oilfl83); loilf84=log(oilfl84); loilf85=log(oilfl85); loilf86=log(oilfl86); loilf87=log(oilfl87); loilf88=log(oilfl88); loilf89=log(oilfl89); loilf90=log(oilfl90); loilf91=log(oilfl91); loilf92=log(oilfl92); loilf93=log(oilfl93); loilf94=log(oilfl94); loilf95=log(oilfl95); loilf96=log(oilfl96); loilf97=log(oilfl97); loilf98=log(oilfl98); loilf99=log(oilfl99); loilf100=log(oilfl100); loilf101=log(oilfl101); loilf102=log(oilfl102); loilf103=log(oilfl103); loilf104=log(oilfl104); loilf105=log(oilfl105); loilf106=log(oilfl106); loilf107=log(oilfl107); loilf108=log(oilfl108); loilf109=log(oilfl109); loilf110=log(oilfl110); loilf111=log(oilfl111); loilf112=log(oilfl112); loilf113=log(oilfl113); loilf114=log(oilfl114); loilf115=log(oilfl115); loilf116=log(oilfl116); loilf117=log(oilfl117); loilf118=log(oilfl118); loilf119=log(oilfl119); loilf120=log(oilfl120); loilf121=log(oilfl121); loilf122=log(oilfl122); loilf123=log(oilfl123); loilf124=log(oilfl124); loilf125=log(oilfl125); loilf126=log(oilfl126); loilf127=log(oilfl127); loilf128=log(oilfl128); loilf129=log(oilfl129); loilf130=log(oilfl130); loilf131=log(oilfl131); loilf132=log(oilfl132); loilf133=log(oilfl133); loilf134=log(oilfl134); loilf135=log(oilfl135); loilf136=log(oilfl136); loilf137=log(oilfl137); loilf138=log(oilfl138); loilf139=log(oilfl139); loilf140=log(oilfl140); loilf141=log(oilfl141); loilf142=log(oilfl142); loilf143=log(oilfl143); lgassprl0=log(gassprl0); lgassprl1=log(gassprl1); lgassprl2=log(gassprl2); lgassprl3=log(gassprl3); lgassprl4=log(gassprl4); lgassprl5=log(gassprl5); lgassprl6=log(gassprl6); lgassprl7=log(gassprl7); lgassprl8=log(gassprl8); lgassprl9=log(gassprl9); lgassprl10=log(gassprl10); lgassprl11=log(gassprl11); lgassprl12=log(gassprl12); lgassprl13=log(gassprl13); lgassprl14=log(gassprl14); lgassprl15=log(gassprl15); lgassprl16=log(gassprl16); lgassprl17=log(gassprl17); lgassprl18=log(gassprl18); lgassprl19=log(gassprl19); lgassprl20=log(gassprl20); lgassprl21=log(gassprl21); lgassprl22=log(gassprl22); lgassprl23=log(gassprl23); lgassprl24=log(gassprl24); lgassprl25=log(gassprl25); lgassprl26=log(gassprl26); lgassprl27=log(gassprl27); lgassprl28=log(gassprl28); lgassprl29=log(gassprl29); lgassprl30=log(gassprl30); lgassprl31=log(gassprl31); lgassprl32=log(gassprl32); lgassprl33=log(gassprl33); lgassprl34=log(gassprl34); lgassprl35=log(gassprl35); lgassprl36=log(gassprl36); lgassprl37=log(gassprl37); lgassprl38=log(gassprl38); lgassprl39=log(gassprl39); lgassprl40=log(gassprl40); lgassprl41=log(gassprl41); lgassprl42=log(gassprl42); lgassprl43=log(gassprl43); lgassprl44=log(gassprl44); lgassprl45=log(gassprl45); lgassprl46=log(gassprl46); lgassprl47=log(gassprl47); lgassprl48=log(gassprl48); lgassprl49=log(gassprl49); lgassprl50=log(gassprl50); lgassprl51=log(gassprl51); lgassprl52=log(gassprl52); lgassprl53=log(gassprl53); lgassprl54=log(gassprl54); lgassprl55=log(gassprl55); lgassprl56=log(gassprl56); lgassprl57=log(gassprl57); lgassprl58=log(gassprl58); lgassprl59=log(gassprl59); lgassprl60=log(gassprl60); lgassprl61=log(gassprl61); lgassprl62=log(gassprl62); lgassprl63=log(gassprl63); lgassprl64=log(gassprl64); lgassprl65=log(gassprl65); lgassprl66=log(gassprl66); lgassprl67=log(gassprl67); lgassprl68=log(gassprl68); lgassprl69=log(gassprl69); lgassprl70=log(gassprl70); lgassprl71=log(gassprl71); lgassprl72=log(gassprl72); lgassprl73=log(gassprl73); lgassprl74=log(gassprl74); lgassprl75=log(gassprl75); lgassprl76=log(gassprl76); lgassprl77=log(gassprl77); lgassprl78=log(gassprl78); lgassprl79=log(gassprl79); lgassprl80=log(gassprl80); lgassprl81=log(gassprl81); lgassprl82=log(gassprl82); lgassprl83=log(gassprl83); lgassprl84=log(gassprl84); lgassprl85=log(gassprl85); lgassprl86=log(gassprl86); lgassprl87=log(gassprl87); lgassprl88=log(gassprl88); lgassprl89=log(gassprl89); lgassprl90=log(gassprl90); lgassprl91=log(gassprl91); lgassprl92=log(gassprl92); lgassprl93=log(gassprl93); lgassprl94=log(gassprl94); lgassprl95=log(gassprl95); lgassprl96=log(gassprl96); lgassprl97=log(gassprl97); lgassprl98=log(gassprl98); lgassprl99=log(gassprl99); lgassprl100=log(gassprl100); lgassprl101=log(gassprl101); lgassprl102=log(gassprl102); lgassprl103=log(gassprl103); lgassprl104=log(gassprl104); lgassprl105=log(gassprl105); lgassprl106=log(gassprl106); lgassprl107=log(gassprl107); lgassprl108=log(gassprl108); lgassprl109=log(gassprl109); lgassprl110=log(gassprl110); lgassprl111=log(gassprl111); lgassprl112=log(gassprl112); lgassprl113=log(gassprl113); lgassprl114=log(gassprl114); lgassprl115=log(gassprl115); lgassprl116=log(gassprl116); lgassprl117=log(gassprl117); lgassprl118=log(gassprl118); lgassprl119=log(gassprl119); lgassprl120=log(gassprl120); lgassprl121=log(gassprl121); lgassprl122=log(gassprl122); lgassprl123=log(gassprl123); lgassprl124=log(gassprl124); lgassprl125=log(gassprl125); lgassprl126=log(gassprl126); lgassprl127=log(gassprl127); lgassprl128=log(gassprl128); lgassprl129=log(gassprl129); lgassprl130=log(gassprl130); lgassprl131=log(gassprl131); lgassprl132=log(gassprl132); lgassprl133=log(gassprl133); lgassprl134=log(gassprl134); lgassprl135=log(gassprl135); lgassprl136=log(gassprl136); lgassprl137=log(gassprl137); lgassprl138=log(gassprl138); lgassprl139=log(gassprl139); lgassprl140=log(gassprl140); lgassprl141=log(gassprl141); lgassprl142=log(gassprl142); lgassprl143=log(gassprl143); msd mainarea; smplif regno^=123 & regno^=220 & regno^=242 & regno^=245 & regno^=321 & regno^=323 & regno^=523 & regno^=746 & regno^=791; msd dumcat12 dumcat13 dumcat14 dumcat15; msd depth leadtime build; leadtime=leadtime/1000; build=build/100; pertes=1; msd pertes; mform(NROW=2,NCOL=1) ll=0; smplif year>=1990; msd contes; msd year1999; mmake dt regno regnony contr categ mainarea lratef ldepth conlength leadtime build year month loep riggser9l1tr highutil oilfl0 oilfl1-oilfl143 gassprl0 gassprl1-gassprl143 year1991-year1999 year2000-year2009 dumcat13 dumcat14 dumcat15; deftes=1; write(file='dataall.asc',format='(325F20.6)') dt; smplif year=1990 & month=2; smpl 144; print gassprl0 gassprl1-gassprl143; print oilfl0 oilfl1-oilfl143; smpl 1 7192; smplif deftes=1; DD1=0; DD2=0; DD3=0; DD4=0; DD5=0; DD6=0; DD7=0; DD9=0; DD10=0; DD11=0; DD12=0; DD13=0; DD14=0; DD15=0; DD16=0; DD17=0; DD18=0; DD19=0; DD20=0; DD21=0; DD22=0; DD23=0; DD24=0; DD25=0; DD26=0; DD27=0; DD29=0; DD30=0; DD31=0; DD32=0; DD35=0; DD36=0; DD37=0; DD38=0; DD39=0; DD40=0; DD41=0; DD42=0; DD43=0; DD44=0; DD45=0; DD46=0; DD47=0; DD48=0; DD49=0; DD50=0; DD52=0; DD54=0; DD55=0; DD56=0; DD57=0; DD58=0; DD59=0; DD60=0; DD61=0; DD62=0; DD63=0; DD64=0; DD65=0; DD66=0; DD67=0; DD68=0; DD69=0; DD70=0; DD71=0; DD72=0; DD73=0; DD74=0; DD75=0; DD76=0; DD77=0; DD78=0; DD79=0; DD80=0; DD81=0; DD82=0; DD83=0; DD84=0; DD85=0; DD86=0; DD87=0; DD88=0; DD89=0; DD90=0; DD91=0; DD92=0; DD93=0; DD94=0; DD95=0; DD96=0; DD97=0; DD98=0; DD99=0; DD100=0; DD101=0; DD102=0; DD103=0; DD104=0; DD105=0; DD106=0; DD107=0; DD108=0; DD109=0; DD110=0; DD111=0; DD112=0; DD114=0; DD115=0; DD116=0; DD117=0; DD118=0; DD119=0; DD120=0; DD121=0; DD122=0; DD123=0; DD124=0; DD125=0; DD126=0; DD127=0; DD128=0; DD129=0; DD130=0; DD131=0; DD132=0; DD133=0; DD134=0; DD135=0; DD136=0; DD137=0; DD138=0; DD139=0; DD140=0; DD141=0; DD142=0; DD143=0; DD144=0; DD145=0; DD146=0; DD147=0; DD148=0; DD149=0; DD150=0; DD151=0; DD152=0; DD153=0; DD154=0; DD155=0; DD156=0; DD157=0; DD158=0; DD159=0; DD160=0; DD161=0; DD162=0; DD163=0; DD164=0; DD165=0; DD166=0; DD168=0; DD169=0; DD170=0; DD171=0; DD172=0; DD173=0; DD174=0; DD176=0; DD177=0; DD178=0; DD179=0; DD180=0; DD181=0; DD182=0; DD183=0; DD184=0; DD185=0; DD186=0; DD187=0; DD188=0; DD189=0; DD190=0; DD191=0; DD192=0; DD193=0; DD194=0; DD195=0; DD196=0; DD197=0; DD198=0; DD199=0; DD200=0; DD201=0; DD202=0; DD203=0; DD204=0; DD205=0; DD206=0; DD207=0; DD208=0; DD209=0; DD210=0; DD211=0; DD212=0; DD213=0; smplif regnony=1; DD1=1; smpl 1 7192; smplif deftes=1; smplif regnony=2; DD2=1; smpl 1 7192; smplif deftes=1; smplif regnony=3; DD3=1; smpl 1 7192; smplif deftes=1; smplif regnony=4; DD4=1; smpl 1 7192; smplif deftes=1; smplif regnony=5; DD5=1; smpl 1 7192; smplif deftes=1; smplif regnony=6; DD6=1; smpl 1 7192; smplif deftes=1; smplif regnony=7; DD7=1; smpl 1 7192; smplif deftes=1; smplif regnony=9; DD9=1; smpl 1 7192; smplif deftes=1; smplif regnony=10; DD10=1; smpl 1 7192; smplif deftes=1; smplif regnony=11; DD11=1; smpl 1 7192; smplif deftes=1; smplif regnony=12; DD12=1; smpl 1 7192; smplif deftes=1; smplif regnony=13; DD13=1; smpl 1 7192; smplif deftes=1; smplif regnony=14; DD14=1; smpl 1 7192; smplif deftes=1; smplif regnony=15; DD15=1; smpl 1 7192; smplif deftes=1; smplif regnony=16; DD16=1; smpl 1 7192; smplif deftes=1; smplif regnony=17; DD17=1; smpl 1 7192; smplif deftes=1; smplif regnony=18; DD18=1; smpl 1 7192; smplif deftes=1; smplif regnony=19; DD19=1; smpl 1 7192; smplif deftes=1; smplif regnony=20; DD20=1; smpl 1 7192; smplif deftes=1; smplif regnony=21; DD21=1; smpl 1 7192; smplif deftes=1; smplif regnony=22; DD22=1; smpl 1 7192; smplif deftes=1; smplif regnony=23; DD23=1; smpl 1 7192; smplif deftes=1; smplif regnony=24; DD24=1; smpl 1 7192; smplif deftes=1; smplif regnony=25; DD25=1; smpl 1 7192; smplif deftes=1; smplif regnony=26; DD26=1; smpl 1 7192; smplif deftes=1; smplif regnony=27; DD27=1; smpl 1 7192; smplif deftes=1; smplif regnony=29; DD29=1; smpl 1 7192; smplif deftes=1; smplif regnony=30; DD30=1; smpl 1 7192; smplif deftes=1; smplif regnony=31; DD31=1; smpl 1 7192; smplif deftes=1; smplif regnony=32; DD32=1; smpl 1 7192; smplif deftes=1; smplif regnony=35; DD35=1; smpl 1 7192; smplif deftes=1; smplif regnony=36; DD36=1; smpl 1 7192; smplif deftes=1; smplif regnony=37; DD37=1; smpl 1 7192; smplif deftes=1; smplif regnony=38; DD38=1; smpl 1 7192; smplif deftes=1; smplif regnony=39; DD39=1; smpl 1 7192; smplif deftes=1; smplif regnony=40; DD40=1; smpl 1 7192; smplif deftes=1; smplif regnony=41; DD41=1; smpl 1 7192; smplif deftes=1; smplif regnony=42; DD42=1; smpl 1 7192; smplif deftes=1; smplif regnony=43; DD43=1; smpl 1 7192; smplif deftes=1; smplif regnony=44; DD44=1; smpl 1 7192; smplif deftes=1; smplif regnony=45; DD45=1; smpl 1 7192; smplif deftes=1; smplif regnony=46; DD46=1; smpl 1 7192; smplif deftes=1; smplif regnony=47; DD47=1; smpl 1 7192; smplif deftes=1; smplif regnony=48; DD48=1; smpl 1 7192; smplif deftes=1; smplif regnony=49; DD49=1; smpl 1 7192; smplif deftes=1; smplif regnony=50; DD50=1; smpl 1 7192; smplif deftes=1; smplif regnony=52; DD52=1; smpl 1 7192; smplif deftes=1; smplif regnony=54; DD54=1; smpl 1 7192; smplif deftes=1; smplif regnony=55; DD55=1; smpl 1 7192; smplif deftes=1; smplif regnony=56; DD56=1; smpl 1 7192; smplif deftes=1; smplif regnony=57; DD57=1; smpl 1 7192; smplif deftes=1; smplif regnony=58; DD58=1; smpl 1 7192; smplif deftes=1; smplif regnony=59; DD59=1; smpl 1 7192; smplif deftes=1; smplif regnony=60; DD60=1; smpl 1 7192; smplif deftes=1; smplif regnony=61; DD61=1; smpl 1 7192; smplif deftes=1; smplif regnony=62; DD62=1; smpl 1 7192; smplif deftes=1; smplif regnony=63; DD63=1; smpl 1 7192; smplif deftes=1; smplif regnony=64; DD64=1; smpl 1 7192; smplif deftes=1; smplif regnony=65; DD65=1; smpl 1 7192; smplif deftes=1; smplif regnony=66; DD66=1; smpl 1 7192; smplif deftes=1; smplif regnony=67; DD67=1; smpl 1 7192; smplif deftes=1; smplif regnony=68; DD68=1; smpl 1 7192; smplif deftes=1; smplif regnony=69; DD69=1; smpl 1 7192; smplif deftes=1; smplif regnony=70; DD70=1; smpl 1 7192; smplif deftes=1; smplif regnony=71; DD71=1; smpl 1 7192; smplif deftes=1; smplif regnony=72; DD72=1; smpl 1 7192; smplif deftes=1; smplif regnony=73; DD73=1; smpl 1 7192; smplif deftes=1; smplif regnony=74; DD74=1; smpl 1 7192; smplif deftes=1; smplif regnony=75; DD75=1; smpl 1 7192; smplif deftes=1; smplif regnony=76; DD76=1; smpl 1 7192; smplif deftes=1; smplif regnony=77 DD77=1; smpl 1 7192; smplif deftes=1; smplif regnony=78; DD78=1; smpl 1 7192; smplif deftes=1; smplif regnony=79; DD79=1; smpl 1 7192; smplif deftes=1; smplif regnony=80; DD80=1; smpl 1 7192; smplif deftes=1; smplif regnony=81; DD81=1; smpl 1 7192; smplif deftes=1; smplif regnony=82; DD82=1; smpl 1 7192; smplif deftes=1; smplif regnony=83; DD83=1; smpl 1 7192; smplif deftes=1; smplif regnony=84; DD84=1; smpl 1 7192; smplif deftes=1; smplif regnony=85; DD85=1; smpl 1 7192; smplif deftes=1; smplif regnony=86; DD86=1; smpl 1 7192; smplif deftes=1; smplif regnony=87; DD87=1; smpl 1 7192; smplif deftes=1; smplif regnony=88; DD88=1; smpl 1 7192; smplif deftes=1; smplif regnony=89; DD89=1; smpl 1 7192; smplif deftes=1; smplif regnony=90; DD90=1; smpl 1 7192; smplif deftes=1; smplif regnony=91; DD91=1; smpl 1 7192; smplif deftes=1; smplif regnony=92; DD92=1; smpl 1 7192; smplif deftes=1; smplif regnony=93; DD93=1; smpl 1 7192; smplif deftes=1; smplif regnony=94; DD94=1; smpl 1 7192; smplif deftes=1; smplif regnony=95; DD95=1; smpl 1 7192; smplif deftes=1; smplif regnony=96; DD96=1; smpl 1 7192; smplif deftes=1; smplif regnony=97; DD97=1; smpl 1 7192; smplif deftes=1; smplif regnony=98; DD98=1; smpl 1 7192; smplif deftes=1; smplif regnony=99; DD99=1; smpl 1 7192; smplif deftes=1; smplif regnony=100; DD100=1; smpl 1 7192; smplif deftes=1; smplif regnony=101; DD101=1; smpl 1 7192; smplif deftes=1; smplif regnony=102; DD102=1; smpl 1 7192; smplif deftes=1; smplif regnony=103; DD103=1; smpl 1 7192; smplif deftes=1; smplif regnony=104; DD104=1; smpl 1 7192; smplif deftes=1; smplif regnony=105; DD105=1; smpl 1 7192; smplif deftes=1; smplif regnony=106; DD106=1; smpl 1 7192; smplif deftes=1; smplif regnony=107; DD107=1; smpl 1 7192; smplif deftes=1; smplif regnony=108; DD108=1; smpl 1 7192; smplif deftes=1; smplif regnony=109; DD109=1; smpl 1 7192; smplif deftes=1; smplif regnony=110; DD110=1; smpl 1 7192; smplif deftes=1; smplif regnony=111; DD111=1; smpl 1 7192; smplif deftes=1; smplif regnony=112; DD112=1; smpl 1 7192; smplif deftes=1; smplif regnony=114; DD114=1; smpl 1 7192; smplif deftes=1; smplif regnony=115; DD115=1; smpl 1 7192; smplif deftes=1; smplif regnony=116; DD116=1; smpl 1 7192; smplif deftes=1; smplif regnony=117; DD117=1; smpl 1 7192; smplif deftes=1; smplif regnony=118; DD118=1; smpl 1 7192; smplif deftes=1; smplif regnony=119; DD119=1; smpl 1 7192; smplif deftes=1; smplif regnony=120; DD120=1; smpl 1 7192; smplif deftes=1; smplif regnony=121; DD121=1; smpl 1 7192; smplif deftes=1; smplif regnony=122; DD122=1; smpl 1 7192; smplif deftes=1; smplif regnony=123; DD123=1; smpl 1 7192; smplif deftes=1; smplif regnony=124; DD124=1; smpl 1 7192; smplif deftes=1; smplif regnony=125; DD125=1; smpl 1 7192; smplif deftes=1; smplif regnony=126; DD126=1; smpl 1 7192; smplif deftes=1; smplif regnony=127; DD127=1; smpl 1 7192; smplif deftes=1; smplif regnony=128; DD128=1; smpl 1 7192; smplif deftes=1; smplif regnony=129; DD129=1; smpl 1 7192; smplif deftes=1; smplif regnony=130; DD130=1; smpl 1 7192; smplif deftes=1; smplif regnony=131; DD131=1; smpl 1 7192; smplif deftes=1; smplif regnony=132; DD132=1; smpl 1 7192; smplif deftes=1; smplif regnony=133; DD133=1; smpl 1 7192; smplif deftes=1; smplif regnony=134; DD134=1; smpl 1 7192; smplif deftes=1; smplif regnony=135; DD135=1; smpl 1 7192; smplif deftes=1; smplif regnony=136; DD136=1; smpl 1 7192; smplif deftes=1; smplif regnony=137; DD137=1; smpl 1 7192; smplif deftes=1; smplif regnony=138; DD138=1; smpl 1 7192; smplif deftes=1; smplif regnony=139; DD139=1; smpl 1 7192; smplif deftes=1; smplif regnony=140; DD140=1; smpl 1 7192; smplif deftes=1; smplif regnony=141; DD141=1; smpl 1 7192; smplif deftes=1; smplif regnony=142; DD142=1; smpl 1 7192; smplif deftes=1; smplif regnony=143; DD143=1; smpl 1 7192; smplif deftes=1; smplif regnony=144; DD144=1; smpl 1 7192; smplif deftes=1; smplif regnony=145; DD145=1; smpl 1 7192; smplif deftes=1; smplif regnony=146; DD146=1; smpl 1 7192; smplif deftes=1; smplif regnony=147; DD147=1; smpl 1 7192; smplif deftes=1; smplif regnony=148; DD148=1; smpl 1 7192; smplif deftes=1; smplif regnony=149; DD149=1; smpl 1 7192; smplif deftes=1; smplif regnony=150; DD150=1; smpl 1 7192; smplif deftes=1; smplif regnony=151; DD151=1; smpl 1 7192; smplif deftes=1; smplif regnony=152; DD152=1; smpl 1 7192; smplif deftes=1; smplif regnony=153; DD153=1; smpl 1 7192; smplif deftes=1; smplif regnony=154; DD154=1; smpl 1 7192; smplif deftes=1; smplif regnony=155; DD155=1; smpl 1 7192; smplif deftes=1; smplif regnony=156; DD156=1; smpl 1 7192; smplif deftes=1; smplif regnony=157; DD157=1; smpl 1 7192; smplif deftes=1; smplif regnony=158; DD158=1; smpl 1 7192; smplif deftes=1; smplif regnony=159; DD159=1; smpl 1 7192; smplif deftes=1; smplif regnony=160; DD160=1; smpl 1 7192; smplif deftes=1; smplif regnony=161; DD161=1; smpl 1 7192; smplif deftes=1; smplif regnony=162; DD162=1; smpl 1 7192; smplif deftes=1; smplif regnony=163; DD163=1; smpl 1 7192; smplif deftes=1; smplif regnony=164; DD164=1; smpl 1 7192; smplif deftes=1; smplif regnony=165; DD165=1; smpl 1 7192; smplif deftes=1; smplif regnony=166; DD166=1; smpl 1 7192; smplif deftes=1; smplif regnony=168; DD168=1; smpl 1 7192; smplif deftes=1; smplif regnony=169; DD169=1; smpl 1 7192; smplif deftes=1; smplif regnony=170; DD170=1; smpl 1 7192; smplif deftes=1; smplif regnony=171; DD171=1; smpl 1 7192; smplif deftes=1; smplif regnony=172; DD172=1; smpl 1 7192; smplif deftes=1; smplif regnony=173; DD173=1; smpl 1 7192; smplif deftes=1; smplif regnony=174; DD174=1; smpl 1 7192; smplif deftes=1; smplif regnony=176; DD176=1; smpl 1 7192; smplif deftes=1; smplif regnony=177; DD177=1; smpl 1 7192; smplif deftes=1; smplif regnony=178; DD178=1; smpl 1 7192; smplif deftes=1; smplif regnony=179; DD179=1; smpl 1 7192; smplif deftes=1; smplif regnony=180; DD180=1; smpl 1 7192; smplif deftes=1; smplif regnony=181; DD181=1; smpl 1 7192; smplif deftes=1; smplif regnony=182; DD182=1; smpl 1 7192; smplif deftes=1; smplif regnony=183; DD183=1; smpl 1 7192; smplif deftes=1; smplif regnony=184; DD184=1; smpl 1 7192; smplif deftes=1; smplif regnony=185; DD185=1; smpl 1 7192; smplif deftes=1; smplif regnony=186; DD186=1; smpl 1 7192; smplif deftes=1; smplif regnony=187; DD187=1; smpl 1 7192; smplif deftes=1; smplif regnony=188; DD188=1; smpl 1 7192; smplif deftes=1; smplif regnony=189; DD189=1; smpl 1 7192; smplif deftes=1; smplif regnony=190; DD190=1; smpl 1 7192; smplif deftes=1; smplif regnony=191; DD191=1; smpl 1 7192; smplif deftes=1; smplif regnony=192; DD192=1; smpl 1 7192; smplif deftes=1; smplif regnony=193; DD193=1; smpl 1 7192; smplif deftes=1; smplif regnony=194; DD194=1; smpl 1 7192; smplif deftes=1; smplif regnony=195; DD195=1; smpl 1 7192; smplif deftes=1; smplif regnony=196; DD196=1; smpl 1 7192; smplif deftes=1; smplif regnony=197; DD197=1; smpl 1 7192; smplif deftes=1; smplif regnony=198; DD198=1; smpl 1 7192; smplif deftes=1; smplif regnony=199; DD199=1; smpl 1 7192; smplif deftes=1; smplif regnony=200; DD200=1; smpl 1 7192; smplif deftes=1; smplif regnony=201; DD201=1; smpl 1 7192; smplif deftes=1; smplif regnony=202; DD202=1; smpl 1 7192; smplif deftes=1; smplif regnony=203; DD203=1; smpl 1 7192; smplif deftes=1; smplif regnony=204; DD204=1; smpl 1 7192; smplif deftes=1; smplif regnony=205; DD205=1; smpl 1 7192; smplif deftes=1; smplif regnony=206; DD206=1; smpl 1 7192; smplif deftes=1; smplif regnony=207; DD207=1; smpl 1 7192; smplif deftes=1; smplif regnony=208; DD208=1; smpl 1 7192; smplif deftes=1; smplif regnony=209; DD209=1; smpl 1 7192; smplif deftes=1; smplif regnony=210; DD210=1; smpl 1 7192; smplif deftes=1; smplif regnony=211; DD211=1; smpl 1 7192; smplif deftes=1; smplif regnony=211; DD212=1; smpl 1 7192; smplif deftes=1; smplif regnony=213; DD213=1; smpl 1 7192; smplif deftes=1; smplif DD1=1 | DD2=1 | DD4=1 | DD5=1 | DD7=1 | DD9=1 | DD10=1 | DD11=1 | DD12=1 | DD13=1 | DD14=1 | DD15=1 | DD16=1 | DD17=1 | DD18=1 | DD19=1 | DD21=1 | DD23=1 | DD24=1 | DD25=1 | DD26=1 | DD27=1 | DD29=1 | DD30=1 | DD32=1 | DD36=1 | DD37=1 | DD38=1 | DD39=1 | DD40=1 | DD41=1 | DD42=1 | DD43=1 | DD44=1 | DD45=1 | DD46=1 | DD47=1 | DD48=1 | DD49=1 | DD50=1 | DD52=1 | DD54=1 | DD56=1 | DD57=1 | DD58=1 | DD59=1 | DD60=1 | DD61=1 | DD62=1 | DD64=1 | DD65=1 | DD66=1 | DD67=1 | DD68=1 | DD69=1 | DD70=1 | DD72=1 | DD73=1 | DD74=1 | DD75=1 | DD76=1 | DD78=1 | DD79=1 | DD80=1 | DD81=1 | DD83=1 | DD84=1 | DD85=1 | DD86=1 | DD87=1 | DD88=1 | DD89=1 | DD90=1 | DD91=1 | DD92=1 | DD93=1 | DD95=1 | DD96=1 | DD97=1 | DD98=1 | DD99=1 | DD100=1 | DD101=1 | DD103=1 | DD104=1 | DD105=1 | DD106=1 | DD107=1 | DD109=1 | DD110=1 | DD111=1 | DD112=1 | DD114=1 | DD115=1 | DD116=1 | DD117=1 | DD118=1 | DD119=1 | DD121=1 | DD122=1 | DD123=1 | DD124=1 | DD125=1 | DD126=1 | DD127=1 | DD128=1 | DD129=1 | DD131=1 | DD133=1 | DD134=1 | DD135=1 | DD136=1 | DD137=1 | DD138=1 | DD139=1 | DD140=1 | DD141=1 | DD142=1 | DD143=1 | DD144=1 | DD145=1 | DD146=1 | DD147=1 | DD148=1 | DD149=1 | DD150=1 | DD151=1 | DD152=1 | DD153=1 | DD154=1 | DD155=1 | DD156=1 | DD157=1 | DD158=1 | DD159=1 | DD160=1 | DD161=1 | DD162=1 | DD164=1 | DD166=1 | DD169=1 | DD170=1 | DD171=1 | DD172=1 | DD173=1 | DD176=1 | DD177=1 | DD178=1 | DD179=1 | DD180=1 | DD181=1 | DD182=1 | DD183=1 | DD184=1 | DD185=1 | DD186=1 | DD187=1 | DD188=1 | DD189=1 | DD190=1 | DD191=1 | DD192=1 | DD193=1 | DD194=1 | DD195=1 | DD196=1 | DD197=1 | DD198=1 | DD199=1 | DD200=1 | DD201=1 | DD202=1 | DD203=1 | DD205=1; ?msd dd1-dd7 dd9-dd27 dd29-dd32 dd36-dd50 dd52 dd54-dd112 dd114-dd166 dd168-dd174 dd176-dd213; mmake dtt regno regnony contr categ mainarea lratef ldepth conlength leadtime build year month loep riggser9l1tr highutil oilfl0 oilfl1-oilfl143 gassprl0 gassprl1-gassprl143 year1991-year1999 year2000-year2009 dumcat13 dumcat14 dumcat15 DD1 DD2 DD4 DD5 DD7 DD9 DD10 DD11 DD12 DD13 DD14 DD15 DD16 DD17 DD18 DD19 DD21 DD23 DD24 DD25 DD26 DD27 DD29 DD30 DD32 DD36 DD37 DD38 DD39 DD40 DD41 DD42 DD43 DD44 DD45 DD46 DD47 DD48 DD49 DD50 DD52 DD54 DD56 DD57 DD58 DD59 DD60 DD61 DD62 DD64 DD65 DD66 DD67 DD68 DD69 DD70 DD72 DD73 DD74 DD75 DD76 DD78 DD79 DD80 DD81 DD83 DD84 DD85 DD86 DD87 DD88 DD89 DD90 DD91 DD92 DD93 DD95 DD96 DD97 DD98 DD99 DD100 DD101 DD103 DD104 DD105 DD106 DD107 DD109 DD110 DD111 DD112 DD114 DD115 DD116 DD117 DD118 DD119 DD121 DD122 DD123 DD124 DD125 DD126 DD127 DD128 DD129 DD131 DD133 DD134 DD135 DD136 DD137 DD138 DD139 DD140 DD141 DD142 DD143 DD144 DD145 DD146 DD147 DD148 DD149 DD150 DD151 DD152 DD153 DD154 DD155 DD156 DD157 DD158 DD159 DD160 DD161 DD162 DD164 DD166 DD169 DD170 DD171 DD172 DD173 DD176 DD177 DD178 DD179 DD180 DD181 DD182 DD183 DD184 DD185 DD186 DD187 DD188 DD189 DD190 DD191 DD192 DD193 DD194 DD195 DD196 DD197 DD198 DD199 DD200 DD201 DD202 DD203 DD205; param ff1 ff2 ff4 ff5 ff7 ff9 ff10 ff11 ff12 ff13 ff14 ff15 ff16 ff17 ff18 ff19 ff21 ff23 ff24 ff25 ff26 ff27 ff29 ff30 ff32 ff36 ff37 ff38 ff39 ff40 ff41 ff42 ff43 ff44 ff45 ff46 ff47 ff48 ff49 ff50 ff52 ff54 ff56 ff57 ff58 ff59 ff60 ff61 ff62 ff64 ff65 ff66 ff67 ff68 ff69 ff70 ff72 ff73 ff74 ff75 ff76 ff78 ff79 ff80 ff81 ff83 ff84 ff85 ff86 ff87 ff88 ff89 ff90 ff91 ff92 ff93 ff95 ff96 ff97 ff98 ff99 ff100 ff101 ff103 ff104 ff105 ff106 ff107 ff108 ff110 ff111 ff112 ff114 ff115 ff116 ff117 ff118 ff119 ff121 ff122 ff123 ff124 ff125 ff126 ff127 ff128 ff129 ff131 ff133 ff134 ff135 ff136 ff137 ff138 ff139 ff140 ff141 ff142 ff143 ff144 ff145 ff146 ff147 ff148 ff149 ff150 ff151 ff152 ff153 ff154 ff155 ff156 ff157 ff158 ff159 ff160 ff161 ff162 ff164 ff166 ff169 ff170 ff171 ff172 ff173 ff176 ff177 ff178 ff179 ff180 ff181 ff182 ff183 ff184 ff185 ff186 ff187 ff188 ff189 ff190 ff191 ff192 ff193 ff194 ff195 ff196 ff197 ff198 ff199 ff200 ff201 ff202 ff203 ff205 b2 0.2254 b4 -0.0470 b5 0.2473 b6 -0.28 b7 -0.25 b8 -0.22 b9 -0.14 b10 0.1407 b11 0.2527 b12 0.2877 b13 0.1414 b14 0.03905 b15 -0.1709 b16 -0.41 b17 -0.4590 b18 -0.692 b19 -0.717 b20 -0.626 b21 -0.3741 b22 -0.57 b23 -0.69 b24 -0.85 b25 1.30204 s2e 0.05 g 0.06 d 0.8 s 5 m 0.016 ; frml eq11 lpcesg= g*gassprl0+ g*(1-g)*gassprl1+ g*(1-g)**2*gassprl2+ g*(1-g)**3*gassprl3+ g*(1-g)**4*gassprl4+ g*(1-g)**5*gassprl5+ g*(1-g)**6*gassprl6+ g*(1-g)**7*gassprl7+ g*(1-g)**8*gassprl8+ g*(1-g)**9*gassprl9 +g*(1-g)**10*gassprl10 +g*(1-g)**11*gassprl11 +g*(1-g)**12*gassprl12 +g*(1-g)**13*gassprl13 +g*(1-g)**14*gassprl14 +g*(1-g)**15*gassprl15 +g*(1-g)**16*gassprl16 +g*(1-g)**17*gassprl17 +g*(1-g)**18*gassprl18 +g*(1-g)**19*gassprl19 +g*(1-g)**20*gassprl20 +g*(1-g)**21*gassprl21 +g*(1-g)**22*gassprl22 +g*(1-g)**23*gassprl23 +g*(1-g)**24*gassprl24 +g*(1-g)**25*gassprl25 +g*(1-g)**26*gassprl26 +g*(1-g)**27*gassprl27 +g*(1-g)**28*gassprl28 +g*(1-g)**29*gassprl29 +g*(1-g)**30*gassprl30 +g*(1-g)**31*gassprl31 +g*(1-g)**32*gassprl32 +g*(1-g)**33*gassprl33 +g*(1-g)**34*gassprl34 +g*(1-g)**35*gassprl35 +g*(1-g)**36*gassprl36 +g*(1-g)**37*gassprl37 +g*(1-g)**38*gassprl38 +g*(1-g)**39*gassprl39 +g*(1-g)**40*gassprl40 +g*(1-g)**41*gassprl41 +g*(1-g)**42*gassprl42 +g*(1-g)**43*gassprl43 +g*(1-g)**44*gassprl44 +g*(1-g)**45*gassprl45 +g*(1-g)**46*gassprl46 +g*(1-g)**47*gassprl47 +g*(1-g)**48*gassprl48 +g*(1-g)**49*gassprl49 +g*(1-g)**50*gassprl50 +g*(1-g)**51*gassprl51 +g*(1-g)**52*gassprl52 +g*(1-g)**53*gassprl53 +g*(1-g)**54*gassprl54 +g*(1-g)**55*gassprl55 +g*(1-g)**56*gassprl56 +g*(1-g)**57*gassprl57 +g*(1-g)**58*gassprl58 +g*(1-g)**59*gassprl59 +g*(1-g)**60*gassprl60 +g*(1-g)**61*gassprl61 +g*(1-g)**62*gassprl62 +g*(1-g)**63*gassprl63 +g*(1-g)**64*gassprl64 +g*(1-g)**65*gassprl65 +g*(1-g)**66*gassprl66 +g*(1-g)**67*gassprl67 +g*(1-g)**68*gassprl68 +g*(1-g)**69*gassprl69 +g*(1-g)**70*gassprl70 +g*(1-g)**71*gassprl71 +g*(1-g)**72*gassprl72 +g*(1-g)**73*gassprl73 +g*(1-g)**74*gassprl74 +g*(1-g)**75*gassprl75 +g*(1-g)**76*gassprl76 +g*(1-g)**77*gassprl77 +g*(1-g)**78*gassprl78 +g*(1-g)**79*gassprl79 +g*(1-g)**80*gassprl80 +g*(1-g)**81*gassprl81 +g*(1-g)**82*gassprl82 +g*(1-g)**83*gassprl83 +g*(1-g)**84*gassprl84 +g*(1-g)**85*gassprl85 +g*(1-g)**86*gassprl86 +g*(1-g)**87*gassprl87 +g*(1-g)**88*gassprl88 +g*(1-g)**89*gassprl89 +g*(1-g)**90*gassprl90 +g*(1-g)**91*gassprl91 +g*(1-g)**92*gassprl92 +g*(1-g)**93*gassprl93 +g*(1-g)**94*gassprl94 +g*(1-g)**95*gassprl95 +g*(1-g)**96*gassprl96 +g*(1-g)**97*gassprl97 +g*(1-g)**98*gassprl98 +g*(1-g)**99*gassprl99 +g*(1-g)**100*gassprl100 +g*(1-g)**101*gassprl101 +g*(1-g)**102*gassprl102 +g*(1-g)**103*gassprl103 +g*(1-g)**104*gassprl104 +g*(1-g)**105*gassprl105 +g*(1-g)**106*gassprl106 +g*(1-g)**107*gassprl107 +g*(1-g)**108*gassprl108 +g*(1-g)**109*gassprl109 +g*(1-g)**110*gassprl110 +g*(1-g)**111*gassprl111 +g*(1-g)**112*gassprl112 +g*(1-g)**113*gassprl113 +g*(1-g)**114*gassprl114 +g*(1-g)**115*gassprl115 +g*(1-g)**116*gassprl116 +g*(1-g)**117*gassprl117 +g*(1-g)**118*gassprl118 +g*(1-g)**119*gassprl119 +g*(1-g)**120*gassprl120 +g*(1-g)**121*gassprl121 +g*(1-g)**122*gassprl122 +g*(1-g)**123*gassprl123 +g*(1-g)**124*gassprl124 +g*(1-g)**125*gassprl125 +g*(1-g)**126*gassprl126 +g*(1-g)**127*gassprl127 +g*(1-g)**128*gassprl128 +g*(1-g)**129*gassprl129 +g*(1-g)**130*gassprl130 +g*(1-g)**131*gassprl131 +g*(1-g)**132*gassprl132 +g*(1-g)**133*gassprl133 +g*(1-g)**134*gassprl134 +g*(1-g)**135*gassprl135 +g*(1-g)**136*gassprl136 +g*(1-g)**137*gassprl137 +g*(1-g)**138*gassprl138 +g*(1-g)**139*gassprl139 +g*(1-g)**140*gassprl140 +g*(1-g)**141*gassprl141 +g*(1-g)**142*gassprl142 +g*(1-g)**143*gassprl143 ; frml eq12 lpceso= g*oilfl0+ g*(1-g)*oilfl1+ g*(1-g)**2*oilfl2+ g*(1-g)**3*oilfl3+ g*(1-g)**4*oilfl4+ g*(1-g)**5*oilfl5+ g*(1-g)**6*oilfl6+ g*(1-g)**7*oilfl7+ g*(1-g)**8*oilfl8+ g*(1-g)**9*oilfl9 +g*(1-g)**10*oilfl10 +g*(1-g)**11*oilfl11 +g*(1-g)**12*oilfl12 +g*(1-g)**13*oilfl13 +g*(1-g)**14*oilfl14 +g*(1-g)**15*oilfl15 +g*(1-g)**16*oilfl16 +g*(1-g)**17*oilfl17 +g*(1-g)**18*oilfl18 +g*(1-g)**19*oilfl19 +g*(1-g)**20*oilfl20 +g*(1-g)**21*oilfl21 +g*(1-g)**22*oilfl22 +g*(1-g)**23*oilfl23 +g*(1-g)**24*oilfl24 +g*(1-g)**25*oilfl25 +g*(1-g)**26*oilfl26 +g*(1-g)**27*oilfl27 +g*(1-g)**28*oilfl28 +g*(1-g)**29*oilfl29 +g*(1-g)**30*oilfl30 +g*(1-g)**31*oilfl31 +g*(1-g)**32*oilfl32 +g*(1-g)**33*oilfl33 +g*(1-g)**34*oilfl34 +g*(1-g)**35*oilfl35 +g*(1-g)**36*oilfl36 +g*(1-g)**37*oilfl37 +g*(1-g)**38*oilfl38 +g*(1-g)**39*oilfl39 +g*(1-g)**40*oilfl40 +g*(1-g)**41*oilfl41 +g*(1-g)**42*oilfl42 +g*(1-g)**43*oilfl43 +g*(1-g)**44*oilfl44 +g*(1-g)**45*oilfl45 +g*(1-g)**46*oilfl46 +g*(1-g)**47*oilfl47 +g*(1-g)**48*oilfl48 +g*(1-g)**49*oilfl49 +g*(1-g)**50*oilfl50 +g*(1-g)**51*oilfl51 +g*(1-g)**52*oilfl52 +g*(1-g)**53*oilfl53 +g*(1-g)**54*oilfl54 +g*(1-g)**55*oilfl55 +g*(1-g)**56*oilfl56 +g*(1-g)**57*oilfl57 +g*(1-g)**58*oilfl58 +g*(1-g)**59*oilfl59 +g*(1-g)**60*oilfl60 +g*(1-g)**61*oilfl61 +g*(1-g)**62*oilfl62 +g*(1-g)**63*oilfl63 +g*(1-g)**64*oilfl64 +g*(1-g)**65*oilfl65 +g*(1-g)**66*oilfl66 +g*(1-g)**67*oilfl67 +g*(1-g)**68*oilfl68 +g*(1-g)**69*oilfl69 +g*(1-g)**70*oilfl70 +g*(1-g)**71*oilfl71 +g*(1-g)**72*oilfl72 +g*(1-g)**73*oilfl73 +g*(1-g)**74*oilfl74 +g*(1-g)**75*oilfl75 +g*(1-g)**76*oilfl76 +g*(1-g)**77*oilfl77 +g*(1-g)**78*oilfl78 +g*(1-g)**79*oilfl79 +g*(1-g)**80*oilfl80 +g*(1-g)**81*oilfl81 +g*(1-g)**82*oilfl82 +g*(1-g)**83*oilfl83 +g*(1-g)**84*oilfl84 +g*(1-g)**85*oilfl85 +g*(1-g)**86*oilfl86 +g*(1-g)**87*oilfl87 +g*(1-g)**88*oilfl88 +g*(1-g)**89*oilfl89 +g*(1-g)**90*oilfl90 +g*(1-g)**91*oilfl91 +g*(1-g)**92*oilfl92 +g*(1-g)**93*oilfl93 +g*(1-g)**94*oilfl94 +g*(1-g)**95*oilfl95 +g*(1-g)**96*oilfl96 +g*(1-g)**97*oilfl97 +g*(1-g)**98*oilfl98 +g*(1-g)**99*oilfl99 +g*(1-g)**100*oilfl100 +g*(1-g)**101*oilfl101 +g*(1-g)**102*oilfl102 +g*(1-g)**103*oilfl103 +g*(1-g)**104*oilfl104 +g*(1-g)**105*oilfl105 +g*(1-g)**106*oilfl106 +g*(1-g)**107*oilfl107 +g*(1-g)**108*oilfl108 +g*(1-g)**109*oilfl109 +g*(1-g)**110*oilfl110 +g*(1-g)**111*oilfl111 +g*(1-g)**112*oilfl112 +g*(1-g)**113*oilfl113 +g*(1-g)**114*oilfl114 +g*(1-g)**115*oilfl115 +g*(1-g)**116*oilfl116 +g*(1-g)**117*oilfl117 +g*(1-g)**118*oilfl118 +g*(1-g)**119*oilfl119 +g*(1-g)**120*oilfl120 +g*(1-g)**121*oilfl121 +g*(1-g)**122*oilfl122 +g*(1-g)**123*oilfl123 +g*(1-g)**124*oilfl124 +g*(1-g)**125*oilfl125 +g*(1-g)**126*oilfl126 +g*(1-g)**127*oilfl127 +g*(1-g)**128*oilfl128 +g*(1-g)**129*oilfl129 +g*(1-g)**130*oilfl130 +g*(1-g)**131*oilfl131 +g*(1-g)**132*oilfl132 +g*(1-g)**133*oilfl133 +g*(1-g)**134*oilfl134 +g*(1-g)**135*oilfl135 +g*(1-g)**136*oilfl136 +g*(1-g)**137*oilfl137 +g*(1-g)**138*oilfl138 +g*(1-g)**139*oilfl139 +g*(1-g)**140*oilfl140 +g*(1-g)**141*oilfl141 +g*(1-g)**142*oilfl142 +g*(1-g)**143*oilfl143 ; frml eq13 lpcesg1= 2.72964*g+ 3.29301*g*(1-g)+ 2.84685*g*(1-g)**2+ 2.46134*g*(1-g)**3+ 2.33889*g*(1-g)**4+ 2.30769*g*(1-g)**5+ 2.39899*g*(1-g)**6+ 2.46062*g*(1-g)**7+ 2.45657*g*(1-g)**8+ 2.39899*g*(1-g)**9+ 2.32832*g*(1-g)**10+ 2.52860*g*(1-g)**11+ 2.71487*g*(1-g)**12+ 2.99479*g*(1-g)**13+ 2.86814*g*(1-g)**14+ 2.68210*g*(1-g)**15+ 2.57101*g*(1-g)**16+ 2.35940*g*(1-g)**17+ 2.50245*g*(1-g)**18+ 2.41802*g*(1-g)**19+ 2.38582*g*(1-g)**20+ 2.37637*g*(1-g)**21+ 2.49442*g*(1-g)**22+ 2.68093*g*(1-g)**23+ 2.91186*g*(1-g)**24+ 3.11265*g*(1-g)**25+ 2.72605*g*(1-g)**26+ 2.63684*g*(1-g)**27+ 2.53329*g*(1-g)**28+ 2.53979*g*(1-g)**29+ 2.66094*g*(1-g)**30+ 2.71726*g*(1-g)**31+ 2.71070*g*(1-g)**32+ 2.71810*g*(1-g)**33+ 2.78906*g*(1-g)**34+ 2.86553*g*(1-g)**35+ 2.87078*g*(1-g)**36+ 2.89268*g*(1-g)**37+ 2.93942*g*(1-g)**38+ 2.95612*g*(1-g)**39+ 2.88932*g*(1-g)**40+ 2.99771*g*(1-g)**41+ 2.97534*g*(1-g)**42+ 3.02577*g*(1-g)**43+ 3.11270*g*(1-g)**44+ 3.30085*g*(1-g)**45+ 3.53334*g*(1-g)**46+ 3.64106*g*(1-g)**47+ 3.81318*g*(1-g)**48+ 3.84692*g*(1-g)**49+ 3.85411*g*(1-g)**50+ 3.98934*g*(1-g)**51+ 4.00624*g*(1-g)**52+ 4.11768*g*(1-g)**53+ 4.17919*g*(1-g)**54+ 4.24686*g*(1-g)**55+ 4.36939*g*(1-g)**56+ 4.28471*g*(1-g)**57+ 4.47519*g*(1-g)**58+ 4.44962*g*(1-g)**59+ 4.60679*g*(1-g)**60+ 4.49624*g*(1-g)**61+ 4.39355*g*(1-g)**62+ 4.45772*g*(1-g)**63+ 4.49612*g*(1-g)**64+ 4.50028*g*(1-g)**65+ 4.60740*g*(1-g)**66+ 4.59027*g*(1-g)**67+ 4.62891*g*(1-g)**68+ 4.58181*g*(1-g)**69+ 4.53033*g*(1-g)**70+ 4.59925*g*(1-g)**71+ 4.68599*g*(1-g)**72+ 4.63898*g*(1-g)**73+ 4.55659*g*(1-g)**74+ 4.54791*g*(1-g)**75+ 4.52166*g*(1-g)**76+ 4.70673*g*(1-g)**77+ 4.54365*g*(1-g)**78+ 4.44664*g*(1-g)**79+ 4.58804*g*(1-g)**80+ 4.49054*g*(1-g)**81+ 4.49936*g*(1-g)**82+ 4.60185*g*(1-g)**83+ 4.74454*g*(1-g)**84+ 4.76803*g*(1-g)**85+ 4.70565*g*(1-g)**86+ 4.71500*g*(1-g)**87+ 4.67901*g*(1-g)**88+ 4.63931*g*(1-g)**89+ 4.58494*g*(1-g)**90+ 4.47620*g*(1-g)**91+ 4.43996*g*(1-g)**92+ 4.43552*g*(1-g)**93+ 4.35371*g*(1-g)**94+ 4.27156*g*(1-g)**95+ 4.18488*g*(1-g)**96+ 4.06160*g*(1-g)**97+ 3.95800*g*(1-g)**98+ 3.94766*g*(1-g)**99+ 3.87815*g*(1-g)**100+ 3.85432*g*(1-g)**101+ 3.74697*g*(1-g)**102+ 3.74376*g*(1-g)**103+ 3.65457*g*(1-g)**104+ 3.66215*g*(1-g)**105+ 3.63969*g*(1-g)**106+ 3.55208*g*(1-g)**107+ 3.47492*g*(1-g)**108+ 3.41624*g*(1-g)**109+ 3.40554*g*(1-g)**110+ 3.47095*g*(1-g)**111+ 3.38340*g*(1-g)**112+ 3.39202*g*(1-g)**113+ 3.29898*g*(1-g)**114+ 3.32112*g*(1-g)**115+ 3.20446*g*(1-g)**116+ 3.20210*g*(1-g)**117+ 3.09947*g*(1-g)**118+ 3.01736*g*(1-g)**119+ 2.94478*g*(1-g)**120+ 2.88404*g*(1-g)**121+ 2.83013*g*(1-g)**122+ 2.80702*g*(1-g)**123+ 2.79871*g*(1-g)**124+ 2.75788*g*(1-g)**125+ 2.78181*g*(1-g)**126+ 2.70548*g*(1-g)**127+ 2.66104*g*(1-g)**128+ 2.62876*g*(1-g)**129+ 2.55345*g*(1-g)**130+ 2.56301*g*(1-g)**131+ 2.47518*g*(1-g)**132+ 2.42977*g*(1-g)**133+ 2.31733*g*(1-g)**134+ 2.25393*g*(1-g)**135+ 2.24472*g*(1-g)**136+ 2.26919*g*(1-g)**137+ 2.25682*g*(1-g)**138+ 2.22244*g*(1-g)**139+ 2.28817*g*(1-g)**140+ 2.30095*g*(1-g)**141+ 2.23759*g*(1-g)**142+ 2.28874*g*(1-g)**143 ; frml eq14 lpceso1= 29.58342*g+ 31.48294*g*(1-g)+ 28.73536*g*(1-g)**2+ 27.91986*g*(1-g)**3+ 27.63735*g*(1-g)**4+ 26.26295*g*(1-g)**5+ 24.82431*g*(1-g)**6+ 26.39570*g*(1-g)**7+ 25.99496*g*(1-g)**8+ 27.83422*g*(1-g)**9+ 29.44725*g*(1-g)**10+ 27.38068*g*(1-g)**11+ 24.88386*g*(1-g)**12+ 25.29767*g*(1-g)**13+ 22.55054*g*(1-g)**14+ 19.44522*g*(1-g)**15+ 19.06829*g*(1-g)**16+ 21.06498*g*(1-g)**17+ 23.35623*g*(1-g)**18+ 22.75421*g*(1-g)**19+ 24.65348*g*(1-g)**20+ 25.87426*g*(1-g)**21+ 25.41486*g*(1-g)**22+ 22.80368*g*(1-g)**23+ 25.17808*g*(1-g)**24+ 27.17213*g*(1-g)**25+ 27.46894*g*(1-g)**26+ 29.00525*g*(1-g)**27+ 30.17361*g*(1-g)**28+ 29.81001*g*(1-g)**29+ 31.21296*g*(1-g)**30+ 32.34524*g*(1-g)**31+ 30.81980*g*(1-g)**32+ 30.60741*g*(1-g)**33+ 30.44867*g*(1-g)**34+ 28.95341*g*(1-g)**35+ 29.52091*g*(1-g)**36+ 30.22349*g*(1-g)**37+ 25.45272*g*(1-g)**38+ 23.59888*g*(1-g)**39+ 23.04774*g*(1-g)**40+ 23.66173*g*(1-g)**41+ 21.18036*g*(1-g)**42+ 16.23828*g*(1-g)**43+ 20.72894*g*(1-g)**44+ 23.51016*g*(1-g)**45+ 20.82989*g*(1-g)**46+ 24.77944*g*(1-g)**47+ 30.43795*g*(1-g)**48+ 39.12727*g*(1-g)**49+ 45.09983*g*(1-g)**50+ 50.40768*g*(1-g)**51+ 48.24397*g*(1-g)**52+ 47.16611*g*(1-g)**53+ 46.25869*g*(1-g)**54+ 45.36193*g*(1-g)**55+ 44.77775*g*(1-g)**56+ 45.16727*g*(1-g)**57+ 47.43020*g*(1-g)**58+ 47.92679*g*(1-g)**59+ 47.88680*g*(1-g)**60+ 45.88206*g*(1-g)**61+ 46.27760*g*(1-g)**62+ 47.17331*g*(1-g)**63+ 47.81620*g*(1-g)**64+ 48.54115*g*(1-g)**65+ 48.26636*g*(1-g)**66+ 47.63263*g*(1-g)**67+ 49.76932*g*(1-g)**68+ 50.96623*g*(1-g)**69+ 51.49820*g*(1-g)**70+ 51.72865*g*(1-g)**71+ 51.40753*g*(1-g)**72+ 51.27197*g*(1-g)**73+ 50.36695*g*(1-g)**74+ 50.63919*g*(1-g)**75+ 51.91146*g*(1-g)**76+ 53.51923*g*(1-g)**77+ 54.85844*g*(1-g)**78+ 54.33015*g*(1-g)**79+ 53.44449*g*(1-g)**80+ 52.53758*g*(1-g)**81+ 52.90756*g*(1-g)**82+ 50.29926*g*(1-g)**83+ 51.90456*g*(1-g)**84+ 55.29837*g*(1-g)**85+ 57.02453*g*(1-g)**86+ 60.82712*g*(1-g)**87+ 62.98666*g*(1-g)**88+ 62.06884*g*(1-g)**89+ 59.53172*g*(1-g)**90+ 61.07211*g*(1-g)**91+ 62.70311*g*(1-g)**92+ 63.69046*g*(1-g)**93+ 60.08116*g*(1-g)**94+ 53.25824*g*(1-g)**95+ 57.53304*g*(1-g)**96+ 65.29520*g*(1-g)**97+ 67.24928*g*(1-g)**98+ 68.15686*g*(1-g)**99+ 67.16003*g*(1-g)**100+ 66.54254*g*(1-g)**101+ 66.40667*g*(1-g)**102+ 65.30166*g*(1-g)**103+ 62.31507*g*(1-g)**104+ 66.08772*g*(1-g)**105+ 70.13478*g*(1-g)**106+ 73.23784*g*(1-g)**107+ 74.29811*g*(1-g)**108+ 77.78230*g*(1-g)**109+ 78.58188*g*(1-g)**110+ 80.56157*g*(1-g)**111+ 74.98885*g*(1-g)**112+ 67.03761*g*(1-g)**113+ 67.58889*g*(1-g)**114+ 75.73779*g*(1-g)**115+ 77.56013*g*(1-g)**116+ 79.02613*g*(1-g)**117+ 78.30771*g*(1-g)**118+ 79.05073*g*(1-g)**119+ 79.84082*g*(1-g)**120+ 84.20542*g*(1-g)**121+ 87.49627*g*(1-g)**122+ 91.39123*g*(1-g)**123+ 84.18000*g*(1-g)**124+ 81.17961*g*(1-g)**125+ 80.20124*g*(1-g)**126+ 78.90981*g*(1-g)**127+ 83.01534*g*(1-g)**128+ 74.86256*g*(1-g)**129+ 54.40452*g*(1-g)**130+ 54.17270*g*(1-g)**131+ 51.86089*g*(1-g)**132+ 45.14131*g*(1-g)**133+ 39.22570*g*(1-g)**134+ 36.11139*g*(1-g)**135+ 35.01278*g*(1-g)**136+ 34.75309*g*(1-g)**137+ 34.62112*g*(1-g)**138+ 34.68509*g*(1-g)**139+ 35.07682*g*(1-g)**140+ 35.14635*g*(1-g)**141+ 35.16582*g*(1-g)**142+ 35.87402*g*(1-g)**143; frml eq15 lratef= m*conlength+ b2*leadtime+b4*riggser9l1tr+b5*highutil+ b6*year1991+ b7*year1992+ b8*year1993+ b9*year1994+ b10*year1995+ b11*year1996+ b12*year1997+ b13*year1998+ b14*year1999+ b15*year2000+ b16*year2001+ b17*year2002+ b18*year2003+ b19*year2004+ b20*year2005+ b21*year2006+ b22*year2007+ b23*year2008+ b24*year2009+ b25*log(((1-d)*(lpceso/lpceso1)**(1-s)+d*(lpcesg/lpcesg1)**(1-s))**(1/(1-s))) + ff1*DD1+ ff2*DD2+ ff4*DD4+ ff5*DD5+ ff7*DD7+ ff9*DD9+ ff10*DD10+ ff11*DD11+ ff12*DD12+ ff13*DD13+ ff14*DD14+ ff15*DD15+ ff16*DD16+ ff17*DD17+ ff18*DD18+ ff19*DD19+ ff21*DD21+ ff23*DD23+ ff24*DD24+ ff25*DD25+ ff26*DD26+ ff27*DD27+ ff29*DD29+ ff30*DD30+ ff32*DD32+ ff36*DD36+ ff37*DD37+ ff38*DD38+ ff39*DD39+ ff40*DD40+ ff41*DD41+ ff42*DD42+ ff43*DD43+ ff44*DD44+ ff45*DD45+ ff46*DD46+ ff47*DD47+ ff48*DD48+ ff49*DD49+ ff50*DD50+ ff52*DD52+ ff54*DD54+ ff56*DD56+ ff57*DD57+ ff58*DD58+ ff59*DD59+ ff60*DD60+ ff61*DD61+ ff62*DD62+ ff64*DD64+ ff65*DD65+ ff66*DD66+ ff67*DD67+ ff68*DD68+ ff69*DD69+ ff70*DD70+ ff72*DD72+ ff73*DD73+ ff74*DD74+ ff75*DD75+ ff76*DD76+ ff78*DD78+ ff79*DD79+ ff80*DD80+ ff81*DD81+ ff83*DD83+ ff84*DD84+ ff85*DD85+ ff86*DD86+ ff87*DD87+ ff88*DD88+ ff89*DD89+ ff90*DD90+ ff91*DD91+ ff92*DD92+ ff93*DD93+ ff95*DD95+ ff96*DD96+ ff97*DD97+ ff98*DD98+ ff99*DD99+ ff100*DD100+ ff101*DD101+ ff103*DD103+ ff104*DD104+ ff105*DD105+ ff106*DD106+ ff107*DD107+ ff108*DD109+ ff110*DD110+ ff111*DD111+ ff112*DD112+ ff114*DD114+ ff115*DD115+ ff116*DD116+ ff117*DD117+ ff118*DD118+ ff119*DD119+ ff121*DD121+ ff122*DD122+ ff123*DD123+ ff124*DD124+ ff125*DD125+ ff126*DD126+ ff127*DD127+ ff128*DD128+ ff129*DD129+ ff131*DD131+ ff133*DD133+ ff134*DD134+ ff135*DD135+ ff136*DD136+ ff137*DD137+ ff138*DD138+ ff139*DD139+ ff140*DD140+ ff141*DD141+ ff142*DD142+ ff143*DD143+ ff144*DD144+ ff145*DD145+ ff146*DD146+ ff147*DD147+ ff148*DD148+ ff149*DD149+ ff150*DD150+ ff151*DD151+ ff152*DD152+ ff153*DD153+ ff154*DD154+ ff155*DD155+ ff156*DD156+ ff157*DD157+ ff158*DD158+ ff159*DD159+ ff160*DD160+ ff161*DD161+ ff162*DD162+ ff164*DD164+ ff166*DD166+ ff169*DD169+ ff170*DD170+ ff171*DD171+ ff172*DD172+ ff173*DD173+ ff176*DD176+ ff177*DD177+ ff178*DD178+ ff179*DD179+ ff180*DD180+ ff181*DD181+ ff182*DD182+ ff183*DD183+ ff184*DD184+ ff185*DD185+ ff186*DD186+ ff187*DD187+ ff188*DD188+ ff189*DD189+ ff190*DD190+ ff191*DD191+ ff192*DD192+ ff193*DD193+ ff194*DD194+ ff195*DD195+ ff196*DD196+ ff197*DD197+ ff198*DD198+ ff199*DD199+ ff200*DD200+ ff201*DD201+ ff202*DD202+ ff203*DD203+ ff205*DD205; eqsub eq15 eq11; eqsub eq15 eq12; eqsub eq15 eq13; eqsub eq15 eq14; lsq eq15; mmake trauma @coef @ses; write(file='des301213.asc', format='(2F10.3)') trauma; msd DD1 DD2 DD4 DD5 DD7 DD9 DD10 DD11 DD12 DD13 DD14 DD15 DD16 DD17 DD18 DD19 DD21 DD23 DD24 DD25 DD26 DD27 DD29 DD30 DD32 DD36 DD37 DD38 DD39 DD40 DD41 DD42 DD43 DD44 DD45 DD46 DD47 DD48 DD49 DD50 DD52 DD54 DD56 DD57 DD58 DD59 DD60 DD61 DD62 DD64 DD65 DD66 DD67 DD68 DD69 DD70 DD72 DD73 DD74 DD75 DD76 DD78 DD79 DD80 DD81 DD83 DD84 DD85 DD86 DD87 DD88 DD89 DD90 DD91 DD92 DD93 DD95 DD96 DD97 DD98 DD99 DD100 DD101 DD103 DD104 DD105 DD106 DD107 DD109 DD110 DD111 DD112 DD114 DD115 DD116 DD117 DD118 DD119 DD121 DD122 DD123 DD124 DD125 DD126 DD127 DD128 DD129 DD131 DD133 DD134 DD135 DD136 DD137 DD138 DD139 DD140 DD141 DD142 DD143 DD144 DD145 DD146 DD147 DD148 DD149 DD150 DD151 DD152 DD153 DD154 DD155 DD156 DD157 DD158 DD159 DD160 DD161 DD162 DD164 DD166 DD169 DD170 DD171 DD172 DD173 DD176 DD177 DD178 DD179 DD180 DD181 DD182 DD183 DD184 DD185 DD186 DD187 DD188 DD189 DD190 DD191 DD192 DD193 DD194 DD195 DD196 DD197 DD198 DD199 DD200 DD201 DD202 DD203 DD205; print @mean;