No subject


Thu Oct 30 11:44:21 EDT 2014


&gt;&gt;&gt; &gt; what I understand is that<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; tend_heat_surf =3D=C2=A0 tend_advection +<br>
&gt;&gt;&gt; &gt;=C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =
=C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 tend_diffusion +<br>
&gt;&gt;&gt; &gt;=C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =
=C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 tend_kpp +<br>
&gt;&gt;&gt; &gt;=C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =
=C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 (TFLUX-oceQsw)/(rhoConst*Cp*DRF(1=
)*hFacC(i,j,1))<br>
&gt;&gt;&gt; &gt; +<br>
&gt;&gt;&gt; &gt; =E2=80=8B<br>
&gt;&gt;&gt; &gt;=C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =
=C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=
=A0 =C2=A0tend_Qsw -<br>
&gt;&gt;&gt; &gt;=C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =
=C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0(WTHMASS(i,j,1) -<br=
>
&gt;&gt;&gt; &gt; TsurfCor)/(DRF(1)*hFacC(i,j,1))<br>
&gt;&gt;&gt; &gt; where<br>
&gt;&gt;&gt; &gt; tend_Qsw =3D oceQsw/(rhoConst*Cp)/(DRF(1)*hFacC(i,j,1))*(=
swfrac(1)-swfrac(2))<br>
&gt;&gt;&gt; &gt; and<br>
&gt;&gt;&gt; &gt; TsurfCor =3D SUM( WTHMASS(i,j,1)*RAC(i,j) ) /=C2=A0 globa=
lArea<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; and &quot;swfrac&quot; is already defined in our earlier =
communications.<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; =E2=80=8BAm I correct?<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; Thanking in advance,<br>
&gt;&gt;&gt; &gt; Abhisek=E2=80=8B<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; *--------------------------------------------------------=
------------------------ABHISEK<br>
&gt;&gt;&gt; &gt; CHAKRABORTY*<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; *Scientist - &quot;SD&quot;Oceanic Sciences Division (*<b=
r>
&gt;&gt;&gt; &gt; *AOSG/EPSA)*<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; *Space Applications Centre*<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; *Indian Space Research Organization*<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; *Ahmedabad - 380 015*<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; *Gujarat, INDIA*<br>
&gt;&gt;&gt; &gt; *Contact: +91-79-2691-6054 (O), +91-79-2686-1929(R)*<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; On Tue, Nov 25, 2014 at 2:21 AM, Jean-Michel Campin &lt;<=
a href=3D"mailto:jmc at ocean.mit.edu">jmc at ocean.mit.edu</a>&gt;<br>
&gt;&gt;&gt; &gt; wrote:<br>
&gt;&gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; Hi Abhisek,<br>
&gt;&gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; Following your suggestions, I think I am now ab=
le to close the heat<br>
&gt;&gt;&gt; &gt; &gt; budget<br>
&gt;&gt;&gt; &gt; &gt; &gt; for the subsurface levels.<br>
&gt;&gt;&gt; &gt; &gt; This is good news.<br>
&gt;&gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; I have attached two figures (one for subsurface=
 at 20 m depth and another<br>
&gt;&gt;&gt; &gt; &gt; &gt; for surface).<br>
&gt;&gt;&gt; &gt; &gt; Could you clarify what is represented on these 2 fig=
s: it is for a<br>
&gt;&gt;&gt; &gt; &gt; single grid point (but where ?) or the horizontal me=
an ?<br>
&gt;&gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; Regarding the surface forcing, which diagnostics are=
 you using:<br>
&gt;&gt;&gt; &gt; &gt;=C2=A0 surForcT ? TFLUX ? or oceQnet and TRELAX ?<br>
&gt;&gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; If you are using TFLUX (which contains the verticall=
y integrated SW<br>
&gt;&gt;&gt; &gt; &gt; heating)<br>
&gt;&gt;&gt; &gt; &gt; for the surface budget, you may want to remove oceQs=
w and then<br>
&gt;&gt;&gt; &gt; &gt; add back tend_Qsw(k=3D1) (as you wrote below).<br>
&gt;&gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; Cheers,<br>
&gt;&gt;&gt; &gt; &gt; Jean-Michel<br>
&gt;&gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; On Fri, Nov 21, 2014 at 07:13:31PM +0530, Abhisek Ch=
akraborty wrote:<br>
&gt;&gt;&gt; &gt; &gt; &gt; Dear Jean-Michel,<br>
&gt;&gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; To close the surface heat budget, I followed yo=
ur suggestion:<br>
&gt;&gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; =E2=80=8B=E2=80=8B<br>
&gt;&gt;&gt; &gt; &gt; TsurfCor =3D SUM( WTHMASS(i,j,1)*RAC(i,j) ) /=C2=A0 =
globalArea<br>
&gt;&gt;&gt; &gt; &gt; &gt;=C2=A0 and finally:<br>
&gt;&gt;&gt; &gt; &gt; &gt;=C2=A0 tend_T =3D (...) - ( WTHMASS(i,j,1) - Tsu=
rfCor )/(drF(1)*hFacC(i,j,1))<br>
&gt;&gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; Also, the shortwave tendency terms is computed =
as:<br>
&gt;&gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; tend_Qsw =3D oceQsw/(rhoConst*Cp)/(drF(k)*hFacC=
(I,j,k))<br>
&gt;&gt;&gt; &gt; &gt; &gt;=C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0=
* ( swfrac(k) - swfrac(k+1) )<br>
&gt;&gt;&gt; &gt; &gt; &gt; where for surface, k=3D1<br>
&gt;&gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; But in this way, the surface tendency term is b=
ecoming one order higher<br>
&gt;&gt;&gt; &gt; &gt; in<br>
&gt;&gt;&gt; &gt; &gt; &gt; magnitude than other terms (advection, diffusio=
n, surface correction.<br>
&gt;&gt;&gt; &gt; &gt; Pls.<br>
&gt;&gt;&gt; &gt; &gt; &gt; see the figure in my earlier mail).<br>
&gt;&gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; So, please guide me how to close the budget at =
surface.<br>
&gt;&gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; regards,<br>
&gt;&gt;&gt; &gt; &gt; &gt; Abhisek<br>
&gt;&gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; *---------------------------------------------------=
-----------------------------ABHISEK<br>
&gt;&gt;&gt; &gt; &gt; &gt; CHAKRABORTY*<br>
&gt;&gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; *Scientist - &quot;SD&quot;Oceanic Sciences Div=
ision (*<br>
&gt;&gt;&gt; &gt; &gt; &gt; *AOSG/EPSA)*<br>
&gt;&gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; *Space Applications Centre*<br>
&gt;&gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; *Indian Space Research Organization*<br>
&gt;&gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; *Ahmedabad - 380 015*<br>
&gt;&gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; *Gujarat, INDIA*<br>
&gt;&gt;&gt; &gt; &gt; &gt; *Contact: +91-79-2691-6054 (O), +91-79-2686-192=
9(R)*<br>
&gt;&gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; On Thu, Nov 20, 2014 at 6:59 PM, Abhisek Chakra=
borty &lt;<br>
&gt;&gt;&gt; &gt; &gt; <a href=3D"mailto:abhisek.sac at gmail.com">abhisek.sac=
@gmail.com</a>&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; wrote:<br>
&gt;&gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt; Hi Jean-Michel,<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt; Sorry for delayed reply.<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt; Following your suggestions, I think I am n=
ow able to close the heat<br>
&gt;&gt;&gt; &gt; &gt; budget<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt; for the subsurface levels.<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt; However, in spite of my different attempts=
 I am not able to close the<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt; budget at surface where I am finding the s=
urface correction term (as<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt; specified in your mail) along with other a=
dvection-diffusion terms are<br>
&gt;&gt;&gt; &gt; &gt; too<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt; small to balance the shortwave heating ter=
m.<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt; I have attached two figures (one for subsu=
rface at 20 m depth and<br>
&gt;&gt;&gt; &gt; &gt; another<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt; for surface).<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt; Pls. let me know if I am missing something=
.<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt; Thanking in advance,<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt; Abhisek<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; *---------------------------------------------------=
-----------------------------ABHISEK<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt; CHAKRABORTY*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt; *Scientist - &quot;SD&quot;Oceanic Science=
s Division (*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt; *AOSG/EPSA)*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt; *Space Applications Centre*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt; *Indian Space Research Organization*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt; *Ahmedabad - 380 015*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt; *Gujarat, INDIA*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt; *Contact: +91-79-2691-6054 (O), +91-79-268=
6-1929(R)*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt; On Wed, Nov 19, 2014 at 1:19 AM, Jean-Mich=
el Campin &lt;<a href=3D"mailto:jmc at ocean.mit.edu">jmc at ocean.mit.edu</a><br=
>
&gt;&gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt; wrote:<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; Hi Abhisek,<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; This budget looks better. Just few rem=
arks (see below).<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; On Tue, Nov 18, 2014 at 11:20:46PM +05=
30, Abhisek Chakraborty wrote:<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; Dear Jean-Michel,<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; Thanks for your prompt suggestion=
.<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; (1) the advection term for k&gt;1=
 will be<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;=C2=A0 =C2=A0 - [ (ADVx_TH(i+1,j,k=
) =E2=80=93 ADVx_TH(I,j,k))/CV +<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;=C2=A0 =C2=A0 =C2=A0(ADVy_TH(I,j+1=
,k) =E2=80=93 ADVy_TH(I,j,k))/CV +<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;=C2=A0 =C2=A0 (ADVr_TH(I,j,k+1) =
=E2=80=93 ADVr_TH(I,j,k))/CV ]<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;=C2=A0 The vertical advective &amp; dif=
fusive transport/flux are (I think)<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;=C2=A0 positive in the direction of the=
 vertical coordinate, i.e.,<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;=C2=A0 when using z-coords, positive up=
ward. So you will have to reverse<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;=C2=A0 (or add a minus sign) in front o=
f:<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;=C2=A0 =C2=A0 =C2=A0 (ADVr_TH(I,j,k+1) =
=E2=80=93 ADVr_TH(I,j,k))/CV<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;=C2=A0 Where CV=3DAc * del RF * Hc=
 =3D RAC(I,j) * (RF(K+1)-RF(K))=C2=A0 *<br>
&gt;&gt;&gt; &gt; &gt; hFacC(I,j,k)<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;=C2=A0 (<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; <a href=3D"http://mitgcm.org/public/r2_manual/latest=
/online_documents/node73.html">http://mitgcm.org/public/r2_manual/latest/on=
line_documents/node73.html</a>)<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; And similarly the diffusion and K=
PP terms. Correct?<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; (2) The source term (shortwave) f=
or k&gt;1 will be the following<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;=C2=A0 =C2=A0 depth=3DRF(k)<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;=C2=A0 =C2=A0 swfrac=3D0.62*exp(de=
pth/0.6) + (1.0-0.62) * exp(depth/20)<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;=C2=A0 =C2=A0if(depth &lt; -200.0)=
 swfrac=3D0.0<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;=C2=A0 =C2=A0source term =3D=C2=A0=
 =C2=A0(oceQsw/(rhoConst*Cp)) * swfrac * hFacC(I,j,k)<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;=C2=A0 =C2=A0Is this correct for k=
&gt;1 ?<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;=C2=A0 This does not look right:<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;=C2=A0 =C2=A0 from depth=3DRF(k)=C2=A0 =
, we compute swfrac(k) as above<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;=C2=A0 =C2=A0 from depth=3DRF(k+1), we =
compute swfrac(k+1) as above except that<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;=C2=A0 =C2=A0 =C2=A0at the bottom (deep=
est level of this column) swfrac(k+1)=3D0.<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;=C2=A0 and the tendency contribution fr=
om shortwave heating should be<br>
&gt;&gt;&gt; &gt; &gt; something<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; like:<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;=C2=A0 =C2=A0 tend_Qsw =3D oceQsw/(rhoC=
onst*Cp)/(drF(k)*hFacC(I,j,k))<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;=C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=
=A0 =C2=A0* ( swfrac(k) - swfrac(k+1) )<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; (3) I am using the tendency as TO=
TTTEND/86400. Is it correct?<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;=C2=A0 yes, providing you are using lin=
ear free-surface.<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; (4) I am using linFSConserveTr=3D=
TRUE. Yet I have to use<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; WTHMASS(i,j,1)*RAC(i,j) to correc=
t the heat budget at the surface?<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;=C2=A0 Yes, but since you are using lin=
FSConserveTr=3DTRUE., the global mean<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;=C2=A0 =C2=A0value of WTHMASS has to be=
 removed in the budget (like it is in the<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; model):<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;=C2=A0 =C2=A0 TsurfCor =3D SUM( WTHMASS=
(i,j,1)*RAC(i,j) ) /=C2=A0 globalArea<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;=C2=A0 and finally:<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;=C2=A0 tend_T =3D (...) - ( WTHMASS(i,j=
,1) - TsurfCor )/(drF(1)*hFacC(i,j,1))<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; Cheers,<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; Jean-Michel<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; Thanking in advance,<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; Abhisek<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; *---------------------------------------------------=
-----------------------------ABHISEK<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; CHAKRABORTY*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; *Scientist - &quot;SD&quot;Oceani=
c Sciences Division (*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; *AOSG/EPSA)*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; *Space Applications Centre*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; *Indian Space Research Organizati=
on*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; *Ahmedabad - 380 015*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; *Gujarat, INDIA*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; *Contact: +91-79-2691-6054 (O), +=
91-79-2686-1929(R)*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; On Tue, Nov 18, 2014 at 10:42 PM,=
 Jean-Michel Campin &lt;<br>
&gt;&gt;&gt; &gt; &gt; <a href=3D"mailto:jmc at ocean.mit.edu">jmc at ocean.mit.e=
du</a><br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; wrote:<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; Hi Abhisek,<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; I would recommand to concent=
rate first on any level except the<br>
&gt;&gt;&gt; &gt; &gt; surface<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; (k=3D1),<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; and try to close the budget =
; then, once you get a clean closed<br>
&gt;&gt;&gt; &gt; &gt; budget<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; for levels at depth, you can=
 revisit the surface level problem.<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; at depth (k&gt;1):<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; the way you estimate the ten=
dency is not what is done in the<br>
&gt;&gt;&gt; &gt; &gt; model,<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; see e.g.:<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; <a href=3D"http://mitgcm.org/public/r2_manual/latest=
/online_documents/node73.html">http://mitgcm.org/public/r2_manual/latest/on=
line_documents/node73.html</a><br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; The diagnostics for the KPP =
non-local term is like a vertical flux<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; (added to fVerT), so it can =
be treated (for budget purpose) the<br>
&gt;&gt;&gt; &gt; &gt; same<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; way<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; as DIFrI_TH.<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; at the surface, you can sear=
ch the mitgcm-support archive and find<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; some<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; information there, e.g.:<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; <a href=3D"http://mitgcm.org=
/pipermail/mitgcm-support/2014-April/009093.html">http://mitgcm.org/piperma=
il/mitgcm-support/2014-April/009093.html</a><br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; but it&#39;s likely that you=
 would need additional diagnostics<br>
&gt;&gt;&gt; &gt; &gt; output,i<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; e.g. WTHMASS &amp; WSLTMASS =
at k=3D1.<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; Cheers,<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; Jean-Michel<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; On Tue, Nov 18, 2014 at 09:4=
9:34PM +0530, Abhisek Chakraborty<br>
&gt;&gt;&gt; &gt; &gt; wrote:<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; Many thanks Ryan, Dimit=
ris, John and Gael for your suggestions.<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; I am not much familiar =
with MATLAB, so for me it was too much<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; difficult<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; to<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; understand the code sha=
red by Gael, though I tried to look into<br>
&gt;&gt;&gt; &gt; &gt; it.<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; As a matter of fact is =
that I am using is the ECCO1 version (1<br>
&gt;&gt;&gt; &gt; &gt; deg<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; global,<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; 50 levels, excluding po=
les; i=3D1:360; j=3D1:160; k=3D1:50) of MITGCM<br>
&gt;&gt;&gt; &gt; &gt; with<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; linear<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; free surface and linFSC=
onserveTR=3DTrue. The diagnostic<br>
&gt;&gt;&gt; &gt; &gt; parameters are<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; written in every day wi=
th diagnostic timephase=3D0. Model<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; parameterizations<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; are KPP and GM-Redi.<br=
>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; I=E2=80=99d like to ana=
lyze the heat budget at a particular location,<br>
&gt;&gt;&gt; &gt; &gt; say,<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; (i,j,k).<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; So, from the suggestion=
s by Ryan, Dimitris and looking into the<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; Gael=E2=80=99s<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; code, I am trying the f=
ollowing:<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; !!Cell volumes (center =
difference)<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; Cvx =3D (RAC(i+1,j) - R=
AC(i-1,j)) * RF(k) * hFacC(i,j,k)<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; Cvy =3D (RAC(i,j+1) - R=
AC(i,j-1)) * RF(k) * hFacC(i,j,k)<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; Cvz =3D RAC(i,j) * (RF(=
k+1) =E2=80=93 RF(k)) * hFacC(i,j,k)<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; !!Divergence of advecti=
on &amp; diffusion by center difference<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; if (k&gt;1) then<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;=C2=A0 =C2=A0 =C2=A0 =C2=
=A0 =C2=A0 adc_dif_vert_k =3D (ADVr_TH(i,j,k+1) =E2=80=93<br>
&gt;&gt;&gt; &gt; &gt; ADVr_TH(I,j,k-1))/cvz<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; +<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;=C2=A0 =C2=A0 =C2=A0 =C2=
=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =
=C2=A0 =C2=A0(DIFrE_TH(I,j,k+1) =E2=80=93<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; DIFrE_TH(I,j,k-1))/cvz +<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;=C2=A0 =C2=A0 =C2=A0 =C2=
=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =
=C2=A0 =C2=A0(DIFrI_TH(I,j,k+1) =E2=80=93<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; DIFrI_TH(I,j,k-1))/cvz<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; else<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;=C2=A0 =C2=A0 =C2=A0 =C2=
=A0 =C2=A0 adc_dif_vert_k =3D 0.0<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; endif<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; adc_dif_hori_k =3D (ADV=
x_TH(i+1,j,k) =E2=80=93 ADVx_TH(i-1,j,k))/cvx +<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;=C2=A0 =C2=A0 =C2=A0 =C2=
=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 (ADVy_TH=
(i,j+1,k) =E2=80=93<br>
&gt;&gt;&gt; &gt; &gt; ADVy_TH(i,j-1,k))/cvy +<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;=C2=A0 =C2=A0 =C2=A0 =C2=
=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 (DIFxE_T=
H(i+1,j,k) =E2=80=93<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; DIFxE_TH(i-1,j,k))/cvx +<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;=C2=A0 =C2=A0 =C2=A0 =C2=
=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 (DIFyE_T=
H(i,j+1,k) =E2=80=93<br>
&gt;&gt;&gt; &gt; &gt; DIFyE_TH(i,j-1,k))/cvy<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; !!To account for the=C2=
=A0 shortwave heating<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; depth=3DRF(k)<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; swfrac=3D0.62*exp(depth=
/0.6) + (1.0-0.62) * exp(depth/20)<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; if(depth &lt; -200.0) t=
hen<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; swfrac=3D0.0<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; endif<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; if (k=3D1)then<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;=C2=A0 =C2=A0 =C2=A0 =C2=
=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 swfrac=3D1.0<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; endif<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; !!Finally balance<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; TOTTTEND/86400.0 =3D - =
adc_dif_vert_k -=C2=A0 adc_dif_hori_k +<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; (oceQsw/(rhoConst*Cp)) =
* swfrac * hFacC(I,j,k) + KPPg_TH(i,j,k)<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; Please let me know whet=
her I am doing properly or not.<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; Thanking in advance,<br=
>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; Abhisek<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; *---------------------------------------------------=
-----------------------------ABHISEK<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; CHAKRABORTY*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; *Scientist - &quot;SD&q=
uot;Oceanic Sciences Division (*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; *AOSG/EPSA)*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; *Space Applications Cen=
tre*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; *Indian Space Research =
Organization*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; *Ahmedabad - 380 015*<b=
r>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; *Gujarat, INDIA*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; *Contact: +91-79-2691-6=
054 (O), +91-79-2686-1929(R)*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; On Tue, Nov 18, 2014 at=
 5:06 AM, gael forget &lt;<a href=3D"mailto:gforget at mit.edu">gforget at mit.ed=
u</a>&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; wrote:<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; Hi Abhisek, et al,=
<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; the following matl=
ab code may provide some guidance :<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; <a href=3D"http://mitgcm.org/viewvc/MITgcm/MITgcm_co=
ntrib/gael/matlab_class/gcmfaces_diags/diags_set_D.m">http://mitgcm.org/vie=
wvc/MITgcm/MITgcm_contrib/gael/matlab_class/gcmfaces_diags/diags_set_D.m</a=
><br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; It does not cover =
every possible combination of model options<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; (e.g., no kpp or t=
hsice) of course, but treats a few. In that<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; code,<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; it is assumed that=
<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; 1) the setting for=
 the diagnostics package was according to<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; <a href=3D"http://mitgcm.org/viewvc/MITgcm/MITgcm_co=
ntrib/gael/verification/global_oce_llc90/input/data.diagnostics">http://mit=
gcm.org/viewvc/MITgcm/MITgcm_contrib/gael/verification/global_oce_llc90/inp=
ut/data.diagnostics</a><br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;=C2=A0 =C2=A0 =C2=
=A0 which includes vertical integration. The relevant<br>
&gt;&gt;&gt; &gt; &gt; diagnostics<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;=C2=A0 =C2=A0 =C2=
=A0 for budget are those with filename like<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &#39;diags/budg2d_snap_set1&#39;,<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; 2) tendencies have=
 been pre-computed by differencing between<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;=C2=A0 =C2=A0 =C2=
=A0consecutive (monthly) snap shots =E2=80=94 so that ETAN<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;=C2=A0 =C2=A0 =C2=
=A0really means dETAN/dt in this context, and similarly for<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;=C2=A0 =C2=A0 =C2=
=A0THETA, SALT, etc. The other variable names used in<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;=C2=A0 =C2=A0 =C2=
=A0diags_set_D.m<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; &lt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; <a href=3D"http://mitgcm.org/viewvc/MITgcm/MITgcm_co=
ntrib/gael/matlab_class/gcmfaces_diags/diags_set_D.m">http://mitgcm.org/vie=
wvc/MITgcm/MITgcm_contrib/gael/matlab_class/gcmfaces_diags/diags_set_D.m</a=
><br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; match<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; those used in pkg/=
diagnostics<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; I should also ment=
ion that diags_set_D.m<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; &lt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; <a href=3D"http://mitgcm.org/viewvc/MITgcm/MITgcm_co=
ntrib/gael/matlab_class/gcmfaces_diags/diags_set_D.m">http://mitgcm.org/vie=
wvc/MITgcm/MITgcm_contrib/gael/matlab_class/gcmfaces_diags/diags_set_D.m</a=
><br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; is<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; part of a<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; broader framework =
(gcmfaces) which is documented in<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; <a href=3D"http://mitgcm.org/viewvc/*checkout*/MITgc=
m/MITgcm_contrib/gael/matlab_class/gcmfaces.pd">http://mitgcm.org/viewvc/*c=
heckout*/MITgcm/MITgcm_contrib/gael/matlab_class/gcmfaces.pd</a><br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; &lt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; <a href=3D"http://mitgcm.org/viewvc/MITgcm/MITgcm_co=
ntrib/gael/matlab_class/gcmfaces.pd">http://mitgcm.org/viewvc/MITgcm/MITgcm=
_contrib/gael/matlab_class/gcmfaces.pd</a><br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; f<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; =E2=80=A6 although=
 I am afraid this documentation is also<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; lacking wrt the bu=
dget diagnostic codes.<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; Cheers,<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; Gael<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; On Nov 17, 2014, a=
t 5:53 PM, Dimitris Menemenlis &lt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; <a href=3D"mailto:dmenemenli=
s at gmail.com">dmenemenlis at gmail.com</a>&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; wrote:<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; =E2=80=A6 and you =
probably already know this, but not clear from your<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; email:<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; shortwave is not d=
umped in surface level only, but rather is<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; distributed<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; with exponentially=
<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; decaying profile i=
n top 200-m, with subtle modifications if<br>
&gt;&gt;&gt; &gt; &gt; this<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; heating<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; is occurring withi=
n<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; the KPP mixing lay=
er depth.<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; On Nov 17, 2014, a=
t 2:44 PM, Ryan Abernathey &lt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; <a href=3D"mailto:ryan.abern=
athey at gmail.com">ryan.abernathey at gmail.com</a>&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; wrote:<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; Abhisek,<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; This is a very com=
mon question that unfortunately is not<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; documented<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; very<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; well. A couple of =
suggestions...<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; For the transport =
terms, (ADV* and DF*), keep in mind that<br>
&gt;&gt;&gt; &gt; &gt; these<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; are<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; FLUX<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; terms, defined on =
the cell boundaries. In order to calculate a<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; heat<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; budget,<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; you need to calcul=
ate the DIVERGENCE of those fluxes.<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; <a href=3D"http://mitgcm.org/public/r2_manual/latest=
/online_documents/node71.html">http://mitgcm.org/public/r2_manual/latest/on=
line_documents/node71.html</a><br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; Maybe this is obvi=
ous, but your email did not make it clear<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; whether you<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; were doing this or=
 not.<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; Likewise, TFLUX is=
 the downward flux at the surface, in<br>
&gt;&gt;&gt; &gt; &gt; W/m^2. In<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; order to<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; convert to tempera=
ture units, you need to divide by<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; HeatCapacity_Cp*rU=
nit2mass. For z coordinates, rUnit2mass is<br>
&gt;&gt;&gt; &gt; &gt; equal<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; to rhoConst.<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; There are some sub=
tleties regarding tracer budgets under<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; different free<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; surface treatments=
. Maybe someone else on the list understands<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; this<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; issue<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; better and is will=
ing to explain. Anyway, those should be<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; second-order<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; effects.<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; Good luck with you=
r analysis.<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; -Ryan<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; On Sat, Nov 15, 20=
14 at 2:29 AM, Abhisek Chakraborty &lt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; <a href=3D"mailto:=
abhisek.sac at gmail.com">abhisek.sac at gmail.com</a>&gt; wrote:<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; Dear Users,<br=
>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; I am trying to=
 analyze heat budget in MITGCM. From the<br>
&gt;&gt;&gt; &gt; &gt; diagnostic<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; outputs I have=
 the following terms<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; TOTTTEND , ADV=
r_TH, ADVx_TH, ADVy_TH, DFrE_TH, DFxE_TH,<br>
&gt;&gt;&gt; &gt; &gt; DFyE_TH,<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; DFrI_TH, KPPg_=
TH, TFLUX (i.e. all terms corresponding to<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; tendency,<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; advection, dif=
fusion, KPP and total heat flux).<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; Apart from the=
se, I have other usual outputs like THETA, SALT<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; etc.<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; The tendency t=
erm is in degC/day, the other terms are in<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; degC*m^3/s<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; and TFLUX is i=
n W/m^2. Thus I have converted the tendency<br>
&gt;&gt;&gt; &gt; &gt; term<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; into<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; degC/s by divi=
ding by 86400.0 and the advection-diffusion<br>
&gt;&gt;&gt; &gt; &gt; terms<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; are<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; divided by cel=
l volume (=3DRAC*RF*HFACC).<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; I am using &qu=
ot;linFSconserveTr=3DTrue&quot;.<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; For the surfac=
e layer I have to consider TFLUX (watt/m^2),<br>
&gt;&gt;&gt; &gt; &gt; but<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; how to<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; convert it int=
o degC/s unit?<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; For subsurface=
 levels, I am trying to equate tendency term<br>
&gt;&gt;&gt; &gt; &gt; with<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; (advection + d=
iffusion + KPP term). But there is a mismatch.<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; Can somebody p=
lease guide me how to achieve the exact heat<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; budget for<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; both surface a=
nd subsurface levels ?<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; Thanking in ad=
vance,<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; Abhisek<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; --<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; *ABHISEK CHAKR=
ABORTY*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; *Scientist - &=
quot;SD&quot;Oceanic Sciences Division (*AOSG/EPSA)*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; *Space Applica=
tions Centre*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; *Indian Space =
Research Organization*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; *Ahmedabad - 3=
80 015*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; *Gujarat, INDI=
A*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt; *Contact: +91-=
79-2691-6054 (O), +91-79-2686-1929(R)*<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; __________________=
_____________________________<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; MITgcm-support mai=
ling list<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; <a href=3D"mailto:=
MITgcm-support at mitgcm.org">MITgcm-support at mitgcm.org</a><br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; <a href=3D"http://=
mitgcm.org/mailman/listinfo/mitgcm-support">http://mitgcm.org/mailman/listi=
nfo/mitgcm-support</a><br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; __________________=
_____________________________<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; MITgcm-support mai=
ling list<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; <a href=3D"mailto:=
MITgcm-support at mitgcm.org">MITgcm-support at mitgcm.org</a><br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt; <a href=3D"http://=
mitgcm.org/mailman/listinfo/mitgcm-support">http://mitgcm.org/mailman/listi=
nfo/mitgcm-support</a><br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; _______________________=
________________________<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; MITgcm-support mailing =
list<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; <a href=3D"mailto:MITgc=
m-support at mitgcm.org">MITgcm-support at mitgcm.org</a><br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; &gt; <a href=3D"http://mitgc=
m.org/mailman/listinfo/mitgcm-support">http://mitgcm.org/mailman/listinfo/m=
itgcm-support</a><br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; ____________________________=
___________________<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; MITgcm-support mailing list<=
br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; <a href=3D"mailto:MITgcm-sup=
port at mitgcm.org">MITgcm-support at mitgcm.org</a><br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt; <a href=3D"http://mitgcm.org=
/mailman/listinfo/mitgcm-support">http://mitgcm.org/mailman/listinfo/mitgcm=
-support</a><br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; _________________________________=
______________<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; MITgcm-support mailing list<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; <a href=3D"mailto:MITgcm-support@=
mitgcm.org">MITgcm-support at mitgcm.org</a><br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; &gt; <a href=3D"http://mitgcm.org/mail=
man/listinfo/mitgcm-support">http://mitgcm.org/mailman/listinfo/mitgcm-supp=
ort</a><br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; ______________________________________=
_________<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; MITgcm-support mailing list<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; <a href=3D"mailto:MITgcm-support at mitgc=
m.org">MITgcm-support at mitgcm.org</a><br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt; <a href=3D"http://mitgcm.org/mailman/l=
istinfo/mitgcm-support">http://mitgcm.org/mailman/listinfo/mitgcm-support</=
a><br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;&gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; &gt; _______________________________________________=
<br>
&gt;&gt;&gt; &gt; &gt; &gt; MITgcm-support mailing list<br>
&gt;&gt;&gt; &gt; &gt; &gt; <a href=3D"mailto:MITgcm-support at mitgcm.org">MI=
Tgcm-support at mitgcm.org</a><br>
&gt;&gt;&gt; &gt; &gt; &gt; <a href=3D"http://mitgcm.org/mailman/listinfo/m=
itgcm-support">http://mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
&gt;&gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt; &gt; &gt; _______________________________________________<br>
&gt;&gt;&gt; &gt; &gt; MITgcm-support mailing list<br>
&gt;&gt;&gt; &gt; &gt; <a href=3D"mailto:MITgcm-support at mitgcm.org">MITgcm-=
support at mitgcm.org</a><br>
&gt;&gt;&gt; &gt; &gt; <a href=3D"http://mitgcm.org/mailman/listinfo/mitgcm=
-support">http://mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
&gt;&gt;&gt; &gt; &gt;<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; &gt; _______________________________________________<br>
&gt;&gt;&gt; &gt; MITgcm-support mailing list<br>
&gt;&gt;&gt; &gt; <a href=3D"mailto:MITgcm-support at mitgcm.org">MITgcm-suppo=
rt at mitgcm.org</a><br>
&gt;&gt;&gt; &gt; <a href=3D"http://mitgcm.org/mailman/listinfo/mitgcm-supp=
ort">http://mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; _______________________________________________<br>
&gt;&gt;&gt; MITgcm-support mailing list<br>
&gt;&gt;&gt; <a href=3D"mailto:MITgcm-support at mitgcm.org">MITgcm-support at mi=
tgcm.org</a><br>
&gt;&gt;&gt; <a href=3D"http://mitgcm.org/mailman/listinfo/mitgcm-support">=
http://mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; _______________________________________________<br>
&gt;&gt; MITgcm-support mailing list<br>
&gt;&gt; <a href=3D"mailto:MITgcm-support at mitgcm.org">MITgcm-support at mitgcm=
.org</a><br>
&gt;&gt; <a href=3D"http://mitgcm.org/mailman/listinfo/mitgcm-support">http=
://mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
&gt;&gt;<br>
&gt;<br>
&gt;<br>
&gt; _______________________________________________<br>
&gt; MITgcm-support mailing list<br>
&gt; <a href=3D"mailto:MITgcm-support at mitgcm.org">MITgcm-support at mitgcm.org=
</a><br>
&gt; <a href=3D"http://mitgcm.org/mailman/listinfo/mitgcm-support">http://m=
itgcm.org/mailman/listinfo/mitgcm-support</a><br>
&gt;<br>
</p>

--089e013c6a7e5c5e6d0508d95e43--



More information about the MITgcm-support mailing list