[cig-commits] r6036 - geodyn/3D/MAG/trunk/idl

[wei at geodynamics.org]

Author: wei
Date: 2007-02-15 21:47:35 -0800 (Thu, 15 Feb 2007)
New Revision: 6036

Modified:
   geodyn/3D/MAG/trunk/idl/magmovieCIG.pro
Log:
Update IDL move procedure for a full-sphere view, added continents on map, changed display window with 2x1 aspect ratio


Modified: geodyn/3D/MAG/trunk/idl/magmovieCIG.pro
===================================================================
--- geodyn/3D/MAG/trunk/idl/magmovieCIG.pro	2007-02-16 02:16:14 UTC (rev 6035)
+++ geodyn/3D/MAG/trunk/idl/magmovieCIG.pro	2007-02-16 05:47:35 UTC (rev 6036)
@@ -1,7 +1,7 @@
 ;A PROCEDURE TO CREATE A BYTE ARRAY FOR MAKING A MOVIE OF
 ;MAGNETOCONVECTION OR DYNAMOS USING RESULTS FROM PROGRAM mag
 
-;THIS PROCEDURE IS MEANT TO DISPLAY RESULTS AT MID-DEPTH OR THE OUTER SURFACE
+;THIS PROCEDURE IS MEANT TO DISPLAY RESULTS AT THE OUTER SURFACE
 ;IN MAP VIEW, ASSUMING SOME LONGITUDINAL SYMMETRY.
 
 ;THIS PROCEDURE ALSO CREATES JPG FILES OF THE IMAGES, FOR WEB MOVIES.
@@ -59,8 +59,8 @@
 ;CHOOSE A VARIABLE
          CHOICE=1
 ;          IF IOP GT 1 THEN BEGIN
-         PRINT,'SELECT DISPLAY (Br surface=1, Br middepth=2):'
-	READ,CHOICE
+;         PRINT,'SELECT DISPLAY (Br surface=1, Br middepth=2):'
+;	READ,CHOICE
          PRINT,'Select Polarity (Reverse=1; Normal=-1):'
 	READ,POLAR
 ;           ENDIF
@@ -89,13 +89,13 @@
 	 WXSIZE=300
 	PRINT,'WINDOW SIZE IN PIXEL?  DEFAULT=300'
          READ,WXSIZE
-         WYSIZE=WXSIZE
+         WYSIZE=0.5*WXSIZE
 
 ;DIMENSION THE BYTE ARRAY
 	FRAMES=BYTARR(WXSIZE,WYSIZE,NFR)
 
 ;CREATE A WINDOW
-	WINDOW,2,XSIZE=WXSIZE,YSIZE=2*WYSIZE
+	WINDOW,2,XSIZE=WXSIZE,YSIZE=WYSIZE
 
 ;READ IN COLATITUDES
 	READF,1,THET
@@ -245,8 +245,12 @@
           ENDIF
 
 ; DEFINE MAPPING PROJECTIONS
+  LONGSHIFT=0
           CASE CHOICE OF
-          1: MAP_SET,INCL,LONSHIFT,0,/ORTHOGRAPHIC,TITLE='CMB RADIAL FIELD'
+          1: BEGIN
+              MAP_SET,INCL,LONGSHIFT,0,/AITOFF, /ADVANCE, $
+			TITLE='CMB RADIAL FIELD',/NOBORDER,COLOR=TEXTCOLOR
+             END
           2: MAP_SET,INCL,LONSHIFT,0,/ORTHOGRAPHIC,TITLE='MIDDEPTH RADIAL FIELD'
           3: MAP_SET,INCL,LONSHIFT,0,/ORTHOGRAPHIC,TITLE='MIDDEPTH RAD. VELOCITY'
           ENDCASE
@@ -259,7 +263,8 @@
             ERASE
 
  
-CONTOUR,POLAR*VI,LON,LAT,LEVELS=VV,/OVERPLOT,C_COLORS=VCOLORS,/CELL_FILL
+
+          CONTOUR,POLAR*VI,LON,LAT,LEVELS=VV,/OVERPLOT,C_COLORS=VCOLORS,/CELL_FILL
             MAP_GRID,COLOR=0,GLINESTYLE=1,GLINETHICK=1
 ; PLOT TANGENT CYLINDER INTERSECTION
 	 IF ICTC EQ 1 THEN BEGIN
@@ -293,6 +298,7 @@
           ERASE
           CONTOUR,POLAR*V,LON,LAT,LEVELS=VV,/OVERPLOT,C_COLORS=VCOLORS,/CELL_FILL
           MAP_GRID,COLOR=0,GLINESTYLE=1,GLINETHICK=1
+	  MAP_CONTINENTS, MLINETHICK=2,COLOR=0  
 
 ; PLOT TANGENT CYLINDER INTERSECTION
 	 IF ICTC EQ 1 THEN BEGIN