Hello,

I am working with OpenIFS CY43R3 in its aqua planet configuration and trying to add a stationary localized heat perturbation which would mimick convective heating in the tropical mid-troposphere.

My intention is to get an ensemble of short-time simulations (order of few days) to analyse the short-time, near-field response.

I have a few ideas how to implement this, but I am unsure which one would be the easiest. First of all, I was thinking of nudging a temperature field which would contain this localized heating. 

Secondly, I believe it is also possible to add an additional temperature tendency to the righthand side of the temperature equation, however this seems quite technical.

Lastly, the upper air initialization file could also be formed in such way that it contains a heat source, however it should be there throughout the whole simulation not only in the initial time.

Which way would be the easiest to proceed?



Best Regards,

Sándor Mahó

4 Comments

  1. It would be great if some OpenIFS users who have actively used the aquaplanet idealised configuration would be able to assist Sándor with his query?  It is well possible that the confluence page on the aquaplanet configuration needs updating.  Any input that we can receive would be very welcome!

    Best regards,  Marcus

  2. Hi Sándor,

    If your intention is to change temperature tendencies in the model then this could be fairly easily achieved within the subroutine callpar.F90.  Within the ECMWF training course programme we have computer practicals with OpenIFS which attempt to achieve exactly that – simulating a local modification of temperature tendencies from physical parameterizations through applying a constant scaling factor to a regional selection of gridboxes. The simplicity of it is that the change is fairly self-contained within callpar.  I believe that changing forcing files used for nudging is much more work intensive.

    Cheers,

    Marcus

  3. Hi Marcus,

    Thank you! I have found the description of NAWDEX perturbation experiments of the OpenIFS workshop from 2019. I think you may have referred to this computer practical. The described approach there is similar to that I intend to achieve. Is the code containing the modifications in callpar.F90 available somewhere?


    Cheers,

    Sándor

    1. Hi Sándor – I will email you directly.

      Marcus