[MITgcm-support] Wind stress without interpolation

Martin Losch Martin.Losch at awi.de
Mon Feb 4 03:41:25 EST 2013


Hi Massimo,

I don't think that there are run time flags for this. The only way I see how you can do it is by modifying model/src/external_fields_load.F
Please have a look at the file. In this file the forcing fields are read and interpolated, e.g., for wind stress:
            fu(i,j,bi,bj)    = bWght*taux0(i,j,bi,bj)
     &                       + aWght*taux1(i,j,bi,bj)
For "fu" to be constant over an "externForcingPeriod", you need to set either aWght or bWght to zero, or remove one of the terms in this sum.

Martin


On Feb 3, 2013, at 4:21 AM, massimo novello <massimonovello at yahoo.com.br> wrote:

> Please, is it possible to use a periodic wind force without interpolate the values? I need that the wind change suddently according to the files (wind_x.bin and wind_y.bin).
> 
> I started the simulation using the exp4 sample with some changes on its configuration (bathymetry, coast line, wind stress, open boundaries,...). 
> 
> Many thanks,
> Massimo Novello
> UERJ - Rio de Janeiro
> 
> From MITgcm manual:
> "There is also the possibility of prescribing time-dependent periodic forcing. To do this, concatenate the successive time records into a single file (for each stress component) ordered in a (x,y,t) fashion and set the following variables: periodicExternalForcing to ’.TRUE.’, externForcingPeriod to the period (in s) of which the forcing varies (typically 1 month), and externForcingCycle to the repeat time (in s) of the forcing (typically 1 year – note: externForcingCycle must be a multiple of externForcingPeriod). With these variables set up, the model will INTERPOLATE the forcing linearly at each iteration."
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support




More information about the MITgcm-support mailing list