[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