Make a coupling file

Prepare nesting data using the output of Grid A. In /make_nest_file/, use conver.f to generate coupling.txt. The format of the coupling file is described below

coupling data

boundary info: num of points (negative means no point), start point

EAST

-1, 1 (no data on east side)

WEST

5, 1 (5 points at west boundary, start from 1)

SOUTH

-1, 0 ( no data on south side)

NORTH

-1, 0 (no data on north side)

TIME SERIES

68.443001 (first step, NOTE: model clock will be initialized as this time)

EAST SIDE (no nesting on east side)

WEST SIDE

0.141220E-02 0.141220E-02 0.141220E-02 0.141220E-02 0.141220E-02 (u)

-0.119260E-10 -0.119260E-10 -0.667390E-10 -0.667390E-10 -0.219240E-10 (v)

0.141890E-02 0.141890E-02 0.141890E-02 0.141890E-02 0.141890E-02 (z)

SOUTH SIDE (no nesting on south side)

NORTH SIDE (no nesting on north side)

68.641998 (next time step)

…

The example above is a case that a model nesting takes place at the WEST (left) boundary of a small domain. Boundaries are defined with the order: EAST, WEST, SOUTH, and NORTH. If the num of points of a boundary is larger than zero, the program will read a time series of (u, v, eta) below ‘XXXX SIDE’. The read format is

READ(11,*) TIME_COUPLING_2

READ(11,119)(U_COUPLING_EAST(I,2),I=1,N_COUPLING_EAST)

READ(11,119)(V_COUPLING_EAST(I,2),I=1,N_COUPLING_EAST)

READ(11,119)(Z_COUPLING_EAST(I,2),I=1,N_COUPLING_EAST)

ENDDO

119 FORMAT(5E16.6)

where N_COUPLING_EAST is Num of points at the EAST boundary.