Re: hydro problem

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Mon Feb 01 2010 - 09:55:33 MST

> I would calculate geopotential height from my model data by hydro NCL
> function.
> I test it, but the resulting field has a pattern equal to surface
> geopotential, also if the geopotential height numerical values seems to
> be corrected.
>
> Could you help me?

I tried the script/data you sent and it works fine.

Note: it looks like you added p(time, lev, lat, lon) to your netCDF.
This is not necessary. They are all repeats of constant pressure levels.
You can create this dynamically:

pi =(/ 1000,950, 925, 900, 850., 800, 700 ,600 ,500 ,400 ,300 ,250 , \
           200 ,150 ,100 ,70 ,50/)

If T(time,lev,lat,lon)

P = conform(T,pi,1) ; P(ntim,klev,nlat,mlon)
http://www.ncl.ucar.edu/Document/Functions/Built-in/conform.shtml
http://www.ncl.ucar.edu/Document/Functions/Built-in/conform_dims.shtml
=================================
Variable: geopoth
Type: double
Total Size: 2506752 bytes
             313344 values
Number of Dimensions: 4
Dimensions and sizes: [time | 1] x [lat | 96] x [lon | 192] x [lev | 17]
Coordinates:
             time: [ 0.. 0]
             lat: [88.57216851400727..-88.57216851400727]
             lon: [ 0..358.125]
             lev: [100000..5000]
Number Of Attributes: 2
   long_name : Geopotential Height
   units : m
(0) P Tv Z
(0) =================
(0) 1000.0 276.5 15.2
(1) 950.0 272.4 427.2
(2) 925.0 270.5 639.1
(3) 900.0 268.5 855.2
(4) 850.0 265.4 1301.9
(5) 800.0 264.7 1772.3
(6) 700.0 261.5 2800.6
(7) 600.0 253.9 3963.6
(8) 500.0 243.5 5291.1
(9) 400.0 231.1 6841.6
(10) 300.0 218.3 8735.0
(11) 250.0 215.0 9891.3
(12) 200.0 217.8 11304.6
(13) 150.0 219.5 13145.7
(14) 100.0 218.4 15744.2
(15) 70.0 218.6 18024.8
(16) 50.0 220.7 20187.1

Good luck
D

>
> Thank you in advance!
> Cheers
> Andrea
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

-- 
======================================================
Dennis J. Shea                  tel: 303-497-1361    |
P.O. Box 3000                   fax: 303-497-1333    |
Climate Analysis Section                             |
Climate & Global Dynamics Div.                       |
National Center for Atmospheric Research             |
Boulder, CO  80307                                   |
USA                        email: shea 'at' ucar.edu |
======================================================


_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk

Received on Mon Feb 1 09:55:39 2010

This archive was generated by hypermail 2.1.8 : Sun Feb 07 2010 - 15:37:25 MST