[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