Hello!
I'm attempting to conduct a simulation which uses climatological (100 year averaged) prescribed wind-stress forcing from the piControl run, but also imposes a wind stress anomaly (from a regression of wind stress on the Nino 3.4 index in the piControl run) scaled by the magnitude of the Nino 3.4 index at each integration. I see that from the tavg.F90 file in the bld/ocn/source directory can output Nino indices, but I'm wondering if it's possible to call these indices into the forcing_ws.f90 file to perform some arithmetic on the prescribed wind stress data that changes at each time step.
I have daily climatological wind stress computed from a piControl simulation, the Nino34 anomaly pattern for wind stress - I just need to scale the anomaly by time varying quantities of SST during this proposed simulation.
As a function:
daily_climatological_wind_stress (x, y, t) + wind_stress_nino34_regression_pattern(x,y) * nino34(t) = prescribed_wind_stress
Thanks in advance for any help ~
I'm attempting to conduct a simulation which uses climatological (100 year averaged) prescribed wind-stress forcing from the piControl run, but also imposes a wind stress anomaly (from a regression of wind stress on the Nino 3.4 index in the piControl run) scaled by the magnitude of the Nino 3.4 index at each integration. I see that from the tavg.F90 file in the bld/ocn/source directory can output Nino indices, but I'm wondering if it's possible to call these indices into the forcing_ws.f90 file to perform some arithmetic on the prescribed wind stress data that changes at each time step.
I have daily climatological wind stress computed from a piControl simulation, the Nino34 anomaly pattern for wind stress - I just need to scale the anomaly by time varying quantities of SST during this proposed simulation.
As a function:
daily_climatological_wind_stress (x, y, t) + wind_stress_nino34_regression_pattern(x,y) * nino34(t) = prescribed_wind_stress
Thanks in advance for any help ~