[gradsusr] How to highlight multiple model domains in one figure

[Jennifer M Adams]

Hi, Mano —
If you want to fill the area of the grid domain with a color of your choice, then you can do this:
‘set rgb 20 100 100 100 100’  ;* a transparent grey
‘set ccolor 20’
‘d const(var,0)’

If you want to draw an outline, then here is an algorithm to try:
1. Draw a blank plot covering the largest domain to set up the scaling environment
2. Construct a ‘draw line’ command by looping over the grid coordinates at the boundaries (e.g. x=1,y=1:YDEF, y=1,x=1:XDEF, etc.) and using the 'q this2that’ commands (either 'q gr2xy’ or ‘q w2xy’) to get the location in X,Y space of each grid point.

cmd=‘draw line ‘
loop over edges
   cmd=cmd%xpos’ ‘ypos’ '
endloop
cmd

—Jennifer

On Oct 15, 2016, at 10:41 AM, Jeff Duda <jeffduda319 at gmail.com<mailto:jeffduda319 at gmail.com>> wrote:

There's no simple or easy way to do it. Here are two methods I have used in the past:

1) Generate a netcdf grid using whatever library you're familiar with (if any). Make the grid the same size as your model domain. You can make the values inside the grid essentially whatever you want so long as you can bound the outer edge of the domain in a single contour. I suggest making all values in the domain either 0.0 or 1.0. Then write a control file to display that grid and make sure the value of UNDEF is set to something other than the value on the grid. Then just display the field using gxout contour and with contour labeling set to off.

2) Use gxout fwrite and define any field on that domain using const like 'define field = const(tmp,1,-a)' which essentially does the same thing as in option 1, except it will write the field to a binary file instead of a netcdf file. Then you can just modify the control file you used for the original data set (assuming you used one), just changing the DSET, DTYPE, and VARS entries to suit your new data. Whatever the name of the field you wrote (in this example, it was defined using the name 'field'), the VARS entry will become
VARS 1
field 0 99 field
ENDVARS

Also make sure to set UNDEF to some value other than the constant value (1 in this example). Then, as in the first suggestion, display the field using gxout contour and with contour labeling off.

Once you have the general idea here, you can manipulate the steps to make small alterations to improve the output image.

Jeff Duda

On Fri, Oct 14, 2016 at 7:59 PM, lpasmanoranjan <lpasmanoranjan at gmail.com<mailto:lpasmanoranjan at gmail.com>> wrote:
Dear GrADS users,
I would like to ask you a favour.

I would like to highlight the following three model domains in one figure around my study region (India as center). But I can't able to draw the closed polygons as in the attachment.

I have three model domains in Lambert conformal conic projection:
MODEL 1: (8S,50E), (8S,90E), (45N,95E), (45N,45E).
MODEL 2: (7S,42E), (7S,92E), (47N,97E), (47N,41E).
MODEL 3:

Would you please suggest how to do it in grads!!

Thank you very much for any kind of help or suggestion.


--
Kind Regards,
Mano


_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org<mailto:gradsusr at gradsusr.org>
http://gradsusr.org/mailman/listinfo/gradsusr




--
Jeff Duda
Post-doctoral research associate
University of Oklahoma School of Meteorology
_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org<mailto:gradsusr at gradsusr.org>
http://gradsusr.org/mailman/listinfo/gradsusr

--
Jennifer Miletta Adams
Center for Ocean-Land-Atmosphere Studies (COLA)
George Mason University



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20161021/36c33e43/attachment.html