Physics Template HOWTO

1. Monthly Data Fields Requested

1. Temperature and Salinity are requested so that the velocity fields from this submitted data does not need to be coordinated with T and S fields from particular OCMIP runs.
2. Heat and Freshwater fluxes are total combined fluxes. For instance, if the model includes climatological values and restoring, both terms should be included.
3. The adjectives Zonal and Meridional are NOT to be taken literally. Velocity and Wind Stress fields should be aligned with the model's local grid. The direction of the vectors is specified in the grid description.
4. Not all models will have eddy induced velocities.

2. Grids

• Eulerian Zonal Velocity
• Eulerian Meridional Velocity
• Eddy Induced Zonal Velocity
• Eddy Induced Meridional Velocity

For each grid, the following quantities are specified :

quantity	dimensions	description
imt	scalar	'Zonal' dimension
jmt	scalar	'Meridional' dimension
lon	(imt,jmt)	longitude of cell center
bounds_lon	(imt,jmt,2,2)	longitude boundaries of cell
lat	(imt,jmt)	latitude of cell center
bounds_lat	(imt,jmt,2,2)	latitude boundaries of cell
area	(imt,jmt)	surface area of cell
vector_angle	(imt,jmt)	Angle that the Local Vector makes with a Latitude Line

The depth axes for all of the velocity grids are assumed to coincide, so they are only specified once.

The bounds quantities follow the same conventions as the tracer grid cell bounds :

• bounds_lon(i,j,1,1) : lower left cell corner
• bounds_lon(i,j,2,1) : lower right cell corner
• bounds_lon(i,j,1,2) : upper left cell corner
• bounds_lon(i,j,2,2) : upper right cell corner

The quantity vector_angle is the counterclockwise angle that a vector on the grid makes with a latitude line. Examples are

Vector Direction	Corresponding Angle
Northeast	pi/4
East	0
South	-pi/2

3. Conventions

1. Longitudes must be expressed in degrees as eastward positive (e.g., 120W is to be expressed as -120; 120E is to be expressed as +120).
2. Latitudes must be expressed in degrees as northward positive. (e.g., 90S is to be expressed as -90; 90N is to be expressed as +90).
3. Depth must be expressed in meters as positive downward (e.g., the depth of 1000 m is to be expressed as 1000).
4. For irregular model grids which must be stored as a 1-D vector instead of a 2-D array (e.g., model AWI), you must set jmt=1, and not imt=1.
5. You must provide model output on a grid without overlapping boxes. For example, if for longitudes of your model
   • i = 1 is the same as i = imt-1 and
   • i = 2 is the same as i = imt,
   then you must only provide model output for values from i=2 to i=imt-1. Future reference to i during OCMIP-2 analysis will thus be shifted by one, relative to that used in your model.
6. As opposed to specifying a mask for each grid, values at interior land points should be reported as the missing_value -1.E34.

4. Template Assumptions

1. Surface Heat and Freshwater fluxes are on the tracer grid specified by write_nc_MaskAreaBathy.f.
2. All velocity fields, Eulerian and Eddy Induced, are on the same vertical grid.
3. Zonal and Meridional Wind Stresses are on the same horizontal grids as the Eulerian Zonal and Meridional Velocities respectively.

5. Downloading the output routines

• write_nc_phys.f
• write_nc_vel_grid.f
• nc_util.f

6. Transfer of output

If this method of transfer is not convient, please notify me.