<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN">
<html><body style='font-size: 10pt; font-family: Tahoma,Arial,Helvetica,sans-serif'>
<p>Hi all, </p>
<p>Using the cfsv2 45 daily, how would one go about breaking down the pratesfc in to 6hr charts? Sample of scrips and ctr file below:</p>
<p> </p>
<p>CTR File:</p>
<p>title CFS_grib.grib2<br />* produced by g2ctl v0.0.8.8<br />* command line options: -verf CFS_grib.grib2<br />* griddef=1:0:(360 x 181):grid_template=0:winds(N/S): lat-lon grid:(360 x 181) units 1e-06 input WE:NS output WE:SN res 48 lat 90.000000 to -90.000000 by 1.000000 lon 0.000000 to 359.000000 by 1.000000 #points=65160:winds(N/S)</p>
<p>dtype grib2<br />ydef 181 linear -90.000000 1<br />xdef 360 linear 0.000000 1.000000<br />tdef 180 linear 06Z10may2016 6hr<br />* PROFILE hPa<br />zdef 3 levels 85000 50000 25000<br />options pascals<br />vars 12<br />HGT500mb 0,100,50000 0,3,5 ** 500 mb Geopotential Height [gpm]<br />PRATEsfc 0,1,0 0,1,7 ** surface Precipitation Rate [kg/m^2/s]<br />PRMSLmsl 0,101,0 0,3,1 ** mean sea level Pressure Reduced to</p>
<p> </p>
<p>Script:</p>
<p>function main(args)</p>
<p>* Parse the arguments: date, hour<br /> if (args = '')<br /> prompt 'Enter forecast date (example, 20110717) --> ' <br /> pull date<br /> prompt 'Enter forecast hour (example, 00 or 06 or 12 or 18) --> ' <br /> pull hour<br /> else <br /> date = subwrd(args,1)<br /> hour = subwrd(args,2)<br /> endif</p>
<p>'reinit'<br />'open CFS_grib.ctl'</p>
<p> </p>
<p><br />* Get info from the descriptor file<br />'q ctlinfo'<br />_ctl = result<br />_undef = getctl(undef)<br />_tdef = getctl(tdef)<br />_zdef = getctl(zdef)</p>
<p>maps = 181</p>
<p>* Get the Time axis info</p>
<p><br />_t1 = 1 ;* 2nd half of timeseries<br />_t2 = 45<br />tsize = _t2 - _t1 + 1<br />'set t '_t1' '_t2<br />'q dims'<br />times = sublin(result,5)<br />_time1 = subwrd(times,6) <br />_time2 = subwrd(times,8)<br />_tdim = _time1' '_time2</p>
<p>tincr = subwrd(_tdef,5)<br />_tdef = 'tdef 'tsize' linear '_time1' 'tincr</p>
<p>huh = subwrd(_tdef,4)</p>
<p>****************************************<br />* PRECIP & SLP<br />****************************************<br /> i = 2<br /> while ( i<maps )</p>
<p>* Set up a few preliminary characteristics<br />setmap(0)<br />setPRECIPcols(1)</p>
<p>'set t 'i</p>
<p>'set gxout shaded'<br />'set csmooth on'<br />'set clevs 0 0.05 0.25 0.5 0.75 1.5 2 2.5 3 4 5 6 7 8 10 12 14 16 18 20 25 30 35 40 45 50'<br />'set ccols 95 95 95 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85'</p>
<p>'d PRATEsfc'</p>
<div>-- <br />
<p>Regards,</p>
<p>Stuart Markham</p>
<p> </p>
</div>
</body></html>