Merapi workup for air saturated water (ASW). Oxygen part already done in Solveq.
Now add N2

<  erpc  ><   ph   >< pfluid ><  temp  ><  tempc ><  volbox-1  ><rhofresh>< rhoroc >
 .1000E-11    .00000   1.00000  25.00000   915.000     1.00000      .00000   2.50000

< step increm >< step limit  >< total mixer >
     .00001         .0001           .0000

<  enth   ><  senth  ><  denth  ><  totwat >< solmin ><   rm   ><  aqgrm >< suprnt >
     0.0000     0.0000     0.0000    90.0000 0.000E+00    0.0000    0.0000 0.100E-19

l/v:0  chg    i    i    n    i  lim         i   it    i    i   in incr  min         a  min
hiP:1  bal frac punc loop step  sol looc enth  ref deal psat crem    P solv neut watr  chg
    1    3    0    1  400    2    1    0    0    0    0    1    0    0    2    0    0    0

 saq> <   name   ><   total moles   >< trial molality  ><       gamma     ><mixer total moles>
    1 H+             1.200e-04        0.100000000000E-07  1.00000000000
    2 H2O            1.0              0.100000000000E+01  1.00000000000
    3 Cl-            0.0000002        0.100000000000E-10  1.00000000000
    4 SO4--          1.200e-04        0.100000001169E-10  1.00000000000
    6 HS-           -1.200e-04        0.100000001169E-56  1.00000000000
   13 Na+            0.0000001        0.100000000000E-10  1.00000000000

<     minerals     >  < min trial moles>

<    reactants     >  < weight percent ><  ppm? >

< mins to suppress >
FMQ_HS-1

< don't fractionate>

<  T-P curve data  >


N2 gas                       100.
   32 NH4+           0.0000001        0.100000000000E-10  1.00000000000

N2,aq  +  1.25 H+  +  3.0 H2O  +  0.75 HS-
     =  2.0 NH4+  +  0.75 SO4--