- Timestamp:
- 2017-09-27T11:11:01+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icewri.F90
r8563 r8564 121 121 CALL iom_put( "iceconc" , at_i * zswi ) ! ice concentration 122 122 CALL iom_put( "icevolu" , vt_i * zswi ) ! ice volume = mean ice thickness over the cell 123 CALL iom_put( "icethick" , h tm_i* zswi ) ! ice thickness123 CALL iom_put( "icethick" , hm_i * zswi ) ! ice thickness 124 124 CALL iom_put( "icehc" , et_i * zswi ) ! ice total heat content 125 125 CALL iom_put( "isnowhc" , et_s * zswi ) ! snow total heat content 126 126 CALL iom_put( "ibrinv" , bvm_i * zswi * 100. ) ! brine volume 127 CALL iom_put( "micesalt" , sm t_i* zswi ) ! mean ice salinity128 CALL iom_put( "snowvol" , vt_s 127 CALL iom_put( "micesalt" , sm_i * zswi ) ! mean ice salinity 128 CALL iom_put( "snowvol" , vt_s * zswi ) ! snow volume 129 129 130 130 IF ( ln_pnd ) THEN … … 136 136 ! Output category-dependent fields 137 137 !---------------------------------- 138 IF ( iom_use('iceconc_cat' ) ) CALL iom_put( "iceconc_cat" , a_i 138 IF ( iom_use('iceconc_cat' ) ) CALL iom_put( "iceconc_cat" , a_i * zswi2 ) ! area for categories 139 139 IF ( iom_use('icethic_cat' ) ) CALL iom_put( "icethic_cat" , h_i * zswi2 ) ! thickness for categories 140 140 IF ( iom_use('snowthic_cat') ) CALL iom_put( "snowthic_cat" , h_s * zswi2 ) ! snow depth for categories 141 IF ( iom_use('salinity_cat') ) CALL iom_put( "salinity_cat" , s m_i * zswi2 ) ! salinity for categories141 IF ( iom_use('salinity_cat') ) CALL iom_put( "salinity_cat" , s_i * zswi2 ) ! salinity for categories 142 142 IF ( iom_use('icetemp_cat' ) ) CALL iom_put( "icetemp_cat", ( SUM( t_i(:,:,:,:), dim=3 ) * r1_nlay_i - rt0 ) * zswi2 ) ! ice temperature 143 143 IF ( iom_use('snwtemp_cat' ) ) CALL iom_put( "snwtemp_cat", ( SUM( t_s(:,:,:,:), dim=3 ) * r1_nlay_s - rt0 ) * zswi2 ) ! snow temperature … … 159 159 IF ( iom_use( "icepres" ) ) CALL iom_put( "icepres" , zswi(:,:) ) ! Ice presence (1 or 0) 160 160 IF ( iom_use( "icemass" ) ) CALL iom_put( "icemass" , rhoic * vt_i(:,:) * zswi(:,:) ) ! Ice mass per cell area 161 IF ( iom_use( "icethic" ) ) CALL iom_put( "icethic" , h tm_i(:,:) * zamask(:,:) + ( 1. - zamask(:,:) ) * zmiss_val ) ! Ice thickness161 IF ( iom_use( "icethic" ) ) CALL iom_put( "icethic" , hm_i(:,:) * zamask(:,:) + ( 1. - zamask(:,:) ) * zmiss_val ) ! Ice thickness 162 162 IF ( iom_use( "snomass" ) ) CALL iom_put( "snomass" , rhosn * vt_s(:,:) * zswi(:,:) + zmiss(:,:) ) ! Snow mass per cell area 163 IF ( iom_use( "snothic" ) ) CALL iom_put( "snothic" , h tm_s(:,:) * zamask(:,:) + ( 1. - zamask(:,:) ) * zmiss_val ) ! Snow thickness163 IF ( iom_use( "snothic" ) ) CALL iom_put( "snothic" , hm_s(:,:) * zamask(:,:) + ( 1. - zamask(:,:) ) * zmiss_val ) ! Snow thickness 164 164 165 165 IF ( iom_use( "iceconc_cat_mv" ) ) CALL iom_put( "iceconc_cat_mv" , a_i(:,:,:) * zswi2(:,:,:) + zmiss2(:,:,:) ) ! Area for categories … … 171 171 IF ( iom_use( "icebotK" ) ) CALL iom_put( "icebotK" , t_bo(:,:) * zswi(:,:) + zmiss(:,:) ) ! Ice bottom temperature 172 172 IF ( iom_use( "iceage" ) ) CALL iom_put( "iceage" , om_i(:,:) * zamask15(:,:) + ( 1. - zamask15(:,:) ) * zmiss_val ) ! Ice age 173 IF ( iom_use( "icesmass" ) ) CALL iom_put( "icesmass" , SUM( s mv_i, DIM = 3 ) * rhoic * 1.0e-3 * zswi(:,:) ) ! Mass of salt in sea ice per cell area174 IF ( iom_use( "icesal" ) ) CALL iom_put( "icesal" , sm t_i(:,:) * zswi(:,:) + zmiss(:,:) ) ! Ice salinity173 IF ( iom_use( "icesmass" ) ) CALL iom_put( "icesmass" , SUM( sv_i, DIM = 3 ) * rhoic * 1.0e-3 * zswi(:,:) ) ! Mass of salt in sea ice per cell area 174 IF ( iom_use( "icesal" ) ) CALL iom_put( "icesal" , sm_i(:,:) * zswi(:,:) + zmiss(:,:) ) ! Ice salinity 175 175 176 176 IF ( iom_use( "icefb" ) ) THEN 177 zfb(:,:) = ( zrho1 * h tm_i(:,:) - zrho2 * htm_s(:,:) )177 zfb(:,:) = ( zrho1 * hm_i(:,:) - zrho2 * hm_s(:,:) ) 178 178 WHERE( zfb < 0._wp ) ; zfb = 0._wp ; END WHERE 179 179 CALL iom_put( "icefb" , zfb(:,:) * zswi(:,:) + zmiss(:,:) ) ! Ice freeboard … … 356 356 CALL histend( kid, snc4set ) ! end of the file definition 357 357 358 CALL histwrite( kid, "sithic", kt, h tm_i, jpi*jpj, (/1/) )359 CALL histwrite( kid, "siconc", kt, at_i , jpi*jpj, (/1/) )360 CALL histwrite( kid, "sitemp", kt, tm_i - rt0 , jpi*jpj, (/1/) )358 CALL histwrite( kid, "sithic", kt, hm_i , jpi*jpj, (/1/) ) 359 CALL histwrite( kid, "siconc", kt, at_i , jpi*jpj, (/1/) ) 360 CALL histwrite( kid, "sitemp", kt, tm_i - rt0 , jpi*jpj, (/1/) ) 361 361 CALL histwrite( kid, "sivelu", kt, u_ice , jpi*jpj, (/1/) ) 362 362 CALL histwrite( kid, "sivelv", kt, v_ice , jpi*jpj, (/1/) ) … … 366 366 CALL histwrite( kid, "sinflx", kt, qns , jpi*jpj, (/1/) ) 367 367 CALL histwrite( kid, "isnowpre", kt, sprecip , jpi*jpj, (/1/) ) 368 CALL histwrite( kid, "sisali", kt, sm t_i, jpi*jpj, (/1/) )368 CALL histwrite( kid, "sisali", kt, sm_i , jpi*jpj, (/1/) ) 369 369 CALL histwrite( kid, "sivolu", kt, vt_i , jpi*jpj, (/1/) ) 370 370 CALL histwrite( kid, "sidive", kt, divu_i*1.0e8 , jpi*jpj, (/1/) ) … … 389 389 CALL histwrite( kid, "sithicat", kt, h_i , jpi*jpj*jpl, (/1/) ) 390 390 CALL histwrite( kid, "siconcat", kt, a_i , jpi*jpj*jpl, (/1/) ) 391 CALL histwrite( kid, "sisalcat", kt, s m_i, jpi*jpj*jpl, (/1/) )391 CALL histwrite( kid, "sisalcat", kt, s_i , jpi*jpj*jpl, (/1/) ) 392 392 CALL histwrite( kid, "sitemcat", kt, tm_i - rt0 , jpi*jpj*jpl, (/1/) ) 393 393 CALL histwrite( kid, "snthicat", kt, h_s , jpi*jpj*jpl, (/1/) )
Note: See TracChangeset
for help on using the changeset viewer.