Ray: I have extracted all source code from twodsodium1 into *.h and *.f files, and made a Makefile, thus eliminating update from the procedure. I have compiled and made an executable, so its ready to test. However, as you suspected, the "rhist" mss files no longer exist. This is the file acquired from mss by the twodsodium script, and read by restart.f. One of these will have to be made, I guess by glbmeans. Restart does the following read: read(lu) zp,tno,shto,zpht,sht, + eox,eo,eo3,ehox,eh,eho2,eoh,enoz,eno,eno2,en4s,enas,enasp, + ene,etion,eo2,en2,eo1d,ete,eti,etn,eiop,en2d,enh2o,enh2, + ench4,enh2o2,enco,enco2,enax,enao,enao2,enaoh where: PARAMETER (MX=21,MZ=96,KMX=96) COMMON/PERTRF/TNO(KMX),TNP(KMX),FDIFP(KMX) common/ehist/shto(mz),eox(mx,mz),eo(mx,mz),eo3(mx,mz),ehox(mx,mz), + eh(mx,mz),eho2(mx,mz),eoh(mx,mz),enoz(mx,mz),eno(mx,mz), + eno2(mx,mz),en4s(mx,mz),enas(mx,mz),enasp(mx,mz),ene(mx,mz), + etion(mx,mz),eo2(mx,mz),en2(mx,mz),eo1d(mx,mz),ete(mx,mz), + eti(mx,mz),etn(mx,mz),eiop(mx,mz),en2d(mx,mz),enh2o(mx,mz), + enh2(mx,mz),ench4(mx,mz),enh2o2(mx,mz),enco(mx,mz),enco2(mx,mz) + ,enax(mx,mz),enao(mx,mz),enao2(mx,mz),enaoh(mx,mz) c c Blank common: c COMMON ZP(KMX),ZPHT(KMX),TNMS(KMX),XOM(KMX),XN2M(KMX),XO2M(KMX), 1ZPMHT(KMX),SHTMS(KMX),AMASS(KMX),RHOMS(KMX),TE(KMX),TI(KMX), 2TN(KMX),XNE(KMX),XNO(KMX),XNO2(KMX),XNN2(KMX),XN4S(KMX),XNNO(KMX), 3XN2D(KMX),RHO(KMX),SHT(KMX),XIN2P(KMX),XIOP(KMX),XINOP(KMX), 4XINP(KMX),AMAS(KMX),XCHUE(KMX),XIO2P(KMX),XCO2(KMX),SHCO2(KMX), 5XNH2O(KMX),XNHOX(KMX),XNHO2(KMX),XNOH(KMX),XNHM(KMX),XNO1D(KMX), 6XNO3(KMX),XNH2(KMX),XIHP(KMX),XIOHP(KMX),XIH2OP(KMX),XNNO2(KMX), 7XNCH4(KMX),XNH2O2(KMX),XNHTOT(KMX),XNCO2(KMX),XNCO(KMX),XNH(KMX), 8O1DEX(KMX),XNHE(KMX),XNHEM(KMX),XNARG(KMX),XNARGM(KMX),XN4SM(KMX), 9ZPMS(KMX),GZ(KMX),QHOXI(KMX),ZQHT(KMX),F107,F107A,CO2MIX,DAY, 1GLAT,GLONG,STL,ZI,SCHAP(KMX),RE,BOLTZ