;
; PURPOSE: Calculate g & h good for r >= Rcp
; INPUT: g & h obtained from running gh_mbrsc
; OUTPUT: gc & hc
; NOTE: gc and hc in mT for calculating Br beyond Rcp
; WRITTEN: 13JAN2003 Xuepu Zhao
; MODIFIED: 08DEC2005 Xuepu Zhao
pro pghfcsss_gaha,ga,ha,gca,hca,ghcf,NMAX=nmax,RCP=rcp,RSS=rss,IPH=iph,JTH=jth

    if N_params() lt 1 then begin
       print,'pghfcsss,ga,ha,gca,hca,ghcf,nmax=,rcp=,rss=,iph=,jth='
       print,'DEFAULT: nmax=15, apar=0.2,rcp=2.5, rss=15, iph=72, jth=30'
       return
    endif
    if not keyword_set(nmax) then nmax=15
    if not keyword_set(apar) then apar=0.2
    if not keyword_set(rcp) then rcp=2.5
    if not keyword_set(rss) then rss=15.0
    if not keyword_set(iph) then iph=72
    if not keyword_set(jth) then jth=30
; Calculation
    pzgrid,iph,jth,phd,thd,lad,cth,sth
    csmph,phd,nmax,cmph,smph
    pdp,cth,nmax,P,dP
    g=100*ga 
    h=100*ha                         ; from G to mT
    psb_hc,nmax,apar,rcp,jth-1,iph-1,sth,cmph,smph,g,h,P,dP,brcp,btcp,bpcp
    ghv_csss,apar,nmax,rcp,rss,jth-1,iph-1,sth,cmph,smph,P,dP,$
             brcp,btcp,bpcp,ghcf,gca,hca

end