[ncl-talk] Interpolate WRF eta levels to pressure levels
Anahita Amiri Farahani
aamir003 at ucr.edu
Wed Jan 8 09:05:58 MST 2020
So you mean this:
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
ana=addfile("wrfout_d01_2015-08-13","r")
;ana=addfile("wrfout_d01_2015-08-13_00:00:00","r")
lat1=ana->XLAT(0,:,0)
lon1=ana->XLONG(0,0,:)
P=ana->P
pm10=ana->PM10
pp = wrf_user_getvar(ana, "pressure",-1)
delete(ana)
printVarSummary(pm10)
printMinMax(pm10,0)
print("====")
printVarSummary(P)
printMinMax(P,0)
print("====")
po=(/1000,925,850,700,600,500,400,300,250,200,150,100/);*100.0
po!0 = "po"
po at long_name = "Pressure Level"
po at units = "hPa"
To = int2p_n_Wrap (pp,pm10,po,1,1)
printVarSummary(To)
printMinMax(To,0)
print("====")
and I have this:
Variable: To
Type: float
Total Size: 4678560 bytes
1169640 values
Number of Dimensions: 4
Dimensions and sizes: [Time | 9] x [po | 12] x [south_north | 95] x
[west_east | 114]
Coordinates:
po: [1000..100]
Number Of Attributes: 7
_FillValue : 9.96921e+36
FieldType : 104
MemoryOrder : XYZ
description : pm10 dry mass
units : ug m^-3
stagger :
coordinates : XLONG XLAT XTIME
ncl 32> printMinMax(To,0)
(0) pm10 dry mass (ug m^-3) : min=0 max=32.0972
ncl 33> print("====")
(0) ====
The max size is much smaller than original pm10.
Thanks,
Ana
On Wed, Jan 8, 2020 at 10:59 AM Dennis Shea <shea at ucar.edu> wrote:
> Yes.
>
> Sent from my iPhone
>
> On Jan 8, 2020, at 8:32 AM, Anahita Amiri Farahani <aamir003 at ucr.edu>
> wrote:
>
>
> Hi Dennis,
>
> I tried this, but it gives me missing values:
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
>
> ana=addfile("wrfout_d01_2015-08-13","r")
>
> lat1=ana->XLAT(0,:,0)
> lon1=ana->XLONG(0,0,:)
> P=ana->P
> pm10=ana->PM10
> delete(ana)
>
> printVarSummary(pm10)
> printMinMax(pm10,0)
> print("====")
> printVarSummary(P)
> printMinMax(P,0)
> print("====")
>
>
> po=(/1000,925,850,700,600,500,400,300,250,200,150,100/)*100.0
> po!0 = "po"
> po at long_name = "Pressure Level"
> po at units = "Pa"
>
> To = int2p_n_Wrap (P,pm10,po,1,1)
> printVarSummary(To)
> printMinMax(To,0)
> print("====")
>
>
>
>
>
> Variable: pm10
>
> Type: float
>
> Total Size: 13255920 bytes
>
> 3313980 values
>
> Number of Dimensions: 4
>
> Dimensions and sizes: [Time | 9] x [bottom_top | 34] x [south_north | 95]
> x [west_east | 114]
>
> Coordinates:
>
> Number Of Attributes: 6
>
> FieldType : 104
>
> MemoryOrder : XYZ
>
> description : pm10 dry mass
>
> units : ug m^-3
>
> stagger :
>
> coordinates : XLONG XLAT XTIME
>
> (0) pm10 dry mass (ug m^-3) : min=0 max=183.558
>
> (0) ====
>
>
> Variable: P
>
> Type: float
>
> Total Size: 13255920 bytes
>
> 3313980 values
>
> Number of Dimensions: 4
>
> Dimensions and sizes: [Time | 9] x [bottom_top | 34] x [south_north | 95]
> x [west_east | 114]
>
> Coordinates:
>
> Number Of Attributes: 6
>
> FieldType : 104
>
> MemoryOrder : XYZ
>
> description : perturbation pressure
>
> units : Pa
>
> stagger :
>
> coordinates : XLONG XLAT XTIME
>
> (0) perturbation pressure (Pa) : min=0.702637 max=2079.95
>
> (0) ====
>
>
> Variable: To
>
> Type: float
>
> Total Size: 4678560 bytes
>
> 1169640 values
>
> Number of Dimensions: 4
>
> Dimensions and sizes: [Time | 9] x [po | 12] x [south_north | 95] x
> [west_east | 114]
>
> Coordinates:
>
> po: [100000..10000]
>
> Number Of Attributes: 7
>
> _FillValue : 9.96921e+36
>
> FieldType : 104
>
> MemoryOrder : XYZ
>
> description : pm10 dry mass
>
> units : ug m^-3
>
> stagger :
>
> coordinates : XLONG XLAT XTIME
>
> (0) pm10 dry mass (ug m^-3) : min=9.96921e+36 max=9.96921e+36
>
> (0) ====
>
>
> P is the perturbation pressure, should I use pressure instead?
>
> Thanks,
> Ana
>
>
>
> On Tue, Jan 7, 2020 at 11:14 PM Dennis Shea <shea at ucar.edu> wrote:
>
>> ***If*** the WRF file has a variable
>> P(Time,bottom_top,south_north,west_east)
>> then those are the pressure levels at each 'eta' coordinate at each grid
>> point and time.
>>
>> If say you want to interpolate the variable
>> T(Time,bottom_top,south_north,west_east)
>> then something like:
>> T = f->T
>> P = f->P
>> printVarSummary(T)
>> printMinMax(T,0)
>> print("====")
>> printVarSummary(P)
>> printMinMax(P,0)
>> print("====")
>>
>> po = (/1000,900,850,700,600,500,400,300,200,100,50,10,5/)*100.0
>> po!0 = "po"
>> po at long_name = "Pressure Level"
>> po at units = "Pa"
>>
>> To =* int2p_n_Wrap*
>> <http://www.ncl.ucar.edu/Document/Functions/Contributed/int2p_n_Wrap.shtml>
>> (P,T,po,1,1)
>> printVarSummary(To)
>> printMinMax(To,0)
>> print("====")
>> ====
>> If P is not on the WRF file, then you will have to get help from WRF
>> support.
>>
>>
>>
>> On Tue, Jan 7, 2020 at 9:12 AM Anahita Amiri Farahani via ncl-talk <
>> ncl-talk at ucar.edu> wrote:
>>
>>> Hi all,
>>>
>>> Is there any function/command in NCL to convert WRF eta levels to
>>> pressure levels?
>>>
>>> Thanks,
>>> Ana
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200108/bf785e78/attachment.html>
More information about the ncl-talk
mailing list