Changeset 7536
- Timestamp:
- 2017-01-10T10:59:45+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/v3_6_CMIP6_ice_diagnostics/NEMOGCM/NEMO/LIM_SRC_3/limwri.F90
r7517 r7536 233 233 zrho1 = ( rau0 - rhoic ) / rau0; zrho2 = rhosn / rau0 234 234 235 CALL iom_put( "icethic" , htm_i * zswi ) ! ice thickness (SIMIP 2016)236 CALL iom_put( "icepres" , zswi ) ! ice presence (1 or 0) (SIMIP 2016)237 CALL iom_put( "snowthic" , htm_s * zswi ) ! snow thickness (SIMIP 2016)238 CALL iom_put( "icemass" , rhoic * vt_i(:,:) * zswi ) ! ice mass per cell area (SIMIP 2016)239 CALL iom_put( "snomass" , rhosn * vt_s(:,:) * zswi ) ! snow mass per cell area (SIMIP 2016)240 CALL iom_put( "icesnt" , ( tm_si - rt0 ) * zswi ) ! snow-ice interface temperature (SIMIP 2016)241 CALL iom_put( "icebot" , ( t_bo - rt0 ) * zswi ) ! ice bottom temperature (SIMIP 2016)242 CALL iom_put( "icesmass" , smt_i * vt_i * rhoic / 1000. * zswi ) ! mass of salt in sea ice per cell area (SIMIP 2016)243 CALL iom_put( "icefb" , ( zrho1 * htm_i(:,:) - zrho2 * htm_s(:,:) ) * zswi ) ! mass of salt in sea ice per cell area (SIMIP 2016)244 245 CALL iom_put( "wfxsum" , wfx_sum ) ! Freshwater flux from sea-ice surface246 CALL iom_put( "dmithd" , - wfx_bog - wfx_bom - wfx_sum & ! Sea-ice mass change from thermodynamics235 IF ( iom_use( "icethic" ) ) CALL iom_put( "icethic" , htm_i * zswi ) ! ice thickness 236 IF ( iom_use( "icepres" ) ) CALL iom_put( "icepres" , zswi ) ! ice presence (1 or 0) 237 IF ( iom_use( "snowthic" ) ) CALL iom_put( "snowthic" , htm_s * zswi ) ! snow thickness 238 IF ( iom_use( "icemass" ) ) CALL iom_put( "icemass" , rhoic * vt_i(:,:) * zswi ) ! ice mass per cell area 239 IF ( iom_use( "snomass" ) ) CALL iom_put( "snomass" , rhosn * vt_s(:,:) * zswi ) ! snow mass per cell area 240 IF ( iom_use( "icesnt" ) ) CALL iom_put( "icesnt" , ( tm_si - rt0 ) * zswi ) ! snow-ice interface temperature 241 IF ( iom_use( "icebot" ) ) CALL iom_put( "icebot" , ( t_bo - rt0 ) * zswi ) ! ice bottom temperature 242 IF ( iom_use( "icesmass" ) ) CALL iom_put( "icesmass" , smt_i * vt_i * rhoic / 1000. * zswi ) ! mass of salt in sea ice per cell area 243 IF ( iom_use( "icefb" ) ) CALL iom_put( "icefb" , ( zrho1 * htm_i(:,:) - zrho2 * htm_s(:,:) ) * zswi ) ! mass of salt in sea ice per cell area 244 245 IF ( iom_use( "wfxsum" ) ) CALL iom_put( "wfxsum" , wfx_sum ) ! Freshwater flux from sea-ice surface 246 IF ( iom_use( "dmithd" ) ) CALL iom_put( "dmithd" , - wfx_bog - wfx_bom - wfx_sum & ! Sea-ice mass change from thermodynamics 247 247 & - wfx_sni - wfx_opw - wfx_res ) 248 CALL iom_put( "dmidyn" , diag_dmi_dyn ) ! Sea-ice mass change from dynamics249 CALL iom_put( "dmiopw" , - wfx_opw ) ! Sea-ice mass change through growth in open water250 CALL iom_put( "dmibog" , - wfx_bog ) ! Sea-ice mass change through basal growth251 CALL iom_put( "dmisni" , - wfx_sni ) ! Sea-ice mass change through snow-to-ice conversion252 CALL iom_put( "dmisum" , - wfx_sum ) ! Sea-ice mass change through surface melting253 CALL iom_put( "dmibom" , - wfx_bom ) ! Sea-ice mass change through bottom melting254 CALL iom_put( "dmtsub" , - wfx_sub ) ! Sea-ice mass change through evaporation and sublimation255 CALL iom_put( "dmsspr" , - wfx_spr ) ! snow mass change through snow fall256 CALL iom_put( "dmsssi" , wfx_sni*rhosn/rhoic ) ! snow mass change through snow-to-ice conversion257 258 CALL iom_put( "dmsmel" , diag_dms_mel ) ! snow mass change through melt259 CALL iom_put( "dmsdyn" , diag_dms_dyn ) ! snow mass change through dynamics260 261 CALL iom_put( "hfxconbo" , diag_fc_bo ) ! bottom conduction flux262 CALL iom_put( "hfxconsu" , diag_fc_su ) ! surface conduction flux263 264 CALL iom_put( "wfxtot" , -wfx_ice ) ! total freshwater flux from sea ice265 266 CALL iom_put( "dmtxdyn" , diag_dmtx_dyn ) ! X-component of sea-ice mass transport267 CALL iom_put( "dmtydyn" , diag_dmty_dyn ) ! Y-component of sea-ice mass transport268 269 CALL iom_put( "utau_oce" , diag_utau_oi*zswi ) ! X-component of ocean stress on sea ice270 CALL iom_put( "vtau_oce" , diag_vtau_oi*zswi ) ! Y-component of ocean stress on sea ice271 272 CALL iom_put( "dssh_dx" , diag_dssh_dx*zswi ) ! Sea-surface tilt term in force balance (x-component)273 CALL iom_put( "dssh_dy" , diag_dssh_dy*zswi ) ! Sea-surface tilt term in force balance (y-component)274 275 CALL iom_put( "corstrx" , diag_corstrx*zswi ) ! Coriolis force term in force balance (x-component)276 CALL iom_put( "corstry" , diag_corstry*zswi ) ! Coriolis force term in force balance (y-component)277 278 CALL iom_put( "intstrx" , diag_intstrx*zswi ) ! Internal force term in force balance (x-component)279 CALL iom_put( "intstry" , diag_intstry*zswi ) ! Internal force term in force balance (y-component)280 281 CALL iom_put( "normstr" , diag_sig1 *zswi ) ! Normal stress282 CALL iom_put( "sheastr" , diag_sig2 *zswi ) ! Shear stress248 IF ( iom_use( "dmidyn" ) ) CALL iom_put( "dmidyn" , diag_dmi_dyn ) ! Sea-ice mass change from dynamics 249 IF ( iom_use( "dmiopw" ) ) CALL iom_put( "dmiopw" , - wfx_opw ) ! Sea-ice mass change through growth in open water 250 IF ( iom_use( "dmibog" ) ) CALL iom_put( "dmibog" , - wfx_bog ) ! Sea-ice mass change through basal growth 251 IF ( iom_use( "dmisni" ) ) CALL iom_put( "dmisni" , - wfx_sni ) ! Sea-ice mass change through snow-to-ice conversion 252 IF ( iom_use( "dmisum" ) ) CALL iom_put( "dmisum" , - wfx_sum ) ! Sea-ice mass change through surface melting 253 IF ( iom_use( "dmibom" ) ) CALL iom_put( "dmibom" , - wfx_bom ) ! Sea-ice mass change through bottom melting 254 IF ( iom_use( "dmtsub" ) ) CALL iom_put( "dmtsub" , - wfx_sub ) ! Sea-ice mass change through evaporation and sublimation 255 IF ( iom_use( "dmsspr" ) ) CALL iom_put( "dmsspr" , - wfx_spr ) ! snow mass change through snow fall 256 IF ( iom_use( "dmsssi" ) ) CALL iom_put( "dmsssi" , wfx_sni*rhosn/rhoic ) ! snow mass change through snow-to-ice conversion 257 258 IF ( iom_use( "dmsmel" ) ) CALL iom_put( "dmsmel" , diag_dms_mel ) ! snow mass change through melt 259 IF ( iom_use( "dmsdyn" ) ) CALL iom_put( "dmsdyn" , diag_dms_dyn ) ! snow mass change through dynamics 260 261 IF ( iom_use( "hfxconbo" ) ) CALL iom_put( "hfxconbo" , diag_fc_bo ) ! bottom conduction flux 262 IF ( iom_use( "hfxconsu" ) ) CALL iom_put( "hfxconsu" , diag_fc_su ) ! surface conduction flux 263 264 IF ( iom_use( "wfxtot" ) ) CALL iom_put( "wfxtot" , -wfx_ice ) ! total freshwater flux from sea ice 265 266 IF ( iom_use( "dmtxdyn" ) ) CALL iom_put( "dmtxdyn" , diag_dmtx_dyn ) ! X-component of sea-ice mass transport 267 IF ( iom_use( "dmtydyn" ) ) CALL iom_put( "dmtydyn" , diag_dmty_dyn ) ! Y-component of sea-ice mass transport 268 269 IF ( iom_use( "utau_oce" ) ) CALL iom_put( "utau_oce" , diag_utau_oi*zswi ) ! X-component of ocean stress on sea ice 270 IF ( iom_use( "vtau_oce" ) ) CALL iom_put( "vtau_oce" , diag_vtau_oi*zswi ) ! Y-component of ocean stress on sea ice 271 272 IF ( iom_use( "dssh_dx" ) ) CALL iom_put( "dssh_dx" , diag_dssh_dx*zswi ) ! Sea-surface tilt term in force balance (x-component) 273 IF ( iom_use( "dssh_dy" ) ) CALL iom_put( "dssh_dy" , diag_dssh_dy*zswi ) ! Sea-surface tilt term in force balance (y-component) 274 275 IF ( iom_use( "corstrx" ) ) CALL iom_put( "corstrx" , diag_corstrx*zswi ) ! Coriolis force term in force balance (x-component) 276 IF ( iom_use( "corstry" ) ) CALL iom_put( "corstry" , diag_corstry*zswi ) ! Coriolis force term in force balance (y-component) 277 278 IF ( iom_use( "intstrx" ) ) CALL iom_put( "intstrx" , diag_intstrx*zswi ) ! Internal force term in force balance (x-component) 279 IF ( iom_use( "intstry" ) ) CALL iom_put( "intstry" , diag_intstry*zswi ) ! Internal force term in force balance (y-component) 280 281 IF ( iom_use( "normstr" ) ) CALL iom_put( "normstr" , diag_sig1 *zswi ) ! Normal stress 282 IF ( iom_use( "sheastr" ) ) CALL iom_put( "sheastr" , diag_sig2 *zswi ) ! Shear stress 283 283 284 284 !--------------------------------
Note: See TracChangeset
for help on using the changeset viewer.