[MITgcm-support] format of input files wrong in regional model with open boundary conditions?

Jonny Williams Jonny.Williams at bristol.ac.uk
Fri Dec 13 06:00:19 EST 2013


Thanks for that

For your information, currently I am specifying the external forcing in the
following way in the EXF_NML_04 namelist in data.exf, following numerous
examples in the verification examples

#
 uwind_lon0        = -30.0D0,
 uwind_lon_inc     = 0.1D0,
 uwind_lat0        = -70.0D0,
 uwind_lat_inc     = 0.1D0,
 uwind_nlon        = 450,
 uwind_nlat        = 800,
#

Thank you!

Jonny







On 13 December 2013 03:46, Michael Schaferkotter <schaferk at bellsouth.net>wrote:

> Jonny:
>
> This 'new' issue has nothing to do with loading the input grid values via
> namelist or file.
>
> So you don't need my suggestion at this time.
>
> Clearly the ATM forcing is valued at points outside the domain that you
> are specifying to the
> inline interpolation code, ie.,exf_interp, which is enabled as a package.
> I would check the namelist parameters associated with the inline
> interpolation of the atmospheric forcing.
>
>
>
>
>
>
> Sent from an iOS device.
>
> On Dec 12, 2013, at 12:36, Jonny Williams <Jonny.Williams at bristol.ac.uk>
> wrote:
>
> Thanks very much for that Michael!
>
> I have now changed my SIZE.h file and I am now getting a different error,
> which is good I suppose!..
>
> *(PID.TID 0000.0001) *** ERROR *** EXF_INTERP_UV:
> fileU="u_992hPa.bin_1979", rec=     1*
> *(PID.TID 0000.0001) *** ERROR *** EXF_INTERP_UV: input grid must
> encompass output grid.*
> *(PID.TID 0000.0001) *** ERROR *** i,j,bi,bj=       1       1     1     1
> , xG,yG= -2.995000E+01 -6.995000E+01*
> *(PID.TID 0000.0001) *** ERROR *** nxIn=      450 , x_in(0,nxIn+1)=
> -3.010000E+01  1.500000E+01*
> *(PID.TID 0000.0001) *** ERROR *** nyIn=      800 , y_in(0,nyIn+1)=
> -7.010000E+01 -1.903000E+02*
>
> The key point here is that the x_in vales are seemingly correct (-30 to
> +15 degrees) but the y_in values are completely wrong (-70 to -190 degrees
> when they should be -70 to +10 degrees). This is not at all what I want and
> is surely unphysical. I will try your idea of using a delXfile etc.
>
> Many thanks indeed for this. I feel I am getting somewhere now with all
> this help which is hugely appreciated.
>
> Jonny
>
>
>
>
>
> On 12 December 2013 18:11, michael schaferkotter <schaferk at bellsouth.net>wrote:
>
>> you/ll need to have sNy to be 80 rather than 10.
>>
>> as SIZE.h indicates, the model thinks the domain is 450 x 100 x 50.
>>
>> also namelist parameter lists can be problematic for different platforms,
>> especially when listing the elements, ie *delX=450*0.1, or **delZ   =
>> 10.00, 10.00, 10.00, 1 ....*
>>
>> *i would run into issues with that often. *
>>
>> *so to avoid those types of issues (platform), i would use *
>>
>> *delXfile='dx.bin',*
>> *delYfile='dy.bin',*
>> *delZfile='dr.bin',  *
>>
>> *where the elements are written into the file using matlab or an f90
>> program.*
>>
>>  *thus avoiding the persnickety namelists with single entries (you have
>> to get the count right, which you don't have with your current SIZE.h)*
>>
>> *here/s a MATLAB snippet (note the precision!)*
>>
>>
>>
>>
>>
>>
>>
>> *prec = 'real*8';ieee = 'b';ny=800;dy=zeros(ny,1);for j=1:ny  dy(j)=0.1;
>> endfid=fopen('dy.bin','w',ieee); fwrite(fid,dy,prec); fclose(fid);michael*
>>
>> Begin forwarded message:
>>
>> Thanks for your suggestion Jody
>>
>> I think I have ironed out the SIZE.h aspect now, I am now using the
>> following...
>>
>> *     &           sNx =  45,*
>> *     &           sNy =  10,*
>> *     &           OLx =   4,*
>> *     &           OLy =   4,*
>> *     &           nSx =   1,*
>> *     &           nSy =   1,*
>> *     &           nPx =  10,*
>> *     &           nPy =  10,*
>> *     &           Nx  = sNx*nSx*nPx,*
>> *     &           Ny  = sNy*nSy*nPy,*
>> *     &           Nr  =   50)*
>>
>> However I still cannot get the model running so I have turned off the
>> OBCS package for now only to encounter another problem, noting that my
>> domain is 450 x 800 in total
>>
>>
>>
>
>
> --
> Dr Jonny Williams
> School of Geographical Sciences
> University of Bristol
> University Road
> BS8 1SS
>
> +44 (0)117 3318352
> jonny.williams at bristol.ac.uk
> bit.ly/jonnywilliams
>
>


-- 
Dr Jonny Williams
School of Geographical Sciences
University of Bristol
University Road
BS8 1SS

+44 (0)117 3318352
jonny.williams at bristol.ac.uk
bit.ly/jonnywilliams
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20131213/99670a42/attachment-0001.htm>


More information about the MITgcm-support mailing list