[MITgcm-support] sea ice experiment

Amitabh Mitra amitabhmitra1 at rediffmail.com
Wed Oct 31 04:54:33 EDT 2007


  
 Sir,

  I am trying to run a sea ice experiment for

 the Antarctic region -70 to -40 South and 

 -178 to 178 East with 1x1 resolution. The 
 
 corresponding 'data' file is attached.

 It runs for 10 days but when I extend it 

 to 20 days the model blows with the following error:

 SGI_2100 11% mpirun -np 1 ./mitgcmuv
 **WARNNING** CONFIG_CHECK: useMNC is TRUE and #undef  

 ALLOW_MNC
 STOP NORMAL END
 STOP
 SGI_2100 12% mpirun -np 1 ./mitgcmuv
 **WARNNING** CONFIG_CHECK: useMNC is TRUE and #undef  

 ALLOW_MNC
 SOLUTION IS HEADING OUT OF BOUNDS: tMin,tMax=   

 -2.478E+03   1.998E+03
 MON_SOLUTION: STOPPING CALCULATION
 S/R EEDIE: Only     0 threads have completed,    
 1  are expected for this configuration!



  Possibly you have different setenv PARALLEL and   

 nThreads?
 STOP MON_SOLUTION: STOPPED DUE TO EXTREME VALUES OF   

 SOLUTION
 STOP
 SGI_2100 13%
------------------------------------------------------


 Please suggest suitable values for these coefficients

 viscAz, viscAh, diffKhT, diffKzT, diffKhS, diffKzS

 and deltaT so that the model runs efficiently.

 I have read the manual but it is not helpful

 in determining these values. Do you have any method

 to calculate these values or it is just trial and 
 
 error method ?

 Thanking you,

 With Regards,

 Amitabh Mitra
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20071031/8495abcf/attachment.htm>
-------------- next part --------------
# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
#
#   tRef            - Reference vertical potential temperature          (deg C)
#   sRef            - Reference vertical salinity                         (PSU)
#   viscAh          - Horizontal eddy viscosity coefficient             (m^2/s)
#   viscAz          - Vertical eddy viscosity coefficient               (m^2/s)
#   diffKhT         - Horizontal temperature diffusivity                (m^2/s)
#   diffKzT         - Vertical temperature diffusivity                  (m^2/s)
#   diffKhS         - Horizontal salt diffusivity                       (m^2/s)
#   diffKzS         - Vertical salt diffusivity                         (m^2/s)
#   f0              - Reference coriolis parameter,
#                     south edge of f on beta plane                       (1/s)
#   beta            - df/dy                                         (s^-1.m^-1)
#   tAlpha          - Linear EOS thermal expansion coefficient           (1/oC)
#   sBeta           - Linear EOS haline contraction coefficient         (1/ppt)
#   gravity         - Acceleration due to gravity                       (m/s^2)
#   gBaro           - Accel. due to gravity used in barotropic equation (m/s^2)
#   rigidLid        - Set to true to use rigid lid
#   implicitFreeSurface - Set to true to use implicit free surface
#   eosType         - Flag for linear or polynomial equation of state
#   momAdvection    - On/Off flag for momentum self transport
#   momViscosity    - On/Off flag for momentum mixing
#
 &PARM01
 tRef= 5.18, 5.13, 5.08, 5.04, 4.86,
       4.62, 4.4, 4.29, 4.25, 4.3,
       4.32, 4.35, 4.17, 4.0, 3.85,
 sRef= 34.13, 34.14, 34.14, 34.17, 34.2,
       34.23, 34.27, 34.31, 34.34, 34.4,
       34.44, 34.46, 34.47, 34.48, 34.48,
 no_slip_sides=.FALSE.,
 no_slip_bottom=.TRUE.,
 viscAz=8.33E-3,
 viscAh=6.51E4,
 diffKhT=39.0625,
 diffKzT=7.5E-5,
 diffKhS=39.0625,
 diffKzS=7.5E-5,
 beta=1.E-11,
 tAlpha=2.E-4,
 sBeta =7.4E-4,
 rigidLid=.FALSE.,
 implicitFreeSurface=.TRUE.,
 eosType='JMD95Z',
 ivdc_kappa=0.,
 implicitDiffusion=.TRUE.,
 useOldFreezing=.TRUE.,
 useRealFreshWaterFlux=.TRUE.,
 useCDscheme=.TRUE.,
 useNHMTerms=.TRUE.,
# turn on looped cells
 hFacMin=.05,
 hFacMindz=50.,
# set precision of data files
 readBinaryPrec=32,
 writeBinaryPrec=32,
 saltStepping=.TRUE.,
 tempStepping=.TRUE.,
 momStepping=.TRUE.,
# implicitDiffusion=.TRUE.,
 implicitViscosity=.TRUE.,
#globalFiles=.TRUE.,
#- not safe to use globalFiles in multi-processors runs; set instead useSingleCpuIO
 useSingleCpuIO=.TRUE.,
 allowFreezing=.FALSE.,
 HeatCapacity_Cp = 3986.D0
 gravity         = 9.8156,
 gBaro           = 9.8156,
 rhoNil          = 1027.D0
# useCDscheme=.TRUE.,
# useNHMTerms=.TRUE.,
 &

# Elliptic solver parameters
#
#   cg2dMaxIters       - Maximum number of 2d solver iterations
#   cg2dTargetResidual - Solver target residual
#
 &PARM02
 cg2dMaxIters=500,
 cg2dTargetResidual=1.D-8,
 &

# Time stepping parameters
#
#   startTime         - Integration starting time                (s)
#   endTime           - Integration ending time                  (s)
#   tauCD             - CD scheme coupling timescale             (s)
#   deltaTMom         - Timestep for momemtum equations          (s)
#   deltaTtracer      - Tracer timestep                          (s)
#   deltaTClock       - Timestep used as model "clock"           (s)
#   abEps             - Adams-Bashforth stabilising factor
#   pChkPtFreq        - Frequency of permanent check pointing    (s)
#   chkPtFreq         - Frequency of rolling check pointing      (s)
#   dumpFreq          - Frequency at which model state is stored (s)
#   tauThetaClimRelax - Relaxation to climatology time scale     (s)
#   tauSaltClimRelax  - Relaxation to climatology time scale     (s)
#
 &PARM03
 tauCD=172800.,
 startTime=0.0,
 endTime=864000.,
 deltaTmom=300.0,
 deltaTtracer=1800.0,
# deltaT=1200.0,
 deltaTClock =1800.0,
 cAdjFreq=0.,
 abEps=0.1,
 pChkptFreq=864000.,
 chkptFreq= 864000.,
 dumpFreq = 864000.,
 tavefreq = 864000.,
# 6 months restoring timescale for salinity
 tauSaltClimRelax = 15552000.0,
# 2 months restoring timescale for temperature
 tauThetaClimRelax = 5184000.0,
 monitorFreq=1.,
 periodicExternalForcing=.TRUE.,
 externForcingPeriod=2592000.,
 externForcingCycle=31104000.,
 forcing_In_AB = .FALSE.
 &

# Gridding parameters
#
#   usingSphericalPolarGrid - On/Off flag for spherical polar coordinates
#   usingCartesianGrid      - On/Off flag for selecting cartesian coordinates
#   delX                    - Zonal grid spacing         (degrees)
#   delY                    - Meridional grid spacing    (degrees)
#   delZ                    - Vertical grid spacing      (m)
#   phiMin                  - Southern boundary latitude (degrees)
#
 &PARM04
 usingCartesianGrid=.FALSE.,
 usingSphericalPolarGrid=.TRUE.,
 delX=356*1.E0,
 delY=30*1.E0,
 delZ=50., 70., 100., 140., 190.,
      240., 290., 340., 390., 440.,
      490., 540., 590., 640., 690., 
 phimin=-70.,
 thetamin=-178.,
 rSphere         = 6371.D3
 &

# Input datasets
#
#   bathyFile       - File containing bathymetry
#   hydrogThetaFile - File containing initial potential temperature data
#   hydrogSaltFile  - File containing initial salinity data
#   zonalWindFile   - File containing zonal wind data
#   meridWindFile   - File containing meridional wind data
#   thetaClimFile   - File containing theta climatology used for relaxation
#   saltClimFile    - File containing salt climatology used for relaxation
#
 &PARM05
 bathyFile       = 'bathymetry.bin',
 hydrogThetaFile = 'lev_monthly_ss_am_temp.bin',
 hydrogSaltFile  = 'lev_monthly_ss_am_salt.bin',
 zonalWindFile =   'ncep_taux.bin',
 meridWindFile =   'ncep_tauy.bin',
 thetaClimFile =   'lev_monthly_sstemp.bin', 
 saltClimFile =    'lev_monthly_sssalt.bin',
 surfQFile =       'ncep_qnet.bin',
## fresh water flux is turned off, uncomment next line to turn on
# (not recommend together with surface salinity restoring)
 EmPmRFile =      'ncep_emp.bin',
&


More information about the MITgcm-support mailing list