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

Thu Jun 4 02:57:39 EDT 2009

Nico,

going back to your initial message, if think that this is the clue:
> EXF WARNING: vstress out of range for bi,bj,i,j,it=       1       
> 1     90    267      0   0.260417995106866E+31
exf_check_range often complains about too high heat fluxes (and then  
they are too high by maybe a factor of 2 or so), but I have never seen  
it complain about u or vstress, and 1e30 points to a single/double or  
little/big endian issue. I think it is the former, because I think  
that you are using the default value of  exf_iprec = 32  
(exf_readparms.F), but you said that your input fields are real*8  
(which would correspond to exf_iprec=64, which should show up in your  
data.exf file).

Martin

On Jun 3, 2009, at 6:29 PM, Nico Wienders wrote:

> 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.
>
>
>
>
> <4x3.jpg><10x1.jpg><1x9.jpg>
>
>
> 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
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support