[ncl-talk] Derivatives in NCL
Gus Correa
gus at ldeo.columbia.edu
Mon Feb 6 12:17:47 MST 2017
Hi Forrest
Using finite differences, perhaps?
http://www.ncl.ucar.edu/Document/Functions/Built-in/center_finite_diff.shtml
http://www.ncl.ucar.edu/Document/Functions/Built-in/center_finite_diff_n.shtml
That should be OK for the vertical direction
(your lv_ISBL0,I presume).
I hope this helps,
Gus Correa
On 02/06/2017 01:40 PM, Black, Forrest (LARC-D318)[UNIVERSITIES SPACE
RESEARCH ASSOCIATION] wrote:
> It appears my question was unclear.
>
>
>
> I would like to know if there are any suggestions for how to find the
> derivative of (Re-h) as mentioned below.
>
>
>
> Thank you,
>
>
>
> - Forrest
>
>
>
> *From:* ncl-talk-bounces at ucar.edu [mailto:ncl-talk-bounces at ucar.edu] *On
> Behalf Of *Black, Forrest (LARC-D318)[UNIVERSITIES SPACE RESEARCH
> ASSOCIATION]
> *Sent:* Monday, February 06, 2017 11:04 AM
> *To:* ncl-talk at ucar.edu
> *Subject:* [ncl-talk] Derivatives in NCL
>
>
>
> Hello all,
>
>
>
> I am trying to solve the following equation in NCL:
>
>
>
> Z = (h*Re)/(Re-h)’
>
>
>
> Z: Geometric Height
>
> h: Geopotential height (3D Array)
>
> Re: Radius of Earth, assumed constant (6356766m)
>
>
>
> I have looked at some old threads, but am still unsure of what functions
> I would use to calculate (Re-h)’.
>
>
>
> Can anyone give me some suggestions on how to tackle this?
>
>
>
> This is the Attribute data for Geopotential Height:
>
>
>
> Variable: Geopotential Height
>
> Type: float
>
> Total Size: 2525324 bytes
>
> 631331 values
>
> Number of Dimensions: 3
>
> Dimensions and sizes: [lv_ISBL0 | 37] x [ygrid_0 | 113] x [xgrid_0 | 151]
>
> Coordinates:
>
> lv_ISBL0: [10000..100000]
>
> Number Of Attributes: 13
>
> center : US National Weather Service - NCEP (WMC)
>
> production_status : Operational products
>
> long_name : Geopotential height
>
> units : gpm
>
> _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, Mass
>
> parameter_template_discipline_category_number : ( 0, 0, 3, 5 )
>
> level_type : Isobaric surface (Pa)
>
> forecast_time : 0
>
> forecast_time_units : hours
>
> initial_time : 01/31/2017 (00:00)
>
>
>
>
>
> Thanks for any help.
>
>
>
> - Forrest
>
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
More information about the ncl-talk
mailing list