<div dir="ltr"><div class="gmail_default" style="font-size:large">Hi Manfredi,</div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">Sorry for the delay in answering you, but these days were quite intense.<br>I tried to NOT specify the two boolean variables in the data.darwin namefile, so now darwin_useEXFwind" is FALSE for default. Then I've left there only my wind speed file for the gas exchange as Christoph suggested, while I left the physics forcing wind in the data namefile.<br></div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">Now it seems to work and it doesn't return the warning related to the constant wind.</div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">I want to do a double-check by looking at the diagnostic output. I will let you know if all works properly.</div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">Thank you very much for the help.</div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">Best regards</div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large">Francesco</div><div class="gmail_default" style="font-size:large"><br></div><div class="gmail_default" style="font-size:large"><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno gio 6 lug 2023 alle ore 20:12 Manfredi Manizza <<a href="mailto:mmanizza@ucsd.edu" target="_blank">mmanizza@ucsd.edu</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><span style="font-size:large">Hi Francesco,</span></div><div><span style="font-size:large"><br></span></div><div><span style="font-size:large">I can't look at the code now because the machine is currently down but I think that </span></div><div><span style="font-size:large">if you read the wind file from "data" and NOT from "data.exf" then </span><span style="font-size:large"> "darwin_useEXFwind = .TRUE." </span></div><div><span style="font-size:large">is NOT correct. What if you set both FALSE? </span></div><div><span style="font-size:large">Will it compile? What are the other options in the code.</span></div><div><span style="font-size:large"><br></span></div><div><font size="4">As Christoph rightly suggested in the various bgc packages of MITgcm </font></div><div><font size="4">the wind speed is used only for the gas exchange coefficient calculation.</font></div><div><font size="4">The wind for ocean physics calculation (as MLD) will be read in other parts of the code.</font></div><div><font size="4"><br></font></div><div><font size="4">If it compiles and you have a wind file already in place you can also</font></div><div><font size="4">write a printout statement in the code to see which number you get</font></div><div><font size="4">in order to debug...</font></div><div><font size="4"><br></font></div><div><font size="4">A warning statement normally does not prevent the code compilation.</font></div><div><br></div><div><font size="4">Good luck.</font></div><div><font size="4">M.</font></div><div><font size="4"><br></font></div><div><font size="4"><br></font></div><div><span style="font-size:large"><br></span></div><div><span style="font-size:large"><br></span></div><span style="font-size:large"><div><span style="font-size:large"><br></span></div></span><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Jul 6, 2023 at 4:11 AM Christoph Voelker <<a href="mailto:christoph.voelker@awi.de" target="_blank">christoph.voelker@awi.de</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear all,<br>
<br>
I don't know how it is handled exactly in Darwin, but I think it is <br>
similar to the DIC package.<br>
<br>
The point here is that biogeochemical packages need wind speed for <br>
calculating the gas exchange rate of carbon (possibly also for other <br>
gases, like O2). That is why you have the possibility to read in a field <br>
of wind speed (note, the speed, not the u and v components). This field <br>
is then solely used for gas exchange, not for physics.<br>
<br>
The physics forcing by wind is either described in data (if you use wind <br>
stress files), or data.exf (if you use 10m winds).<br>
<br>
Why is this treated searately? One could also use the components, to <br>
calculate wind speed, of course. But especially if one uses monthly <br>
averaged winds, the wind speed calculated from these averages <br>
underestimates the averaged true wind speed, which also contains a <br>
contribution from the deviations from the average (u'*u' and v'*v').<br>
<br>
Cheers, Christoph<br>
<br>
<br>
On 06.07.23 12:26, Dattilo, Francesco wrote:<br>
> Thank you very much for the answer Manfredi,<br>
> <br>
> However, still having set the same options you provided to me, if I run <br>
> the model like this I receive this warning:<br>
> <br>
> "** WARNING ** DARWIN_EXF_READPARMS: default constant wind speed <br>
> ',windconst <<a href="https://urldefense.com/v3/__https://gud.mit.edu/darwin3/ident?v=darwin3;i=windconst__;!!Mih3wA!AnSU1oWcLN8nr6iiJ0kauXvdiWCntoLiyc9wc8P9El-B7SjcmoU4PW6-R4VMgMuyZd0hx2Q8v3lr1FY0ZEu3PNqDRSbU$" rel="noreferrer" target="_blank">https://urldefense.com/v3/__https://gud.mit.edu/darwin3/ident?v=darwin3;i=windconst__;!!Mih3wA!AnSU1oWcLN8nr6iiJ0kauXvdiWCntoLiyc9wc8P9El-B7SjcmoU4PW6-R4VMgMuyZd0hx2Q8v3lr1FY0ZEu3PNqDRSbU$</a> >,' <br>
> m/s is used'"<br>
> <br>
> This seems because I have to set "darwin_useEXFwind = .TRUE." and <br>
> "DARWIN_haveWind = .FALSE" in data.darwin.<br>
> <br>
> My confusion comes from this, because in the example mentioned before <br>
> there is no reference to these two boolean values.<br>
> <br>
> Do you have any ideas about it?<br>
> <br>
> Thank you very much<br>
> <br>
> Francesco<br>
> <br>
> Il giorno mer 5 lug 2023 alle ore 17:42 Manfredi Manizza <br>
> <<a href="mailto:mmanizza@ucsd.edu" target="_blank">mmanizza@ucsd.edu</a> <mailto:<a href="mailto:mmanizza@ucsd.edu" target="_blank">mmanizza@ucsd.edu</a>>> ha scritto:<br>
> <br>
> Hi Francesco,<br>
> <br>
> it depends on the files and the set-up.<br>
> <br>
> Normally on a 3D domain set-up you would use "data.exf" where you<br>
> specify the input file of wind speed if you use interannual winds<br>
> as in ECCO.<br>
> So you do not read that from the "data" file but it seems not to be<br>
> your case....<br>
> <br>
> You might want to check this 1-D set-up of the ross sea:<br>
> <br>
> <a href="https://urldefense.com/v3/__https://github.com/MITgcm-contrib/ecco_darwin/tree/master/idealized/1D_darwin/ross_sea/run__;!!Mih3wA!AnSU1oWcLN8nr6iiJ0kauXvdiWCntoLiyc9wc8P9El-B7SjcmoU4PW6-R4VMgMuyZd0hx2Q8v3lr1FY0ZEu3PB2vygDK$" rel="noreferrer" target="_blank">https://urldefense.com/v3/__https://github.com/MITgcm-contrib/ecco_darwin/tree/master/idealized/1D_darwin/ross_sea/run__;!!Mih3wA!AnSU1oWcLN8nr6iiJ0kauXvdiWCntoLiyc9wc8P9El-B7SjcmoU4PW6-R4VMgMuyZd0hx2Q8v3lr1FY0ZEu3PB2vygDK$</a> <<a href="https://urldefense.com/v3/__https://github.com/MITgcm-contrib/ecco_darwin/tree/master/idealized/1D_darwin/ross_sea/run__;!!Mih3wA!AnSU1oWcLN8nr6iiJ0kauXvdiWCntoLiyc9wc8P9El-B7SjcmoU4PW6-R4VMgMuyZd0hx2Q8v3lr1FY0ZEu3PB2vygDK$" rel="noreferrer" target="_blank">https://urldefense.com/v3/__https://github.com/MITgcm-contrib/ecco_darwin/tree/master/idealized/1D_darwin/ross_sea/run__;!!Mih3wA!AnSU1oWcLN8nr6iiJ0kauXvdiWCntoLiyc9wc8P9El-B7SjcmoU4PW6-R4VMgMuyZd0hx2Q8v3lr1FY0ZEu3PB2vygDK$</a> ><br>
> <br>
> It might give you a hint about how to read the wind in your 1-D case.<br>
> <br>
> In this experiment the wind forcing is read in the "data" file.<br>
> <br>
> Here a file where to follow the correct options:<br>
> <br>
> <a href="https://urldefense.com/v3/__https://github.com/MITgcm-contrib/ecco_darwin/blob/master/idealized/1D_darwin/ross_sea/code/EXF_OPTIONS.h__;!!Mih3wA!AnSU1oWcLN8nr6iiJ0kauXvdiWCntoLiyc9wc8P9El-B7SjcmoU4PW6-R4VMgMuyZd0hx2Q8v3lr1FY0ZEu3PFXBQ-1c$" rel="noreferrer" target="_blank">https://urldefense.com/v3/__https://github.com/MITgcm-contrib/ecco_darwin/blob/master/idealized/1D_darwin/ross_sea/code/EXF_OPTIONS.h__;!!Mih3wA!AnSU1oWcLN8nr6iiJ0kauXvdiWCntoLiyc9wc8P9El-B7SjcmoU4PW6-R4VMgMuyZd0hx2Q8v3lr1FY0ZEu3PFXBQ-1c$</a> <<a href="https://urldefense.com/v3/__https://github.com/MITgcm-contrib/ecco_darwin/blob/master/idealized/1D_darwin/ross_sea/code/EXF_OPTIONS.h__;!!Mih3wA!AnSU1oWcLN8nr6iiJ0kauXvdiWCntoLiyc9wc8P9El-B7SjcmoU4PW6-R4VMgMuyZd0hx2Q8v3lr1FY0ZEu3PFXBQ-1c$" rel="noreferrer" target="_blank">https://urldefense.com/v3/__https://github.com/MITgcm-contrib/ecco_darwin/blob/master/idealized/1D_darwin/ross_sea/code/EXF_OPTIONS.h__;!!Mih3wA!AnSU1oWcLN8nr6iiJ0kauXvdiWCntoLiyc9wc8P9El-B7SjcmoU4PW6-R4VMgMuyZd0hx2Q8v3lr1FY0ZEu3PFXBQ-1c$</a> ><br>
> <br>
> I never used it but I hope this can be useful for you to set-up<br>
> yours correctly.<br>
> <br>
> Good luck.<br>
> <br>
> Manfredi<br>
> <br>
> <br>
> <br>
> <br>
> <br>
> On Wed, Jul 5, 2023 at 2:15 AM Dattilo, Francesco <<a href="mailto:fdattilo@ogs.it" target="_blank">fdattilo@ogs.it</a><br>
> <mailto:<a href="mailto:fdattilo@ogs.it" target="_blank">fdattilo@ogs.it</a>>> wrote:<br>
> <br>
> Hi everyone,<br>
> <br>
> I'm working with the Darwin package. I'm a bit confused about<br>
> how to set wind as external forcing.<br>
> I saw that it is possible to specify the zonal wind and the<br>
> meridional wind in the "data" file, and this seems to be<br>
> consistent with the example<br>
> <<a href="https://urldefense.com/v3/__https://mitgcm.readthedocs.io/en/latest/examples/barotropic_gyre/barotropic_gyre.html*sec-eg-baro__;Iw!!Mih3wA!F0NpPv7kBp5gSoyJ2jM9Iva1V7k3g7RPj8LMYuP_v6xziAU6aZe1H0hJG2K2LEr1zkauGJQzZoIqHFirQA$" rel="noreferrer" target="_blank">https://urldefense.com/v3/__https://mitgcm.readthedocs.io/en/latest/examples/barotropic_gyre/barotropic_gyre.html*sec-eg-baro__;Iw!!Mih3wA!F0NpPv7kBp5gSoyJ2jM9Iva1V7k3g7RPj8LMYuP_v6xziAU6aZe1H0hJG2K2LEr1zkauGJQzZoIqHFirQA$</a>> . However, it seems that I have also another way, by modifying the "data.darwin" file, specifying the windFile and "darwin_useEXFwind = .TRUE." and "DARWIN_haveWind = .FALSE."<br>
> but in this case it is not clear to me:<br>
> <br>
> 1) which kind of file has to be the "windFile"<br>
> 2) if it is correct to set "darwin_useEXFwind = .TRUE." and<br>
> "DARWIN_haveWind = .FALSE."<br>
> <br>
> Using a 1D model I also don't know what's the proper choice.<br>
> <br>
> Thank you very much for every suggestion<br>
> <br>
> Best regards<br>
> <br>
> Francesco Dattilo<br>
> _______________________________________________<br>
> MITgcm-support mailing list<br>
> <a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a> <mailto:<a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a>><br>
> <a href="https://urldefense.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!F0NpPv7kBp5gSoyJ2jM9Iva1V7k3g7RPj8LMYuP_v6xziAU6aZe1H0hJG2K2LEr1zkauGJQzZoIJYEEeUQ$" rel="noreferrer" target="_blank">https://urldefense.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!F0NpPv7kBp5gSoyJ2jM9Iva1V7k3g7RPj8LMYuP_v6xziAU6aZe1H0hJG2K2LEr1zkauGJQzZoIJYEEeUQ$</a> <<a href="https://urldefense.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!F0NpPv7kBp5gSoyJ2jM9Iva1V7k3g7RPj8LMYuP_v6xziAU6aZe1H0hJG2K2LEr1zkauGJQzZoIJYEEeUQ$" rel="noreferrer" target="_blank">https://urldefense.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!F0NpPv7kBp5gSoyJ2jM9Iva1V7k3g7RPj8LMYuP_v6xziAU6aZe1H0hJG2K2LEr1zkauGJQzZoIJYEEeUQ$</a>><br>
> <br>
> <br>
> <br>
> -- <br>
> <br>
> Dr. Manfredi Manizza<br>
> Geosciences Research Division<br>
> Scripps Institution of Oceanography<br>
> University of California San Diego<br>
> 9500 Gilman Drive La Jolla, CA 92093-0244<br>
> email :<a href="mailto:mmanizza@ucsd.edu" target="_blank">mmanizza@ucsd.edu</a> <mailto:<a href="mailto:mmanizza@ucsd.edu" target="_blank">mmanizza@ucsd.edu</a>><br>
> tel : +1-858-534-7094<br>
> web :<a href="http://bluemoon.ucsd.edu/mmanizza/" rel="noreferrer" target="_blank">http://bluemoon.ucsd.edu/mmanizza/</a> <<a href="http://bluemoon.ucsd.edu/mmanizza/" rel="noreferrer" target="_blank">http://bluemoon.ucsd.edu/mmanizza/</a>> <br>
> <br>
> web :<a href="https://mmanizza.scrippsprofiles.ucsd.edu/" rel="noreferrer" target="_blank">https://mmanizza.scrippsprofiles.ucsd.edu/</a> <<a href="https://mmanizza.scrippsprofiles.ucsd.edu/" rel="noreferrer" target="_blank">https://mmanizza.scrippsprofiles.ucsd.edu/</a>><br>
> <br>
> _______________________________________________<br>
> MITgcm-support mailing list<br>
> <a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a> <mailto:<a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a>><br>
> <a href="https://urldefense.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!AnSU1oWcLN8nr6iiJ0kauXvdiWCntoLiyc9wc8P9El-B7SjcmoU4PW6-R4VMgMuyZd0hx2Q8v3lr1FY0ZEu3PHkhECh9$" rel="noreferrer" target="_blank">https://urldefense.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!AnSU1oWcLN8nr6iiJ0kauXvdiWCntoLiyc9wc8P9El-B7SjcmoU4PW6-R4VMgMuyZd0hx2Q8v3lr1FY0ZEu3PHkhECh9$</a> <br>
> <<a href="https://urldefense.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!AnSU1oWcLN8nr6iiJ0kauXvdiWCntoLiyc9wc8P9El-B7SjcmoU4PW6-R4VMgMuyZd0hx2Q8v3lr1FY0ZEu3PHkhECh9$" rel="noreferrer" target="_blank">https://urldefense.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!AnSU1oWcLN8nr6iiJ0kauXvdiWCntoLiyc9wc8P9El-B7SjcmoU4PW6-R4VMgMuyZd0hx2Q8v3lr1FY0ZEu3PHkhECh9$</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.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!AnSU1oWcLN8nr6iiJ0kauXvdiWCntoLiyc9wc8P9El-B7SjcmoU4PW6-R4VMgMuyZd0hx2Q8v3lr1FY0ZEu3PHkhECh9$" rel="noreferrer" target="_blank">https://urldefense.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!AnSU1oWcLN8nr6iiJ0kauXvdiWCntoLiyc9wc8P9El-B7SjcmoU4PW6-R4VMgMuyZd0hx2Q8v3lr1FY0ZEu3PHkhECh9$</a> <br>
<br>
-- <br>
Christoph Völker<br>
Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research<br>
Am Handelshafen 12<br>
27570 Bremerhaven<br>
<br>
+49 (0)471-4831-1848<br>
<a href="mailto:christoph.voelker@awi.de" target="_blank">christoph.voelker@awi.de</a><br>
<a href="https://urldefense.com/v3/__https://orcid.org/0000-0003-3032-114X__;!!Mih3wA!AnSU1oWcLN8nr6iiJ0kauXvdiWCntoLiyc9wc8P9El-B7SjcmoU4PW6-R4VMgMuyZd0hx2Q8v3lr1FY0ZEu3PAGXMfHj$" rel="noreferrer" target="_blank">https://urldefense.com/v3/__https://orcid.org/0000-0003-3032-114X__;!!Mih3wA!AnSU1oWcLN8nr6iiJ0kauXvdiWCntoLiyc9wc8P9El-B7SjcmoU4PW6-R4VMgMuyZd0hx2Q8v3lr1FY0ZEu3PAGXMfHj$</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.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!AnSU1oWcLN8nr6iiJ0kauXvdiWCntoLiyc9wc8P9El-B7SjcmoU4PW6-R4VMgMuyZd0hx2Q8v3lr1FY0ZEu3PHkhECh9$" rel="noreferrer" target="_blank">https://urldefense.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!AnSU1oWcLN8nr6iiJ0kauXvdiWCntoLiyc9wc8P9El-B7SjcmoU4PW6-R4VMgMuyZd0hx2Q8v3lr1FY0ZEu3PHkhECh9$</a> <br>
</blockquote></div><br clear="all"><div><br></div><span class="gmail_signature_prefix">-- </span><br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><pre cols="72">Dr. Manfredi Manizza
Geosciences Research Division
Scripps Institution of Oceanography
University of California San Diego
9500 Gilman Drive La Jolla, CA 92093-0244
email : <a href="mailto:mmanizza@ucsd.edu" target="_blank">mmanizza@ucsd.edu</a>
tel : +1-858-534-7094
web : <a href="http://bluemoon.ucsd.edu/mmanizza/" target="_blank">http://bluemoon.ucsd.edu/mmanizza/</a> </pre><pre cols="72">web : <a href="https://mmanizza.scrippsprofiles.ucsd.edu/" style="font-family:Arial,Helvetica,sans-serif" target="_blank">https://mmanizza.scrippsprofiles.ucsd.edu/</a></pre></div></div></div></div>
_______________________________________________<br>
MITgcm-support mailing list<br>
<a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
</blockquote></div>