** Calcualte g & h based on WSO and MDI Br
gh_wbrsc.pro		  wbrsc => ghf 
gh_mbrsc.pro		  mbrsc => ghf
 zgrid.pro		Asign grid for longitude and latitude
 csmph.pro < csmphj.pro (m,phj => cmpj,smpj)
 pdp.pro   < pdpi.pro   (n,thi => Pi,dPi)
IDL> gh_mbrsc,'MDI.1968.fits',9,g9,h9,ghmf='mdigh_1968.asc',/ifl

** Calculate Br Bt Bp at a specfic point with r < Rcp
zgrid.pro	jph,ith => phd,thd,lad,cth,sth,cmp,scs
pb_hc.pro       g,h => signed pb  Calculate B between 1 and Rcp
 hchf
 csmphj		m,phj => cmpj,smpj
 gfnj		g,h,cmpj,smpj,Nmax => gfj,dgfj
 pdpi		n,thi => Pi,dPi
 brtp
IDL> zgrid,72,30,phd,thd,lad,cth,sth,cmp,scs
IDL> ppb_hc,g9,h9,2.0,cth(10),sth(10),phd(20)*!DTOR,Br,Bt,Bp

** Calculate gca & hca good for r > Rcp
pghfcsss_gaha.pro           gaha => gca,hca,ghfcsss
 zgrid.pro
 csmph.pro < csmphj.pro (m,phj => cmpj,smpj)
 pdp.pro   < pdpi.pro   (n,thi => Pi,dPi)
 sb_hc		
  pb_hc
 ghv_csss
IDL> ghfcsss_gaha,g9,h9,gca,hca,'ghfcsss_1968.asc',nmax=9

** Calculate BrBtBp for each point with  Rcp < r <= Rss
ppb_csss.pro            g,h,gc,hc => signed pb between rcp and rss
 pb_hc                    g,h => signed pb
 pb0_csss.pro             gc,hc => unsigned pb between rcp and rss
  cssshf
  csmphj
  gfnj
  pdpi
  brtp_csss
IDL> ppb_csss,g9,h9,gca,hca,5.,thd(10)*!DTOR,phd(20)*!DTOR,nmax=9

**Calculate BrBtBp on a spherical surface with Rcp < r <= Rss
psb_csss.pro               g,h;gc,hc => sb between rcp and rss
 csmph.pro
 pdp.pro
 gfn.pro
 ppb_csss
IDL> psb_csss,g9,h9,gca,hca,5,phd,thd,9,0.2,2.5,15.0,bbr,bbt,bbp