1 |
C $Header: /u/gcmpack/MITgcm_contrib/rpa_layers/layers/layers_mnc_init.F,v 1.2 2009/09/16 18:04:49 jmc Exp $ |
2 |
C $Name: $ |
3 |
|
4 |
#include "LAYERS_OPTIONS.h" |
5 |
|
6 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
7 |
|
8 |
SUBROUTINE LAYERS_MNC_INIT( myThid ) |
9 |
|
10 |
C Initialize MNC lookup tables for LAYERS. |
11 |
|
12 |
IMPLICIT NONE |
13 |
#include "SIZE.h" |
14 |
#include "EEPARAMS.h" |
15 |
#include "PARAMS.h" |
16 |
#include "LAYERS_SIZE.h" |
17 |
#include "LAYERS.h" |
18 |
|
19 |
C !INPUT PARAMETERS: |
20 |
INTEGER myThid |
21 |
CEOP |
22 |
|
23 |
#ifdef ALLOW_LAYERS |
24 |
#ifdef ALLOW_MNC |
25 |
|
26 |
IF (useMNC) THEN |
27 |
|
28 |
C-- Need to figure out how to create new grids. We need |
29 |
C 'U_xy_Hn__G__t' |
30 |
C 'V_xy_Hn__G__t' |
31 |
|
32 |
CALL MNC_CW_ADD_VNAME('UFlux','U_xy_Hn__G__t', 4,5,myThid) |
33 |
CALL MNC_CW_ADD_VATTR_TEXT('UFlux','units','m^2/s', myThid) |
34 |
CALL MNC_CW_ADD_VATTR_TEXT('UFlux','description', |
35 |
& 'Eastward Velocity Integrated Over Isopycnal Layer', |
36 |
& myThid) |
37 |
|
38 |
CALL MNC_CW_ADD_VNAME('VFlux','V_xy_Hn__G__t', 4,5,myThid) |
39 |
CALL MNC_CW_ADD_VATTR_TEXT('VFlux','units','m^2/s', myThid) |
40 |
CALL MNC_CW_ADD_VATTR_TEXT('VFlux','description', |
41 |
& 'Northward Velocity Integrated Over Isopycnal Layer', |
42 |
& myThid) |
43 |
|
44 |
|
45 |
ENDIF |
46 |
|
47 |
#endif /* ALLOW_MNC */ |
48 |
#endif /* ALLOW_GMREDI */ |
49 |
|
50 |
RETURN |
51 |
END |
52 |
|
53 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |