[gradsusr] ctl development

[James T. Potemra]

Stephen,

You seem to be mixing the syntax of GrADS commands, e.g., "set lon 80 
180", and descriptor file syntax.  These are commands that you can enter 
via the command line or in a script:

> clear
> set lon 80 180
> set lat -70 0
> set levels 1000  950 900 850 800 750 700 650 600 550 500 450 400 350 
> 300 250 200 150 100 75 50 25 10
> set cint30
> set gxout shaded
> set csmooth on
> set clevs 4920 4960 5000 5040 5080 5120 5160 5200 5240 5280
>  5320  5360 5400 5440 5480 5520 5560 5600 5640 5680 5720 5760 5800 
> 5840 5880 5920 5960 6000 6040 6080 6120'
> set ccols 61  63  65  67  69  71  73  75  77  79
> 81  83  85  87  89  18  20  22  24  26  28  30  32  34
> 36  38  41  43  45  47  49  51
> d hcurl(u,v)
While these are descriptor file contents:
> dset ^model.grb
> title "thickness 1000 500"
> vars 8
> ps        0    1,  1,  0,  0 Surface pressure [hPa]
> ua        7   33,100 Eastward wind [m/s]
> va        7   34,100 Northward wind [m/s]
> zg        7    7,100 Geopotential height [m]
> ta        7   11,100 Air Temperature [K]
> hus       7   51,100 Specific humidity [kg/kg]
> ts        0   11,105,  2 Surface (2m) air temperature [K]
> pr        0   59,  1,  0,  0 Total precipitation rate [kg/(m^2*s)]
> endvars
Jim


On 4/27/15 8:58 PM, Steve Sgalowski wrote:
>
> australia.ctl
> looking to create a australia.ctl file ,
> so that i an run scripts on top level, but pull out data from ncep 
>  useing the australia.ctl file as the main descriptor .
>
> wondering if i have the setup corectly worded, and coded right
>
> Stephen
>
> ***************************************************************
> clear
> dset ^model.grb
> title "thickness 1000 500"
> set lon 80 180
> set lat -70 0
> set levels 1000  950 900 850 800 750 700 650 600 550 500 450 400 350 
> 300 250 200 150 100 75 50 25 10
> set cint30
> set gxout shaded
> set csmooth on
> set clevs 4920 4960 5000 5040 5080 5120 5160 5200 5240 5280
>  5320  5360 5400 5440 5480 5520 5560 5600 5640 5680 5720 5760 5800 
> 5840 5880 5920 5960 6000 6040 6080 6120'
> set ccols 61  63  65  67  69  71  73  75  77  79
> 81  83  85  87  89  18  20  22  24  26  28  30  32  34
> 36  38  41  43  45  47  49  51
> d hcurl(u,v)
> vars 8
> ps        0    1,  1,  0,  0 Surface pressure [hPa]
> ua        7   33,100 Eastward wind [m/s]
> va        7   34,100 Northward wind [m/s]
> zg        7    7,100 Geopotential height [m]
> ta        7   11,100 Air Temperature [K]
> hus       7   51,100 Specific humidity [kg/kg]
> ts        0   11,105,  2 Surface (2m) air temperature [K]
> pr        0   59,  1,  0,  0 Total precipitation rate [kg/(m^2*s)]
> endvars
> ************************************************************************************
>
>
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20150428/37e04896/attachment.html