[cig-commits] r20556 - mc/1D/hc/trunk

[becker at geodynamics.org]

Author: becker
Date: 2012-08-05 06:40:23 -0700 (Sun, 05 Aug 2012)
New Revision: 20556

Modified:
   mc/1D/hc/trunk/hc_propagator.c
   mc/1D/hc/trunk/sh_syn.c
Log:
Improved geographic region precisions for sh_syn



Modified: mc/1D/hc/trunk/hc_propagator.c
===================================================================
--- mc/1D/hc/trunk/hc_propagator.c	2012-08-03 22:16:21 UTC (rev 20555)
+++ mc/1D/hc/trunk/hc_propagator.c	2012-08-05 13:40:23 UTC (rev 20556)
@@ -20,8 +20,8 @@
   //
   HC_HIGH_PREC den1,den2,f[4],r,rlm1,rlp1,rmlm2,rml,v2,rs;
   long int np[4][4][4];
-  int lp1,lp2,lp3,lm1,lm2,lpp,lmm,l2p3,l2p1,l2m1,lltp1,lltp2;
-  int i,j,k,os1,os2;
+  long int lp1,lp2,lp3,lm1,lm2,lpp,lmm,l2p3,l2p1,l2m1,lltp1,lltp2;
+  long int i,j,k,os1,os2;
   //fprintf(stderr,"hc_evalpa: %i %g %g %g\n",l,r1,r2,visc);
   //
   //    PASSED PARAMETERS:  L: DEGREE,

Modified: mc/1D/hc/trunk/sh_syn.c
===================================================================
--- mc/1D/hc/trunk/sh_syn.c	2012-08-03 22:16:21 UTC (rev 20555)
+++ mc/1D/hc/trunk/sh_syn.c	2012-08-05 13:40:23 UTC (rev 20556)
@@ -119,9 +119,9 @@
       }
       /*  */
       dphi = DEG2RAD(dx);
-      nphi = DEG2RAD(e-w)/dphi + 1;
+      nphi = ((e-w)/dx) + 1;
       dtheta = DEG2RAD(dy);
-      ntheta = DEG2RAD(n-s)/dtheta + 1;
+      ntheta = ((n-s)/dy) + 1;
       npoints = nphi * ntheta;
 
       /*  */