load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
 load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
 load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;  Model Ouputs ;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

  DATADir = "./"
  FILES = systemfunc (" ls -1 " + DATADir + "wrfout_* ")
  numFILES = dimsizes(FILES)
  print("numFILES = " + numFILES)


  ifil = 0
  a = addfile(FILES(ifil)+".nc","r")

  lat = wrf_user_getvar(a,"XLAT",0)
  lon = wrf_user_getvar(a,"XLONG",0)

 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
 ;;;Here we calculate the desired variables for every time available do loop;;;;;;
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

      it = 0 
      FirstTime = True

      do ifil = 0,numFILES-1

      a = addfile(FILES(ifil)+".nc","r")

     pres = wrf_user_getvar(a,"pres",it) 
     tk = wrf_user_getvar(a,"tk",it) 

     Rspe=287.058
     density=pres/(Rspe*tk)

     lat    = wrf_user_getvar(a,"XLAT",it)
     lon    = wrf_user_getvar(a,"XLONG",it)  
     dust1  = wrf_user_getvar(a,"DUST_1",it)        
     dust2  = wrf_user_getvar(a,"DUST_2",it)        
     dust3  = wrf_user_getvar(a,"DUST_3",it)        
     dust4  = wrf_user_getvar(a,"DUST_4",it) 
     dust5  = wrf_user_getvar(a,"DUST_5",it)   

     level=0

     dust_tot=(dust1(level,:,:)+dust2(level,:,:)+dust3(level,:,:)+dust4(level,:,:)+dust5(level,:,:))*density(level,:,:)
     dustnew = avg(dust_tot)

     dimx = dimsizes(dust_tot) 
     nlat = dimx(1)  
     mlon = dimx(0)
     npts = nlat*mlon



     if (FirstTime) then         ; NCL
     fname = "Dust_NEW.txt"
     data = new(npts,"string")
     end if                      ; NCL

     npt=-1

     do ml =0,mlon-1
      do nl=0,nlat-1
       npt = npt+1
      
     if (FirstTime) then
     data(npt) = sprinti("%05i ",(npt+1))
     data(npt) = data(npt) + sprintf("%7.3f ", lat(ml,nl))
     data(npt) = data(npt) + sprintf("%7.3f ", lon(ml,nl))
     end if

   data(npt) = data(npt) + sprintf("%7.3f ", dustnew(ml,nl))

      end do              ; do-loop lat
   end do                 ; do-loop lon

   FirstTime=False
   end do
   

  asciiwrite (fname, data)