[cig-commits] r6793 - geodyn/3D/MAG/trunk/src

[wei at geodynamics.org]

Author: wei
Date: 2007-05-07 12:32:23 -0700 (Mon, 07 May 2007)
New Revision: 6793

Modified:
   geodyn/3D/MAG/trunk/src/cmbcoeff.f
Log:
cleanup the test code in cmbcoeff

Modified: geodyn/3D/MAG/trunk/src/cmbcoeff.f
===================================================================
--- geodyn/3D/MAG/trunk/src/cmbcoeff.f	2007-05-07 19:26:51 UTC (rev 6792)
+++ geodyn/3D/MAG/trunk/src/cmbcoeff.f	2007-05-07 19:32:23 UTC (rev 6793)
@@ -150,8 +150,8 @@
 c   define MAG harmonic-Gauss harmonic conversion factors
 c   conalm,conblm
 c
-      conalm=fact1
-      conblm=-fact1
+      conalm=real(fact1)
+      conblm=-real(fact1)
 
 c      if (l .le. 0 .or. m .lt. 0 .or. m .gt. l) then
 c       write(6,'(''bad l or m in getgauss'')')
@@ -166,11 +166,11 @@
 
 c           return
 c      else ! form dimensionless fully normalized potential coeffs
-        aalm(l,m)=glm(l,m)/(anano*escale*fact2*conalm)
-        ablm(l,m)=hlm(l,m)/(anano*escale*fact2*conblm)
+c        aalm(l,m)=glm(l,m)/(anano*escale*fact2*conalm)
+c        ablm(l,m)=hlm(l,m)/(anano*escale*fact2*conblm)
 
-      write(21,2106) l,m,aalm(l,m),ablm(l,m),glm(l,m),hlm(l,m)
- 2206 format(/,2x,2i3,2x,2(f15.5),2x,2(f15.5))
+c      write(21,2106) l,m,aalm(l,m),ablm(l,m),glm(l,m),hlm(l,m)
+c 2206 format(/,2x,2i3,2x,2(f15.5),2x,2(f15.5))
 
 c      return
 c      endif
@@ -189,8 +189,8 @@
       do 41 j=0,i
          write(22,2201) glm(i,j)
          write(22,2202) hlm(i,j)
- 2201 format(2x,f15.5)
- 2202 format(4x,f15.5)
+ 2201 format(2x,f10.5)
+ 2202 format(4x,f10.5)
    41 continue
    43 continue