<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Hello<div class=""><br class=""></div><div class=""><br class=""></div><div class="">You turn the fluxes into tendency by, for example, </div><div class=""><br class=""></div><div class=""><div class=""><div class="">for k = 1:NZ</div><div class=""> VOL(:,:,k) = hFacC(:,:,k).*RAC.*DRF(k);</div><div class=""> AREAWEST(:,:,k) = DYG.*DRF(k).*hFacW(:,:,k);</div><div class=""> AREASOUTH(:,:,k) = DXG.*DRF(k).*hFacS(:,:,k);</div><div class=""> AREACELL(:,:,k) = RAC;</div><div class="">end</div><div class="">RecipVol = VOL;</div></div></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><div class=""><div class="">%zonal</div><div class="">Q = ADVx;</div><div class="">tm1 = Q(ILON+1,ILAT,IDEP);% 'ADVx_TH '</div><div class="">tm2 = Q(ILON ,ILAT,IDEP);% 'ADVx_TH '</div><div class="">Ax = (tm1 - tm2);</div><div class="">%meridional</div><div class="">Q = ADVy;</div><div class="">tm1 = Q(ILON,ILAT+1,IDEP);% 'ADVy_TH '</div><div class="">tm2 = Q(ILON,ILAT ,IDEP);% 'ADVy_TH '</div><div class="">Ay = (tm1 - tm2);</div><div class="">%vertical</div><div class="">Q1 = zeros(NX,NY,NZ+1,'single');%extra z lvl for vert der.</div><div class="">Q1(:,:,1:NZ) = ADVr;</div><div class="">tm1 = Q1(ILON,ILAT,IDEP );% 'ADVr_TH '</div><div class="">tm2 = Q1(ILON,ILAT,IDEP+1);% 'ADVr_TH '</div><div class="">Az = (tm1 - tm2);</div></div></div><div class=""><br class=""></div><div class="">%Also need this to close budget</div><div class=""><div class="">%zonal</div><div class="">Q = UVEL;</div><div class="">tm1 = Q(ILON+1,ILAT,IDEP).*AREAWEST(ILON+1,ILAT,IDEP);%UVEL</div><div class="">tm2 = Q(ILON ,ILAT,IDEP).*AREAWEST(ILON,ILAT,IDEP);%UVEL</div><div class="">DU = (tm1 - tm2);</div><div class="">%meridional</div><div class="">Q = VVEL;</div><div class="">tm1 = Q(ILON,ILAT+1,IDEP).*AREASOUTH(ILON,ILAT+1,IDEP);%VVEL</div><div class="">tm2 = Q(ILON,ILAT ,IDEP).*AREASOUTH(ILON,ILAT,IDEP);%VVEL</div><div class="">DV = (tm1 - tm2);</div><div class="">%vertical</div><div class="">Q1(:,:,1:NZ) = WVEL;</div><div class="">Q1(:,:,1)=0; %NEED TO SET SURFACE W*T = 0</div><div class="">tm1 = Q1(ILON,ILAT,IDEP ).*AREACELL(ILON,ILAT,IDEP);%WVEL</div><div class="">tm2 = Q1(ILON,ILAT,IDEP+1).*AREACELL(ILON,ILAT,IDEP+1);%WVEL</div><div class="">DW = (tm1 - tm2);</div></div><div class=""><br class=""></div><div class=""><div class=""><br class=""></div></div><div class=""><div class="">gTracer = -RecipVol(ILON,ILAT,IDEP) ...</div><div class=""> .*( (Ax+Ay+Az);</div></div><div class=""> - TH.*(DU+DV+DW) );</div><div class=""><br class=""></div><div class="">Diffusive are handled like advective, so would just add them in, i.e.: </div><div class=""><div class="">gTracer = -RecipVol(ILON,ILAT,IDEP) ...</div><div class=""> .*( (Ax+Ay+Az+Dx+Dy+Dz);</div></div><div class=""> - TH.*(DU+DV+DW) );</div><div class=""><br class=""></div><div class="">Hope this helps</div><div class=""><br class=""></div><div class="">Matt</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""><div><blockquote type="cite" class=""><div class="">On Sep 24, 2017, at 6:33 AM, liuyi315 <<a href="mailto:liuyi315@mails.ucas.ac.cn" class="">liuyi315@mails.ucas.ac.cn</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; font-family: 新宋体; font-size: 14px; line-height: 1.7;" class=""><div class="">Hi, all users,</div><div class="">I am trying to do some diagnostics analysis of heat transport and vertical redistribution in the Indonesian Seas using MITgcm. I have output ADVx/y/r_TH, DFr/x/yE_TH and DFrI_TH using "data.diagnostics". I found the unit of these variables is "degC.m^3/s", and my understanding is that the advection term of heat flux (i.e. ADVx/y/r_TH) is the product of potential temperature and volume transport. However, after I checked the results, I found the magnitude of heat transport is about one to two orders larger than observed data, although the simulated temperature and velocity compares well with observed data. Thus <span style="line-height: 1.7;" class="">I believe that my understanding of the calculation of heat flux advection term has something wrong. </span></div><div class=""><span style="line-height: 1.7;" class="">Could anyone tell me how the model computes the advection and diffusion term of heat flux?</span></div><div class=""><span style="line-height: 1.7;" class="">Thanks very much.</span></div><div class=""><span style="line-height: 1.7;" class=""><br class=""></span></div><div class=""><span style="line-height: 1.7;" class="">Best regards, </span></div><div class=""><span class=""><br class=""></span></div><div class=""><span class="">Yi Liu</span></div><div id="ntes-pcmail-signature" style="font-family: 新宋体;" class=""></div><br class=""></div><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">_______________________________________________</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; float: none; display: inline !important;" class="">MITgcm-support mailing list</span><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><a href="mailto:MITgcm-support@mitgcm.org" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">MITgcm-support@mitgcm.org</a><br style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a></div></blockquote></div><br class=""></div></body></html>