load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

begin

a = addfile("./Apparent_heat_sink.nc","r")

q1 = a->q1

Q1 = a->Q1

lat = a->lat

lon = a->lon

level = a->level

time = a->time

Q1!0 = "time"

Q1&time = time

Q1!1 = "lat"

Q1&lat = lat

Q1&lat@units = "degrees_north"

Q1!2 = "lon"

Q1&lon = lon

Q1&lon@units = "degrees_east"

Q1mo6 = Q1(5::12,:,:)

Q1mo7 = Q1(6::12,:,:)

Q1mo8 = Q1(7::12,:,:)

Q2avg = (Q1mo6+Q1mo7+Q1mo8)/3.0

QQ = dim_avg_n_Wrap(Q1,0)

Q1avg = dim_avg_n_Wrap(Q2avg,0)

copy_VarCoords(QQ,Q1avg)

printVarSummary(Q1avg)

wks = gsn_open_wks("png","Q1")

res = True

res@gsnFrame = False

res@gsnDraw = False

mpres = res

mpres@mpFillOn = False 

mpres@mpOutlineOn = True

mpres@mpDataSetName = "Earth..4"

mpres@mpDataBaseVersion = "MediumRes"

; mpres@mpOutlineSpecifiers = "China:states"

mpres@mpGeophysicalLineThicknessF = 0.5

mpres@mpNationalLineThicknessF = 0.5

; mpres@mpNationalLineColor = "transparent"

mpres@mpProvincialLineColor = "transparent"

mpres@mpProvincialLineThicknessF = 0.001

; mpres@mpGeophysicalLineColor = "transparent"

mpres@mpOutlineSpecifiers = (/"China:states","Taiwan","Disputed area between India and China","India:Arunachal Pradesh"/) 

mpres@mpFillDrawOrder = "PostDraw"

        mpres@mpMaskAreaSpecifiers = (/"China:states","Taiwan","Disputed area between India and China","India:Arunachal Pradesh"/)

        mpres@mpOceanFillColor = 0

        mpres@mpInlandWaterFillColor = 0

; mpres@mpProvincialLineColor = "White"

mpres@mpProvincialLineThicknessF = 1

; mpres@mpMinLatF = 15

; mpres@mpMaxLatF = 40

; mpres@mpMinLonF = 90

; mpres@mpMaxLonF = 125 

mpres@tmXTOn = False

mpres@tmYROn = False

mpres@tmXBLabelFontHeightF = 0.015

mpres@tmYLLabelFontHeightF = 0.015

mpres@pmTickMarkDisplayMode = "Always"

; mpres@mpLandFillColor = "white"

cnres = res

cnres@cnFillOn = True

cnres@cnLinesOn = False

map = gsn_csm_map(wks,mpres)

contour = gsn_csm_contour(wks,Q1avg,cnres)

overlay(map,contour)

draw(map)

frame(wks)

end