Re: Plot from WRF Grib files / renaming variables

Hello Dennis,

this was exactly the information I was looking for - this was missing
in my script:

wind_speed@lat2d = lat2d
wind_speed@lon2d = lon2d

FYI, I have certain lat/lon boundaries that I need to use in
mpLimitMode = "LatLon". Now I am able to set these without errors and
create my plots as I need.

Thank you for your help!

Best regards,

Jonas Kaufmann

On Tue, Jul 3, 2012 at 3:58 PM, Dennis Shea <shea@ucar.edu> wrote:
> Based upon
>
> GridType : Lambert Conformal Secant or Tangent, Conical or bipolar
> long_name : latitude
> corners : ( 31.373, 31.37483, 57.37227, 57.36925 )
>
> GridType : Lambert Conformal Secant or Tangent, Conical or bipolar
> long_name : longitude
> corners : ( -10.177, 29.7688, 42.64273, -23.05522 )
>
>
>
>> res1@mpProjection = "Mercator"
>> res1@mpLimitMode = "Corners"
>> res1@mpLeftCornerLatF = lat2d@corners(0) ; 31.373
>> res1@mpLeftCornerLonF = lon2d@corners(0) ; -10.177
>> res1@mpRightCornerLatF = lat2d@corners(2) ; 57.37227
>> res1@mpRightCornerLonF = lon2d@corners(2) ; 42.6273
>
> should be all that you need. There is no need to set
> mpLimitMode = "LatLon" or mpMinLatF, etc
> In fact, I don't think you want to.
>
> Also, it is assumed you set
>
> wind_speed@lat2d = lat2d
> wind_speed@lon2d = lon2d
>
>
> On 7/3/12 5:30 AM, Jonas Kaufmann wrote:
>>
>> Hello,
>>
>> I am trying to create a plot based on my WRF run. I am not using the
>> direct WRF output files, but the converted Grib files instead and
>> convert them again to NetCDF. My WRF domain is set up over Europe and
>> I want to plot the area in Mercator projection.
>>
>> I can successfully plot the whole area with (example only):
>>
>> fin = addfile(grib_file,"r")
>> wind_u = fin->U_GRD_GDS3_HTGL(:,:)
>> wind_v = fin->V_GRD_GDS3_HTGL(:,:)
>> wind_speed = sqrt(wind_u^2+wind_v^2)
>> lat2d = fin->g3_lat_0
>> lon2d = fin->g3_lon_1
>> [...]
>> res1@mpProjection = "Mercator"
>> res1@mpLimitMode = "Corners"
>> res1@mpLeftCornerLatF = lat2d@corners(0)
>> res1@mpLeftCornerLonF = lon2d@corners(0)
>> res1@mpRightCornerLatF = lat2d@corners(2)
>> res1@mpRightCornerLonF = lon2d@corners(2)
>> [...]
>> plot = gsn_csm_vector_scalar_map(wks, wind_u, wind_v, wind_speed, res1)
>>
>> ...however I cannot limit the range of the plot using mpLimitMode =
>> "LatLon" and setting mpMinLatF etc. This always gives me the error:
>> check_for_y_lat_coord: Warning: Data either does not contain a valid
>> latitude coordinate array or doesn't contain one at all.
>> [...]
>>
>> According to the error message FAQ
>> (http://www.ncl.ucar.edu/Document/Graphics/error_msg.shtml), I should
>> rename the variable names to the automatically found variable names
>> like "lat" and "lon". However I do not manage to get this to run
>> correctly.
>>
>> What is the correct syntax to rename the variables and display my plot
>> with mpLimitMode = "LatLon"? As an attachment you can find a
>> ncl_filedump of the WRF file. Thank you for your help!
>>
>>
>> Best regards,
>>
>>
>> Jonas Kaufmann
>>
>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Jul 3 13:34:26 2012