parameters from /raid22/jmuglia/c13/newcodenewfactors, except fetopsed=0.005, wd0=20, kfemax=0.2e-3.
ahisop = 1.2e7, jdiar = 0.2

In order to increase wcdeni, addisop = 2.e8 in Pacific, kappa_h=0.35 

Uses egbert tidal dissipiation energy!

Set
               diff_cbt(i,k,j) = max(kappa_h, min(100.
     &,             zkappa + kappa_h))

at all latitudes, because otherwise I get too low LGM AMOC. 
 


ligand concentrations 1.5 nM (lig = 1.5e-3)

Increase scavenging to counter the higher ligand concentration and have reasonable tracer values: 
      kfeorg = 1.0*(1./86400.),          kfecol = 0.005/86400.

LGM run. 

PMIP3 ice sheets. 
LGM wind stress anomalies are an average among PMIP3 models' anomalies.
A_diffqY multiplied by 0 in SH.
Dust fields from Lambert, LGM, variable sol. (/raid21/jmuglia/iron/lambertdata/O_feflux_lambert_lgm_new.nc).
LGM sub grid bathymetry. 


Transient atmospheric d14c from intcal09 (/home/jmuglia/UVic_2.9/intcal09/A_dc14_intcal.nc).
Start at year -24400 for the ocean to be closer to equlibrium witht he atmosphere in the initial phase.  

Use /raid23/jmuglia/ironruns3/lgm-default-diffqy025/rest.-00026500.01.01.nc as restart file.