Back to home page

MITgcm

 
 

    


File indexing completed on 2021-08-12 05:12:27 UTC

view on githubraw file Latest commit 0320e252 on 2021-08-11 16:08:52 UTC
de859229b1 Ed H*0001 #include "SEAICE_OPTIONS.h"
                0002 
                0003 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0004 CBOP
                0005 C     !ROUTINE: SEAICE_MNC_INIT
0320e25227 Mart*0006 
de859229b1 Ed H*0007 C     !INTERFACE:
                0008       SUBROUTINE SEAICE_MNC_INIT( myThid )
0320e25227 Mart*0009 
de859229b1 Ed H*0010 C     !DESCRIPTION:
                0011 C     Initialize SEAICE variables within MNC
                0012 
                0013 C     !USES:
                0014       IMPLICIT NONE
                0015 #include "SIZE.h"
                0016 #include "EEPARAMS.h"
                0017 #include "PARAMS.h"
6e2f4e58fa Mart*0018 #ifdef SEAICE_ALLOW_EVP
7303eab4f2 Patr*0019 #include "SEAICE_SIZE.h"
6e2f4e58fa Mart*0020 #include "SEAICE_PARAMS.h"
                0021 #endif /* SEAICE_ALLOW_EVP */
de859229b1 Ed H*0022 
                0023 C     !INPUT PARAMETERS:
                0024       INTEGER myThid
                0025 CEOP
                0026 
                0027 #ifdef ALLOW_MNC
                0028       IF (useMNC) THEN
                0029 
6e2f4e58fa Mart*0030 #ifdef SEAICE_ALLOW_EVP
                0031        IF ( SEAICEuseEVP ) THEN
                0032         CALL MNC_CW_ADD_VNAME('si_sigma1','Cen_xy_Hn__-__t',3,4,myThid)
                0033         CALL MNC_CW_ADD_VATTR_TEXT('si_sigma1','units','N/m^2', myThid)
                0034         CALL MNC_CW_ADD_VATTR_TEXT('si_sigma1','description',
                0035      &       'divergence of stress tensor', myThid)
                0036         CALL MNC_CW_ADD_VNAME('si_sigma2','Cen_xy_Hn__-__t',3,4,myThid)
                0037         CALL MNC_CW_ADD_VATTR_TEXT('si_sigma2','units','N/m^2', myThid)
                0038         CALL MNC_CW_ADD_VATTR_TEXT('si_sigma2','description',
                0039      &       'tension of stress tensor', myThid)
                0040         CALL MNC_CW_ADD_VNAME('si_sigma12','Cor_xy_Hn__-__t',3,4,myThid)
                0041         CALL MNC_CW_ADD_VATTR_TEXT('si_sigma12','units','N/m^2', myThid)
                0042         CALL MNC_CW_ADD_VATTR_TEXT('si_sigma12','description',
                0043      &       'shear of stress tensor', myThid)
                0044        ENDIF
                0045 #endif /* SEAICE_ALLOW_EVP */
                0046 
b726014783 Ed H*0047         CALL MNC_CW_ADD_VNAME('si_UWIND','U_xy_Hn__-__t',3,4,myThid)
                0048         CALL MNC_CW_ADD_VATTR_TEXT('si_UWIND','units','m/s', myThid)
                0049         CALL MNC_CW_ADD_VATTR_TEXT('si_UWIND','description',
de859229b1 Ed H*0050      &       'Surface (10m) zonal wind velocity', myThid)
                0051 
b726014783 Ed H*0052         CALL MNC_CW_ADD_VNAME('si_VWIND','V_xy_Hn__-__t',3,4,myThid)
                0053         CALL MNC_CW_ADD_VATTR_TEXT('si_VWIND','units','m/s', myThid)
                0054         CALL MNC_CW_ADD_VATTR_TEXT('si_VWIND','description',
de859229b1 Ed H*0055      &       'Surface (10m) meridional wind velocity', myThid)
                0056 
b726014783 Ed H*0057         CALL MNC_CW_ADD_VNAME('si_EVAP','Cen_xy_Hn__-__t',3,4,myThid)
                0058         CALL MNC_CW_ADD_VATTR_TEXT('si_EVAP','units','m/s', myThid)
                0059         CALL MNC_CW_ADD_VATTR_TEXT('si_EVAP','description',
de859229b1 Ed H*0060      &       'Evaporation', myThid)
                0061 
b726014783 Ed H*0062         CALL MNC_CW_ADD_VNAME('si_PRECIP','Cen_xy_Hn__-__t',3,4,myThid)
                0063         CALL MNC_CW_ADD_VATTR_TEXT('si_PRECIP','units','m/s', myThid)
                0064         CALL MNC_CW_ADD_VATTR_TEXT('si_PRECIP','description',
de859229b1 Ed H*0065      &       'Precipitation', myThid)
                0066 
b726014783 Ed H*0067         CALL MNC_CW_ADD_VNAME('si_RUNOFF','Cen_xy_Hn__-__t',3,4,myThid)
                0068         CALL MNC_CW_ADD_VATTR_TEXT('si_RUNOFF','units','m/s', myThid)
                0069         CALL MNC_CW_ADD_VATTR_TEXT('si_RUNOFF','description',
de859229b1 Ed H*0070      &       'River and glacier runoff', myThid)
                0071 
b726014783 Ed H*0072         CALL MNC_CW_ADD_VNAME('si_ATEMP','Cen_xy_Hn__-__t',3,4,myThid)
                0073         CALL MNC_CW_ADD_VATTR_TEXT('si_ATEMP','units','deg K', myThid)
                0074         CALL MNC_CW_ADD_VATTR_TEXT('si_ATEMP','description',
de859229b1 Ed H*0075      &       'Surface (2m) air temperature', myThid)
                0076 
b726014783 Ed H*0077         CALL MNC_CW_ADD_VNAME('si_AQH','Cen_xy_Hn__-__t',3,4,myThid)
                0078         CALL MNC_CW_ADD_VATTR_TEXT('si_AQH','units','kg/kg', myThid)
                0079         CALL MNC_CW_ADD_VATTR_TEXT('si_AQH','description',
de859229b1 Ed H*0080      &       'Surface (2m) specific humidity', myThid)
                0081 
b726014783 Ed H*0082         CALL MNC_CW_ADD_VNAME('si_SWDOWN','Cen_xy_Hn__-__t',3,4,myThid)
                0083         CALL MNC_CW_ADD_VATTR_TEXT('si_SWDOWN','units','W/m^2', myThid)
                0084         CALL MNC_CW_ADD_VATTR_TEXT('si_SWDOWN','description',
de859229b1 Ed H*0085      &       'Downward shortwave radiation', myThid)
                0086 
b726014783 Ed H*0087         CALL MNC_CW_ADD_VNAME('si_LWDOWN','Cen_xy_Hn__-__t',3,4,myThid)
                0088         CALL MNC_CW_ADD_VATTR_TEXT('si_LWDOWN','units','W/m^2', myThid)
                0089         CALL MNC_CW_ADD_VATTR_TEXT('si_LWDOWN','description',
de859229b1 Ed H*0090      &       'Downward longwave radiation', myThid)
                0091 
ec1f837e9d Ed H*0092         CALL MNC_CW_ADD_VNAME('si_UICE','U_xy_Hn__-__t',3,4,myThid)
                0093         CALL MNC_CW_ADD_VATTR_TEXT('si_UICE','units','m/s', myThid)
                0094         CALL MNC_CW_ADD_VATTR_TEXT('si_UICE','description',
                0095      &       'Ice U velocity', myThid)
                0096 
                0097         CALL MNC_CW_ADD_VNAME('si_VICE','V_xy_Hn__-__t',3,4,myThid)
                0098         CALL MNC_CW_ADD_VATTR_TEXT('si_VICE','units','m/s', myThid)
                0099         CALL MNC_CW_ADD_VATTR_TEXT('si_VICE','description',
                0100      &       'Ice V velocity', myThid)
                0101 
                0102         CALL MNC_CW_ADD_VNAME('si_HEFF','Cen_xy_Hn__-__t',3,4,myThid)
                0103         CALL MNC_CW_ADD_VATTR_TEXT('si_HEFF','units','m', myThid)
                0104         CALL MNC_CW_ADD_VATTR_TEXT('si_HEFF','description',
                0105      &       'Effective ice thickness',myThid)
                0106 
                0107         CALL MNC_CW_ADD_VNAME('si_AREA','Cen_xy_Hn__-__t',3,4,myThid)
                0108         CALL MNC_CW_ADD_VATTR_TEXT('si_AREA','units','m^2/m^2',myThid)
                0109         CALL MNC_CW_ADD_VATTR_TEXT('si_AREA','description',
                0110      &       'Fractional ice-covered area',myThid)
                0111 
de859229b1 Ed H*0112 C       Time averages
b726014783 Ed H*0113         CALL MNC_CW_ADD_VNAME('si_UICEtave','U_xy_Hn__-__t',3,4,myThid)
                0114         CALL MNC_CW_ADD_VATTR_TEXT('si_UICEtave','units','m/s', myThid)
                0115         CALL MNC_CW_ADD_VATTR_TEXT('si_UICEtave','description',
de859229b1 Ed H*0116      &       'Time-averaged ', myThid)
                0117 
b726014783 Ed H*0118         CALL MNC_CW_ADD_VNAME('si_VICEtave','V_xy_Hn__-__t',3,4,myThid)
                0119         CALL MNC_CW_ADD_VATTR_TEXT('si_VICEtave','units','m/s', myThid)
                0120         CALL MNC_CW_ADD_VATTR_TEXT('si_VICEtave','description',
de859229b1 Ed H*0121      &       'Time-averaged ', myThid)
                0122 
b726014783 Ed H*0123         CALL MNC_CW_ADD_VNAME('si_FUtave','U_xy_Hn__-__t',3,4,myThid)
                0124         CALL MNC_CW_ADD_VATTR_TEXT('si_FUtave','units','N/m^2', myThid)
                0125         CALL MNC_CW_ADD_VATTR_TEXT('si_FUtave','description',
de859229b1 Ed H*0126      &       'Time-averaged zonal surface wind stress', myThid)
                0127 
b726014783 Ed H*0128         CALL MNC_CW_ADD_VNAME('si_FVtave','V_xy_Hn__-__t',3,4,myThid)
                0129         CALL MNC_CW_ADD_VATTR_TEXT('si_FVtave','units','N/m^2',myThid)
                0130         CALL MNC_CW_ADD_VATTR_TEXT('si_FVtave','description',
de859229b1 Ed H*0131      &       'Time-averaged Meridional surface wind stress', myThid)
                0132 
b726014783 Ed H*0133         CALL MNC_CW_ADD_VNAME(
                0134      &       'si_EmPmRtave','Cen_xy_Hn__-__t',3,4,myThid)
                0135         CALL MNC_CW_ADD_VATTR_TEXT(
e31dd4223a Mart*0136      &       'si_EmPmRtave','units','kg/m^2/s',myThid)
b726014783 Ed H*0137         CALL MNC_CW_ADD_VATTR_TEXT('si_EmPmRtave','description',
de859229b1 Ed H*0138      &       'Time-averaged net upward freshwater flux', myThid)
                0139 
b726014783 Ed H*0140         CALL MNC_CW_ADD_VNAME(
                0141      &       'si_QNETtave','Cen_xy_Hn__-__t',3,4,myThid)
                0142         CALL MNC_CW_ADD_VATTR_TEXT(
                0143      &       'si_QNETtave','units','W/m^2', myThid)
                0144         CALL MNC_CW_ADD_VATTR_TEXT('si_QNETtave','description',
de859229b1 Ed H*0145      &       'Time-averaged net upward surface heat flux', myThid)
                0146 
b726014783 Ed H*0147         CALL MNC_CW_ADD_VNAME(
                0148      &       'si_QSWtave','Cen_xy_Hn__-__t',3,4,myThid)
                0149         CALL MNC_CW_ADD_VATTR_TEXT(
                0150      &       'si_QSWtave','units','W/m^2',myThid)
                0151         CALL MNC_CW_ADD_VATTR_TEXT('si_QSWtave','description',
                0152      &       'Time-averaged net upward shortwave radiation',myThid)
                0153 
                0154         CALL MNC_CW_ADD_VNAME(
                0155      &       'si_HEFFtave','Cen_xy_Hn__-__t',3,4,myThid)
                0156         CALL MNC_CW_ADD_VATTR_TEXT(
                0157      &       'si_HEFFtave','units','m', myThid)
                0158         CALL MNC_CW_ADD_VATTR_TEXT('si_HEFFtave','description',
                0159      &       'Time-averaged effective ice thickness',myThid)
                0160 
                0161         CALL MNC_CW_ADD_VNAME(
                0162      &       'si_AREAtave','Cen_xy_Hn__-__t',3,4,myThid)
                0163         CALL MNC_CW_ADD_VATTR_TEXT(
                0164      &       'si_AREAtave','units','m^2/m^2',myThid)
                0165         CALL MNC_CW_ADD_VATTR_TEXT('si_AREAtave','description',
                0166      &       'Time-averaged fractional ice-covered area',myThid)
de859229b1 Ed H*0167 
                0168       ENDIF
                0169 #endif /*  ALLOW_MNC     */
                0170 
                0171       RETURN
                0172       END
                0173