 azimuthal velocity to represent the zonal flow generated in the overlying envelope. The shell rotates around the $z$-axis at the imposed rotation rate $\Omega$. The aspect ratio is $\gamma=r_i/r_o$ where $r_i$ is the inner sphere radius, corresponding to a rocky core, and $r_o$ the outer sphere radius, corresponding to the top of the fully conducting region. The fluid is assumed incompressible with constant density $\rho$ and constant temperature, that is, no convective motions are computed. The assumption of incompressibility is made for simplicity, although the pressure scale height at the depths of the conducting layer is roughly $8000$km [@Gui04], that is, about $1/5$ of the thickness of the layer. The effects of compressibility may well play a role in the dynamics of the conducting regions [see for instance @Evo04].

For simplicity we model the angular momentum coupling with the external zonal flow as a rigid boundary condition for the velocity at the outer boundary, rather than as a shear stress condition. The flow is driven through a boundary forcing rather than a volume forcing to avoid directly imposing bidimensionality to the velocity field. As we are interested in the bulk magnetohydrodynamical process, the exact nature of the coupling (electromagnetic or viscous, shear stress or rigid) with the upper molecular hydrogen layer is not crucial for our study. We discuss the implication of the choice of the rigid boundary condition in section \[sec:hydro\_jupiter\]. The radial profile of electrical conductivity is not well constrained in the gas giants. In particular the existence of a first order or continuous transition between the molecular and metallic hydrogen phase is still an open question, although high-pressure experiments are in favor of a continuous transition [@Nel99]. We choose to model the outer boundary as electrically insulating to simplify the coupling between the layers. The conductivity is assumed constant throughout the whole modeled conducting layer. As we do not model the molecular hydrogen layer, we assume zonal geostrophic balance within this envelope for simplicity. The amplitude of the zonal motions at the outer boundary of our model is therefore the same as the surface winds. This idealized representation of the dynamics of the molecular hydrogen layer would be altered if the magnetic forces upset the zonal geostrophic balance. Depending on the magnitude and radial profile of the electrical conductivity, the amplitude of the zonal motions might be reduced, and the zonal flow contours would tend
azimuthal velocity to represent the zonal flow generated in the overlying envelope. The carapace rotate around the $ z$-axis at the imposed rotation rate $ \Omega$. The aspect proportion is $ \gamma = r_i / r_o$ where $ r_i$ is the inner sphere radius, match to a rocky core, and $ r_o$ the out sphere radius, represent to the top of the fully conduct area. The fluid is assumed incompressible with constant concentration $ \rho$ and constant temperature, that is, no convective apparent motion are calculate. The assumption of incompressibility is made for simplicity, although the imperativeness scale height at the astuteness of the conducting layer is approximately $ 8000$km [ @Gui04 ], that is, about $ 1/5 $ of the thickness of the level. The effects of compressibility may well play a role in the dynamics of the conducting area [ see for instance @Evo04 ]. 

 For simplicity we model the angular momentum coupling with the external zonal flow as a inflexible boundary condition for the velocity at the outer boundary, rather than as a shear stress condition. The flow is driven through a boundary storm rather than a volume force to invalidate directly imposing bidimensionality to the velocity discipline. As we are concerned in the bulk magnetohydrodynamical process, the exact nature of the coupling (electromagnetic or syrupy, shear stress or rigid) with the upper molecular hydrogen layer is not crucial for our study. We discuss the implication of the choice of the rigid boundary condition in section   \[sec: hydro\_jupiter\ ]. The radial visibility of electrical conduction is not well constrained in the gas giants. In particular the existence of a first holy order or continuous transition between the molecular and metallic hydrogen phase is still an open question, although high - pressure experiments are in favor of a continuous conversion [ @Nel99 ]. We choose to model the outer limit as electrically insulating to simplify the coupling between the layers. The conductivity is assumed changeless throughout the whole model conducting layer. As we do not model the molecular hydrogen level, we assume zonal geostrophic balance within this envelope for simplicity. The amplitude of the zonal motion at the outer limit of our model is therefore the same as the airfoil winds. This idealized representation of the dynamics of the molecular hydrogen layer would be change if the magnetic forces disturb the zonal geostrophic balance. Depending on the magnitude and radial profile of the electrical conductivity, the amplitude of the zonal motion might be reduced, and the zonal flow contours would tend
 azlmuthal velocity to reprtsent the zonal flow genxrated jn the oxerlying envelope. The shell cotares aeound the $z$-axis at the imposed gotation eate $\Imega$. The eapect ratio ia $\gamke=r_i/r_o$ where $r_i$ ls the innes sphere radiuv, zoxresponding to a rocky core, and $r_o$ tre outet dphere radius, sorrtspjndihg to the top of the fully conductjng regpon. The fluid is sssumed incompressible witj cojstant density $\rho$ and constabt tqnperature, thxt is, no convective mojions are computed. The assumptiov of nncompressivioitj is made foc simpjicity, althombh the pressute scale heighb at vhe eepths of the conductmng layer is roughly $8000$km [@Gui04], tvac is, about $1/5$ of the thucjness of dhe uqyef. Tge erfects of compressigility may qell play a role in trv dynamics of the cjnqucting regions [see for instance @Evo04].

For sijplicity we model the abgular momentum coupllng with ehe external zonal flow as a rigid boundary condidion hof tkc xwllcity at the outer boundary, rather than as a fgesr stress conditlon. The flow is drovfn jhrough a bounaary forcjng rather than a golume sorcibg to avoyd dorectly imposing bidimensiobality to thv veoocity field. As we are intererted in tne bulk magnetohydrodyncmical process, thf exact nzgure of the coupuinn (enectromagnetic or viscous, fhear strxss ox rigid) dith the u[per molecklar mfdrogen layer is nlt crocial xor our stkdy. We discuss the implication of the choice pf thv rigid bjundavy condition in section \[sec:hydto\_jupiter\]. The fadial prorile of electrical conductivity ls not well ronstrainqd ib thw gas gkxnts. In particilar the vxnstence od a first order or covfinuous transitnub between the mplezulwr aid meedllic hydrogan pfasd is sgill an opek qjestoon, although high-prevsurs experiments are on favor od a contynuous transiyion [@Nel99]. We choose to mmdeu the oujer boundary as electrically inaulating ho fimplify the couilind between che layers. The conductivity is assumed cinstant throughouj the whole modeled cpkducting layxr. As re do not model the molecular hydrogen layer, wt assume zonal geostropgic banance within this envelope for simplicity. The amplitude of the zonal motions ar the outer boundady og our mjbel is eherxfure the same as nhe surface winds. This idealized represettction of the dynamics of the kouecular hydroeen layer would be altersd if thr magnetic forces upset the zonsl geostrophic balance. Dcpending in ths magnpjudx and radial profile of thz elecrrical xondmctivity, the aoplotmde of the zonel jotions might be reduced, and the zonal fpow contours wound tend
azimuthal velocity to represent the zonal flow the envelope. The rotates around the rate The aspect ratio $\gamma=r_i/r_o$ where $r_i$ the inner sphere radius, corresponding to rocky core, and $r_o$ the outer sphere radius, corresponding to the top of fully conducting region. The fluid is assumed incompressible with constant density $\rho$ and temperature, is, convective are computed. The assumption of incompressibility is made for simplicity, although the pressure scale height at depths of the conducting layer is roughly $8000$km that is, about $1/5$ the thickness of the layer. effects compressibility may play role the dynamics of conducting regions [see for instance @Evo04]. For simplicity we model the angular momentum coupling with the external flow as boundary condition the at outer boundary, rather a shear stress condition. The flow a boundary forcing rather than a volume forcing avoid directly bidimensionality to the velocity field. As are interested in the bulk magnetohydrodynamical process, the nature of the coupling (electromagnetic or viscous, shear stress or rigid) with the upper molecular is not crucial for study. We discuss implication the of rigid boundary in section \[sec:hydro\_jupiter\]. The radial profile of electrical conductivity is not constrained in the gas giants. In particular the existence of order continuous transition between molecular and metallic hydrogen is an open question, although are favor transition We to model the outer as electrically insulating to simplify coupling between the layers. throughout the whole modeled conducting layer. As we not model the molecular hydrogen layer, we zonal geostrophic balance within this envelope for simplicity. The amplitude of the motions at boundary of our model is therefore the same the surface winds. This representation of the dynamics of the molecular hydrogen layer be if the forces upset the geostrophic balance. Depending the magnitude and of the conductivity, amplitude motions might be reduced, and the flow would tend
 azimuthal velocity to represEnt the zonaL flow GenEraTeD in tHe ovErlying envelopE. the sHell rotates around the $z$-aXis at ThE ImpoSEd RotatIon rate $\oMeGA$. the AsPeCt rAtIO iS $\gammA=r_i/R_o$ where $R_i$ is the innEr sPhEre radius, corREsPonding to a RocKy core, and $r_o$ tHe oUter spHeRe rADius, cOrrEsponDing to THe top oF the fully CoNDuctinG Region. THE FlUid iS assumed incompresSIbLE with constant dEnsity $\RhO$ AnD COnsTanT temperatuRe, That iS, No conveCTiVE MOtiONs are computed. the assumptiON of IncompReSsiBIlity iS made FoR SimPlicity, althOugh The pressuRe scalE Height aT The deptHs of thE coNduCtinG LaYeR is RoUGhlY $8000$Km [@gui04], THat Is, about $1/5$ oF tHe ThickNess OF THE layEr. THe efFects Of compressibiLitY may WEll Play a Role iN the DyNamicS of the ConduCtIng regions [see foR insTance @Evo04].

FOr sImPliCiTy we mODel the AngUlaR momentUm couplINg wItH THE eXternal zonal flow as A rIGId Boundary ConditIOn FoR The velocItY at The oUTEr bouNdarY, RaTher than As a sheAR sTrEss condItIon. The FlOw iS drIven tHRougH a bounDary forcIng raTHer than a volume FOrcing to avoid DIrECTlY ImpoSinG bidimensioNaliTY to tHe veLOcIty FIeld. AS we arE iNTeREsted in the bulk magneToHydrodYnamiCal process, the Exact naturE OF The couplIng (eLEcTRomagnetic or viScous, Shear stresS Or rigid) wIth thE upper moLecular hyDROgen layeR is Not CruCiaL FOr Our study. We disCUSs thE iMplicatIon Of the chOicE of The RigId Boundary cOndition In SeCtIoN \[seC:hydrO\_Jupiter\]. THe RadIaL prOfile OF electRical CondUcTiVIty Is not weLL cONStraInEd In thE gaS gIants. in paRTicUlar the Existence Of a FIrst OrDeR or contInuous transitIoN between thE mOleCular aND Metallic Hydrogen phase is still an oPEn questIon, AlthoUgh hIgh-pressuRe eXperimEntS Are in fAvor of A contInUouS TRansiTIOn [@nel99]. we Choose to moDEL thE outeR bOundAry as elEctrically insulatiNG to Simplify the coUplIng bETWeEn tHE lAYerS. THE coNDUctivity is assumEd constant ThROuGhout the whOLe mOdEled conDucting Layer. aS we do noT model the Molecular HyDrogEN LayEr, we assume Zonal geoStrophic bALance WItHin thIs eNvelopE fOr sImpliCity. ThE AmpLitudE of the ZoNal motIons aT tHe outer bOundary of our model is therEfore tHe samE as The surfacE wiNDs. THis idealiZed rEpresentatIon Of tHe dynAmiCS of thE molECuLar HYdrogEn laYEr would be ALtEreD IF tHe magnetic fORCEs uPset tHe zONal geoStroPhic balance. DependINg on the magnituDe anD RAdiAl pROfilE oF the electrical ConDuCTIvity, the AmPlitude of thE zonal moTiONs migHt be reDuced, aNd the zoNAL fLOw contOurs WouLd tend
 azimuthal velocity to rep resent the  zona l f low  g ener ated  in the overly i ng e nvelope. The shell rot atesar o undt he  $z$- axis at th e imp os ed  ro ta t io n rat e $ \Omega$ . The aspe ctra tio is $\gam m a= r_i/r_o$ w her e $r_i$ is t heinnersp her e  radi us,  corr espond i ng toa rocky c or e , and$ r_o$ th e ou tersphere radius, co r re s ponding to the  top o ft he f ull y c onductingre gion. The flu i di s ass u med incompres sible withc ons tant d en sit y  $\rho $ and  c o nst ant tempera ture , that is , no c o nvectiv e  motion s arecom put ed.T he  a ssu mp t ion of  in c omp ressibil it yis ma de f o r s impl ici ty,altho ugh the press ure  sca l e h eight  at t he d ep ths o f thecondu ct ing layer is ro ughl y $8000$k m [ @G ui0 4] , tha t  is, a bou t $ 1/5$ of  the th i ckn es s o fthe layer. The eff ec t s  o f compre ssibil i ty  m a y well p la y a  rol e in th e dy n am ics of t he con d uc ti ng regi on s [see  f orins tance @Evo 04].
For simp licit y  we model thea ngular moment u mc o up l ingwit h the exter nalz onal  flo w  a s a rigid  boun da r yc ondition for the ve lo city a t the  outer bounda ry, rather t h an as ashea r  s t ress condition . The  flow is d r iven thr ougha bounda ry forcin g rather t han  avol ume f or cing to avoid d irec tl y impos ing  bidime nsi ona lit y t othe veloc ity fiel d.  A swe  ar e int e rested i nthe  b ulk  magn e tohydr odyna mica lpr o ces s, thee xa c t  nat ur eof t heco uplin g (e l ect romagne tic or vi sco u s, s he ar  stress  or rigid) wi th  the upper  m ole cularh y drogen l ayer is not crucial for our stu dy.  We d iscu ss the im pli cation  of the ch oice o f the  r igi d bound a r ycon di tion in se c t ion  \[se c: hydr o\_jupi ter\]. The radialp rof ile of electr ica l co n d uc tiv i ty isno t  we l l  constrained in  the gas g ia n ts . In parti c ula rthe exi stenceof af irst or der or co ntinuoustr ansi t i onbetween th e molecu lar and m e talli c  h ydrog enphaseis  st ill a n open que stion , alth ou gh hig h-pre ss ure expe riments are in favor of  a con tinuo ustransitio n [ @ Nel 99]. We c hoos e to model  th e o uterbou n daryas e l ec tri c allyinsu l ating tos im pli f y  t he coupling b e twe en th e l a yers.Theconductivity is a s sumed constant  thr o u gho utt he w ho le modeled con duc ti n g  layer.As  we do notmodel th em olecu lar hy drogen  layer, w ea ssumezona l g eostrophi c b al a nce wit hi nt his en velo pe  for s implic i ty.T h e amplitude of t he zo n a l mot i ons  at t he  outerb ound ary of our  model is t herefo re t he sa me as t he  surfa cewi nds. Thisi dealizedrepre sentati on  ofthe  dynam icso f  themole cu lar  hydrogen l ay e rwo u ldbe a ltere dif t he magnet i c forces  up s et thezo nal g eostro p hi c balance. D epe nding o n the magn i tude an d  radi al pro file o f the e l ect ri cal con duc t i vity, the  amplitud e  of  t he z onal mot io ns m ig htbe  reduced, and t he zonal  flow co n t ours w oul d  ten d 
 azimuthal_velocity to_represent the zonal flow_generated in_the_overlying envelope._The_shell rotates around_the $z$-axis at_the imposed rotation rate_$\Omega$. The aspect_ratio_is $\gamma=r_i/r_o$ where $r_i$ is the inner sphere radius, corresponding to a rocky core,_and_$r_o$ the_outer_sphere_radius, corresponding to the top_of the fully conducting region._The fluid_is assumed incompressible with constant density $\rho$ and_constant_temperature, that is,_no convective motions are computed. The assumption of incompressibility_is made for simplicity, although the_pressure scale height_at_the_depths of the conducting_layer is roughly $8000$km [@Gui04], that_is, about $1/5$ of the thickness_of the layer. The effects of compressibility_may well play a role in_the dynamics of the conducting_regions [see_for instance @Evo04].

For simplicity we_model the angular_momentum coupling_with the external_zonal flow as a rigid boundary_condition for the_velocity at the outer boundary, rather_than_as a shear_stress_condition._The flow_is driven through_a_boundary forcing_rather_than a volume forcing to avoid_directly_imposing bidimensionality to the velocity field. As_we are interested in_the_bulk magnetohydrodynamical process, the_exact nature of the coupling_(electromagnetic or viscous, shear stress or_rigid) with_the upper_molecular hydrogen layer is not crucial for our study. We discuss_the implication of the choice of_the rigid boundary condition_in section \[sec:hydro\_jupiter\]._The_radial profile of_electrical_conductivity is_not well constrained in the gas giants._In particular_the existence of a first order_or continuous transition between_the_molecular and metallic hydrogen phase is_still an open question, although high-pressure_experiments are in favor of_a_continuous_transition [@Nel99]. We choose to_model the outer boundary as electrically_insulating to simplify_the coupling between the layers. The conductivity_is_assumed constant throughout the whole modeled_conducting_layer. As we do not model_the_molecular_hydrogen layer, we assume zonal_geostrophic balance within this envelope for_simplicity. The amplitude of the zonal motions at the_outer boundary of_our model is therefore the_same_as_the surface winds. This idealized representation of the dynamics of_the molecular_hydrogen layer would_be altered if the magnetic forces upset the zonal geostrophic_balance. Depending on the magnitude and radial_profile of the electrical conductivity, the amplitude of the zonal motions_might be reduced, and the zonal flow contours_would tend