[ncl-talk] solar radiation
Barry Lynn
barry.h.lynn at gmail.com
Thu Jan 24 10:58:29 MST 2019
Please try taking SWDOWN by itself.
If you need direct normal then write again.
Barry
On Thu, Jan 24, 2019 at 2:15 PM Komal Shukla <komalshukla1992 at gmail.com>
wrote:
> 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
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
--
Barry H. Lynn, Ph.D
Senior Associate Scientist, Lecturer,
The Institute of the Earth Science,
The Hebrew University of Jerusalem,
Givat Ram, Jerusalem 91904, Israel
Tel: 972 547 231 170
Fax: (972)-25662581
C.E.O, Weather It Is, LTD
Weather and Climate Focus
http://weather-it-is.com
Jerusalem, Israel
Local: 02 930 9525
Cell: 054 7 231 170
Int-IS: x972 2 930 9525
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190124/dc4a702d/attachment.html>
More information about the ncl-talk
mailing list