<div dir="ltr">Thanks a lot Jean-Michel and Martin. I think these suggestions will help me do what I want. <div><br></div><div>Dhruv</div></div><br><div class="gmail_quote"><div dir="ltr">On Mon, Oct 22, 2018 at 11:34 AM Jean-Michel Campin <<a href="mailto:jmc@mit.edu">jmc@mit.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Dhruv and Martin,<br>
<br>
This is right: the full 3-D specification of GM and Redi coefficient is currently only <br>
available with pkg/ctrl. <br>
<br>
However, there is an option to combine a 1-D vertical profile with a 2-D horizontal map<br>
for both GM and Redi coefficient, without the need to use pkg/ctrl:<br>
 K_GM(i,j,k)   = GM_background_K * GM_bolFac1d(k) * GM_bolFac2d(i,j)<br>
 K_Redi(i,j,k) = GM_isopycK      * GM_isoFac1d(k) * GM_isoFac2d(i,j)<br>
and fields:<br>
 GM_bolFac1d  loaded from:  GM_bol1dFile<br>
 GM_bolFac2d  loaded from:  GM_bol2dFile<br>
 GM_isoFac1d  loaded from:  GM_iso1dFile<br>
 GM_isoFac2d  loaded from:  GM_iso2dFile<br>
<br>
Cheers,<br>
Jean-Michel<br>
<br>
On Mon, Oct 22, 2018 at 11:43:45AM +0200, Martin Losch wrote:<br>
> Hi Dhruv,<br>
> <br>
> Hong pointed out to me, that I was wrong: it is possible to specify 3D GM and Redi diffusivities. There are parameters in the gmredi namelist (GM_background_K3dFile, GM_isopycK3dFile), but the code for reading the fields is connected to the ctrl-package. The corresponding 5d-fields KapGM and and KapRedi are defined in pkg/ctrl/CTRL_FIELDS.h, so that they are only available when the control package is compiled (set in pacackage.config) and the CPP-flags ALLOW_KAPGM_CONTROL, ALLOW_KAPREDI_CONTROL are defined.<br>
> <br>
> There is probably a good reason for doing this in this way, but I am not quite sure if I can suggest how you want use this information. Either you compile the ctrl-package (which will require that you add a definition of nobcs into SIZE.h, please example SIZE.h in the verification/*/code_ad/ directories) and define the cpp flags and keep your fingers crossed, or you take the definition of fields out of the ctrl package and put them into the GMREDI-header file and adjust ini_mixing and/or gmredi_init_varia.F (where I would have put reading the fields in the first place, but then again, since they have been meant for ecco-type optimization, there is probably a good reason for doing it they way it has been done).<br>
> <br>
> Maybe others have an opinion on this?<br>
> <br>
> Martin<br>
> <br>
> > On 19. Oct 2018, at 15:55, Martin Losch <<a href="mailto:Martin.Losch@awi.de" target="_blank">Martin.Losch@awi.de</a>> wrote:<br>
> > <br>
> > Hi Dhruv,<br>
> > <br>
> >> On 19. Oct 2018, at 03:21, Dhruv Balwada <<a href="mailto:db194@nyu.edu" target="_blank">db194@nyu.edu</a>> wrote:<br>
> >> <br>
> >> Hi, <br>
> >> <br>
> >> I am running a MITgcm channel simulation in coarse resolution with the GMRedi package, forced with winds and buoyancy at the surface. I also advect tracers using the ptracers package. <br>
> >> I have 2 questions:<br>
> >> a) Is it possible to provide MITgcm with an arbitrary 3D GM diffusivity and a 3D Redi diffusivity field as input, in the current MITgcm release? If 3D fields cannot be specified, then is it possible to provide a vector with the depth variation? I am interested in prescribing arbitrary fields, rather than a particular diffusivity parameterization scheme.<br>
> > <br>
> > I don???t think so, you can only specify 3D fields of vertical diffusivies, afaik.<br>
> >> <br>
> >> b) When the model is being run with both GMRedi and ptracers packages, what is the  flux diagnostic for the ptracer tracer that has the influence of the GMRedi diffusion tensor? <br>
> > Don???t the normal diffusive flux diagnostics should contain this (DFxETr01, DFyETr01, DFrETr01, DFrITr01, etc.)?<br>
> > <br>
> > Martin<br>
> >> <br>
> >> Thank you,<br>
> >> Dhruv<br>
> >> _______________________________________________<br>
> >> MITgcm-support mailing list<br>
> >> <a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
> >> <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.mitgcm.org_mailman_listinfo_mitgcm-2Dsupport&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=3huS5z96YHlUbTqcqC5pEg&m=1mJdyRCUuog3AXxsKEs6r5oHRi4CU4_EK0cWY-kQAn0&s=3CKiZEh9yDizkNEIUfxMmq9KoZfg8TCU51sj1asXZ2A&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.mitgcm.org_mailman_listinfo_mitgcm-2Dsupport&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=3huS5z96YHlUbTqcqC5pEg&m=1mJdyRCUuog3AXxsKEs6r5oHRi4CU4_EK0cWY-kQAn0&s=3CKiZEh9yDizkNEIUfxMmq9KoZfg8TCU51sj1asXZ2A&e=</a><br>
> > <br>
> <br>
> _______________________________________________<br>
> MITgcm-support mailing list<br>
> <a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
> <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.mitgcm.org_mailman_listinfo_mitgcm-2Dsupport&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=3huS5z96YHlUbTqcqC5pEg&m=1mJdyRCUuog3AXxsKEs6r5oHRi4CU4_EK0cWY-kQAn0&s=3CKiZEh9yDizkNEIUfxMmq9KoZfg8TCU51sj1asXZ2A&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.mitgcm.org_mailman_listinfo_mitgcm-2Dsupport&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=3huS5z96YHlUbTqcqC5pEg&m=1mJdyRCUuog3AXxsKEs6r5oHRi4CU4_EK0cWY-kQAn0&s=3CKiZEh9yDizkNEIUfxMmq9KoZfg8TCU51sj1asXZ2A&e=</a><br>
_______________________________________________<br>
MITgcm-support mailing list<br>
<a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.mitgcm.org_mailman_listinfo_mitgcm-2Dsupport&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=3huS5z96YHlUbTqcqC5pEg&m=1mJdyRCUuog3AXxsKEs6r5oHRi4CU4_EK0cWY-kQAn0&s=3CKiZEh9yDizkNEIUfxMmq9KoZfg8TCU51sj1asXZ2A&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.mitgcm.org_mailman_listinfo_mitgcm-2Dsupport&d=DwICAg&c=slrrB7dE8n7gBJbeO0g-IQ&r=3huS5z96YHlUbTqcqC5pEg&m=1mJdyRCUuog3AXxsKEs6r5oHRi4CU4_EK0cWY-kQAn0&s=3CKiZEh9yDizkNEIUfxMmq9KoZfg8TCU51sj1asXZ2A&e=</a><br>
</blockquote></div>