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" load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl" begin ;fg_file="/data/alim/projects/osse/geo_hyper/WRF_fc/bkg/wrfout_d01_2008052400_bkg_nmm.grb" fg_file="wrfout_d01_2008052400_bkg_nmm.grb" fg =addfile(fg_file,"r") fg_T =fg->TMP_GDS0_HYBL printVarSummary(fg_T) printMinMax(fg_T,0) ; all levels ;;; lat2d=fg->g0_lat_1 ; No! g0_lat_1 is one dimensional ;;; lon2d=fg->g0_lon_2 ; " g0_lon_2 " " " inc_T = fg_T ;CTRL inc_T@long_name="inc_T" inc_T@units="K" plot=new(1, graphic) i=98 j=i-3 k=i-6 l=i-9 print("===========================") printMinMax(fg_T(i,:,:),0) print("===========================") print(fg_T&g0_lon_2+" "+fg_T(i,122,:)) ; all lon pts near equator print("===========================") nmsg = num(ismissing(fg_T(i,:,:))) print("i="+i+" nmsg="+nmsg) print("fg_T(i,:,:) i="+i+" nmsg="+nmsg) kmsg = num(ismissing(inc_T(i,:,:))) print("i="+i+" kmsg="+kmsg) print("inc_T(i,:,:) i="+i+" kmsg="+kmsg) print("===========================") wks = gsn_open_wks("pdf" ,"OSSE_Global_T_outer_loop_2_"+i+"_"+j+"_"+k+"_"+l+"_tm00") res = True ; plot mods desired ; res@gsnDraw=False ; res@gsnFrame=False res@gsnMaximize = True res@gsnSpreadColors = True res@mpLandFillColor = -1 ; set land to be transperant res@mpOutlineOn = True res@mpOutlineBoundarySets = "AllBoundaries" res@mpCenterLonF = 0.0 ; default is 0.0 ;;; res@tfDoNDCOverlay = True ; not appropriate for _ce res@cnFillOn = True res@cnLinesOn = False res@cnLineLabelsOn = False res@lbBoxLinesOn = True res@gsnAddCyclic = False ; data has cyclic point res@tmXBLabelFontHeightF = 0.01 ; resize tick labels res@tmYLLabelFontHeightF = 0.01 res@pmTickMarkDisplayMode = "Always" res@lbLabelAutoStride = True res@lbLabelBarOn = True res@lbLabelFontHeightF = 0.015 res@cnMissingValFillColor = "background" ; "green", "black" ;;; inc_T@lat2d = lat2d ; No! See above ;;; inc_T@lon2d = lon2d plot=gsn_csm_contour_map_ce(wks,inc_T(i,:,:), res) ; CTRL end