forked from 3rdparty/wrf-python
You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
87 lines
2.8 KiB
87 lines
2.8 KiB
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/wrf/WRFUserARW.ncl" |
|
|
|
in_files = (/"/Users/ladwig/Documents/wrf_files/wrf_vortex_multi/wrfout_d02_2005-08-28_00:00:00", \ |
|
"/Users/ladwig/Documents/wrf_files/wrf_vortex_multi/wrfout_d02_2005-08-28_12:00:00", \ |
|
"/Users/ladwig/Documents/wrf_files/wrf_vortex_multi/wrfout_d02_2005-08-29_00:00:00"/) |
|
|
|
input_file = addfiles(in_files, "r") |
|
|
|
timeidx = 6 |
|
|
|
; 3D interpolation to new vertical coordinates |
|
; interp t to theta |
|
fld1 = wrf_user_getvar(input_file, "tk", timeidx) |
|
vert_coord = "theta" |
|
;interp_levels = (/200,300,500,1000/) |
|
interp_levels = ispan(200,1000,50) |
|
opts = True |
|
opts@extrapolate = True |
|
opts@field_type = "T" |
|
opts@logP = True |
|
opts@time = timeidx |
|
|
|
fld1_intrp = wrf_user_vert_interp(input_file,fld1,vert_coord,interp_levels,opts) |
|
|
|
printVarSummary(fld1_intrp) |
|
|
|
; interp t to theta-e |
|
fld2 = fld1 |
|
vert_coord := "theta-e" |
|
fld2_intrp = wrf_user_vert_interp(input_file,fld2,vert_coord,interp_levels,opts) |
|
|
|
printVarSummary(fld2_intrp) |
|
|
|
; interp t to pressure |
|
fld3 = fld1 |
|
vert_coord := "pressure" |
|
interp_levels := (/850,500/) |
|
fld3_intrp = wrf_user_vert_interp(input_file,fld3,vert_coord,interp_levels,opts) |
|
|
|
printVarSummary(fld3_intrp) |
|
|
|
; interp t to ght_msl |
|
fld4 = fld1 |
|
vert_coord := "ght_msl" |
|
interp_levels := (/1,2/) |
|
fld4_intrp = wrf_user_vert_interp(input_file,fld4,vert_coord,interp_levels,opts) |
|
|
|
printVarSummary(fld4_intrp) |
|
|
|
; interp t to ght_agl |
|
fld5 = fld1 |
|
vert_coord := "ght_agl" |
|
interp_levels := (/1,2/) |
|
fld5_intrp = wrf_user_vert_interp(input_file,fld1,vert_coord,interp_levels,opts) |
|
|
|
printVarSummary(fld5_intrp) |
|
|
|
; interp ht to pres |
|
fld6 = wrf_user_getvar(input_file, "height", timeidx) |
|
vert_coord := "pressure" |
|
opts@field_type = "ght" |
|
interp_levels := (/500,50/) |
|
fld6_intrp = wrf_user_vert_interp(input_file,fld6,vert_coord,interp_levels,opts) |
|
|
|
printVarSummary(fld6_intrp) |
|
|
|
; interp pres to theta |
|
fld7 = wrf_user_getvar(input_file, "pressure", timeidx) |
|
vert_coord := "theta" |
|
opts@field_type = "pressure" |
|
interp_levels := (/200,300,500,1000/) |
|
fld7_intrp = wrf_user_vert_interp(input_file,fld7,vert_coord,interp_levels,opts) |
|
|
|
printVarSummary(fld7_intrp) |
|
|
|
; interp theta-e to pressure |
|
fld8 = wrf_user_getvar(input_file, "eth", timeidx) |
|
vert_coord := "pressure" |
|
opts@field_type = "T" |
|
interp_levels := (/850,500,5/) |
|
fld8_intrp = wrf_user_vert_interp(input_file,fld8,vert_coord,interp_levels,opts) |
|
|
|
printVarSummary(fld8_intrp) |
|
|
|
|