I am new to ROMS, and would like to know whether I am setting the time dimension of the forcing variables correctly. sustr and svstr are given at a single time, and I want them to be applied at every time step. It appears that the routine get_2dfld.F is reading the forcing file properly, but I found the "

**Index:**" part of the text output (see below) suspicious. If sustr and svstr have shape (1,eta_u,xi_u) and (1,eta_v,xi_v), how can get_2dfld.F read two indices? I tried setting the "cycle_length" attribute in the sms_time variable to 0., but then I get the following error from routine get_cycle.F: GET_CYCLE - unable to get value for attribute: cycle_length in variable: sms_time This attribute value is expected to be of the same external type as the variable. What would be the best way to specify this stationary forcing?

#======================================

GET_2DFLD - surface v-momentum stress, t = 0 00:00:00

(Rec=0000001,

**Index=1**, File: msc1_frc.nc)

(Tmin= 0.0000 Tmax= 0.0000)

(Min = -1.23370321E-05 Max = -5.10689844E-06)

297 1 00:45:00 3.560151E-02 1.359857E+04 1.359861E+04 2.012812E+14

(158,071,01) 8.525388E-02 9.717515E-02 1.090288E+01 3.073193E+00

GET_2DFLD - surface u-momentum stress, t = 0 00:00:00

(Rec=0000001,

**Index=2**, File: msc1_frc.nc)

(Tmin= 0.0000 Tmax= 0.0000)

(Min = -1.99545517E-05 Max = -8.68067430E-06)

GET_2DFLD - surface v-momentum stress, t = 0 00:00:00

(Rec=0000001,

**Index=2**, File: msc1_frc.nc)

(Tmin= 0.0000 Tmax= 0.0000)

(Min = -1.23370321E-05 Max = -5.10689844E-06)

298 1 00:50:00 3.566264E-02 1.359856E+04 1.359859E+04 2.012812E+14

(158,071,01) 8.391898E-02 9.654408E-02 1.073583E+01 3.160133E+00

GET_2DFLD - surface u-momentum stress, t = 0 00:00:00

(Rec=0000001,

**Index=1**, File: msc1_frc.nc)

(Tmin= 0.0000 Tmax= 0.0000)

(Min = -1.99545517E-05 Max = -8.68067430E-06)

#======================================

Here is the output of ncdump -h from the forcing file:

dimensions:

xi_u = 160 ;

eta_u = 130 ;

xi_v = 161 ;

eta_v = 129 ;

**sms_time = 1**;

variables:

double sms_time(sms_time) ;

sms_time:field = "surface momentum stress time" ;

sms_time:units = "day" ;

sms_time:

**cycle_length = 1.**;

double sustr(sms_time, eta_u, xi_u) ;

sustr:field = "surface u-momentum stress, scalar, series" ;

sustr:units = "Newton meter-2" ;

sustr:time = "sms_time" ;

double svstr(sms_time, eta_v, xi_v) ;

svstr:field = "surface v-momentum stress, scalar, series" ;

svstr:units = "Newton meter-2" ;

svstr:time = "sms_time" ;

// global attributes:

:type = "ROMS forcing file" ;

:title = "First test, rectangular non-rotated grid." ;

:out_file = "msc1_frc.nc" ;

:grd_file = "msc1_grd.nc" ;

:history = "2014-04-02 17:28:03.206938" ;

}

Thank you for your time.