[MITgcm-support] Problem using partial cells

[Isa Rosso]

Hi all,

I'm trying to use partial cells in my high-res model (1/80 deg resolution, open boundaries, 150 vertical layers, realistic topography). The model was running so far without this capability, with a time step up to 30 seconds. I tried to run it setting the hFacMin = 0.1, reducing the time step to 5 seconds, but the model crashes. 

Is there some particular settings I should consider?

Below the PARAM01 of the data file.


Thanks,
isa


&PARM01
 tRef               	= 150*5.0,
 sRef               = 150*34.5,
 viscAr 		= 5.6614e-04,
 diffKrT 		= 1.e-5,
 diffKrS 		= 1.e-5,
 no_slip_sides 	= .FALSE.,
 no_slip_bottom 	= .TRUE.,
 rhonil 		= 1027.5,
 rhoConstFresh  	= 999.8,
 eosType 	= 'JMD95Z',
 hFacMin 	= 0.1,
 hFacInf 	 	= 0.01,
 hFacSup 	= 2.,
 select_rStar 	= 2,
 nonlinFreeSurf 	= 4,
 implicitDiffusion 	= .TRUE.,
 implicitViscosity 	= .TRUE.,
 viscC4Leith 	 	= 2.0,
 viscC4Leithd	= 2.0,
 viscC2Leith 		= 1.5,
 viscC2Leithd 	= 1.5,
 viscA4GridMax 	= 0.5,
 viscAhGridMax 	= 0.5,
 viscC4Smag 	= 2.2,
 diffK4T 	= 1.E8,
 diffK4S 	= 1.E8,
 useAreaViscLength = .TRUE.,
 sideDragFactor 	= 0.,
 useAbsVorticity 	= .TRUE.,
 selectVortScheme 	= 2,
 selectKEscheme 	= 3,
 bottomDragQuadratic = 0.0034,
 tempAdvScheme 	= 7,
 saltAdvScheme 		= 7,
 StaggerTimeStep 	= .TRUE.,
 multiDimAdvection 	= .TRUE.,
 vectorInvariantMomentum = .TRUE.,
 implicitFreeSurface 	= .TRUE.,
 exactConserv 	= .TRUE.,
 debuglevel 		= -1,
 convertFW2Salt 	= -1
 useRealFreshWaterFlux = .TRUE.,
 useSingleCPUio 	= .TRUE.,
 globalFiles 		= .TRUE.,
 allowFreezing 	= .TRUE.,
 &

***************
Isa Rosso
PhD student

Research School of Earth Sciences
The Australian National University
Canberra ACT 0200
Australia

ph: +61 2 6125 9958

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20120903/eefe22dd/attachment-0001.htm>