[ncl-talk] plotting problem with JRA55
Dennis Shea
shea at ucar.edu
Tue Mar 7 14:17:44 MST 2017
In my opinion, the source nc files are not well constructed.
---
NCL variables are data-structures. They contain the values and meta data:
(a) long_name, (b) units; and (c) coordinate variables such as time,
level/depth, lat/lon (rectilinear grids).
The 'coordinate variables' associated with each variable are assumed to
follow the NUG netCDF rules:
(1) variable name and dimension name are the same;
(2) numeric;
(3) one-dimensional;
(4) monotonically increasing or decreasing
An issue:
Variable: u_wv
Dimensions and sizes: [initial_time0_hours | 648] x [g0_lat_1 | 53] x
[g0_lon_2 | 113]
Coordinates:
initial_time0_hours: [1385736..1858344]
g0_lat_1: [90..25] <==== GOOD monotonically
decreasing
g0_lon_2: [250..30] <==== BAD for lon 'coordinate'
variable
This is going: 250, ...,358.75?, 0, 1.25 ,....) *not*
monotonically increasing
=====
One fix:
g0_lon_2 = a->g0_lon_2
g0_lon_2 = where(g0_lon_2.gt.180, g0_lon_2-360, g0_lon_2)
g0_lon_2&g0_lon_2 = g0_lon_2
printVarSummary(g0_lon_2)
u_wv = a->UWV_GDS0_EATM_S123(1:648,:,:)
Variable: u_wv
Type: float
Total Size: 15523488 bytes
3880872 values
Number of Dimensions: 3
Dimensions and sizes: [initial_time0_hours | 648] x [g0_lat_1 | 53] x
[g0_lon_2 | 113]
Coordinates:
initial_time0_hours: [1385736..1858344]
g0_lat_1: [90..25]
g0_lon_2: [250..30]
===
u_wv&g0_lon_2 = g0_lon_2 ; replace with 'correct' ordering
Variable: u_wv
Coordinates:
initial_time0_hours: [1385736..1858344]
g0_lat_1: [90..25]
g0_lon_2: [-110..30] ,=== monotonically increasing
Same for v_wv
v_wv&g0_lon_2 = g0_lon_2
===
Any other variables also.
Good luck
On Tue, Mar 7, 2017 at 1:33 PM, MIRA BERDAHL <mlosic at scarletmail.rutgers.edu
> wrote:
> Hi Mary,
> Yes that is right, I have changed the script to use gsn_csm_vector_map and
> there is a map displayed now. Thank you! Trouble is, the longitude is still
> wrong, so it is plotting my subregion in the wrong place. The latitudes
> look good, but the range of longitudes should be -110/30 (or 250/30).
> I've attached my new plots so you can see where it is plotting them.
> Any thoughts on why this would be shifted like this?
> Thanks for any ideas,
> Mira
>
> I've copied my script and variable information below again:
>
>
>
>
> 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"
> ;************************************************
> begin
> ;************************************************
> ; read in netCDF files
> ;************************************************
> a = addfile("Uwv.mon.mean.nc","r") ; u moisture transport
> b = addfile("Vwv.mon.mean.nc","r") ; v moisture transport
>
> ;************************************************
> ; read in zonal [u] and meridional [v] moisture transport (water vapor)
> ;************************************************
>
> u_wv = a->UWV_GDS0_EATM_S123(1:648,:,:)
> v_wv = b->VWV_GDS0_EATM_S123(1:648,:,:) ;
>
> lat = a->g0_lat_1
> lon = a->g0_lon_2
>
> printVarSummary(lat)
> printVarSummary(lon)
>
> printVarSummary(u_wv)
> printVarSummary(v_wv)
>
>
> ; Calculate the seasonal averages.
> uwvDJF = month_to_season(u_wv, "DJF")
> vwvDJF = month_to_season(v_wv, "DJF")
>
>
> printVarSummary(uwvDJF)
> printVarSummary(vwvDJF)
>
>
> ; from the matlab script i wrote: findExtremeYrs, i pulled out the extreme
> years (> or < 1std) that i want to average and plot here.
>
>
> ;ans = 7 (longitude of IL)
> ;We find the years where the Icelandic low is in an “extreme east” (high
> Long) position to be: 1966, 1968, 1969, 1974, 1983, 1984, 1994, 1995, 1999
> and 2005.
> ;The years with “extreme west” (low long) position of the Icelandic Low
> are: 1963, 1964, 1980, 1985, 1987, 1991, 1992, 1996, 2003 and 2006.
>
>
> ; this data starts at 1958 (this is index 0), so 1963=5, 1973=15 etc.
>
>
> uDJF_lon_hi = uwvDJF((/8,10,11,16,25,26,36,37,41,47/),:,:)
> uDJF_lon_lo = uwvDJF((/5,6,22,27,29,33,34,38,45,48/),:,:)
>
> vDJF_lon_hi = vwvDJF((/8,10,11,16,25,26,36,37,41,47/),:,:)
> vDJF_lon_lo = vwvDJF((/5,6,22,27,29,33,34,38,45,48/),:,:)
>
> uAvgTime_hi = dim_avg_n_Wrap(uDJF_lon_hi,0)
> uAvgTime_lo = dim_avg_n_Wrap(uDJF_lon_lo,0)
>
> printVarSummary(uAvgTime_hi)
> printVarSummary(uAvgTime_lo)
>
> vAvgTime_hi = dim_avg_n_Wrap(vDJF_lon_hi,0)
> vAvgTime_lo = dim_avg_n_Wrap(vDJF_lon_lo,0)
>
> printVarSummary(vAvgTime_hi)
> printVarSummary(vAvgTime_lo)
>
>
> ; sneaky way to copy metadata over first.
> diff_u = uAvgTime_hi;
> diff_v = vAvgTime_hi;
>
> diff_u = uAvgTime_lo - uAvgTime_hi
> diff_v = vAvgTime_lo - vAvgTime_hi
>
> printVarSummary(diff_u)
> printVarSummary(diff_v)
>
> ;************************************************
> ; create plot
> ;************************************************
> wks = gsn_open_wks("ps","Panel_lon_WV") ; open a ps file
> gsn_define_colormap(wks,"temp1")
>
> plot = new(3,graphic) ; create a plot array
>
> ;---- set common resources for all plots
> res = True
> res at gsnDraw = False ; dont draw
> res at gsnFrame = False ; dont advance frame
> res at cnInfoLabelOn = False ; trn off cn info label
> res at gsnAddCyclic = False ; has to do with wrapping the lonitude at 0/360
> res at cnFillPalette = "matlab_jet"
>
>
> ;***********************************************
> ; ----wind vector plot
> ;***********************************************
> vcres = res
> vcres at vcRefAnnoOrthogonalPosF = -1.0 ; move ref vector up
> vcres at vcRefMagnitudeF = 10.0 ; define vector ref mag
> vcres at vcRefLengthF = 0.045 ; define length of vec ref
> vcres at vcGlyphStyle = "CurlyVector" ; turn on curly vectors
> vcres at vcMinDistanceF = 0.017
> vcres at mpFillOn = False ; turn off gray fill
> vcres at mpOutlineBoundarySets = "National" ; turn on country boundaries
> ;vcres at mpFillBoundarySets = AllBoundaries
> vcres at mpGeophysicalLineColor = "Navy" ; color of cont. outlines
> vcres at mpGeophysicalLineThicknessF = 1.5 ; thickness of outlines
> vcres at gsnRightString = "" ; turn off thue units string
>
> plot(0) = gsn_csm_vector_map(wks,uAvgTime_hi,vAvgTime_hi,vcres)
> plot(1) = gsn_csm_vector_map(wks,uAvgTime_lo,vAvgTime_lo,vcres)
> plot(2) = gsn_csm_vector_map(wks, diff_u, diff_v,vcres)
>
> ;************************************************
> ; create panel
> ;************************************************
> resP = False ; modify the panel plot
> resP at txString = "lon - water vapor transport"
> gsn_panel(wks,(/plot/),(/3,1/),resP) ; now draw as one plot;
>
>
> end
>
>
>
> ************************************************************
> *******************
> Variable: lat
> Type: float
> Total Size: 212 bytes
> 53 values
> Number of Dimensions: 1
> Dimensions and sizes: [g0_lat_1 | 53]
> Coordinates:
> g0_lat_1: [90..25]
> Number Of Attributes: 9
> La1 : 90
> Lo1 : 0
> La2 : -90
> Lo2 : -1.25
> Di : 1.25
> Dj : 1.25
> units : degrees_north
> GridType : Cylindrical Equidistant Projection Grid
> long_name : latitude
>
> Variable: lon
> Type: float
> Total Size: 452 bytes
> 113 values
> Number of Dimensions: 1
> Dimensions and sizes: [g0_lon_2 | 113]
> Coordinates:
> g0_lon_2: [250..30]
> Number Of Attributes: 9
> La1 : 90
> Lo1 : 0
> La2 : -90
> Lo2 : -1.25
> Di : 1.25
> Dj : 1.25
> units : degrees_east
> GridType : Cylindrical Equidistant Projection Grid
> long_name : longitude
>
> Variable: u_wv
> Type: float
> Total Size: 15523488 bytes
> 3880872 values
> Number of Dimensions: 3
> Dimensions and sizes: [initial_time0_hours | 648] x [g0_lat_1 | 53] x
> [g0_lon_2 | 113]
> Coordinates:
> initial_time0_hours: [1385736..1858344]
> g0_lat_1: [90..25]
> g0_lon_2: [250..30]
> Number Of Attributes: 15
> N : <ARRAY of 12 elements>
> statistical_process_duration : instantaneous (beginning at
> reference time at intervals of 6 hours)
> statistical_process_descriptor : average of N uninitialized analyses
> forecast_time_units : hours
> forecast_time : 0
> level : 0
> parameter_number : 157
> parameter_table_version : 200
> gds_grid_type : 0
> level_indicator : 200
> _FillValue : 1e+20
> units : kg/m/s
> long_name : Zonal water vapour flux
> center : Japanese Meteorological Agency - Tokyo (RSMC)
> sub_center : 241
>
> Variable: v_wv
> Type: float
> Total Size: 15523488 bytes
> 3880872 values
> Number of Dimensions: 3
> Dimensions and sizes: [initial_time0_hours | 648] x [g0_lat_1 | 53] x
> [g0_lon_2 | 113]
> Coordinates:
> initial_time0_hours: [1385736..1858344]
> g0_lat_1: [90..25]
> g0_lon_2: [250..30]
> Number Of Attributes: 15
> N : <ARRAY of 12 elements>
> statistical_process_duration : instantaneous (beginning at
> reference time at intervals of 6 hours)
> statistical_process_descriptor : average of N uninitialized analyses
> forecast_time_units : hours
> forecast_time : 0
> level : 0
> parameter_number : 152
> parameter_table_version : 200
> gds_grid_type : 0
> level_indicator : 200
> _FillValue : 1e+20
> units : Kg/m/s
> long_name : Meridional water vapour flux
> center : Japanese Meteorological Agency - Tokyo (RSMC)
> sub_center : 241
>
> Variable: uwvDJF
> Type: float
> Total Size: 1293624 bytes
> 323406 values
> Number of Dimensions: 3
> Dimensions and sizes: [initial_time0_hours | 54] x [g0_lat_1 | 53] x
> [g0_lon_2 | 113]
> Coordinates:
> initial_time0_hours: [1385736..1850328]
> g0_lat_1: [90..25]
> g0_lon_2: [250..30]
> Number Of Attributes: 16
> N : <ARRAY of 12 elements>
> statistical_process_duration : instantaneous (beginning at
> reference time at intervals of 6 hours)
> statistical_process_descriptor : average of N uninitialized analyses
> forecast_time_units : hours
> forecast_time : 0
> level : 0
> parameter_number : 157
> parameter_table_version : 200
> gds_grid_type : 0
> level_indicator : 200
> _FillValue : 1e+20
> units : kg/m/s
> long_name : DJF: Zonal water vapour flux
> center : Japanese Meteorological Agency - Tokyo (RSMC)
> sub_center : 241
> NMO : 0
>
> Variable: vwvDJF
> Type: float
> Total Size: 1293624 bytes
> 323406 values
> Number of Dimensions: 3
> Dimensions and sizes: [initial_time0_hours | 54] x [g0_lat_1 | 53] x
> [g0_lon_2 | 113]
> Coordinates:
> initial_time0_hours: [1385736..1850328]
> g0_lat_1: [90..25]
> g0_lon_2: [250..30]
> Number Of Attributes: 16
> N : <ARRAY of 12 elements>
> statistical_process_duration : instantaneous (beginning at
> reference time at intervals of 6 hours)
> statistical_process_descriptor : average of N uninitialized analyses
> forecast_time_units : hours
> forecast_time : 0
> level : 0
> parameter_number : 152
> parameter_table_version : 200
> gds_grid_type : 0
> level_indicator : 200
> _FillValue : 1e+20
> units : Kg/m/s
> long_name : DJF: Meridional water vapour flux
> center : Japanese Meteorological Agency - Tokyo (RSMC)
> sub_center : 241
> NMO : 0
>
> Variable: uAvgTime_hi
> Type: float
> Total Size: 23956 bytes
> 5989 values
> Number of Dimensions: 2
> Dimensions and sizes: [g0_lat_1 | 53] x [g0_lon_2 | 113]
> Coordinates:
> g0_lat_1: [90..25]
> g0_lon_2: [250..30]
> Number Of Attributes: 17
> NMO : 0
> sub_center : 241
> center : Japanese Meteorological Agency - Tokyo (RSMC)
> long_name : DJF: Zonal water vapour flux
> units : kg/m/s
> _FillValue : 1e+20
> level_indicator : 200
> gds_grid_type : 0
> parameter_table_version : 200
> parameter_number : 157
> level : 0
> forecast_time : 0
> forecast_time_units : hours
> statistical_process_descriptor : average of N uninitialized analyses
> statistical_process_duration : instantaneous (beginning at
> reference time at intervals of 6 hours)
> N : <ARRAY of 12 elements>
> average_op_ncl : dim_avg_n over dimension(s): initial_time0_hours
>
> Variable: uAvgTime_lo
> Type: float
> Total Size: 23956 bytes
> 5989 values
> Number of Dimensions: 2
> Dimensions and sizes: [g0_lat_1 | 53] x [g0_lon_2 | 113]
> Coordinates:
> g0_lat_1: [90..25]
> g0_lon_2: [250..30]
> Number Of Attributes: 17
> NMO : 0
> sub_center : 241
> center : Japanese Meteorological Agency - Tokyo (RSMC)
> long_name : DJF: Zonal water vapour flux
> units : kg/m/s
> _FillValue : 1e+20
> level_indicator : 200
> gds_grid_type : 0
> parameter_table_version : 200
> parameter_number : 157
> level : 0
> forecast_time : 0
> forecast_time_units : hours
> statistical_process_descriptor : average of N uninitialized analyses
> statistical_process_duration : instantaneous (beginning at
> reference time at intervals of 6 hours)
> N : <ARRAY of 12 elements>
> average_op_ncl : dim_avg_n over dimension(s): initial_time0_hours
>
> Variable: vAvgTime_hi
> Type: float
> Total Size: 23956 bytes
> 5989 values
> Number of Dimensions: 2
> Dimensions and sizes: [g0_lat_1 | 53] x [g0_lon_2 | 113]
> Coordinates:
> g0_lat_1: [90..25]
> g0_lon_2: [250..30]
> Number Of Attributes: 17
> NMO : 0
> sub_center : 241
> center : Japanese Meteorological Agency - Tokyo (RSMC)
> long_name : DJF: Meridional water vapour flux
> units : Kg/m/s
> _FillValue : 1e+20
> level_indicator : 200
> gds_grid_type : 0
> parameter_table_version : 200
> parameter_number : 152
> level : 0
> forecast_time : 0
> forecast_time_units : hours
> statistical_process_descriptor : average of N uninitialized analyses
> statistical_process_duration : instantaneous (beginning at
> reference time at intervals of 6 hours)
> N : <ARRAY of 12 elements>
> average_op_ncl : dim_avg_n over dimension(s): initial_time0_hours
>
> Variable: vAvgTime_lo
> Type: float
> Total Size: 23956 bytes
> 5989 values
> Number of Dimensions: 2
> Dimensions and sizes: [g0_lat_1 | 53] x [g0_lon_2 | 113]
> Coordinates:
> g0_lat_1: [90..25]
> g0_lon_2: [250..30]
> Number Of Attributes: 17
> NMO : 0
> sub_center : 241
> center : Japanese Meteorological Agency - Tokyo (RSMC)
> long_name : DJF: Meridional water vapour flux
> units : Kg/m/s
> _FillValue : 1e+20
> level_indicator : 200
> gds_grid_type : 0
> parameter_table_version : 200
> parameter_number : 152
> level : 0
> forecast_time : 0
> forecast_time_units : hours
> statistical_process_descriptor : average of N uninitialized analyses
> statistical_process_duration : instantaneous (beginning at
> reference time at intervals of 6 hours)
> N : <ARRAY of 12 elements>
> average_op_ncl : dim_avg_n over dimension(s): initial_time0_hours
>
> Variable: diff_u
> Type: float
> Total Size: 23956 bytes
> 5989 values
> Number of Dimensions: 2
> Dimensions and sizes: [g0_lat_1 | 53] x [g0_lon_2 | 113]
> Coordinates:
> g0_lat_1: [90..25]
> g0_lon_2: [250..30]
> Number Of Attributes: 17
> average_op_ncl : dim_avg_n over dimension(s): initial_time0_hours
> N : <ARRAY of 12 elements>
> statistical_process_duration : instantaneous (beginning at
> reference time at intervals of 6 hours)
> statistical_process_descriptor : average of N uninitialized analyses
> forecast_time_units : hours
> forecast_time : 0
> level : 0
> parameter_number : 157
> parameter_table_version : 200
> gds_grid_type : 0
> level_indicator : 200
> _FillValue : 1e+20
> units : kg/m/s
> long_name : DJF: Zonal water vapour flux
> center : Japanese Meteorological Agency - Tokyo (RSMC)
> sub_center : 241
> NMO : 0
>
> Variable: diff_v
> Type: float
> Total Size: 23956 bytes
> 5989 values
> Number of Dimensions: 2
> Dimensions and sizes: [g0_lat_1 | 53] x [g0_lon_2 | 113]
> Coordinates:
> g0_lat_1: [90..25]
> g0_lon_2: [250..30]
> Number Of Attributes: 17
> average_op_ncl : dim_avg_n over dimension(s): initial_time0_hours
> N : <ARRAY of 12 elements>
> statistical_process_duration : instantaneous (beginning at
> reference time at intervals of 6 hours)
> statistical_process_descriptor : average of N uninitialized analyses
> forecast_time_units : hours
> forecast_time : 0
> level : 0
> parameter_number : 152
> parameter_table_version : 200
> gds_grid_type : 0
> level_indicator : 200
> _FillValue : 1e+20
> units : Kg/m/s
> long_name : DJF: Meridional water vapour flux
> center : Japanese Meteorological Agency - Tokyo (RSMC)
> sub_center : 241
> NMO : 0
>
>
>
> On Tue, Mar 7, 2017 at 8:16 AM, Mary Haley <haley at ucar.edu> wrote:
>
>> Mira,
>>
>> In order to generate vectors over a map, you need to call
>> gsn_csm_vector_map, and not gsn_csm_vector map.
>>
>> I think that's the only change you need to make to your script.
>>
>> --Mary
>>
>>
>> On Mon, Mar 6, 2017 at 8:18 PM, MIRA BERDAHL <
>> mlosic at scarletmail.rutgers.edu> wrote:
>>
>>> Hi,
>>> I'm trying to plot moisture transport from JRA55 output. I can generate
>>> figures, where the latitude looks correct, but my longitude is not
>>> correctly labelled. Also, there is no map generated, just the vector
>>> field.
>>> Can anyone see what the problem is with my code? I also get the
>>> following error:
>>> warning:VectorFieldSetValues: irregular coordinate array vfXArray
>>> non-monotonic: defaulting vfXArray
>>> I've copied my code and the variable information below. I've attached
>>> the figure I can generate.
>>> Thank you in advance,
>>> Mira
>>>
>>> ************************************************************
>>> ******************
>>>
>>> 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"
>>> ;************************************************
>>> begin
>>> ;************************************************
>>> ; read in netCDF files
>>> ;************************************************
>>> a = addfile("Uwv.mon.mean.nc","r") ; u moisture transport
>>> b = addfile("Vwv.mon.mean.nc","r") ; v moisture transport
>>>
>>> ;************************************************
>>> ; read in zonal [u] and meridional [v] moisture transport (water vapor)
>>> ;************************************************
>>>
>>> u_wv = a->UWV_GDS0_EATM_S123(1:648,:,:)
>>> v_wv = b->VWV_GDS0_EATM_S123(1:648,:,:) ;
>>>
>>> lat = a->g0_lat_1
>>> lon = a->g0_lon_2
>>>
>>> printVarSummary(lat)
>>> printVarSummary(lon)
>>>
>>> printVarSummary(u_wv)
>>> printVarSummary(v_wv)
>>>
>>>
>>> ; Calculate the seasonal averages.
>>> uwvDJF = month_to_season(u_wv, "DJF")
>>> vwvDJF = month_to_season(v_wv, "DJF")
>>>
>>>
>>> printVarSummary(uwvDJF)
>>> printVarSummary(vwvDJF)
>>>
>>>
>>> ; from the matlab script i wrote: findExtremeYrs, i pulled out the
>>> extreme years (> or < 1std) that i want to average and plot here.
>>>
>>>
>>> ;ans = 7 (longitude of IL)
>>> ;We find the years where the Icelandic low is in an “extreme east” (high
>>> Long) position to be: 1966, 1968, 1969, 1974, 1983, 1984, 1994, 1995, 1999
>>> and 2005.
>>> ;The years with “extreme west” (low long) position of the Icelandic Low
>>> are: 1963, 1964, 1980, 1985, 1987, 1991, 1992, 1996, 2003 and 2006.
>>>
>>>
>>> ; this data starts at 1958 (this is index 0), so 1963=5, 1973=15 etc.
>>>
>>>
>>> uDJF_lon_hi = uwvDJF((/8,10,11,16,25,26,36,37,41,47/),:,:)
>>> uDJF_lon_lo = uwvDJF((/5,6,22,27,29,33,34,38,45,48/),:,:)
>>>
>>> vDJF_lon_hi = vwvDJF((/8,10,11,16,25,26,36,37,41,47/),:,:)
>>> vDJF_lon_lo = vwvDJF((/5,6,22,27,29,33,34,38,45,48/),:,:)
>>>
>>> uAvgTime_hi = dim_avg_n_Wrap(uDJF_lon_hi,0)
>>> uAvgTime_lo = dim_avg_n_Wrap(uDJF_lon_lo,0)
>>>
>>> printVarSummary(uAvgTime_hi)
>>> printVarSummary(uAvgTime_lo)
>>>
>>> vAvgTime_hi = dim_avg_n_Wrap(vDJF_lon_hi,0)
>>> vAvgTime_lo = dim_avg_n_Wrap(vDJF_lon_lo,0)
>>>
>>> printVarSummary(vAvgTime_hi)
>>> printVarSummary(vAvgTime_lo)
>>>
>>>
>>> ; sneaky way to copy metadata over first.
>>> diff_u = uAvgTime_hi;
>>> diff_v = vAvgTime_hi;
>>>
>>> diff_u = uAvgTime_lo - uAvgTime_hi
>>> diff_v = vAvgTime_lo - vAvgTime_hi
>>>
>>> printVarSummary(diff_u)
>>> printVarSummary(diff_v)
>>>
>>> ;************************************************
>>> ; create plot
>>> ;************************************************
>>> wks = gsn_open_wks("ps","Panel_lon_WV") ; open a ps file
>>> gsn_define_colormap(wks,"temp1")
>>>
>>> plot = new(3,graphic) ; create a plot array
>>>
>>> ;---- set common resources for all plots
>>> res = True
>>> res at gsnDraw = False ; dont draw
>>> res at gsnFrame = False ; dont advance frame
>>> res at cnInfoLabelOn = False ; trn off cn info label
>>> res at gsnAddCyclic = False ; has to do with wrapping the lonitude at 0/360
>>> res at cnFillPalette = "matlab_jet"
>>>
>>> ;***********************************************
>>> ; ----wind vector plot
>>> ;***********************************************
>>> vcres = res
>>> vcres at vcRefAnnoOrthogonalPosF = -1.0 ; move ref vector up
>>> vcres at vcRefMagnitudeF = 10.0 ; define vector ref mag
>>> vcres at vcRefLengthF = 0.045 ; define length of vec ref
>>> vcres at vcGlyphStyle = "CurlyVector" ; turn on curly vectors
>>> vcres at vcMinDistanceF = 0.017
>>> vcres at mpFillOn = False ; turn off gray fill
>>> vcres at mpOutlineBoundarySets = "National" ; turn on country boundaries
>>> ;vcres at mpFillBoundarySets = AllBoundaries
>>> vcres at mpGeophysicalLineColor = "Navy" ; color of cont. outlines
>>> vcres at mpGeophysicalLineThicknessF = 1.5 ; thickness of outlines
>>> vcres at gsnRightString = "" ; turn off thue units string
>>>
>>> plot(0) = gsn_csm_vector(wks,uAvgTime_hi,vAvgTime_hi,vcres)
>>> plot(1) = gsn_csm_vector(wks,uAvgTime_lo,vAvgTime_lo,vcres)
>>> plot(2) = gsn_csm_vector(wks, diff_u, diff_v,vcres)
>>>
>>> ;************************************************
>>> ; create panel
>>> ;************************************************
>>> resP = False ; modify the panel plot
>>> resP at txString = "lon - water vapor transport"
>>> gsn_panel(wks,(/plot/),(/3,1/),resP) ; now draw as one plot;
>>>
>>>
>>> end
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>> ************************************************************
>>> *******************
>>> Variable: lat
>>> Type: float
>>> Total Size: 212 bytes
>>> 53 values
>>> Number of Dimensions: 1
>>> Dimensions and sizes: [g0_lat_1 | 53]
>>> Coordinates:
>>> g0_lat_1: [90..25]
>>> Number Of Attributes: 9
>>> La1 : 90
>>> Lo1 : 0
>>> La2 : -90
>>> Lo2 : -1.25
>>> Di : 1.25
>>> Dj : 1.25
>>> units : degrees_north
>>> GridType : Cylindrical Equidistant Projection Grid
>>> long_name : latitude
>>>
>>> Variable: lon
>>> Type: float
>>> Total Size: 452 bytes
>>> 113 values
>>> Number of Dimensions: 1
>>> Dimensions and sizes: [g0_lon_2 | 113]
>>> Coordinates:
>>> g0_lon_2: [250..30]
>>> Number Of Attributes: 9
>>> La1 : 90
>>> Lo1 : 0
>>> La2 : -90
>>> Lo2 : -1.25
>>> Di : 1.25
>>> Dj : 1.25
>>> units : degrees_east
>>> GridType : Cylindrical Equidistant Projection Grid
>>> long_name : longitude
>>>
>>> Variable: u_wv
>>> Type: float
>>> Total Size: 15523488 bytes
>>> 3880872 values
>>> Number of Dimensions: 3
>>> Dimensions and sizes: [initial_time0_hours | 648] x [g0_lat_1 | 53] x
>>> [g0_lon_2 | 113]
>>> Coordinates:
>>> initial_time0_hours: [1385736..1858344]
>>> g0_lat_1: [90..25]
>>> g0_lon_2: [250..30]
>>> Number Of Attributes: 15
>>> N : <ARRAY of 12 elements>
>>> statistical_process_duration : instantaneous (beginning at
>>> reference time at intervals of 6 hours)
>>> statistical_process_descriptor : average of N uninitialized
>>> analyses
>>> forecast_time_units : hours
>>> forecast_time : 0
>>> level : 0
>>> parameter_number : 157
>>> parameter_table_version : 200
>>> gds_grid_type : 0
>>> level_indicator : 200
>>> _FillValue : 1e+20
>>> units : kg/m/s
>>> long_name : Zonal water vapour flux
>>> center : Japanese Meteorological Agency - Tokyo (RSMC)
>>> sub_center : 241
>>>
>>> Variable: v_wv
>>> Type: float
>>> Total Size: 15523488 bytes
>>> 3880872 values
>>> Number of Dimensions: 3
>>> Dimensions and sizes: [initial_time0_hours | 648] x [g0_lat_1 | 53] x
>>> [g0_lon_2 | 113]
>>> Coordinates:
>>> initial_time0_hours: [1385736..1858344]
>>> g0_lat_1: [90..25]
>>> g0_lon_2: [250..30]
>>> Number Of Attributes: 15
>>> N : <ARRAY of 12 elements>
>>> statistical_process_duration : instantaneous (beginning at
>>> reference time at intervals of 6 hours)
>>> statistical_process_descriptor : average of N uninitialized
>>> analyses
>>> forecast_time_units : hours
>>> forecast_time : 0
>>> level : 0
>>> parameter_number : 152
>>> parameter_table_version : 200
>>> gds_grid_type : 0
>>> level_indicator : 200
>>> _FillValue : 1e+20
>>> units : Kg/m/s
>>> long_name : Meridional water vapour flux
>>> center : Japanese Meteorological Agency - Tokyo (RSMC)
>>> sub_center : 241
>>>
>>> Variable: uwvDJF
>>> Type: float
>>> Total Size: 1293624 bytes
>>> 323406 values
>>> Number of Dimensions: 3
>>> Dimensions and sizes: [initial_time0_hours | 54] x [g0_lat_1 | 53] x
>>> [g0_lon_2 | 113]
>>> Coordinates:
>>> initial_time0_hours: [1385736..1850328]
>>> g0_lat_1: [90..25]
>>> g0_lon_2: [250..30]
>>> Number Of Attributes: 16
>>> N : <ARRAY of 12 elements>
>>> statistical_process_duration : instantaneous (beginning at
>>> reference time at intervals of 6 hours)
>>> statistical_process_descriptor : average of N uninitialized
>>> analyses
>>> forecast_time_units : hours
>>> forecast_time : 0
>>> level : 0
>>> parameter_number : 157
>>> parameter_table_version : 200
>>> gds_grid_type : 0
>>> level_indicator : 200
>>> _FillValue : 1e+20
>>> units : kg/m/s
>>> long_name : DJF: Zonal water vapour flux
>>> center : Japanese Meteorological Agency - Tokyo (RSMC)
>>> sub_center : 241
>>> NMO : 0
>>>
>>> Variable: vwvDJF
>>> Type: float
>>> Total Size: 1293624 bytes
>>> 323406 values
>>> Number of Dimensions: 3
>>> Dimensions and sizes: [initial_time0_hours | 54] x [g0_lat_1 | 53] x
>>> [g0_lon_2 | 113]
>>> Coordinates:
>>> initial_time0_hours: [1385736..1850328]
>>> g0_lat_1: [90..25]
>>> g0_lon_2: [250..30]
>>> Number Of Attributes: 16
>>> N : <ARRAY of 12 elements>
>>> statistical_process_duration : instantaneous (beginning at
>>> reference time at intervals of 6 hours)
>>> statistical_process_descriptor : average of N uninitialized
>>> analyses
>>> forecast_time_units : hours
>>> forecast_time : 0
>>> level : 0
>>> parameter_number : 152
>>> parameter_table_version : 200
>>> gds_grid_type : 0
>>> level_indicator : 200
>>> _FillValue : 1e+20
>>> units : Kg/m/s
>>> long_name : DJF: Meridional water vapour flux
>>> center : Japanese Meteorological Agency - Tokyo (RSMC)
>>> sub_center : 241
>>> NMO : 0
>>>
>>> Variable: uAvgTime_hi
>>> Type: float
>>> Total Size: 23956 bytes
>>> 5989 values
>>> Number of Dimensions: 2
>>> Dimensions and sizes: [g0_lat_1 | 53] x [g0_lon_2 | 113]
>>> Coordinates:
>>> g0_lat_1: [90..25]
>>> g0_lon_2: [250..30]
>>> Number Of Attributes: 17
>>> NMO : 0
>>> sub_center : 241
>>> center : Japanese Meteorological Agency - Tokyo (RSMC)
>>> long_name : DJF: Zonal water vapour flux
>>> units : kg/m/s
>>> _FillValue : 1e+20
>>> level_indicator : 200
>>> gds_grid_type : 0
>>> parameter_table_version : 200
>>> parameter_number : 157
>>> level : 0
>>> forecast_time : 0
>>> forecast_time_units : hours
>>> statistical_process_descriptor : average of N uninitialized
>>> analyses
>>> statistical_process_duration : instantaneous (beginning at
>>> reference time at intervals of 6 hours)
>>> N : <ARRAY of 12 elements>
>>> average_op_ncl : dim_avg_n over dimension(s): initial_time0_hours
>>>
>>> Variable: uAvgTime_lo
>>> Type: float
>>> Total Size: 23956 bytes
>>> 5989 values
>>> Number of Dimensions: 2
>>> Dimensions and sizes: [g0_lat_1 | 53] x [g0_lon_2 | 113]
>>> Coordinates:
>>> g0_lat_1: [90..25]
>>> g0_lon_2: [250..30]
>>> Number Of Attributes: 17
>>> NMO : 0
>>> sub_center : 241
>>> center : Japanese Meteorological Agency - Tokyo (RSMC)
>>> long_name : DJF: Zonal water vapour flux
>>> units : kg/m/s
>>> _FillValue : 1e+20
>>> level_indicator : 200
>>> gds_grid_type : 0
>>> parameter_table_version : 200
>>> parameter_number : 157
>>> level : 0
>>> forecast_time : 0
>>> forecast_time_units : hours
>>> statistical_process_descriptor : average of N uninitialized
>>> analyses
>>> statistical_process_duration : instantaneous (beginning at
>>> reference time at intervals of 6 hours)
>>> N : <ARRAY of 12 elements>
>>> average_op_ncl : dim_avg_n over dimension(s): initial_time0_hours
>>>
>>> Variable: vAvgTime_hi
>>> Type: float
>>> Total Size: 23956 bytes
>>> 5989 values
>>> Number of Dimensions: 2
>>> Dimensions and sizes: [g0_lat_1 | 53] x [g0_lon_2 | 113]
>>> Coordinates:
>>> g0_lat_1: [90..25]
>>> g0_lon_2: [250..30]
>>> Number Of Attributes: 17
>>> NMO : 0
>>> sub_center : 241
>>> center : Japanese Meteorological Agency - Tokyo (RSMC)
>>> long_name : DJF: Meridional water vapour flux
>>> units : Kg/m/s
>>> _FillValue : 1e+20
>>> level_indicator : 200
>>> gds_grid_type : 0
>>> parameter_table_version : 200
>>> parameter_number : 152
>>> level : 0
>>> forecast_time : 0
>>> forecast_time_units : hours
>>> statistical_process_descriptor : average of N uninitialized
>>> analyses
>>> statistical_process_duration : instantaneous (beginning at
>>> reference time at intervals of 6 hours)
>>> N : <ARRAY of 12 elements>
>>> average_op_ncl : dim_avg_n over dimension(s): initial_time0_hours
>>>
>>> Variable: vAvgTime_lo
>>> Type: float
>>> Total Size: 23956 bytes
>>> 5989 values
>>> Number of Dimensions: 2
>>> Dimensions and sizes: [g0_lat_1 | 53] x [g0_lon_2 | 113]
>>> Coordinates:
>>> g0_lat_1: [90..25]
>>> g0_lon_2: [250..30]
>>> Number Of Attributes: 17
>>> NMO : 0
>>> sub_center : 241
>>> center : Japanese Meteorological Agency - Tokyo (RSMC)
>>> long_name : DJF: Meridional water vapour flux
>>> units : Kg/m/s
>>> _FillValue : 1e+20
>>> level_indicator : 200
>>> gds_grid_type : 0
>>> parameter_table_version : 200
>>> parameter_number : 152
>>> level : 0
>>> forecast_time : 0
>>> forecast_time_units : hours
>>> statistical_process_descriptor : average of N uninitialized
>>> analyses
>>> statistical_process_duration : instantaneous (beginning at
>>> reference time at intervals of 6 hours)
>>> N : <ARRAY of 12 elements>
>>> average_op_ncl : dim_avg_n over dimension(s): initial_time0_hours
>>>
>>> Variable: diff_u
>>> Type: float
>>> Total Size: 23956 bytes
>>> 5989 values
>>> Number of Dimensions: 2
>>> Dimensions and sizes: [g0_lat_1 | 53] x [g0_lon_2 | 113]
>>> Coordinates:
>>> g0_lat_1: [90..25]
>>> g0_lon_2: [250..30]
>>> Number Of Attributes: 17
>>> average_op_ncl : dim_avg_n over dimension(s): initial_time0_hours
>>> N : <ARRAY of 12 elements>
>>> statistical_process_duration : instantaneous (beginning at
>>> reference time at intervals of 6 hours)
>>> statistical_process_descriptor : average of N uninitialized
>>> analyses
>>> forecast_time_units : hours
>>> forecast_time : 0
>>> level : 0
>>> parameter_number : 157
>>> parameter_table_version : 200
>>> gds_grid_type : 0
>>> level_indicator : 200
>>> _FillValue : 1e+20
>>> units : kg/m/s
>>> long_name : DJF: Zonal water vapour flux
>>> center : Japanese Meteorological Agency - Tokyo (RSMC)
>>> sub_center : 241
>>> NMO : 0
>>>
>>> Variable: diff_v
>>> Type: float
>>> Total Size: 23956 bytes
>>> 5989 values
>>> Number of Dimensions: 2
>>> Dimensions and sizes: [g0_lat_1 | 53] x [g0_lon_2 | 113]
>>> Coordinates:
>>> g0_lat_1: [90..25]
>>> g0_lon_2: [250..30]
>>> Number Of Attributes: 17
>>> average_op_ncl : dim_avg_n over dimension(s): initial_time0_hours
>>> N : <ARRAY of 12 elements>
>>> statistical_process_duration : instantaneous (beginning at
>>> reference time at intervals of 6 hours)
>>> statistical_process_descriptor : average of N uninitialized
>>> analyses
>>> forecast_time_units : hours
>>> forecast_time : 0
>>> level : 0
>>> parameter_number : 152
>>> parameter_table_version : 200
>>> gds_grid_type : 0
>>> level_indicator : 200
>>> _FillValue : 1e+20
>>> units : Kg/m/s
>>> long_name : DJF: Meridional water vapour flux
>>> center : Japanese Meteorological Agency - Tokyo (RSMC)
>>> sub_center : 241
>>> NMO : 0
>>>
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>>
>
> _______________________________________________
> 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/20170307/89b7a9a8/attachment-0001.html
More information about the ncl-talk
mailing list