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