[MITgcm-support] add stochastic forcing in surface fluxes

Martin Losch Martin.Losch at awi.de
Mon Jan 8 11:29:18 EST 2018


Hi Olivier,

I would try the end of model/src/external_forcing_surf.F where the fields surfaceForcingT and surfaceForcingS are computed from various input fields (including restoring). At the end you can add your stochastic forcing.

Alternatively you can modify the forcing fields (Qnet and EmPmR, different units that surfaceForcingT/S) at the beginning of the routine.

Martin

> On 26. Dec 2017, at 17:05, Olivier Arzel <Olivier.Arzel at univ-brest.fr> wrote:
> 
> Hello all
> 
> I would like to add stochastic forcing to prescribed surface fluxes.
> 
> The stochastic forcing is build here from white noise using random numbers picked up from a normal distribution and a specified autocorrelation timescale.
> 
> The easiest way to do it in the model would be to read the forcing files (which contain both the prescribed and stochastic parts) (for example surfQFile for the heat flux) every n timesteps (n small). Unfortunately, the size of these forcing files in this case is too big for simulations of typically > 100 years.
> 
> The other possibility would be to first read the prescribed part using 'surfQFile', then compute the stochastic part online, and finally add it to the prescribed flux. What is the routine I need to modify to do this, and has this been attempted before ?
> 
> Thank you,
> 
> Olivier
> 
> <Mail Attachment.png>
> http://stockage.univ-brest.fr/~oarzel/
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support



More information about the MITgcm-support mailing list