[MITgcm-support] cal, exf, obcs: problem at the tile boundaries

[Nico Wienders]

My apologies for insisting on this matter, but i still haven't been  
able to find what
is the problem here. The errors i am getting through exf_check_range.F  
are very localized.

My domain is as said 400*801.

	If i use 12 proc, NPx=4 Npy=3
	The error are localized on the western side of the first tile (1,1),  
at three different locations
	in longitude. (cf. 4x3.jpg).


	if i use 10 procs, NPx=10 Npy=3
	The error is still localized on the western side of the first tile  
(1,1) but now at several locations
	(cf. 10x1.jpg).

	if i use 9 procs, NPx=1 Npy=9
	The error is localized on the western side of the first tile (1,1) at  
one location.
	(cf. 1*9.jpg)
	

Curiously the spreading of this error becomes bigger in longitude if i  
reduce the number of forcing fields.
Thank you.




-------------- next part --------------
A non-text attachment was scrubbed...
Name: 4x3.jpg
Type: image/jpeg
Size: 29920 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20090603/500b272f/attachment.jpg>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 10x1.jpg
Type: image/jpeg
Size: 29070 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20090603/500b272f/attachment-0001.jpg>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 1x9.jpg
Type: image/jpeg
Size: 26182 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20090603/500b272f/attachment-0002.jpg>
-------------- next part --------------



On May 21, 2009, at 8:05 PM, Nico Wienders wrote:
>
> I will bother you once again with one of my stupid questions
> and, in advance, deeply apologize for making you waste your
> precious time. But the truth is that i am actually stuck with an
> error which i haven't been able to solve. It should be something
> quite silly i assume.
>
> In brief:
>
> I am running an experiment with initial conditions, open boundaries
> and climatological forcings. I get an error that my forcing fields  
> have
> incorrect values on the edge of the tiles, which is not true.
> Any idea? Thank you. data.obcs, data.exf and SIZE.h are attached.
>
>
> Error message exemple:
>
> EXF WARNING: sflux out of range for bi,bj,i,j,it=          1       
> 1     89    267      0   0.218961325573595E+00
> EXF WARNING: vstress out of range for bi,bj,i,j,it=       1       
> 1     90    267      0   0.260417995106866E+31
> EXF WARNING: swflux out of range for bi,bj,i,j,it=        1       
> 1     90    267      0   0.481624151346128E+28
> EXF WARNING: sflux out of range for bi,bj,i,j,it=          1       
> 1     91    267      0   0.218993820760348E+00
> EXF WARNING: sflux out of range for bi,bj,i,j,it=          1       
> 1     92    267      0  -0.416280946298594E+33
>
>
> #DATA.OBCS
>
> # Open-boundaries
> &OBCS_PARM01
> OB_Jsouth=400*1,
> OB_Ieast=801*400,
>
>
> # useOBCSprescribe = .TRUE.,
> useOBCSsponge=.TRUE.,
>
>  OBEsFile='input_files/salteast',
>  OBEtFile='input_files/tempeast',
>  OBEuFile='input_files/ueast',
>  OBEvFile='input_files/veast',
>  OBSsFile='input_files/saltsouth',
>  OBStFile='input_files/tempsouth',
>  OBSuFile='input_files/usouth',
>  OBSvFile='input_files/vsouth',
>  OBNsFile='input_files/saltnorth',
>  OBNtFile='input_files/tempnorth',
>  OBNuFile='input_files/unorth',
>  OBNvFile='input_files/vnorth',
> &
>
> # Orlanski parameters
> &OBCS_PARM02
> &
>
> # Sponge Layer Parameters.
> &OBCS_PARM03
> Urelaxobcsinner=10.E0,
> Urelaxobcsbound=5.E0,
> Vrelaxobcsinner=10.E0,
> Vrelaxobcsbound=5.E0,
> spongeThickness=8,
> &
>
>
> ---------------------------------------------------------
>
>
> #DATA.EXF
>
> # External Forcing Data
> # *********************
> #
> &EXF_NML_01
> readStressOnCgrid = .TRUE.,
> &
>
>
> &EXF_NML_02
> hfluxfile='/home/c1/bill/MITgcm/verification/exp4/input/bigjapan/ 
> forcing/seas/shf',
> hfluxstartdate1=19911215,
> hfluxstartdate2=000000,
> hfluxperiod= 7884000.0,
> sfluxfile='/home/c1/bill/MITgcm/verification/exp4/input/bigjapan/ 
> forcing/seas/swf',
> sfluxstartdate1=19911215,
> sfluxstartdate2=000000,
> sfluxperiod= 7884000.0,
> swfluxfile='/home/c1/bill/MITgcm/verification/exp4/input/bigjapan/ 
> forcing/seas/srf',
> swfluxstartdate1=19911215,
> swfluxstartdate2=000000,
> swfluxperiod= 7884000.0,
> ustressfile='/home/c1/bill/MITgcm/verification/exp4/input/bigjapan/ 
> forcing/seas/ustress',
> ustressstartdate1=19911215,
> ustressstartdate2=000000,
> ustressperiod= 7884000.0,
> vstressfile='/home/c1/bill/MITgcm/verification/exp4/input/bigjapan/ 
> forcing/seas/vstress',
> vstressstartdate1=19911215,
> vstressstartdate2=000000,
> vstressperiod= 7884000.0,
> &
>
> &EXF_NML_03
> &
>
> &EXF_NML_04
> &
>
> &EXF_NML_OBCS
> obcsSstartdate1   = 19911215,
> obcsSstartdate2   = 000000,
> obcsSperiod       = 7884000.0,
>
> obcsEstartdate1   = 19911215,
> obcsEstartdate2   = 000000,
> obcsEperiod       = 7884000.0,
>
> obcsWstartdate1   = 19911215,
> obcsWstartdate2   = 000000,
> obcsWperiod       = 7884000.0,
>
> ---------------------------------------------------------
>
> #SIZE.h
>
>     &           sNx = 100,
>     &           sNy = 267,
>     &           OLx =   4,
>     &           OLy =   4,
>     &           nSx =   1,
>     &           nSy =   1,
>     &           nPx =   4,
>     &           nPy =   3,
>     &           Nx  = sNx*nSx*nPx,
>     &           Ny  = sNy*nSy*nPy,
>     &           Nr  =   45)
>
>
>
>
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support