[MITgcm-support] cal, exf, obcs: problem at the tile boundaries
Nico Wienders
wienders at ocean.fsu.edu
Thu May 21 20:05:22 EDT 2009
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)
More information about the MITgcm-support
mailing list