[RegCNET] Model Output CTL files

[XUNQIANG BI]

Hi, Ashley:

If you just want the LANDUSE field, there is a ASCII file (named
LANDUSE) under RegCM/PreProc/Terrain,

And if you also want other fields, I think you'd better use a
FORTRAN code to read them (in direct access big-endian binary format).

Or you can change the DOMAIN.CTL file by getting rid of the pdef and
reset the xdef as 95 and ydef as 70. Since the model grid is actually 
with 95*70 equal-area meshes, so the figure you get would be distorted.
But you can print the real-field on screen (by set gxout grid or fgrid).

On Wed, 26 Sep 2007, Ashley Watson wrote:

>  I am running my model at 95 x 70 grid  (iy=70, jx=95, ds=45km  in
> domain.param). My output DOMAN.CTL says:
>
> xdef  373 linear  -38.34  0.2027
> ydef  159 linear   40.84  0.2027
>
> I am wondering why there are more grid points in CTL files than actual and
> why there is need to do linear interpolation,

You need read the GrADS documentation about pdef to find the reason.

Best regards,
Xunqiang Bi
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   Dr. Xunqiang Bi         email:bixq at ictp.it
   Earth System Physics Group
   The Abdus Salam ICTP
   Strada Costiera, 11
   P.O. BOX 586, 34100 Trieste, ITALY
   Tel: +39-040-2240302  Fax: +39-040-2240449
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~