[ncl-talk] ncl-talk Digest, Vol 216, Issue 16
R. Lucas
rousslucas at gmail.com
Sat Nov 27 02:03:05 MST 2021
Thank you dennis I succeeded, thanks to your advice
Le jeu. 25 nov. 2021 à 20:00, <ncl-talk-request at mailman.ucar.edu> a écrit :
> Send ncl-talk mailing list submissions to
> ncl-talk at mailman.ucar.edu
>
> To subscribe or unsubscribe via the World Wide Web, visit
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
> or, via email, send a message with subject or body 'help' to
> ncl-talk-request at mailman.ucar.edu
>
> You can reach the person managing the list at
> ncl-talk-owner at mailman.ucar.edu
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of ncl-talk digest..."
>
>
> Today's Topics:
>
> 1. Variable size (R. Lucas)
> 2. Re: Variable size (Dennis Shea)
> 3. Re: NCL function wrf_user_getvar (Israt Jahan)
> 4. Re: NCL function wrf_user_getvar (Dennis Shea)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Thu, 25 Nov 2021 10:26:02 +0100
> From: "R. Lucas" <rousslucas at gmail.com>
> To: ncl-talk at mailman.ucar.edu
> Subject: [ncl-talk] Variable size
> Message-ID:
> <
> CAEcwEprvOmxHLTF55Ki4UsD85FgxNkgMbmTcpWNcdJvJHSpv5g at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Dear NCL users
> i have 2 problems :
>
> 1) Is this varibable level looks like 10m above ground ?
> 2) This variable is 3d. How can I remove first dimension, somthing like
> (0,:,:)
> because I have a fatal error :
> fatal:["Execute.c":7925]:Number of subscripts on right-hand-side do not
> match
> number of dimensions of variable: (3), Subscripts used: (1)
>
> Thanks in advance.
> Rouss.
>
>
> float UGRD_P0_L103_GLC0 ( lv_HTGL3, ygrid_0, xgrid_0 )
> center : US National Weather Service - NCEP (WMC)
> production_status : Operational products
> long_name : U-component of wind
> units : m s-1
> _FillValue : 1e+20
> coordinates : gridlat_0 gridlon_0
> grid_type : Lambert Conformal can be secant or tangent, conical
> or bipolar
> parameter_discipline_and_category : Meteorological products,
> Momentum
> parameter_template_discipline_category_number : ( 0, 0, 2,
> 2 )
> level_type : Specified height level above ground (m)
> forecast_time : 0
> forecast_time_units : hours
> initial_time : 11/23/2021 (06:00)
>
> grib2 (WRF ARW)
> NCL 6.6.2
> UBUNTU 20.04
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <
> https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20211125/bc5ab177/attachment-0001.html
> >
>
> ------------------------------
>
> Message: 2
> Date: Thu, 25 Nov 2021 10:21:45 -0700
> From: Dennis Shea <shea at ucar.edu>
> To: "R. Lucas" <rousslucas at gmail.com>
> Cc: ncl-talk <ncl-talk at mailman.ucar.edu>
> Subject: Re: [ncl-talk] Variable size
> Message-ID:
> <CAOF1d_49Fz0Qb051ku3aOKBC=jxHksHpnnu==MywF=
> X9bUBCjg at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> UGRD_P0_L103_GLC0 (* lv_HTGL3*, ygrid_0, xgrid_0 )
> level_type : Specified height level above ground (m)
> --------------------------------------------------------------------------
>
> Likely, * lv_HTGL3, *is a coordinate variable.
>
> print( UGRD_P0_L103_GLC0*&lv_HTGL3 *)
>
> will tell you the levels (m)
>
> So ...
>
> u2d = UGRD_P0_L103_GLC0 (0,:,: )
>
> will yield a 2-dimensional variable corresponding to the 0-th dimension.
>
> Let's say that * lv_HTGL3 *contained (/0, 2, 10/)
>
> You could also use NCL's coordinate subscripting via the *{...}* syntax
>
> u2d_0 = UGRD_P0_L103_GLC0 (*{0}*,:,: )
> u2d_2 = UGRD_P0_L103_GLC0 (*{2}*,:,: )
> u2d_10 = UGRD_P0_L103_GLC0 (*1**0}*,:,: )
>
>
> On Thu, Nov 25, 2021 at 2:26 AM R. Lucas via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
> > Dear NCL users
> > i have 2 problems :
> >
> > 1) Is this varibable level looks like 10m above ground ?
> > 2) This variable is 3d. How can I remove first dimension, somthing like
> > (0,:,:)
> > because I have a fatal error :
> > fatal:["Execute.c":7925]:Number of subscripts on right-hand-side do not
> > match
> > number of dimensions of variable: (3), Subscripts used: (1)
> >
> > Thanks in advance.
> > Rouss.
> >
> >
> > float UGRD_P0_L103_GLC0 ( lv_HTGL3, ygrid_0, xgrid_0 )
> > center : US National Weather Service - NCEP (WMC)
> > production_status : Operational products
> > long_name : U-component of wind
> > units : m s-1
> > _FillValue : 1e+20
> > coordinates : gridlat_0 gridlon_0
> > grid_type : Lambert Conformal can be secant or tangent,
> > conical or bipolar
> > parameter_discipline_and_category : Meteorological products,
> > Momentum
> > parameter_template_discipline_category_number : ( 0, 0,
> 2,
> > 2 )
> > level_type : Specified height level above ground (m)
> > forecast_time : 0
> > forecast_time_units : hours
> > initial_time : 11/23/2021 (06:00)
> >
> > grib2 (WRF ARW)
> > NCL 6.6.2
> > UBUNTU 20.04
> > _______________________________________________
> > ncl-talk mailing list
> > ncl-talk at mailman.ucar.edu
> > List instructions, subscriber options, unsubscribe:
> > https://mailman.ucar.edu/mailman/listinfo/ncl-talk
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <
> https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20211125/5a551a41/attachment-0001.html
> >
>
> ------------------------------
>
> Message: 3
> Date: Thu, 25 Nov 2021 12:27:26 -0500
> From: Israt Jahan <israt.jahan at uconn.edu>
> To: ncl-talk at mailman.ucar.edu
> Subject: Re: [ncl-talk] NCL function wrf_user_getvar
> Message-ID:
> <
> CANJTbXaN4fN31Zp-DRLM_Od-ZRyUvkDJUOV868mYRxrnBDDhsQ at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Hello,
> It would be much appreciated if someone could shed some light on this. Does
> wrf_user_getvar(NCDF_File,"z",0) measure z from ground level or sea level?
> I need to find temperature at 1000 m above ground level. I am doing the
> following:
> 1. ht=wrf_user_getvar(NCDF_File,"z",0) ;
> where NCDF_File is a list of WRF netCDF files
> 2. tk=wrf_user_getvar(NCDF_File,"tk",0) ;get temperature in K
> 2. temp_1km=wrf_user_interp_level(tk,ht,1000,0) ;
> get temperature at 1 km above ground level
> Is my approach correct?
> Best,
> Israt
>
>
> On Tue, Nov 16, 2021 at 4:26 PM Israt Jahan <israt.jahan at uconn.edu> wrote:
>
> > Dear NCL users,
> > I have a question about the NCL function "wrf_user_getvar". I see this
> > function can extract data from WRF files and calculate the requested
> > diagnostic. One of the available diagnostics is "z/height" which is
> labeled
> > as the "full model height in m" in this site wrf_user_getvar (ucar.edu)
> > <
> https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml
> >.
> > From WRF output files, I need to find temperature at 1000 m above ground
> > level. I am doing the following:
> > 1. ht=wrf_user_getvar(NCDF_File,"z",0) where NCDF_File is a list of WRF
> > netCDF files
> > 2. tk=wrf_user_getvar(NCDF_File,"tk",0) to get temperature in K
> > 2. temp_1km=wrf_user_interp_level(tk,ht,1000,0) to get temperature at 1
> km
> >
> > My question is: is the diagnostic "z/height" that we get from
> > "wrf_user_getvar" measured from ground level or sea level?
> > If it is measured from sea level, what should I do to get z from ground
> > level?
> >
> > Thank you so much.
> > Best,
> > Israt
> >
> >
> > --
> > Regards,
> > Israt Jahan
> > Graduate Assistant
> > Department of Civil and Environmental Engineering
> > School of Engineering | University of Connecticut
> > Email: israt.jahan at uconn.edu
> >
> --
> Regards,
> Israt Jahan
> Graduate Assistant
> Department of Civil and Environmental Engineering
> School of Engineering | University of Connecticut
> Email: israt.jahan at uconn.edu
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <
> https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20211125/a7f2ce89/attachment-0001.html
> >
>
> ------------------------------
>
> Message: 4
> Date: Thu, 25 Nov 2021 11:01:09 -0700
> From: Dennis Shea <shea at ucar.edu>
> To: Israt Jahan <israt.jahan at uconn.edu>
> Cc: ncl-talk <ncl-talk at mailman.ucar.edu>
> Subject: Re: [ncl-talk] NCL function wrf_user_getvar
> Message-ID:
> <CAOF1d_5xiRj_ck=V1mZney7Rkk2fiXu8Quzupkr5-=mYziJ=
> LA at mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> The **wrf_user* *functions and documentation were created by WRF
> developers not NCL developers.
> NCL provides an interface to to provided derivation code [usually fortran].
> Hence, we can not support the underlying codes. Please contact WRF.
> ========
> That said:
>
> NCDF = addfile("wrfout_d02_2017-10-21_00_00_00.nc","r")
> ht = wrf_user_getvar(NCDF,"z",0)
> printVarSummary(ht) ; ?? [bottom_top | ..] x
> [south_north | ..] x [west_east | ..]?? height
> tk = wrf_user_getvar(NCDF,"tk",0)
> printVarSummary(tk) ; ; ?? [bottom_top | ..] x
> [south_north | ..] x [west_east | ..]?? temperature
>
> J = 45 ; arbitrary latitude index
> I = 50 ; longitude
>
> print( ht(:,J,I)+" "+tk(:,J,I)) ; print all vertical values at a
> specified location
>
> ===
> *wrf_user_interp_level*
> <
> https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_interp_level.shtml
> >
> was introduced in the NCL 6.6 release. I do not have NCL 6.6 on my MAC
>
> temp_1km=*wrf_user_interp_level*(tk,ht,1000,0)
> printVarSummary( temp_1km )
>
> However, you could print the derived 1000m temperature at the J,K
>
> print( temp_1km(J,K) )
>
> This should indicate to you if the derived value is correct.
>
> ===
>
>
>
>
>
> On Thu, Nov 25, 2021 at 10:27 AM Israt Jahan via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
> >
> > Hello,
> > It would be much appreciated if someone could shed some light on this.
> > Does wrf_user_getvar(NCDF_File,"z",0) measure z from ground level or sea
> > level?
> > I need to find temperature at 1000 m above ground level. I am doing the
> > following:
> > 1. ht=wrf_user_getvar(NCDF_File,"z",0) ;
> > where NCDF_File is a list of WRF netCDF files
> > 2. tk=wrf_user_getvar(NCDF_File,"tk",0) ;get temperature in K
> > 2. temp_1km=wrf_user_interp_level(tk,ht,1000,0) ;
> > get temperature at 1 km above ground level
> > Is my approach correct?
> > Best,
> > Israt
> >
> >
> > On Tue, Nov 16, 2021 at 4:26 PM Israt Jahan <israt.jahan at uconn.edu>
> wrote:
> >
> >> Dear NCL users,
> >> I have a question about the NCL function "wrf_user_getvar". I see this
> >> function can extract data from WRF files and calculate the requested
> >> diagnostic. One of the available diagnostics is "z/height" which is
> labeled
> >> as the "full model height in m" in this site wrf_user_getvar (ucar.edu)
> >> <
> https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_getvar.shtml
> >.
> >> From WRF output files, I need to find temperature at 1000 m above ground
> >> level. I am doing the following:
> >> 1. ht=wrf_user_getvar(NCDF_File,"z",0) where NCDF_File is a list of WRF
> >> netCDF files
> >> 2. tk=wrf_user_getvar(NCDF_File,"tk",0) to get temperature in K
> >> 2. temp_1km=wrf_user_interp_level(tk,ht,1000,0) to get temperature at 1
> >> km
> >>
> >> My question is: is the diagnostic "z/height" that we get from
> >> "wrf_user_getvar" measured from ground level or sea level?
> >> If it is measured from sea level, what should I do to get z from ground
> >> level?
> >>
> >> Thank you so much.
> >> Best,
> >> Israt
> >>
> >>
> >> --
> >> Regards,
> >> Israt Jahan
> >> Graduate Assistant
> >> Department of Civil and Environmental Engineering
> >> School of Engineering | University of Connecticut
> >> Email: israt.jahan at uconn.edu
> >>
> > --
> > Regards,
> > Israt Jahan
> > Graduate Assistant
> > Department of Civil and Environmental Engineering
> > School of Engineering | University of Connecticut
> > Email: israt.jahan at uconn.edu
> > _______________________________________________
> > ncl-talk mailing list
> > ncl-talk at mailman.ucar.edu
> > List instructions, subscriber options, unsubscribe:
> > https://mailman.ucar.edu/mailman/listinfo/ncl-talk
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <
> https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20211125/5ee16f93/attachment-0001.html
> >
>
> ------------------------------
>
> Subject: Digest Footer
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
> ------------------------------
>
> End of ncl-talk Digest, Vol 216, Issue 16
> *****************************************
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20211127/f491eaae/attachment.html>
More information about the ncl-talk
mailing list