[MITgcm-support] Bulk formula

[Martin Losch]

Hi Oceane,

did you check your radiative and sensibe fluxes, and surface wind stress that the bulk formulae computes? (use the diagnostics package to do so) Are they “reasonable” and comparable to the heat flux implied by restoring?

You SSS increases, this may indicate too much evaporation, i.e. too much latent heat loss -> too much cooling? (just speculating)

It would help if you provided your data.exf, data, EXF_OPTIONS.h, CPP_OPTIONS.h along with your question. Also (or alternatively) the STDOUT.0000 would help (to see which parameter values are actually used).

Martin
> On 24. Jan 2020, at 05:45, Richet, Oceane (O&A, Hobart) <Oceane.Richet at csiro.au> wrote:
> 
> MITgcmers,
> 
> Previously I was running my regional simulations with a restoring temperature and salinity at the surface and at the lateral boundaries. I recently moved to a configuration where now I am using the bulk formula (BULK_FORCE package) at the surface for a more natural surface boundary condition without changing the forcing at the lateral boundaries. Both simulations are forced with a seasonal forcing.
> 
> In the case of the restoring forcing at the surface, I used the monthly SST and SSS from the COSIMA modelhttp://cosima.org.au/ (same for the lateral boundaries). For the bulk formula, I use the JRA55 outputs which have been used to force the COSIMA model. Now, the forcing instead to be monthly inputs is 3-hourly inputs - the forcing at the lateral boundaries is still the same as previously but now it has been linearly interpolated in time 3-hourly.
just a comment, not important: You don’t need to do this interpolation yourself. With exf/obcs the monthly values are interpolately linearly in time anyway. It’s enough to provide the original monthly obcs data.
> 
> What I observe between these 2 configurations after 12 years of running is that the SST is much colder and does not match the lateral boundary forcing (see figure SST12 - bulk formula). This behavior does not seem to be observe for the SSS (figure SSS_12 - bulk formula). So I am wondering if it is physical or if there is an issue with the bulk formula like for example a units problem as, for the bulk formula, my input file for the temperature is in Kevin and in Celsius for the lateral open boundaries. Also the wind stress is now calculated through the bulk formula from the U and V wind fields.
units of temperature are OK.
> 
> Has someone already observed this discontinuity between the surface and the lateral boundaries?
> 
> Cheers,
> 
> Océane Richet
> <itf111_SSS.png><itf111_SST.png><itf111b_SSS.png><itf111b_SST.png><SSTxy_cosima.png>_______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support