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"
VALUE=new((/31,17,174,168/),float)
LUE=new((/31,17,174,168/),float)
VALUE2=new((/174,168/),float)
VALUE3=new((/174,168/),float)
MAX1=new((/31/),float)
MAX2=new((/31/),float)
wrffile=addfile("/data1/loucx/COPY/WRF/WRF/run/wrfout_d03_2020-08-01_00:00:00","r")
lat2d=wrffile->XLAT(0,3:176,3:170)
lon2d=wrffile->XLONG(0,3:176,3:170)
dims = dimsizes(lat2d)
DataDir1 = "/data1/loucx/COPY/CNAQ/data1/cctm/"
DataDir2 = "/data1/loucx/COPY/CNAQ/data/cctm/"
flist1 = systemfunc("ls "+DataDir1+"ACONC.*_d03")
flist2 = systemfunc("ls "+DataDir2+"ACONC.*