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

Sergey Vinogradov, PhD svinogra at aer.com
Thu May 21 22:12:32 EDT 2009


Nicolas
First guess: endianness of the files (big endian vs little endian).
Your forcing files endianness may be different from your platform's.
This may lead to this sort of problems.

Sent from my iPhone 3G

On May 21, 2009, at 8:05 PM, Nico Wienders <wienders at ocean.fsu.edu>  
wrote:

>
> Comrades,
>
> 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=801*1,
> OB_Ieast=400*801,
>
>
> # 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



More information about the MITgcm-support mailing list