pkg/gmredi/gmredi_mnc_init.F

3e63c76015 Ed H*0001 #include "GMREDI_OPTIONS.h"
                0002
                0003 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0004 CBOP 0
                0005 C     !ROUTINE: GMREDI_MNC_INIT
                0006
                0007 C     !INTERFACE:
                0008       SUBROUTINE GMREDI_MNC_INIT( myThid )
                0009
                0010 C     !DESCRIPTION:
                0011 C     Initialize MNC lookup tables for GMRedi.
                0012
                0013 C     !USES:
                0014       IMPLICIT NONE
                0015 #include "SIZE.h"
                0016 #include "EEPARAMS.h"
                0017 #include "PARAMS.h"
                0018 #include "GMREDI.h"
                0019
                0020 C     !INPUT PARAMETERS:
                0021       INTEGER myThid
                0022 CEOP
                0023
                0024 #ifdef ALLOW_GMREDI
                0025 #ifdef ALLOW_MNC
                0026
                0027       IF (useMNC) THEN
                0028
                0029         CALL MNC_CW_ADD_VNAME('Kux','U_xy_Hn__C__t', 4,5,myThid)
                0030         CALL MNC_CW_ADD_VATTR_TEXT('Kux','units','m^2/s', myThid)
                0031         CALL MNC_CW_ADD_VATTR_TEXT('Kux','description',
                0032      &       'GM Tensor Kux is K_11 element, X direction at U point',
                0033      &       myThid)
                0034 C        CALL MNC_CW_ADD_VATTR_TEXT('Kux',
                0035 C    &       'coordinates','XC YC RC iter', myThid)
                0036
                0037         CALL MNC_CW_ADD_VNAME('Kvy','V_xy_Hn__C__t', 4,5,myThid)
                0038         CALL MNC_CW_ADD_VATTR_TEXT('Kvy','units','m^2/s', myThid)
                0039         CALL MNC_CW_ADD_VATTR_TEXT('Kvy','description',
                0040      &       'GM Tensor Kvy is K_22 element, Y direction at V point',
                0041      &       myThid)
                0042
                0043         CALL MNC_CW_ADD_VNAME('Kuz','U_xy_Hn__C__t', 4,5,myThid)
                0044         CALL MNC_CW_ADD_VATTR_TEXT('Kuz','units','m^2/s', myThid)
                0045         CALL MNC_CW_ADD_VATTR_TEXT('Kuz','description',
                0046      &       'GM Tensor Kuz is K_13 element, Z direction at U point',
                0047      &       myThid)
                0048
                0049         CALL MNC_CW_ADD_VNAME('Kvz','V_xy_Hn__C__t', 4,5,myThid)
                0050         CALL MNC_CW_ADD_VATTR_TEXT('Kvz','units','m^2/s', myThid)
                0051         CALL MNC_CW_ADD_VATTR_TEXT('Kvz','description',
                0052      &       'GM Tensor Kvz is K_23 element, Z direction at V point',
                0053      &       myThid)
                0054
                0055         CALL MNC_CW_ADD_VNAME('Kwx','Cen_xy_Hn__L__t', 4,5,myThid)
                0056         CALL MNC_CW_ADD_VATTR_TEXT('Kwx','units','m^2/s', myThid)
                0057         CALL MNC_CW_ADD_VATTR_TEXT('Kwx','description',
                0058      &       'GM Tensor Kwx is K_31 element, X direction at W point',
                0059      &       myThid)
                0060
                0061         CALL MNC_CW_ADD_VNAME('Kwy','Cen_xy_Hn__L__t', 4,5,myThid)
                0062         CALL MNC_CW_ADD_VATTR_TEXT('Kwy','units','m^2/s', myThid)
                0063         CALL MNC_CW_ADD_VATTR_TEXT('Kwy','description',
                0064      &       'GM Tensor Kwy is K_32 element, Y direction at W point',
                0065      &       myThid)
                0066
                0067         CALL MNC_CW_ADD_VNAME('Kwz','Cen_xy_Hn__L__t', 4,5,myThid)
                0068         CALL MNC_CW_ADD_VATTR_TEXT('Kwz','units','m^2/s', myThid)
                0069         CALL MNC_CW_ADD_VATTR_TEXT('Kwz','description',
                0070      &       'GM Tensor Kwz is K_33 element, Z direction at W point',
                0071      &       myThid)
                0072
7a34c45a18 Mart*0073         CALL MNC_CW_ADD_VNAME('VisbK','Cen_xy_Hn__-__t',3,4,myThid)
3e63c76015 Ed H*0074         CALL MNC_CW_ADD_VATTR_TEXT('VisbK','units','m^2/s', myThid)
                0075         CALL MNC_CW_ADD_VATTR_TEXT('VisbK','description',
                0076      &       'GM mixing/stirring coef after Visbeck et al.',
                0077      &       myThid)
                0078
                0079         CALL MNC_CW_ADD_VNAME('PsiX','Cen_xy_Hn__C__t',4,5,myThid)
                0080         CALL MNC_CW_ADD_VATTR_TEXT('PsiX','units','m^2/s', myThid)
                0081         CALL MNC_CW_ADD_VATTR_TEXT('PsiX','description',
                0082      &       'GM bolus velocity streamfunction in X',
                0083      &       myThid)
                0084
                0085         CALL MNC_CW_ADD_VNAME('PsiY','Cen_xy_Hn__C__t',4,5,myThid)
                0086         CALL MNC_CW_ADD_VATTR_TEXT('PsiY','units','m^2/s', myThid)
                0087         CALL MNC_CW_ADD_VATTR_TEXT('PsiY','description',
                0088      &       'GM bolus velocity streamfunction in Y',
                0089      &       myThid)
                0090
                0091       ENDIF
                0092
                0093 #endif /* ALLOW_MNC */
                0094 #endif /* ALLOW_GMREDI */
                0095
                0096       RETURN
                0097       END
                0098
                0099 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|