[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