[cig-commits] r4687 - geodyn/3D/MAG/trunk/idl
wei at geodynamics.org
wei at geodynamics.org
Tue Oct 3 17:27:23 PDT 2006
Author: wei
Date: 2006-10-03 17:27:23 -0700 (Tue, 03 Oct 2006)
New Revision: 4687
Modified:
geodyn/3D/MAG/trunk/idl/magsym.pro
Log:
Added more plot options in magsym
Modified: geodyn/3D/MAG/trunk/idl/magsym.pro
===================================================================
--- geodyn/3D/MAG/trunk/idl/magsym.pro 2006-10-04 00:22:06 UTC (rev 4686)
+++ geodyn/3D/MAG/trunk/idl/magsym.pro 2006-10-04 00:27:23 UTC (rev 4687)
@@ -749,9 +749,8 @@
XYOUTS,XT1,YT1,PTITLE,CHARSIZE=CSZ*SZF,ALIGNMENT=0.5,/NORMAL,COLOR=TEXTCOLOR
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;MAKE A COLOR SCALE BAR & PRINT VARIABLE NAME
+;MAKE A COLOR SCALE BAR FOR RADIAL FIELD
IF LCOLOR GE 1 THEN BEGIN
-; XBAR1=0.36 & DXBAR=5/((NLV-1)*XDIM) & YBAR1=0.45
XBAR1=0.115 & DXBAR=5/((NLV-1)*XDIM) & YBAR1=0.46
SCALE1=STRING(BRV(0)/BFAC) & SCALE2=STRING(BRV(NLV-1)/BFAC)
SCALE1=STRMID(SCALE1,2,7) & SCALE2=STRMID(SCALE2,2,7)
@@ -765,8 +764,6 @@
XBAR1=XBAR1 + DXBAR
ENDFOR
XYOUTS,XBAR1-3*DXBAR,YBAR1-DXBAR,SCALE2,COLOR=TEXTCOLOR,SIZE=1.5*SCWINDOW,/NORM
-; XYOUTS,0.5,YBAR1-DXBAR,'mT',COLOR=TEXTCOLOR,/NORM,SIZE=1.5*SCWINDOW,ALIGNMENT=0.5
-; XYOUTS,0.5,YBAR2+DXBAR/2,'mt',COLOR=TEXTCOLOR,/NORM,SIZE=1.25*SCWINDOW,ALIGNMENT=0.5
ENDIF
IF LCOLOR EQ 0 THEN BEGIN
XYOUTS,0.5,.47,'PTITLEB',COLOR=TEXTCOLOR,/NORM,SIZE=1.25*SCWINDOW,$
@@ -796,7 +793,7 @@
ENDIF
XYOUTS,XT2,YT2,PTITLE,CHARSIZE=CSZ*SZF,ALIGNMENT=0.5,/NORMAL,COLOR=TEXTCOLOR
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;MAKE A COLOR SCALE BAR & PRINT VARIABLE NAME
+;MAKE A COLOR SCALE BAR FOR VELOCITY
IF LCOLOR GE 1 THEN BEGIN
XBAR1=0.6 & DXBAR=5/((NLV-1)*XDIM) & YBAR1=0.46
SCALE1=STRING(VV(0)/VFAC) & SCALE2=STRING(VV(NLV-1)/VFAC)
@@ -840,7 +837,7 @@
ENDIF
XYOUTS,XT2,YT2,PTITLE,CHARSIZE=CSZ*SZF,ALIGNMENT=0.5,/NORMAL,COLOR=TEXTCOLOR
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;MAKE A COLOR SCALE BAR & PRINT VARIABLE NAME
+;MAKE A COLOR SCALE BAR FOR UPWELLING
IF LCOLOR GE 1 THEN BEGIN
XBAR1=0.6 & DXBAR=5/((NLV-1)*XDIM) & YBAR1=0.46
SCALE1=STRING(UPV(0)/UPFAC) & SCALE2=STRING(UPV(NLV-1)/UPFAC)
@@ -884,7 +881,7 @@
MAP_GRID,COLOR=0,LATDEL=30,LONDEL=45,GLINETHICK=1
XYOUTS,XT3,YT3,PTITLE,CHARSIZE=CSZ*SZF,ALIGNMENT=0.5,/NORMAL,COLOR=TEXTCOLOR
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;MAKE A COLOR SCALE BAR & PRINT VARIABLE NAME
+;MAKE A COLOR SCALE BAR FOR TEMPERATURE
IF LCOLOR GE 1 THEN BEGIN
XBAR1=0.115 & DXBAR=5/((NLV-1)*XDIM) & YBAR1=0.05
SCALE1=STRING(TV(0)/TFAC) & SCALE2=STRING(TV(NLV-1)/TFAC)
@@ -921,7 +918,7 @@
XYOUTS,XT3,YT3,PTITLE,CHARSIZE=CSZ*SZF,ALIGNMENT=0.5,$
/NORMAL,COLOR=TEXTCOLOR
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;MAKE A COLOR SCALE BAR & PRINT VARIABLE NAME
+;MAKE A COLOR SCALE BAR FOR MAGNETIC FIELD
IF LCOLOR GE 1 THEN BEGIN
XBAR1=0.6 & DXBAR=5/((NLV-1)*XDIM) & YBAR1=0.16
SCALE1=STRING(BRV(0)/BFAC) & SCALE2=STRING(BRV(NLV-1)/BFAC)
@@ -965,7 +962,7 @@
MAP_GRID,COLOR=0,LATDEL=30,LONDEL=45,GLINETHICK=1
XYOUTS,XT4,YT4,'HELICITY',CHARSIZE=CSZ*SZF,ALIGNMENT=0.5,/NORMAL,COLOR=TEXTCOLOR
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;MAKE A COLOR SCALE BAR & PRINT VARIABLE NAME
+;MAKE A COLOR SCALE BAR FOR HELICITY
IF LCOLOR GE 1 THEN BEGIN
XBAR1=0.6 & DXBAR=5/((NLV-1)*XDIM) & YBAR1=0.05
SCALE1=STRING(HV(0)/HFAC) & SCALE2=STRING(HV(NLV-1)/HFAC)
@@ -1000,7 +997,7 @@
XYOUTS,XT4,YT4,'VORTICITY IMAGE',CHARSIZE=CSZ*SZF,ALIGNMENT=0.5,$
/NORMAL,COLOR=TEXTCOLOR
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;MAKE A COLOR SCALE BAR & PRINT VARIABLE NAME
+;MAKE A COLOR SCALE BAR FOR VORTICITY
IF LCOLOR GE 1 THEN BEGIN
XBAR1=0.115 & DXBAR=5/((NLV-1)*XDIM) & YBAR1=0.16
SCALE1=STRING(WV(0)/WFAC) & SCALE2=STRING(WV(NLV-1)/WFAC)
@@ -1087,9 +1084,9 @@
XYOUTS,9.25/XDIM,20.3/YDIM,'VELOCITY r='+SNRA,SIZE=CSZ*SZF,/NORMAL,$
ALIGNMENT=0.5,COLOR=TEXTCOLOR
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;MAKE A COLOR SCALE BAR & PRINT VARIABLE NAME
+;MAKE A COLOR SCALE BAR FOR UPWELLING VELOCITY
IF LCOLOR GE 1 THEN BEGIN
- XBAR1=0.115 & DXBAR=5/((NLV-1)*XDIM) & YBAR1=0.45
+ XBAR1=0.6 & DXBAR=5/((NLV-1)*XDIM) & YBAR1=0.45
SCALE1=STRING(UPV(0)/UPFAC) & SCALE2=STRING(UPV(NLV-1)/UPFAC)
SCALE1=STRMID(SCALE1,2,7) & SCALE2=STRMID(SCALE2,2,7)
XYOUTS,XBAR1-2*DXBAR,YBAR1-DXBAR,SCALE1,COLOR=TEXTCOLOR,SIZE=1.5*SCWINDOW,/NORM
@@ -1137,9 +1134,9 @@
XYOUTS,9.25/XDIM,9.8/YDIM,'MAGNETIC FIELD r='+SNBR,SIZE=CSZ*SZF,$
/NORMAL,ALIGNMENT=0.5,COLOR=TEXTCOLOR
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
-;MAKE A COLOR SCALE BAR & PRINT VARIABLE NAME
+;MAKE A COLOR SCALE BAR FOR RADIAL MAGNETIC FIELD
IF LCOLOR GE 1 THEN BEGIN
- XBAR1=0.115 & DXBAR=5/((NLV-1)*XDIM) & YBAR1=0.025
+ XBAR1=0.6 & DXBAR=5/((NLV-1)*XDIM) & YBAR1=0.025
SCALE1=STRING(BRV(0)/BFAC) & SCALE2=STRING(BRV(NLV-1)/BFAC)
SCALE1=STRMID(SCALE1,2,7) & SCALE2=STRMID(SCALE2,2,7)
XYOUTS,XBAR1-2*DXBAR,YBAR1-DXBAR,SCALE1,COLOR=TEXTCOLOR,SIZE=1.5*SCWINDOW,/NORM
@@ -1205,6 +1202,28 @@
XYOUTS,XT1,YT1,'TEMPERATURE',CHARSIZE=CSZ*SZF,ALIGNMENT=0.5,$
/NORMAL,COLOR=TEXTCOLOR
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;MAKE A COLOR SCALE BAR FOR TEMPERATURE
+ IF LCOLOR GE 1 THEN BEGIN
+ XBAR1=0.115 & DXBAR=5/((NLV-1)*XDIM) & YBAR1=0.46
+ SCALE1=STRING(TV(0)/TFAC) & SCALE2=STRING(TV(NLV-1)/TFAC)
+ SCALE1=STRMID(SCALE1,2,7) & SCALE2=STRMID(SCALE2,2,7)
+ XYOUTS,XBAR1-2*DXBAR,YBAR1-DXBAR,SCALE1,COLOR=TEXTCOLOR,SIZE=1.5*SCWINDOW,/NORM
+
+ FOR L=0, NLV-2 DO BEGIN
+ XBAR2=XBAR1 + DXBAR & YBAR2=YBAR1 + DXBAR
+ XBAR=[XBAR1,XBAR2,XBAR2,XBAR1,XBAR1] & YBAR=[YBAR1,YBAR1,YBAR2,YBAR2,YBAR1]
+ POLYFILL,XBAR,YBAR,COLOR=VCOLORS(L),/NORM
+ ;OPLOT,XBAR,YBAR,COLOR=TEXTCOLOR,LINESTYLE=0,/NOCLIP ;ADD BOUNDARIES
+ XBAR1=XBAR1 + DXBAR
+ ENDFOR
+ XYOUTS,XBAR1-3*DXBAR,YBAR1-DXBAR,SCALE2,COLOR=TEXTCOLOR,SIZE=1.5*SCWINDOW,/NORM
+ ENDIF
+ IF LCOLOR EQ 0 THEN BEGIN
+ XYOUTS,0.5,.47,'PTITLEB',COLOR=TEXTCOLOR,/NORM,SIZE=1.25*SCWINDOW,$
+ ALIGNMENT=0.5
+ ENDIF
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
GOTO, LABELCOV
L2: !P.POSITION=XY2 & IFR=2
@@ -1275,6 +1294,28 @@
XYOUTS,XT3,YT3,'Z-VORTICITY',CHARSIZE=CSZ*SZF,ALIGNMENT=0.5,$
/NORMAL,COLOR=TEXTCOLOR
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;MAKE A COLOR SCALE BAR FOR EQUATORIAL W_Z
+ IF LCOLOR GE 1 THEN BEGIN
+ XBAR1=0.115 & DXBAR=5/((NLV-1)*XDIM) & YBAR1=0.05
+ SCALE1=STRING(WV(0)/WFAC) & SCALE2=STRING(WV(NLV-1)/WFAC)
+ SCALE1=STRMID(SCALE1,2,7) & SCALE2=STRMID(SCALE2,2,7)
+ XYOUTS,XBAR1-2*DXBAR,YBAR1-DXBAR,SCALE1,COLOR=TEXTCOLOR,SIZE=1.5*SCWINDOW,/NORM
+
+ FOR L=0, NLV-2 DO BEGIN
+ XBAR2=XBAR1 + DXBAR & YBAR2=YBAR1 + DXBAR
+ XBAR=[XBAR1,XBAR2,XBAR2,XBAR1,XBAR1] & YBAR=[YBAR1,YBAR1,YBAR2,YBAR2,YBAR1]
+ POLYFILL,XBAR,YBAR,COLOR=VCOLORS(L),/NORM
+ ;OPLOT,XBAR,YBAR,COLOR=TEXTCOLOR,LINESTYLE=0,/NOCLIP ;ADD BOUNDARIES
+ XBAR1=XBAR1 + DXBAR
+ ENDFOR
+ XYOUTS,XBAR1-3*DXBAR,YBAR1-DXBAR,SCALE2,COLOR=TEXTCOLOR,SIZE=1.5*SCWINDOW,/NORM
+ ENDIF
+ IF LCOLOR EQ 0 THEN BEGIN
+ XYOUTS,0.5,.47,'PTITLEB',COLOR=TEXTCOLOR,/NORM,SIZE=1.25*SCWINDOW,$
+ ALIGNMENT=0.5
+ ENDIF
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
GOTO, LABELCOV
@@ -1288,6 +1329,29 @@
XYOUTS,XT4,YT4,'Z-FIELD',CHARSIZE=CSZ*SZF,ALIGNMENT=0.5,$
/NORMAL,COLOR=TEXTCOLOR
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+;MAKE A COLOR SCALE BAR FOR Z_FIELD
+ IF LCOLOR GE 1 THEN BEGIN
+ XBAR1=0.6 & DXBAR=5/((NLV-1)*XDIM) & YBAR1=0.05
+ SCALE1=STRING(BRV(0)/BFAC) & SCALE2=STRING(BRV(NLV-1)/BFAC)
+ SCALE1=STRMID(SCALE1,2,7) & SCALE2=STRMID(SCALE2,2,7)
+ XYOUTS,XBAR1-2*DXBAR,YBAR1-DXBAR,SCALE1,COLOR=TEXTCOLOR,SIZE=1.5*SCWINDOW,/NORM
+
+ FOR L=0, NLV-2 DO BEGIN
+ XBAR2=XBAR1 + DXBAR & YBAR2=YBAR1 + DXBAR
+ XBAR=[XBAR1,XBAR2,XBAR2,XBAR1,XBAR1] & YBAR=[YBAR1,YBAR1,YBAR2,YBAR2,YBAR1]
+ POLYFILL,XBAR,YBAR,COLOR=VCOLORS(L),/NORM
+ ;OPLOT,XBAR,YBAR,COLOR=TEXTCOLOR,LINESTYLE=0,/NOCLIP ;ADD BOUNDARIES
+ XBAR1=XBAR1 + DXBAR
+ ENDFOR
+ XYOUTS,XBAR1-3*DXBAR,YBAR1-DXBAR,SCALE2,COLOR=TEXTCOLOR,SIZE=1.5*SCWINDOW,/NORM
+ ENDIF
+ IF LCOLOR EQ 0 THEN BEGIN
+ XYOUTS,0.5,.47,'PTITLEB',COLOR=TEXTCOLOR,/NORM,SIZE=1.25*SCWINDOW,$
+ ALIGNMENT=0.5
+ ENDIF
+;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
+
; PLOT INNER AND OUTER CIRCLE AND COVER INNER CORE
LABELCOV: OPLOT,/POLAR,REQ(*,0),PHIEQ(*,0),COLOR=TEXTCOLOR ; DRAW OUTER EQUATOR
More information about the cig-commits
mailing list