[Mitgcm-support] Re: about AIM

mitgcm-support at dev.mitgcm.org mitgcm-support at dev.mitgcm.org
Wed Jul 9 15:48:02 EDT 2003


Hi Jean-Michel,

I'm not sure why FWIND0 is set to 0.42 but this is the value I always used. 
I believe the value was originally changed by Christophe Herbaut while he 
was working on AIM (before me) and doing some debugging (see the comment 
line cchdbg).

Daniel

At 05:31 PM 9/22/2002 -0400, Jean-Michel Campin wrote:
>Hi,
>
>I have few questions about AIM:
>(Daniel, if you have some idea about the 1rst one ...Thanks)
>
>1) Does any one remember why we set FWIND0 to .42 instead of
>the value Franco use in his paper (=0.6)
>It's set in INPHYS, line 118:
>
>C
>C---  3. Constants for surface fluxes (common SFLCON):
>C
>       FWIND0 = 0.6 _d 0
>cchdbg *****************************************************
>       FWIND0 = 0.6 _d 0*0.7 _d 0
>cchdbg ******************************************************
>
>Since FWIND0 is use to compute the wind at the surface from
>the wind at the 1rst level, and used after in the stress,
>this can contribute to the strong surface wind we found in the
>runs.
>
>2) with full-cell and no Shapiro on Q, the global average of P-E
>balance the surface correction term. The 2 are different by less
>than 10^-4 mm/day.
>
>3) Concerning Temperature, it does not work:
>the same run with full-cell and no Shapiro on theta :
>net heat flux (top) - (surface) = 0.9 W/m2  => warm the atmosphere
>surface correction on Theta     = 0.5 W/m2  => warm the atmosphere
>latent heat due to exess P-E    = 0.3 W/m2  => warm the atmosphere
>
>My impression is that there is a problem in Heat conservation:
>the dynamics conserves theta but the physics conserves the absolute
>temperature:
>If there is a radiative flux Fr from level 1 --> 2,
>   Cp*dT1 = -Fr ; Cp*dT2 = +Fr   ==> T is conserved.
>but
>   PI_1*dTheta_1= -Fr ; PI_2*dTheta_2 = +Fr  ==> Theta is not conserved
>because PI_1 & PI_2 are differents (PI = Cp*(p/Po)**kappa)
>
>And since the dynamics advects theta, only theta can be conserved,
>its my impression.
>
>4) for Chris: to put partial cell in AIM, can I start with the code in
>the repository or would be better to use the coupled version ?
>
>Jean-Michel




More information about the MITgcm-support mailing list