[ncl-talk] solar radiation

[Komal Shukla]

Dear NCL team,

I am trying to extract solar radiaion (w/m2) from my wrf  outputs, and
compare it with ground observation.

SWDOWN = wrf_user_getvar(aa,"SWDOWN", -1 )      ; 3D tc

GLW = wrf_user_getvar(aa,"GLW", -1 )      ; 3D tc
I am taking total of swdown and glw and comparing it with ground obs.
The total(GLW+SWDOWN = around 1000 w/m2 in summer noon) is way high than
the actual observation (around 340 w/m2 in summer noon).

My ncl script is

rh2 = wrf_user_getvar(aa,"SWDOWN", -1 )      ; 3D tc

 rh2_day=new(24,float)

 rh2_dd=new((/numFILES,24/),float)

rh2_day_ist=new(24,float) ; f

 rh2_t=rh2(:,Y1,X1)

  print(rh2_t)

asciiwrite("rh_pb_720", rh2_t)

  do id=0,(numFILES-1)

  do ih=0,23

      idh=ih+id*24

      rh2_dd(id,ih)=rh2_t(idh)

  end do

  end do

 rh2_day_m=dim_avg_n(rh2_dd,0)


Is there some other way to calculate surface solar radiation?
I'd appreciate any help.

Kim
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190124/b248bce7/attachment.html>