    /*PL1 UNIONIZED AMMONIA ROUTINE*/                                   00000100
1/*121100*/          NH3:    PROC (PH,T,P,IT,IP,Z);                     13450000
 /*121200*/        /* BY AUTH OF RXG*/                                  13460000
 /*121300*/        /* INPUT    PH, T(TEMP), P(PAR VALUE)        */      13470000
 /*121400*/        /*          IT  1 FOR C TEMP     2 FOR F TEMP     */ 13480000
 /*121500*/        /*          IP  1 FOR PAR NO 610(NH3-N)   */         13490000
 /*121600*/        /*          IP 2 FOR PAR NO 71845(NH3-NH4)   */      13500000
 /*121700*/        /* OUTPUT   Z(UN-IONIZED NH3-NH3)    */              13510000
 /*121800*/        DCL (T,PH,P,Z) FLOAT DEC(6), (IT,IP) FIXED BIN(31);  13520000
 /*121900*/                  DCL (Q,CF) FLOAT DEC(6);                   13530000
 /*122000*/        IF P = 0E0 THEN DO; Z =    0.0; RETURN; END;         13540000
 /*122100*/                  IF P< 0E0 THEN GO TO BAD;                  13550000
 /*122200*/        IF IT=2 THEN Q=5E0*(T-32E0)/9E0 /* F TO C */;        13560000
 /*122300*/                  ELSE IF IT=1 THEN Q=T;                     13570000
 /*122400*/                  ELSE GO TO BAD;                            13580000
 /*122500*/        IF IP=1 THEN CF=17E0/14E0/* 17/14 FOR N-TO-NH3  */;  13590000
 /*122600*/ELSE IF IP=2 THEN CF=17E0/18E0;/*     17/18 FOR NH4-TO-NH3 */13600000
 /*122700*/                  ELSE GO TO BAD;                            13610000
 /*122800*/                  Z=9.0182E-2+2729.92E0/(Q+273.18E0);        13620000
 /*122900*/                  Z=1E0/(10E0**(Z-PH)+1E0);                  13630000
 /*123000*/                  Z=P*Z*CF;                                  13640000
 /*123100*/                  GO TO RET;                                 13650000
 /*123200*/          BAD:     Z = .1E-20; RETURN;                       13660000
 /*123300*/          RET:    END;                                       13670000
