[cig-commits] r6301 - mc/3D/CitcomS/trunk/lib

[tan2 at geodynamics.org]

Author: tan2
Date: 2007-03-19 16:29:00 -0700 (Mon, 19 Mar 2007)
New Revision: 6301

Modified:
   mc/3D/CitcomS/trunk/lib/Composition_related.c
Log:
Removing read_comp(), which is not used anymore

Modified: mc/3D/CitcomS/trunk/lib/Composition_related.c
===================================================================
--- mc/3D/CitcomS/trunk/lib/Composition_related.c	2007-03-19 23:25:14 UTC (rev 6300)
+++ mc/3D/CitcomS/trunk/lib/Composition_related.c	2007-03-19 23:29:00 UTC (rev 6301)
@@ -421,60 +421,3 @@
 
     return;
 }
-
-
-
-/******************* READ COMP ***********************************************/
-/*                                                                           */
-/* This function is similar to read_temp. It is used to read the composition */
-/* from file for post-proceesing.                                            */
-
-void read_comp(E)
-     struct All_variables *E;
-{
-    int i,ii,m,mm,ll;
-    char output_file[255],input_s[1000];
-
-    double g;
-    FILE *fp;
-
-
-
-    ii = E->monitor.solution_cycles;
-    sprintf(output_file,"%s.comp.%d.%d",E->control.old_P_file,E->parallel.me,ii);
-
-    if ((fp=fopen(output_file,"r"))==NULL)
-        {
-            fprintf(stderr,"ERROR(read_temp) - %s not found\n",output_file);
-            fflush(stderr);
-            exit(10);
-        }
-
-    fgets(input_s,1000,fp);
-    sscanf(input_s,"%d %d %f",&ll,&mm,&E->monitor.elapsed_time);
-
-    for(m=1;m<=E->sphere.caps_per_proc;m++)
-        {
-            E->composition.comp_node[m]=(double *)malloc((E->lmesh.nno+1)*sizeof(double));
-
-            fgets(input_s,1000,fp);
-            sscanf(input_s,"%d %d",&ll,&mm);
-            for(i=1;i<=E->lmesh.nno;i++)
-                {
-                    if (fgets(input_s,1000,fp)==NULL)
-                        {
-                            fprintf(stderr,"ERROR(read_comp) -data for node %d not found\n",i);
-                            fflush(stderr);
-                            exit(10);
-                        }
-                    sscanf(input_s,"%lf",&g);
-                    E->composition.comp_node[m][i] = g;
-
-                }
-        }
-
-    fclose (fp);
-    return;
-}
-
-