[MITgcm-devel] problem with SEAICEadvSnow

[Martin Losch]

Hi Dimitris,

Do you actually use the 10m cap for heff? This is what I find in the  
current code (and this part has not change for a long time):
> #ifdef DO_WE_NEED_THIS
> c          HEFF(I,J,1,bi,bj)=MIN(MAX_HEFF,HEFF(I,J,1,bi,bj))
> #endif

Here are the parameters that I normally use (in data.seaice):
   SEAICEuseFlooding = .true., ! you may not want to use that, but I  
find it does not have much of an impact, if snow advection is turned  
on by the next line:
   SEAICEadvSnow     = .true.,
   SEAICEadvScheme   = 1, ! this is not the default. The default is 2  
(2nd order central differences with modified euler time step, the  
original code)
   DIFF1 = 0., ! some number that gets non-dimensionalized by the  
grid spacing. DIFF1 is the same for all scalars (HEFF,AREA,HSNOW)

On 6 Feb 2007, at 14:49, Dimitris Menemenlis wrote:

> Martin, welcome back.  Where have you been?  We have missed you!
>
>> I cannot remember having any problems with my 2deg model
>
> Ice and snow seem OK most everywhere.  Only a problem is one region  
> right along
> Antarctic coastline, maybe because of strong katabatic winds?
>
>> However, I normally use advection scheme 1 (1st order upwind) and  
>> not diffusion (DIFF1 = 0.)
>
> Is advection scheme 1 the default?  If not, how would I set it?  Is  
> there a particular advection scheme that would be more stable than  
> the default that I can try?
>
> Do you set DIFF1 to zero both for ice and snow?  Do you do that by  
> commenting lines in the code or is there parameters to set in  
> data.seaice?
>
>> Why should the snow get unstable and the ice not, they usually  
>> have similar
>> distributions (in my case).
>
> I will take closer look at ice too and get back if I see  
> something.  Maybe ice does not blow up because of 10-m cap?  The  
> integration completes normally for "SEAICEadvSnow=.false."
>
> On a different topic, please talk with JM and Chris and let me know  
> how you would like me to respond to Elizabeth.
>
> D.
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel