[cig-commits] r6792 - geodyn/3D/MAG/trunk/idl
wei at geodynamics.org
wei at geodynamics.org
Mon May 7 12:26:51 PDT 2007
Author: wei
Date: 2007-05-07 12:26:51 -0700 (Mon, 07 May 2007)
New Revision: 6792
Modified:
geodyn/3D/MAG/trunk/idl/magsym.pro
Log:
Updated plot option#7 in magsym that includes a global magnetic
field plot and two polar plots with continents.
Modified: geodyn/3D/MAG/trunk/idl/magsym.pro
===================================================================
--- geodyn/3D/MAG/trunk/idl/magsym.pro 2007-05-07 17:39:48 UTC (rev 6791)
+++ geodyn/3D/MAG/trunk/idl/magsym.pro 2007-05-07 19:26:51 UTC (rev 6792)
@@ -1673,103 +1673,96 @@
; *************************************************************************
LABEL7: PRINT, "VIEWS OF RADIAL FIELD ON OUTER BOUNDARY"
- PRINT, "SELECT, Br=0; Br*cos(theta)=1:"
- READ, IBRSEL
- PRINT, "DRAW IC TANGENT CYLINDER? YES==1"
- READ,ICTC
-
- NRB=0
- LONSHIFT=0
+ PRINT, "SELECT, Br=0; Br*cos(theta)=1:"
+ READ, IBRSEL
+ PRINT, "DRAW IC TANGENT CYLINDER? YES=1"
+ READ,ICTC
+ PRINT, 'DRAW CONTINENTS? YES=1'
+ READ,ICONT
+ PRINT, 'ENTER CENTRAL LONGITUDE (-180 to +180):'
+ READ,LONSHIFT
+ PRINT, 'REVERSE POLARITY IN DISPLAY? (YES=-1)'
+ READ,IPOL
+
+ NRB=0
IPAGE=7
- BRTITLE=''
-
+ BRTITLE=''
+
LABEL7A: ERASE
- IF IBRSEL LT 1 THEN BEGIN
- BRTITLE='RADIAL FIELD ON OUTER BOUNDARY'
+ IF IBRSEL LT 1 THEN BEGIN
+ BRTITLE='RADIAL FIELD ON OUTER BOUNDARY'
WRAP(0:NPFULL-1,*)=BR(*,*,0) & WRAP(NPFULL,*)=BR(0,*,0)
ENDIF
- IF IBRSEL GE 1 THEN BEGIN
- BRTITLE='Br*cos(theta) ON OUTER BOUNDARY'
- WRAP(0:NPFULL-1,*)=BR(*,*,0)*CTHET2(*,*)
+ IF IBRSEL GE 1 THEN BEGIN
+ BRTITLE='Br*cos(theta) ON OUTER BOUNDARY'
+ WRAP(0:NPFULL-1,*)=BR(*,*,0)*CTHET2(*,*)
WRAP(NPFULL,*)=BR(0,*,0)*CTHET2(0,*)
ENDIF
-
-;MAP FIRST RADIAL FIELD
-!P.POSITION=XY1
- PTITLE='LONGITUDE 0'
- MAP_SET,0,0,0,/ORTHOGRAPHIC,/ADVANCE,/NOBORDER,/NOERASE
- IF LCOLOR EQ 1 THEN CONTOUR,WRAP,LONWRAP,LAT,LEVELS=BRV/BFAC,/OVERPLOT,$
- /CELL_FILL,C_COLORS=LCINV*VCOLORS,C_LINESTYLE=(BRV LT LIMIT)
- IF LCOLOR EQ 0 THEN CONTOUR,WRAP,LONWRAP,LAT,LEVELS=BRV/BFAC,/OVERPLOT,$
- COLOR=0,C_LINESTYLE=(BRV LT 0)
- MAP_GRID,COLOR=0,LATDEL=30,LONDEL=45,GLINETHICK=1
-; PLOT TANGENT CYLINDER INTERSECTION
- IF ICTC EQ 1 THEN BEGIN
- TCLAT=180*ACOS(R(NR-1)/R(NRB))/PI
- PLOTS,2*FINDGEN(180), REPLICATE(TCLAT,180),PSYM=3,SYMSIZE=.5,COLOR=0
- PLOTS,2*FINDGEN(180),REPLICATE(-TCLAT,180),PSYM=3,SYMSIZE=.5,COLOR=0
- ENDIF
-XYOUTS,XT1,YT1,PTITLE,CHARSIZE=CSZ*SZF,ALIGNMENT=0.5,/NORMAL,COLOR=TEXTCOLOR
+ IF IPOL EQ -1 THEN WRAP=-WRAP ;REVERSE DISPLAY POLARITY
-;MAP SECOND RADIAL FIELD
-!P.POSITION=XY2
- PTITLE='LONGITUDE 180'
- MAP_SET,0,180,0,/ORTHOGRAPHIC,/ADVANCE,/NOBORDER,/NOERASE
- IF LCOLOR EQ 1 THEN CONTOUR,WRAP,LONWRAP,LAT,LEVELS=BRV/BFAC,/OVERPLOT,$
+;GLOBAL MAP RADIAL FIELD
+!P.POSITION=XC1
+ MAP_SET,0,LONSHIFT,0,/AITOFF,/ADVANCE,/NOBORDER,/NOERASE
+ IF LCOLOR EQ 1 THEN CONTOUR,WRAP,LONWRAP,LAT,LEVELS=BRV/BFAC,/OVERPLOT,$
/CELL_FILL,C_COLORS=LCINV*VCOLORS,C_LINESTYLE=(BRV LT LIMIT)
- IF LCOLOR EQ 0 THEN CONTOUR,WRAP,LONWRAP,LAT,LEVELS=BRV/BFAC,/OVERPLOT,$
- COLOR=0,C_LINESTYLE=(BRV LT 0)
- MAP_GRID,COLOR=0,LATDEL=30,LONDEL=45,GLINETHICK=1
+ IF LCOLOR EQ 0 THEN CONTOUR,WRAP,LONWRAP,LAT,LEVELS=BRV/BFAC,/OVERPLOT,$
+ COLOR=0,C_LINESTYLE=(BRV LT 0)
+ IF (ICONT EQ 1) THEN MAP_CONTINENTS, MLINETHICK=1,COLOR=0
+ MAP_GRID,COLOR=0,LATDEL=30,LONDEL=45,GLINETHICK=1
; PLOT TANGENT CYLINDER INTERSECTION
- IF ICTC EQ 1 THEN BEGIN
- TCLAT=180*ACOS(R(NR-1)/R(NRB))/PI
- PLOTS,2*FINDGEN(180), REPLICATE(TCLAT,180),PSYM=3,SYMSIZE=.5,COLOR=0
- PLOTS,2*FINDGEN(180),REPLICATE(-TCLAT,180),PSYM=3,SYMSIZE=.5,COLOR=0
- ENDIF
-XYOUTS,XT2,YT2,PTITLE,CHARSIZE=CSZ*SZF,ALIGNMENT=0.5,/NORMAL,COLOR=TEXTCOLOR
+ IF ICTC EQ 1 THEN BEGIN
+ TCLAT=180*ACOS(R(NR-1)/R(NRB))/PI
+ PLOTS,2*FINDGEN(180), REPLICATE(TCLAT,180),PSYM=3,SYMSIZE=.5,COLOR=0
+ PLOTS,2*FINDGEN(180),REPLICATE(-TCLAT,180),PSYM=3,SYMSIZE=.5,COLOR=0
+ ENDIF
-;MAP THIRD RADIAL FIELD
-!P.POSITION=XY3
+;MAP NORTH POLAR RADIAL FIELD
+!P.POSITION=XY3
PTITLE='NORTH POLAR'
- MAP_SET,90,0,0,/ORTHOGRAPHIC,/ADVANCE,/NOBORDER,/NOERASE
- IF LCOLOR EQ 1 THEN CONTOUR,WRAP,LONWRAP,LAT,LEVELS=BRV/BFAC,/OVERPLOT,$
+ MAP_SET,90,0,0,/ORTHOGRAPHIC,/ADVANCE,/NOBORDER,/NOERASE
+ IF LCOLOR EQ 1 THEN CONTOUR,WRAP,LONWRAP,LAT,LEVELS=BRV/BFAC,/OVERPLOT,$
/CELL_FILL,C_COLORS=LCINV*VCOLORS,C_LINESTYLE=(BRV LT LIMIT)
- IF LCOLOR EQ 0 THEN CONTOUR,WRAP,LONWRAP,LAT,LEVELS=BRV/BFAC,/OVERPLOT,$
- COLOR=0,C_LINESTYLE=(BRV LT 0)
- MAP_GRID,COLOR=0,LATDEL=30,LONDEL=45,GLINETHICK=1
+ IF LCOLOR EQ 0 THEN CONTOUR,WRAP,LONWRAP,LAT,LEVELS=BRV/BFAC,/OVERPLOT,$
+ COLOR=0,C_LINESTYLE=(BRV LT 0)
+ IF (ICONT EQ 1) THEN MAP_CONTINENTS, MLINETHICK=1,COLOR=0
+ MAP_GRID,COLOR=0,LATDEL=30,LONDEL=45,GLINETHICK=1
; PLOT TANGENT CYLINDER INTERSECTION
- IF ICTC EQ 1 THEN BEGIN
- TCLAT=180*ACOS(R(NR-1)/R(NRB))/PI
- PLOTS,2*FINDGEN(180), REPLICATE(TCLAT,180),PSYM=3,SYMSIZE=.5,COLOR=0
- PLOTS,2*FINDGEN(180),REPLICATE(-TCLAT,180),PSYM=3,SYMSIZE=.5,COLOR=0
- ENDIF
+ IF ICTC EQ 1 THEN BEGIN
+ TCLAT=180*ACOS(R(NR-1)/R(NRB))/PI
+ PLOTS,2*FINDGEN(180), REPLICATE(TCLAT,180),PSYM=3,SYMSIZE=.5,COLOR=0
+ PLOTS,2*FINDGEN(180),REPLICATE(-TCLAT,180),PSYM=3,SYMSIZE=.5,COLOR=0
+ ENDIF
XYOUTS,XT3,YT3,PTITLE,CHARSIZE=CSZ*SZF,ALIGNMENT=0.5,/NORMAL,COLOR=TEXTCOLOR
-;MAP FOURTH RADIAL FIELD
-!P.POSITION=XY4
+;MAP SOUTH POLAR RADIAL FIELD
+!P.POSITION=XY4
PTITLE='SOUTH POLAR'
- MAP_SET,-90,180,0,/ORTHOGRAPHIC,/ADVANCE,/NOBORDER,/NOERASE
- IF LCOLOR EQ 1 THEN CONTOUR,WRAP,LONWRAP,LAT,LEVELS=BRV/BFAC,/OVERPLOT,$
+ MAP_SET,-90,180,0,/ORTHOGRAPHIC,/ADVANCE,/NOBORDER,/NOERASE
+ IF LCOLOR EQ 1 THEN CONTOUR,WRAP,LONWRAP,LAT,LEVELS=BRV/BFAC,/OVERPLOT,$
/CELL_FILL,C_COLORS=LCINV*VCOLORS,C_LINESTYLE=(BRV LT LIMIT)
- IF LCOLOR EQ 0 THEN CONTOUR,WRAP,LONWRAP,LAT,LEVELS=BRV/BFAC,/OVERPLOT,$
- COLOR=0,C_LINESTYLE=(BRV LT 0)
- MAP_GRID,COLOR=0,LATDEL=30,LONDEL=45,GLINETHICK=1
+ IF LCOLOR EQ 0 THEN CONTOUR,WRAP,LONWRAP,LAT,LEVELS=BRV/BFAC,/OVERPLOT,$
+ COLOR=0,C_LINESTYLE=(BRV LT 0)
+ IF (ICONT EQ 1) THEN MAP_CONTINENTS, MLINETHICK=1,COLOR=0
+ MAP_GRID,COLOR=0,LATDEL=30,LONDEL=45,GLINETHICK=1
; PLOT TANGENT CYLINDER INTERSECTION
- IF ICTC EQ 1 THEN BEGIN
- TCLAT=180*ACOS(R(NR-1)/R(NRB))/PI
- PLOTS,2*FINDGEN(180), REPLICATE(TCLAT,180),PSYM=3,SYMSIZE=.5,COLOR=0
- PLOTS,2*FINDGEN(180),REPLICATE(-TCLAT,180),PSYM=3,SYMSIZE=.5,COLOR=0
- ENDIF
+ IF ICTC EQ 1 THEN BEGIN
+ TCLAT=180*ACOS(R(NR-1)/R(NRB))/PI
+ PLOTS,2*FINDGEN(180), REPLICATE(TCLAT,180),PSYM=3,SYMSIZE=.5,COLOR=0
+ PLOTS,2*FINDGEN(180),REPLICATE(-TCLAT,180),PSYM=3,SYMSIZE=.5,COLOR=0
+ ENDIF
XYOUTS,XT4,YT4,PTITLE,CHARSIZE=CSZ*SZF,ALIGNMENT=0.5,/NORMAL,COLOR=BWCOL
XYOUTS,0.5,22.7/YDIM,BRTITLE,/NORM,SIZE=CSB*SZF,$
- ALIGNMENT=0.5,COLOR=BWCOL
- XYOUTS,0.5,22.0/YDIM,ANNTEXT3,/NORM,SIZE=CSZ*SZF,ALIGNMENT=0.5,COLOR=TEXTCOLOR
+ ALIGNMENT=0.5,COLOR=BWCOL
- GOTO,LABEL0
+XYOUTS,0.5,22.0/YDIM,ANNTEXT3,/NORM,SIZE=CSZ*SZF,ALIGNMENT=0.5,COLOR=TEXTCOLOR
+ GOTO,LABEL0
+;end of radial cmb field maps option
; *************************************************************************
+
+; *************************************************************************
; *************** PLOTTING OF ZONAL PROFILES *********************
; *************************************************************************
More information about the cig-commits
mailing list