C==============================================================
1. Time series:
      09/09/1999 19:09 Yongli Wang IGPP/UCLA

This is the simplified version of CTIM. The only change to the packeage
is:

1. Unnecessary files are delected.
2. Different files are put to different directories.
3. Several nag functions is replaced by new functions, so that now
   the program does not need nag library.
4. Results are tested, only with several minor warning messages.
5. This program can be used as the base of further CTIM testing!

C================================================================
2. Time series:
      09/10/1999 09:09 Yongli Wang IGPP/UCLA

Add the output code to output data for comparison.

C================================================================
3. Time series:
      09/10/1999 09:50 Yongli Wang IGPP/UCLA

1. Move moddta,startup,hist to upper directory, so that we can save disk
   space when saving them.
2. Finish programming file make_compare.f, which can be used to 
   generate the 'compare_full' and 'compare' to do the comparasion.
3. move data to ../data.
4. Changed tim.f to fit change 1.
5. Set program to stop at nn=288, thus we can only need to get the 
   first 288 step result with compare, while full comparison is only 
   needed when major lot of changings have been made.
6. Changed clean to fit changing in 1.

   Test Run: Results are good till nn=288.

C================================================================
4. Time series:
      09/10/1999 10:20 Yongli Wang IGPP/UCLA

1. Delete most of the commented out part.

   Test Run: Results are good till nn=288.

C================================================================
5. Time series:
      09/10/1999 13:05 Yongli Wang IGPP/UCLA

1. Implement read_hough() in ctim_init.
2. Implement common block /c_hough/ for hough package.

   Test Run: Results are good till nn=288.

C================================================================
6. Time series:
      09/10/1999 14:15 Yongli Wang IGPP/UCLA

1. Implement get_ssa() in ctim_init.
2. Implement common block: /c_ssa/ and /tide / in varcom.h.
3. Implement common blocks: /switch**/ and changed the relative parts
   in the package.

   Test Run: Results are good till nn=288.

C================================================================
7. Time series:
      09/10/1999 15:45 Yongli Wang IGPP/UCLA

1. Implementation of common block:
   /alldata/,/input/ /c_smooth_condition/, /c_dteld/, /c_save_no/,
   /c_loop_no/,/c_Is/,/corr/.

   Test Run: Results are good till nn=288.

C================================================================
8. Time series:
      09/10/1999 16:32 Yongli Wang IGPP/UCLA

1. Implementation of common block:
   /c_sza/,/grid_global/,c_fo2/,/c_facjh/,/c_asymmfac/,/c_kp/,/c_sum/,
   /c_xx/,/c_sda/,/c_c7/,/c_current/,/c_hi/,/param/,/c_aeuv/,
   /c_pflux_char/,/c_common/,/pole/,/basic/,/c_om/.

   Test Run: Results are good till nn=288.

2. Implementation of common block:
   driff,grid_th,grid_ph,grid_mag,c_temp2,c_temp1,c_for_edep,c_elden,
   c_elec,c_vion,q_radiation,c_ml,c_angdif,c_p1233,c_stfac,c_timef,
   c_div.

   Test Run: Results are good till nn=288.

C================================================================
9. Time series:
      09/10/1999 18:00 Yongli Wang IGPP/UCLA

1. Implementation of common block:
   energy,ionden,substorm,glo,eld,delt,c_hemi,sx,c_sth2,c_th,counters,c_diff2,
   e_field,c_magfield.
2. Change the related part in elfmods, rice.f, etc. 

   Test Run: Results are good till nn=288.

3. Implementation of common block:
   c_surround,c_sumenergy,s_plus_min,c_ramin_rvin,blk8,s,newdat,
   c_tempory,rice,blk20,c_qiont

   Test Run: Results are good till nn=288.

4. Implementation of common block:
   a,epgrid,fost,time,c_base_time.

   Test Run: Results are good till nn=288.
   Test Run: Results are good till nn=1440.

C================================================================
10. Time series:
      09/10/1999 21:50 Yongli Wang IGPP/UCLA

1. Implement param.h, and change vardec.h to fit it. 
2. Change test time to nn=288.
3. Add /c_pp/ common block.
4. Change tucan to main program and tim to subroutine program.  

   Test Run: Results are good till nn=288.

C================================================================
11. Time series:
      09/10/1999 22:10 Yongli Wang IGPP/UCLA

1. Implement final version of ctim_init.f into this version of ctim. 
   Related changings are made to finish this job. Several major 
   components is recomment out for compiling and testing.

   Test Run: Results are good till nn=288.

2. Get rid of 444 stuff.
3. Copy ctim_cpcpd.inp, ctim_gwatts.inp, ctim_levpi.inp to the directory.
4. Implement part of the subroutine read_basic_param().
5. Implement set_time().

   Test Run: Results are good till nn=288.

C================================================================
12. Time series:
      09/10/1999 22:10 Yongli Wang IGPP/UCLA

1. Delete sw14 part.
2. Make use sw4=.T. to simply the program.
3. Delete sw5 and sw6 part.
4. Replace aa,bb,yy with the values that they have in the program.
5. Do the similar thing in tim.f.
6. Delete sw21,sw22,sw23 part in tucan.f.
7. Delete fnum which seems useless.
9. Put exsize = 0.0 directly in the program.
10. Put examp = 1. directly in the program.
11. Put fden = 1. directly in the program.
12. delete sw15 from the program.
13. delete sw13 from the program.

   Test Run: Results are good till nn=288.
   Test Run: Results are good till nn=1440.

C================================================================
13. Time series:
      09/11/1999 12:35 Yongli Wang IGPP/UCLA

1. delete sw30 part from the program.
2. delete sw18 part from the program.
3. delete sw19 part from the program.
4. delete sw33 part from the program.
5. delete sw34 part from the program.
6. Change tim.f into two subroutines:
   ctim_open_files() and ctim_read_switchs()
   and change relative part in ctim_init.f.

   Test Run: Results are good till nn=288.

7. delete sw7 part from the program.
8. delete xx from tucan.f.
9. delete several display lines for den, etc. 
10. totally implement inshaun(). 

   Test Run: Results are good till nn=288.

11. delete airglw,glow(). 
12. delete vpole.
13. delete sw26 part in tucan.f.

   Test Run: Results are good till nn=288.

C================================================================
14. Time series:
      09/11/1999 12:35 Yongli Wang IGPP/UCLA

1. Fully implement read_basic_param() in ctim_init.f.

   Test Run: Results are good till nn=288.

2. changing the orders of different groups. 

   Test Run: Results are good till nn=288.

C================================================================
15. Time series:
      09/11/1999 15:40 Yongli Wang IGPP/UCLA

1. Fully implement read_init_field_value().
2. delete setting of psao, etc before this function, which seems
   to be useless in the program.
3. delete iflag part. 
4. Fully implement get_var0().

   Test Run: Results are good till nn=288.

5. replace several variables with numbers in tucan.f.
6. Fully implement init_loop_param().

   Test Run: Results are good till nn=288.

C================================================================
16. Time series:
      09/11/1999 16:40 Yongli Wang IGPP/UCLA

1. Partial get_param1() implementation.

   Test Run: Results are good till nn=288.

C================================================================
17. Time series:
      09/11/1999 16:50 Yongli Wang IGPP/UCLA

1. replace the rest sw34 in tucan.f.
2. delete ex and ey parts, which seems useless in tucan.f. But still
   keep them in varcom.h, which could be used in other routines. 
3. Change certain parts before loop to make them more similar to the 
   standard routine of get_param1().

   Test Run: Results are good till nn=288.

4. Further change certain parts before loop to make them more 
   similar to the standard routine of get_param1(). Get rid of 
   rkm,rkt,which are useless. 

   Test Run: Results are good till nn=288.

C================================================================
18. Time series:
      09/11/1999 17:50 Yongli Wang IGPP/UCLA

1. Unite umut.
2. Delete ssa in this part. 
3. Delete epsp and other pole values definition, because they will be
   set soon after it with subroutine, read_init_field_value().
4. Delete vx and vy, because they will be
   set soon after it with subroutine, read_init_field_value().

   Test Run: Results are good till nn=288.
   Test Run: Results are good till nn=864.

C================================================================
19. Time series:
      09/11/1999 19:50 Yongli Wang IGPP/UCLA

1. Fully implement get_param1().

   Test Run: Results are good till nn=288.


C================================================================
20. Time series:
      09/11/1999 20:00 Yongli Wang IGPP/UCLA

1. Finish most part of the implemetaion of ctim_init, except
   CALL MHD_Efield().
2. Fully shorten the part before loop in tucan.f to be one subroutine:
   call ctim_init().
3. Fully implement subroutine ctim_output() in the program.

   Test Run: Results are good till nn=288.

C================================================================
21. Time series:
      09/11/1999 20:10 Yongli Wang IGPP/UCLA

1. Delete all theh displaying parts related with iwrite6.
2. Restructure the part with iwrite2.
3. Delete some display relating with i9,i10 and nnsolar.

   Test Run: Results are good till nn=288.

C================================================================
22. Time series:
      09/11/1999 20:40 Yongli Wang IGPP/UCLA

1. Implement get_nn(),set_time(),display_time(),get_ut(),counter().

   Test Run: Results are good till nn=288.

2. Implement init_nn2_loop().
3. Implement get_grid_param1(),get_grid_param2().

   Test Run: Results are good till nn=288.
 
4. Delete the rest of sw6 part.
5. Replace sw27 with real value. 
6. Replace sw3 with real value. 
7. Get rid of elx,ely, which seem to be useless.
8. Replace sw14 with real value. Delete sheff method.
9. Implement get_electric_field().

   Test Run: Results are good till nn=288.

C================================================================
23. Time series:
      09/11/1999 22:00 Yongli Wang IGPP/UCLA

1. Replace sw31 with its real value.
2. Fully implement get_magnetic_field().
3. Replace sw2 with its real value.
4. Fully implement get_MAG_coord(nn).

   Test Run: Results are good till nn=288.

5. fully implement get_electron_density().

   Test Run: Results are good till nn=288.
   Test Run: Results are good. But errors less than 1% appears in 
             Vx and Vy. Which might be the changing of the certain parts
	     before. But I do not think they matter much in the whole
	     simulation. At the same time, the Ionosphere error appears at 
	     the same locations as the standard version, which further 
	     proves that we can use this result safely.
   
6. Implement get_om().
7. Implement get_param3().

   Test Run: Results are good till nn=576, except the small errors.

8. Implement ctim_smoothing(). 

   Test Run: Results are good till nn=576, except the small errors.

9. Implement get_cp_eps().

   Test Run: Results are good till nn=576, except the small errors.

10. Implement get_var1(nn).
11.Include 'Hell!' part to test whether a part is used.

   Test Run: Results are good till nn=864, except the small errors.

C================================================================
24. Time series:
      09/12/1999 15:00 Yongli Wang IGPP/UCLA

1. Delete a part that is useless.
2. Replace sw11 with its real value.
3. Implement get_var2(nn).

   Test Run: Results are good till nn=288.

C================================================================
25. Time series:
      09/12/1999 16:20 Yongli Wang IGPP/UCLA
 
   Test Run: Results are good till nn=1440, except some small errors 
             in the middle of vx and vy.

C================================================================
26. Time series:
      09/12/1999 23:30 Yongli Wang IGPP/UCLA
 
1. Add vxn,vyn,epsn,htn and replace it into propriate parts of tucan.f.
2. Implement get_rvin_ramin().
3. Implement get_abc().

   Test Run: Results are good till nn=288.

C================================================================
27. Time series:
      09/13/1999 08:00 Yongli Wang IGPP/UCLA
 
   Test Run: Results are good till nn=1440. Unless four very small errors
   in all the six files.
 
1. Implement ctim_assym().

   Test Run: Results are good till nn=288.

2. change get_surrounding() relative part to make it more like this
   routine.

   Test Run: Results are good till nn=288.

3. Implement get_surrounding().

   Test Run: Results are good till nn=288.

4. Replace sw24 with its real value.
5. Implement get_electron_density1(nn,nn2,ut,alden(1,m,l)).

   Test Run: Results are good till nn=576.

6. Implement get_solar_effect().

   Test Run: Results are good till nn=720.

C================================================================
28. Time series:
      09/13/1999 16:40 Yongli Wang IGPP/UCLA
 
1. Implement get_energy_deposition().
2. Replace sw17 with its real value.
3. Replace sw20 with its real value.

   Test Run: Results are good till nn=288.

4. Get rid of a part between:  get_abc() and get_rvin_ramin().
5. Partial Changing ctim_flux() related part. 

   Test Run: Results are generally good till nn=1440. But two more slight
   errors appear.

6. Change parts of vel1,vel2,vel3 to vx,vy,eps.

7. Implement ctim_flux().

   Test Run: Results are good till nn=432.

8. Implement ctim_integration().
9. Implement get_sum_energy().
10.Implement get_pole_values().
11.Move randgs() and rand() into tool.f from tucan.f.

   Test Run: Results are good till nn=1440, except 5 small errors.

C================================================================
28. Time series:
      09/13/1999 16:40 Yongli Wang IGPP/UCLA
 
1. Get rid of elenh and rht, which are useless in the program.
2. Get rid of nnq.
3. Clearence. 

   Test Run: Results are good till nn=1440, except 5 small errors.

C================================================================
28. Time series:
      09/14/1999 16:40 Yongli Wang IGPP/UCLA
 
1. Implement reduce_time().
2. Implement get_param2().
3. Delete all the format sentences.
4. Implement ctim_loop().

   Test Run: Results are good till nn=1440, except 5 small errors.


C================================================================
28. Time series:
      09/14/1999 16:40 Yongli Wang IGPP/UCLA
 
1. Change displaying.
2. Get rid of vxl,vyl,epsl stuff.
3. Delete vel1,vel2,vel3.
4. Delete vx2,vy2,eps2.
5. Create b1b2.inp.
6. Create switch.inp.

   Test Run: Results are good till nn=1440, though there are some small errors
             between it and the origin_modification results.

   file: ctim_199909151910.tar.gz

C================================================================
29. Time series:
      09/15/1999 19:20 Yongli Wang IGPP/UCLA

1. Move tim.f part into ctim_init.f.
2. Change ctim_open_files() and simplify tim.inp. 

   Test Run: Results are good till nn=180.

   file: ctim_199909160945.tar.gz

C================================================================
30. Time series:
      09/16/1999 09:50 Yongli Wang IGPP/UCLA

1. Delete random part in the electric field.

   Test Run: Results are good till nn=180. Test run shows that there are
   generally the same pattens in old and new results. But in detail, there
   are certain differences which is quite obvious and unneglectable.

2. Change the structure of ctim_init.f.

   Test Run: Results are good till nn=180.

   file: ctim_199909161100.tar.gz

C================================================================
31. Time series:
      09/16/1999 11:00 Yongli Wang IGPP/UCLA

1. Change time to 01/10/1997 12:00.
2. Add Jimmy's electric field.

   Test Run: Results are good till nn=1440.
   results are saved in data_final_MHD.

   file: ctim_199909171645.tar.gz

C================================================================
32. Time series:
      09/17/1999 17:20 Yongli Wang IGPP/UCLA

1. Change to ctim_ name system. 

   Test Run: Results are good till nn=420.

   file: ctim_199909171745.tar.gz

C================================================================
33. Time series:
      09/17/1999 20:00 Yongli Wang IGPP/UCLA

1. Implement geo_sm() to do the coordinate transformation.
2. Implement trans() to output data in SM coordinate system together
   with GEO coordinate system. 
3. Add new make_compare_new.f to produce new compare file.

   Test Run: Results are good till nn=1440.

C================================================================
34. Time series:
      09/20/1999 17:30 Yongli Wang IGPP/UCLA

1. When save to SM, add a very small number 0.00001 to those non-positive
   values, which will be helpful when plotting them in log scale.
2. Change output file name's number to 6 bits long and in seconds.

   Test Run: Results are good till nn=360.(the same as run 33.)

3. Change MHD basetime to 1997 09 22:00 from 1997 10 10:00. 

   Test Run: Results are good till nn=1200. There are no MHD input files
   after that time, so this has to be the end of the CTIM.

4. Add Combine.f, which will produce Combine. Combine will combine the separate
   files of certain time to a unit file, which will be used in later plotting. 
   Test run: Working good!

C================================================================
35. Time series:
      09/26/1999 19:20 Yongli Wang IGPP/UCLA

1. changed MHD_Efield part, which determines the direction of Exns.

   Test run: Working good!

C================================================================
36. Time series:
      10/03/1999 18:40 Yongli Wang IGPP/UCLA

1. Start campaign after travel in Boulder.
2. Change the procedure of getting integrated electron density and include
   in it the volume factor.
3. Include output of temperature and electric field.
4. change electric field to Foster field and test the results of it.
5. before major modification in electric field.

C================================================================
37. Time series:
      10/05/1999 18:40 Yongli Wang IGPP/UCLA

1. change sda producing process with Jimmy's function package.
   result -21.9 =>-22.2
2. Delete corrected magnetic field part in spole.f
3. replace all the transformation in spole.f with Jimmy's 
   coordinate transformation package.

   Results are good in all the time.

C================================================================
38. Time series:
      10/06/1999 10:00 Yongli Wang IGPP/UCLA

1. Test Jimmy's electric field input.
2. Finish efield.f, which can read in MHD simulation potential field
   and get electric field in GEO coordinate system.
   It also provides output for DX drawing of potential field,electric
   field. (Results great!)
3. Finish ctim_elecfield.f, which can get any electric field in GEO
   system at any grid point.
4. Modify ctim_get_electric_field() in ctim_loop.f
5. Change Exns and Eyns to (91,20) from (1,20,20)
6. Finish changing all the relative parts about electric field.

   file: ctim_199910062130.tar.gz

C================================================================
39. Time series:
      10/06/1999 21:30 Yongli Wang IGPP/UCLA

1. Runs greatest without electric field input.


C================================================================
40. Time series:
      10/07/1999 20:30 Yongli Wang IGPP/UCLA

1. Runs good but stop at 1177 because of the errors of molion.f
2. Implement save_potential_Efield(nn) into ctim_output.f

   Runs great until step 1177.

3. Changed with new version of molion.f.
   Runs great until step 1177.

   file: ctim_199910081900.tar.gz

C================================================================
41. Time series:
      10/07/1999 20:30 Yongli Wang IGPP/UCLA

1. Runs very good without electric field input.

  file: ctim_199910101345.tar.gz

C================================================================
42. Time series:
      10/10/1999 15:20 Yongli Wang IGPP/UCLA

1. Delete all vxl,vyl,epsl.
2. Delete all vx2,vy2,eps2.
   
   Runs good at step = 120.

1. Change inshaun.f, use param.h to replace parameter defination.
2. Implement alldat common block in inshaun.f.
3. Change small errors of blk2 common block in varcom.h.
4. Implement epgrid common block in inshaun.f.
5. Change errors in /blk20/ in varcom.h.
6. Implement blk20 common block in inshaun.f. 
7. Implement /blk8/ in varcom.h and inshaun.f
8. Implement /blk12/ in varcom.h and inshaun.f.
9. Implement /blk4/ in varcom.h and inshaun.f.
10. Implement /blk5/ in varcom.h and inshaun.f.
11. Implement /blk7/ in varcom.h and inshaun.f.
12. Implement /blk11/ in varcom.h and inshaun.f.
13. Finalize modifying ishaun.f head.

   Runs good at step = 120.


14. First modernize shajutor.f.

   Runs good at step = 120.

   file: ctim_199910111050.tar.gz

C================================================================
43. Time series:
      10/11/1999 10:50 Yongli Wang IGPP/UCLA

1. Finish changing inshaun.f

   Runs good at step = 120.

   file: ctim_199910111130.tar.gz

   Runs good at step = 1440.

   file: ctim_199910111300.tar.gz

C================================================================
44. Time series:
      10/11/1999 13:00 Yongli Wang IGPP/UCLA

1. Revise specht and the relative part of the program.
2. Major change of ctim_get_param1() in ctim_init.f.

   Runs good at step = 1440.

3. Finish part of unit writting. 
   Runs good at step = 120.

C================================================================
45. Time series:
      10/11/1999 14:40 Yongli Wang IGPP/UCLA

1. Slightly revise spole2().
2. Partially revise elecden.f.
3. Correct a fatal error at ctim_get_param1() in ctim_init.f.

   Runs good at step = 1440.

   file: ctim_199910111820.tar.gz

C================================================================
46. Time series:
      10/11/1999 18:20 Yongli Wang IGPP/UCLA

1. Some major revises the loop structure of tucan.f.

   Runs good at step = 1440.

2. Correct a small error in varcom.h about h.  
3. Partly modification of solar.f.

   Runs good at step = 1440.

   file: ctim_199910121040.tar.gz

C================================================================
47. Time series:
      10/12/1999 10:50 Yongli Wang IGPP/UCLA

1. Use half of the potential to get the electric field. 

   Runs good at step = 1200.

   file: ctim_199910121240.tar.gz

2. Major change from nn to nnloop system.
3. Delete Foster field part.
4. Delete original MHD E field part.
5. comment out MHD electric field part to test nn modification.

   Runs good at step = 1440.

6. furtehr Changing the program to make it nnloop compatible.

   Runs good at step = 780.

   file: ctim_199910121800.tar.gz

7. Add MHD electric field.

   Runs excellent at step = 1440!!!!!!!!!!!!!!

   file: ctim_199910122130.tar.gz

C================================================================
48. Time series:
      10/15/1999 21:50 Yongli Wang IGPP/UCLA

1. First finish implementation of edep.

   Runs ok at step = 1440. But there are some major problems.

   file: ctim_199910161330.tar.gz

2. Correct some errors of the program.

   Runs excellent at step = 1440!!!!!!!!!!!!!!

   file: ctim_199910162030.tar.gz

C================================================================
49. Time series:
      10/16/1999 20:30 Yongli Wang IGPP/UCLA

1. Make the program more beautiful.

   filename: ctim_199910162050.tar.gz


2. Correct a big error which is very difficult to find:
   ht(n,m+1,l) misrepresent by ht(n,m+1,1)

   Runs excellent at step = 1440!!!!!!!!!!!!!!

3. Change several small irrationals in ctim_flux().

   Runs good at step = 120.

   file: ctim_199910172045.tar.gz

4. Modify somepart of ctim_get_pole_values().

   Runs good at step = 900.

   file: ctim_199910172130.tar.gz

5. New version of documentation: ctim_3.tex

   file: ctim_199910172245.tar.gz


C================================================================
50. Time series:
      10/19/1999 16:20 Yongli Wang IGPP/UCLA

1. Delete the effect of levpi. 
2. Huge modification of output routines.

   Runs good at step = 360.

   file: ctim_199910191800.tar.gz
   
3. Deleted elfmods.f.
4. Revise output routines to show data in SM correctly.

   Runs good at step = 1440.

   file: ctim_199910201240.tar.gz

C================================================================
51. Time series:
      10/21/1999 13:15 Yongli Wang IGPP/UCLA

1. Implement a data converting routine at the beginning of tucan.f, 
   with which original CTIM data can be transformed to Jimmy's data
   format for plotting.
 
   file: ctim_199910211320.tar.gz

2. Changing certain parts of the program to output mean energy and 
   energy flux of discrete and diffuse electron, sigh,sigp,pot and 
   pacurr from MHD result. As well as flux and mean energy from CTIM.
3. Correct a problem about file saving of MHD potential and MHD FAC, 
   which include negative values.  

   Runs good at step = 1440.
    
   file: ctim_199910221320.tar.gz
   
C================================================================
52. Time series:
      11/12/1999 11:30 Yongli Wang IGPP/UCLA

    Save for run suggested by Tim!

C================================================================
53. Time series:
      11/21/1999 11:30 Yongli Wang IGPP/UCLA

1. Modernize the program and make it more likely to run with g77.

   file: ctim_199911211130.tar.gz

2. Finally solved all the problems and make it possible for the program
   to run with g77 fully without pgf77.


   Results are extremely good!
   file: ctim_199911212300.tar.gz

C================================================================
54. Time series:
      11/22/1999 10:45 Yongli Wang IGPP/UCLA

1. Include outpur of on2,t,cpcpd,hmf2 and nmf2 
   at certain location (Millstone) in the program.

   file: ctim_199911221045.tar.gz

2. Test new way of getting nmf2 and hmf2!

   Results: The profile is really so crude! This should be the reason
            why we should amend this program later!

3. Setting time related F107 input and fully implement it in the 
   program and realize full time control of this value.

   file: ctim_199911241035.tar.gz

C================================================================
55. Time series:
      11/25/1999 20:45 Yongli Wang IGPP/UCLA

1. Run from midnight Jan. 9, 1997, using initial field data from old CTIM output
   for a whole day, 24 hours. 

   file: ctim_199911252045.tar.gz

2. Change the program to use new MHD run results. Thought the new results are
   in 2 min resolution, this time we only use 4 min data. And later 
   implementation to use all the data will be conducted!

3. Using switch -O2 in Makefile's g77 compiler and hope to achieve 
   5 times faster speed!
4. Change some output file name in tucan.f and output directory for history
   file in ctim_output.f
   
?5. Combine CTIM and CTIM_MHD result to one plot.

?  file: ctim_199911271615.tar.gz

6. run for half potential b017 CTIM-MHD

   file: ctim_199911292145.tar.gz

C================================================================
56. Time series:
      12/10/1999 15:00 Yongli Wang IGPP/UCLA

1.  run for half potential b441 CTIM-MHD.
   
    result: very good!
    file: ctim_199912101500.tar.gz

C================================================================
57. Time series:
      01/05/2000 14:30 Yongli Wang IGPP/UCLA


1. Version after structure modification of the files.

   file: ctim_200001041430.tar.gz

2. Various changing to Shaun.f and molion.f.
    result: very good! But still could not run well with full MHD potential.
    file: ctim_200001081830.tar.gz

3. Massive delete unneeded variables and subroutines.
    result: very good in half potential case!
    file: ctim_200001082330.tar.gz

4. Corrected the interpols.f error. Till now the NAN error has been
   totally eliminated!

    result: very good in full potential case!
    file: ctim_200001161800.tar.gz

5. Stop working to make b441 full potential case run perfect. So far
   we have fully make the program run very well!
   file: ctim_200001171430.tar.gz

C================================================================
58. Time series:
      02/09/2000 14:45 Yongli Wang IGPP/UCLA

1. Add output of alden for DX analysis!
2. Correct the error of getting setting nn in ctim_get_var0().
