/home/tgcm/bin/Diffsrc executing from /home/foster/tiegcm/tiegcm-dev8/modsrc.compox at Tue Sep 5 12:49:14 MDT 2006 Using tgcmroot from TGCMROOT env var: /home/tgcm Source code directory = /home/roble/tiegcm/tiegcmdr10/modsrc.compox ======================================================================== Diff of /home/roble/tiegcm/tiegcmdr10/modsrc.compox/aurora.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/aurora.F /home/roble/tiegcm/tiegcmdr10/modsrc.compox/aurora.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/aurora.F are identical ======================================================================== Diff of /home/roble/tiegcm/tiegcmdr10/modsrc.compox/bndry.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/bndry.F 30a31,33 > ! OX bottom boundary (set in sub comp_ox, and used in sub bndry_comp): > real :: xoxlb(nlonp4,nlat) > ! 578d580 < real :: xoxlb(nlat) 594c596,599 < xoxlb(:) = 2.e-5 ! bottom boundary of ox (zp -10) --- > ! xoxlb(i,j) is calculated in comp_ox (comp_ox must be called > ! before this routine) > ! xoxlb(:) = 3.e-6 ! bottom boundary of ox (zp -10) > ! xoxlb(:) = 2.e-5 ! bottom boundary of ox (zp -10) 616c621 < fb(i,2,lat) = 2.*xoxlb(lat)*rmo3a(i)/barm(1,i) --- > fb(i,2,lat) = 2.*xoxlb(i,lat)*rmo3a(i)/barm(1,i) ======================================================================== >>> WARNING: Cannot find source file /home/roble/tiegcm/tiegcmdr10/modsrc.compox/chemrates.F ======================================================================== Diff of /home/roble/tiegcm/tiegcmdr10/modsrc.compox/comp.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/comp.F /home/roble/tiegcm/tiegcmdr10/modsrc.compox/comp.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/comp.F are identical ======================================================================== Diff of /home/roble/tiegcm/tiegcmdr10/modsrc.compox/comp_major.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/comp_major.F /home/roble/tiegcm/tiegcmdr10/modsrc.compox/comp_major.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/comp_major.F are identical ======================================================================== Diff of /home/roble/tiegcm/tiegcmdr10/modsrc.compox/comp_o1d.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/comp_o1d.F /home/roble/tiegcm/tiegcmdr10/modsrc.compox/comp_o1d.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/comp_o1d.F are identical ======================================================================== Diff of /home/roble/tiegcm/tiegcmdr10/modsrc.compox/comp_ox.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/comp_ox.F 30a31,36 > use init_module,only: istep > ! > ! xoxlb(nlonp4,nlat) bottom boundary of OX (bndry module data in bndry.F), > ! is defined by this routine, then used by sub bndry_comp in bndry.F. > ! > use bndry_module,only: xoxlb ! (nlonp4,nlat) 106a113,114 > call addfld('NO_ox','NO_ox',' ',no(:,i0:i1), > | 'lev',k0,k1,'lon',i0,i1,lat) 161,162c169,170 < | beta9(k,i,lat) *no (k,i)*rmassinv_no+ < | pdo3m(k,i))) --- > | beta9(k,i,lat) *no (k,i)*rmassinv_no)+ > | pdo3m(k,i)) 197,198c205,206 < | rkm37(k,i,lat)*h (k,i)*rmassinv_h *o3_o(k,i)+ < | beta9 (k,i,lat)*no (k,i)*rmassinv_no *o3_o(k,i)+ --- > | (rkm37(k,i,lat)*h (k,i)*rmassinv_h *o3_o(k,i)+ > | beta9 (k,i,lat)*no (k,i)*rmassinv_no *o3_o(k,i))+ 209,211c217,219 < | (rkm37(k,i,lat)*h (k,i)*rmassinv_h *o3_o(k,i)+ < | beta9(k,i,lat)*no (k,i)*rmassinv_no *o3_o(k,i)+ < | pdo3m(k,i)*o3_o(k,i))*o_ox(k,i)*xnmbar(k,i)/barmm(k,i) --- > | ((rkm37(k,i,lat)*h (k,i)*rmassinv_h *o3_o(k,i)+ > | beta9(k,i,lat)*no (k,i)*rmassinv_no *o3_o(k,i))* > | xnmbar(k,i)+pdo3m(k,i)*o3_o(k,i))*o_ox(k,i) 226a235,243 > > ! write(6,"('comp_ox: lat=',i3,' k=',i3,' i=',i3, > ! | ' o3_o=',e12.4,' o_ox=',e12.4,' ox_prod1=',e12.4, > ! | ' ox_prod2=',e12.4,' ox_loss1=',e12.4,' ox_loss2=',e12.4, > ! | /,' o2_prod1,2,3=',3e12.4,' o2_loss=',e12.4,' rmtru=', > ! | e12.4)") lat,k,i,o3_o(k,i),o_ox(k,i),ox_prod1(k,i), > ! | ox_prod2(k,i),ox_loss1(k,i),ox_loss2(k,i),o2_prod1(k,i), > ! | o2_prod2(k,i),o2_prod3(k,i),o2_loss(k,i),rmtru_ox(k,i) > 257c274,275 < --- > ! > ! Upper boundaries (nlevp1) are not calculated: 267a286,288 > ! > ! Solve quadratic for lower boundary of OX: > xoxlb(i,lat) = 2.e-5 268a290 > ======================================================================== Diff of /home/roble/tiegcm/tiegcmdr10/modsrc.compox/compart.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/compart.F 96,101c96,99 < call addfld('O1_OUT',' ',' ',o1_out(:,lon0:lon1), < | 'lev',lev0,lev1,'lon',lon0,lon1,lat) < call addfld('O3_OUT',' ',' ',o3_out(:,lon0:lon1), < | 'lev',lev0,lev1,'lon',lon0,lon1,lat) < call addfld('OX_OUT',' ',' ',ox(:,lon0:lon1), < | 'lev',lev0,lev1,'lon',lon0,lon1,lat) --- > ! call addfld('O1_OUT',' ',' ',o1_out(:,lon0:lon1), > ! | 'lev',lev0,lev1,'lon',lon0,lon1,lat) > ! call addfld('O3_OUT',' ',' ',o3_out(:,lon0:lon1), > ! | 'lev',lev0,lev1,'lon',lon0,lon1,lat) ======================================================================== Diff of /home/roble/tiegcm/tiegcmdr10/modsrc.compox/cons.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/cons.F /home/roble/tiegcm/tiegcmdr10/modsrc.compox/cons.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/cons.F are identical ======================================================================== Diff of /home/roble/tiegcm/tiegcmdr10/modsrc.compox/duv.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/duv.F /home/roble/tiegcm/tiegcmdr10/modsrc.compox/duv.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/duv.F are identical ======================================================================== Diff of /home/roble/tiegcm/tiegcmdr10/modsrc.compox/dynamics.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/dynamics.F 50d49 < use o2o_module,only: comp_o2o 597,624d595 < ! Calculate matrix coefficients for partitioning of OX: < ! subroutine comp_ox(tn,o2,ox,op,o1,o1d,o3,no,n4s,n2d,o2p,n2p,nplus, < ! | nop,ne,xiop2p,xiop2d,barm,xnmbar,lev0,lev1,lon0,lon1,lat) < ! < ! call comp_ox( < ! | tn (levd0,lond0,lat,itp), ! tn (deg K) < ! | o2 (levd0,lond0,lat,itp), ! o2 (mmr) < ! | ox (levd0,lond0,lat,itp), ! ox (mmr) < ! | op (levd0,lond0,lat,itc), ! o+ (mmr) < ! | o1 (levd0,lond0,lat,itp), ! o (mmr) < ! | o1d (levd0,lond0,lat), ! o(1d) (mmr) < ! | o3 (levd0,lond0,lat,itp), ! o3 (mmr) < ! | no (levd0,lond0,lat,itc), ! no (mmr) < ! | n4s (levd0,lond0,lat,itc), ! n4s (mmr) < ! | n2d (levd0,lond0,lat,itc), ! n(2d) (mmr) < ! | o2p (levd0,lond0,lat,itp), ! o2+ < ! | n2p (levd0,lond0,lat), ! n2+ < ! | nplus (levd0,lond0,lat), ! n+ < ! | nop (levd0,lond0,lat), ! no+ < ! | ne (levd0,lond0,lat,itc), ! ne (cm3) < ! | xiop2p(levd0,lond0,lat), ! o+(2p) < ! | xiop2d(levd0,lond0,lat), ! o+(2d) < ! | barm (levd0,lond0,lat,itp), ! mbar < ! | xnmbar(levd0,lond0,lat), ! p0e(-z)... < ! | 1,nlevp1,lon0,lon1,lat) < ! if (debug) write(6,"('dynamics after comp_ox: lat=',i3, < ! | ' istep=',i5)") lat,istep < ! 902,903c873 < ! subroutine comp_ox(tn,o2,ox,op,o1,o1d,o3,no,n4s,n2d,o2p,n2p,nplus, < ! | nop,ne,xiop2p,xiop2d,barm,xnmbar,lev0,lev1,lon0,lon1,lat) --- > ! (also define OX lbc in xoxlb) 929,951d898 < ! New comp_o2o: < ! call comp_o2o( < ! | tn (levd0,lond0,lat,itp), ! neutral temperature < ! | o2 (levd0,lond0,lat,itp), ! molecular oxygen (mmr) < ! | o1 (levd0,lond0,lat,itp), ! atomic oxygen (mmr) < ! | barm (levd0,lond0,lat,itp), ! mean molecular weight < ! | op (levd0,lond0,lat,itp), ! O+ < ! | no (levd0,lond0,lat,itp), ! NO < ! | n4s (levd0,lond0,lat,itp), ! N(4S) < ! | n2d (levd0,lond0,lat,itp), ! N(2D) < ! | o2p (levd0,lond0,lat,itp), ! O2+ < ! | ne (levd0,lond0,lat,itp), ! electron density < ! | o1d (levd0,lond0,lat), ! o1d (mmr) from comp_o1d < ! | o3 (levd0,lond0,lat,itc), ! o3 (mmr) from comp_o1d < ! | n2p (levd0,lond0,lat), ! N2+ < ! | nplus(levd0,lond0,lat), ! N+ < ! | nop (levd0,lond0,lat), ! NO+ < ! | xiop2p(levd0,lond0,lat), ! from oplus < ! | xiop2d(levd0,lond0,lat), ! from oplus < ! | 1,nlevp1,lon0,lon1,lat) < ! if (debug) write(6,"('dynamics after comp_o2o: lat=',i3, < ! | ' istep=',i5)") lat,istep < ! ======================================================================== Diff of /home/roble/tiegcm/tiegcmdr10/modsrc.compox/fields.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/fields.F /home/roble/tiegcm/tiegcmdr10/modsrc.compox/fields.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/fields.F are identical ======================================================================== Diff of /home/roble/tiegcm/tiegcmdr10/modsrc.compox/mgw.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/mgw.F /home/roble/tiegcm/tiegcmdr10/modsrc.compox/mgw.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/mgw.F are identical ======================================================================== Diff of /home/roble/tiegcm/tiegcmdr10/modsrc.compox/minor.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/minor.F /home/roble/tiegcm/tiegcmdr10/modsrc.compox/minor.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/minor.F are identical ======================================================================== Diff of /home/roble/tiegcm/tiegcmdr10/modsrc.compox/nchist.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/nchist.F /home/roble/tiegcm/tiegcmdr10/modsrc.compox/nchist.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/nchist.F are identical ======================================================================== Diff of /home/roble/tiegcm/tiegcmdr10/modsrc.compox/oplus.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/oplus.F /home/roble/tiegcm/tiegcmdr10/modsrc.compox/oplus.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/oplus.F are identical ======================================================================== Diff of /home/roble/tiegcm/tiegcmdr10/modsrc.compox/qinite.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/qinite.F /home/roble/tiegcm/tiegcmdr10/modsrc.compox/qinite.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/qinite.F are identical ======================================================================== Diff of /home/roble/tiegcm/tiegcmdr10/modsrc.compox/qrj.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/qrj.F /home/roble/tiegcm/tiegcmdr10/modsrc.compox/qrj.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/qrj.F are identical ======================================================================== Diff of /home/roble/tiegcm/tiegcmdr10/modsrc.compox/settei.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/settei.F /home/roble/tiegcm/tiegcmdr10/modsrc.compox/settei.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/settei.F are identical ======================================================================== Diff of /home/roble/tiegcm/tiegcmdr10/modsrc.compox/tuvbnd.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/tuvbnd.F /home/roble/tiegcm/tiegcmdr10/modsrc.compox/tuvbnd.F and /home/foster/tiegcm/tiegcm-dev8/modsrc.compox/tuvbnd.F are identical