transf.f: 	12-7-2000      -ifstatement J.EQ.JMAXM/2+1 for equator
				at the end of the j-loop
			       -periodic points for sigma1m and sigma2m
			        at 180 lon and at the poles
			       -second history file for sigma1m and sigma2m
			       -second history file for sigma1 and sigma2
			                                ssigma1 and ssigma2
				for control (local: sig1save,sig2save arrays)

threed.f: 	12-7-2000      -adjust bracket where THETAM(I) is calculated
				for southern hemisphere
				
Version with updates of changes from tgcm12 to Art Richmond version (local in modsrc)
 
threed.f:	12-11-2000	-THETAM(I) = ACOS(SQRT(R0/(R00+ZZM(I,J,K)))*COSLTM)
				 not changed R00->RE at the lower boundary
				 PHIM3D,ED1M3D,ED2M3D not smooth at the equator
		12-19-2000	-THETAM(I) = ACOS(SQRT(R0/(R00+ZZM(I,J,K)))*COSLTM)
		                 changed however due to computer precision 
				 ZZM < z0 at lower boundary (\sum wt \ne 1 in
				 apex_mod.f)
				 therefore zm , zzm set to z0 if zm < z0 in
				 transf.f		 
				 
transf.f	12-11-2000	-SINI(I,J) = ZB(I,J)/BMOD(I,J) ZB same as ZZB
				 defined in fieldz.h for tgcm14
				 PHIM3D,ED1M3D,ED2M3D not smooth at the equator
		12-19-2000	-ZZM < z0 at lower boundary (\sum wt \ne 1 in
				 apex_mod.f)
				 therefore zm , zzm set to z0 if zm < z0 in
				 transf.f		 

transmag.h	12-11-2000	-BE3(IMAXGP,0:JMAXGP) added before in transf.f
				 defined

tgcm.F		12-11-2000	-magnetic data from apxparm instead rdmag


including modified mudpack solver and muhcr2 used as direct solver
  new version of mud.f (version 5.0)

modified files:
adheel.f   (new copied)   04-16-2001
cnmmod.f   (new)          04-16-2001
divide.f   (new copied)   04-16-2001
dyncal.f                  04-16-2001
intgrls.f                 04-16-2001
modmud.f   (new)          04-16-2001
muh2cr.f   (new)          04-16-2001
stencil.f  (new copied)   04-16-2001
stenmd.f   (new copied)   04-16-2001


copy of all files from modsrc14me 
  to include these modification within tgcm15

am 02/19/02:
  - try to run coupled with mtiegcm15.me ( Christphe Peymirat's coupling with
    the magnetosphere).
  - starting point ot include electrodynamo changes in tgcm15
  
am 03/19/02:
  - changes in tail.F because POWER didn't fit if flux was integrated
    (s. ~/Tiegcm.d/work.d)
 
am 03/27/02:
  - this is the code for the new version tgcm16 
    all other files are from /home/tgcm/tgcm15
  - this runs with Christophes model mtiegcm15me.tmp