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
0004
0005
0320e25227 Mart*0006
de859229b1 Ed H*0007
0008 SUBROUTINE SEAICE_MNC_INIT( myThid )
0320e25227 Mart*0009
de859229b1 Ed H*0010
0011
0012
0013
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
0024 INTEGER myThid
0025
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
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