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

[wei at geodynamics.org]

Author: wei
Date: 2007-04-27 12:22:49 -0700 (Fri, 27 Apr 2007)
New Revision: 6710

Modified:
   geodyn/3D/MAG/trunk/src/cmbcoeff.f
Log:
Add some comments in the code.

Modified: geodyn/3D/MAG/trunk/src/cmbcoeff.f
===================================================================
--- geodyn/3D/MAG/trunk/src/cmbcoeff.f	2007-04-27 15:57:49 UTC (rev 6709)
+++ geodyn/3D/MAG/trunk/src/cmbcoeff.f	2007-04-27 19:22:49 UTC (rev 6710)
@@ -20,7 +20,6 @@
       dimension la(nlma),ma(nlma)
       dimension alm(lmax,lmax),blm(lmax,lmax)
       dimension glm(lmax,lmax),hlm(lmax,lmax)
-      dimension cc(2.*nlma)
 c
 c   constants
 c
@@ -31,16 +30,16 @@
 c
 c   nominal Earth parameters
 c
-      sigma=6.0e+05     ! Core electrical conductivity in S/m
-      rho=11.e+03       ! Core mean density in kg m^{-3}
-      omega=7.2924e-05  ! Rotation angular frequency
+      sigma=6.0e+05                 ! Core electrical conductivity in S/m
+      rho=11.e+03                   ! Core mean density in kg m^{-3}
+      omega=7.2924e-05              ! Rotation angular frequency
       escale=sqrt(rho*omega/sigma)  ! "Elsasser number" scaling
-      re=6.371e+06      ! Earth radius in m
-      rc=3.485e+06      ! Core radius in m
-      ri=1.215e+06      ! Inner core radius in m
-      cdepth=rc-ri      ! Depth of outer core
-      rcre=rc/re        ! core/surface radius ratio
-      rcsqinv=(cdepth/rc)**2  ! dim-less inverse squared core radius
+      re=6.371e+06                  ! Earth radius in m
+      rc=3.485e+06                  ! Core radius in m
+      ri=1.215e+06                  ! Inner core radius in m
+      cdepth=rc-ri                  ! Depth of outer core
+      rcre=rc/re                    ! core/surface radius ratio
+      rcsqinv=(cdepth/rc)**2        ! dim-less inverse squared core radius
 
 c
 c  write header to cc file
@@ -89,9 +88,10 @@
    36 continue
 c
 
-c then read the contents of the cc-file block in pairs,
-c and assign the new indices
-
+c separate the complex spherical harmonic coefficients
+c for the poloidal magnetic potential b(lm,1) into real and
+c imaginary part and assign two new indices
+c  
       do 37 lm=1,nlma
        c1=real(b(lm,1))
        c2=aimag(b(lm,1))
@@ -103,7 +103,7 @@
 
 c   Convertion starts here
 c   the following code converts between Gauss coefficients (glm, hlm)
-c       and the spherical harmonic coefficients (alm,blm)
+c   and the spherical harmonic coefficients (alm,blm)
 c   of the magnetic potential at harmonic degree l
 c   and order m used by the numerical dynamo model MAG.
 c