[MITgcm-support] boundary condition on tracer
Joseph LaCasce
josephhl at met.no
Thu Mar 31 03:25:15 EST 2005
Hello!
I'm trying to understand how the boundary conditions are implemented for
a tracer in the model. Take a square basin with a flat bottom. It looks
like---
in GAD_GRAD_X for instance---that no-heat flux is imposed on one side
of the basin (the west?). But where is it imposed for the other boundary?
Thanks!
Joe LaCasce
*DO* j <http://mitgcm.org/sealion/code_reference/vdb/names/8M.htm>=1-Oly <http://mitgcm.org/sealion/code_reference/vdb/names/333.htm>,sNy <http://mitgcm.org/sealion/code_reference/vdb/names/107.htm>+Oly <http://mitgcm.org/sealion/code_reference/vdb/names/333.htm>
dfx <http://mitgcm.org/sealion/code_reference/vdb/names/1OT.htm>*(*1-Olx <http://mitgcm.org/sealion/code_reference/vdb/names/334.htm>,j <http://mitgcm.org/sealion/code_reference/vdb/names/8M.htm>*)*=0.
*DO* i <http://mitgcm.org/sealion/code_reference/vdb/names/8L.htm>=1-Olx <http://mitgcm.org/sealion/code_reference/vdb/names/334.htm>+1,sNx <http://mitgcm.org/sealion/code_reference/vdb/names/106.htm>+Olx <http://mitgcm.org/sealion/code_reference/vdb/names/334.htm>
dfx <http://mitgcm.org/sealion/code_reference/vdb/names/1OT.htm>*(*i <http://mitgcm.org/sealion/code_reference/vdb/names/8L.htm>,j <http://mitgcm.org/sealion/code_reference/vdb/names/8M.htm>*)* = -diffKh <http://mitgcm.org/sealion/code_reference/vdb/names/1PG.htm>*xA <http://mitgcm.org/sealion/code_reference/vdb/names/EF.htm>*(*i <http://mitgcm.org/sealion/code_reference/vdb/names/8L.htm>,j <http://mitgcm.org/sealion/code_reference/vdb/names/8M.htm>*)*
& *_recip_dxC <http://mitgcm.org/sealion/code_reference/vdb/names/3Z1.htm>*(*i <http://mitgcm.org/sealion/code_reference/vdb/names/8L.htm>,j <http://mitgcm.org/sealion/code_reference/vdb/names/8M.htm>,bi <http://mitgcm.org/sealion/code_reference/vdb/names/84.htm>,bj <http://mitgcm.org/sealion/code_reference/vdb/names/85.htm>*)*
& **(*Tracer <http://mitgcm.org/sealion/code_reference/vdb/names/1UR.htm>*(*i <http://mitgcm.org/sealion/code_reference/vdb/names/8L.htm>,j <http://mitgcm.org/sealion/code_reference/vdb/names/8M.htm>*)*-Tracer <http://mitgcm.org/sealion/code_reference/vdb/names/1UR.htm>*(*i <http://mitgcm.org/sealion/code_reference/vdb/names/8L.htm>-1,j <http://mitgcm.org/sealion/code_reference/vdb/names/8M.htm>*)**)*
& *CosFacU <http://mitgcm.org/sealion/code_reference/vdb/names/90L.htm>*(*j <http://mitgcm.org/sealion/code_reference/vdb/names/8M.htm>,bi <http://mitgcm.org/sealion/code_reference/vdb/names/84.htm>,bj <http://mitgcm.org/sealion/code_reference/vdb/names/85.htm>*)*
*ENDDO*
*ENDDO*
--
J.H. LaCasce, Phone: +47 22 96 33 11; FAX: +47 22 69 63 55
Norwegian Meteorological Institute, P.O. Box 43 Blindern, 0313 Oslo
or (UPS) Niels Henrik Abelsvei 40, 0371 Oslo, Norway
http://met.no/english/r_and_d_activities/people/josephhl.html
More information about the MITgcm-support
mailing list