fwrite bends values

[Mary Jo Nath]

Try using "set x 1 360" instead of "set lon -180 180" when you do the
"fwrite". You are probably the equivalent of "set x 1 361" when you use
"set lon -180 180". (Do "q dims" to see what X is set to.)

MJ

Erko Jakobson wrote:
> Hello
>
> I have problem that when I save my results using fwrite, then it bends it,
> so the new picture is a little rotated and stretched out (I use mproj). I
> found it so that the Greenland temperature minimum was moved to the
> Greenland Sea.
> Without 'set gxout fwrite' it showed the minimum directly on the Greenland.
> Have anyone had similar problems and where might be the error source?
>
> ***********************
> *Writing script:
>
> 'open /home/erkok/Desktop/data/Temperature/ta.ctl'
> 'set lon -180 180' ; 'set mproj nps'
> 'set mpdset mres' ; 'set gxout shaded'
>      'set fwrite -le -sq /home/erkok/Desktop/data/qvaver3/temp.dat'
>      'set gxout fwrite'
>           'd thbl'
>      'disable fwrite'
> *************************
>
> *Original ta.ctl
> dset ^ta%y4%m2%d2+%h2.ml
> index ^ta.idx
> options template
> undef 9.999E+20
> title ta19580101+00.ml
> *  produced by grib2ctl v0.9.12.5p32y
> dtype grib 255
> options yrev
> ydef 36 linear 55.000000 1
> xdef 360 linear -180.000000 1.000000
> tdef 64284 linear 00Z01jan1958 6hr
> zdef 1 linear 1 1
> vars 1
> Thbl  0 130,109,60  ** Temperature K
> ENDVARS
> ***********************
>
> *New ctl to open the temp.dat
> dset ^temp.dat
> undef 9.999E+20
> title temp
> options sequential little_endian
> ydef 36 linear 55.000000 1
> xdef 360 linear -180.000000 1.000000
> tdef 1 linear 00Z01jan1958 1mo
> zdef 1 levels 1 1
> vars 1
> temp2  27 122,109,0 ** temperature
> ENDVARS
> *************************
>