- Timestamp:
- 2016-11-21T10:38:43+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_CNRS_2016/NEMOGCM/NEMO/LIM_SRC_3/limwri.F90
r6140 r7278 154 154 ENDIF 155 155 156 IF ( iom_use( "icecolf" ) ) THEN 157 DO jj = 1, jpj 158 DO ji = 1, jpi 159 rswitch = MAX( 0._wp , SIGN( 1._wp , at_i(ji,jj) ) ) 160 z2d(ji,jj) = hicol(ji,jj) * rswitch 161 END DO 162 END DO 163 CALL iom_put( "icecolf" , z2d ) ! frazil ice collection thickness 164 ENDIF 165 156 IF ( iom_use( "icecolf" ) ) CALL iom_put( "icecolf", hicol ) ! frazil ice collection thickness 157 166 158 CALL iom_put( "isst" , sst_m ) ! sea surface temperature 167 159 CALL iom_put( "isss" , sss_m ) ! sea surface salinity … … 187 179 CALL iom_put( "destrp" , diag_trp_es ) ! advected snw enthalpy (W/m2) 188 180 189 CALL iom_put( "sfxbog" , sfx_bog * rday ) ! salt flux from b rines190 CALL iom_put( "sfxbom" , sfx_bom * rday ) ! salt flux from b rines191 CALL iom_put( "sfxsum" , sfx_sum * rday ) ! salt flux from brines192 CALL iom_put( "sfxsni" , sfx_sni * rday ) ! salt flux from brines193 CALL iom_put( "sfxopw" , sfx_opw * rday ) ! salt flux from brines181 CALL iom_put( "sfxbog" , sfx_bog * rday ) ! salt flux from bottom growth 182 CALL iom_put( "sfxbom" , sfx_bom * rday ) ! salt flux from bottom melt 183 CALL iom_put( "sfxsum" , sfx_sum * rday ) ! salt flux from surface melt 184 CALL iom_put( "sfxsni" , sfx_sni * rday ) ! salt flux from snow ice formation 185 CALL iom_put( "sfxopw" , sfx_opw * rday ) ! salt flux from open water formation 194 186 CALL iom_put( "sfxdyn" , sfx_dyn * rday ) ! salt flux from ridging rafting 195 CALL iom_put( "sfxres" , sfx_res * rday ) ! salt flux from limupdate (resultant)187 CALL iom_put( "sfxres" , sfx_res * rday ) ! salt flux from residual 196 188 CALL iom_put( "sfxbri" , sfx_bri * rday ) ! salt flux from brines 189 CALL iom_put( "sfxsub" , sfx_sub * rday ) ! salt flux from sublimation 197 190 CALL iom_put( "sfx" , sfx * rday ) ! total salt flux 198 191 … … 233 226 CALL iom_put ('hfxspr' , hfx_spr(:,:) ) ! Heat content of snow precip 234 227 228 IF ( iom_use( "vfxthin" ) ) THEN ! ice production for open water + thin ice (<20cm) => comparable to observations 229 DO jj = 1, jpj 230 DO ji = 1, jpi 231 z2d(ji,jj) = vt_i(ji,jj) / MAX( at_i(ji,jj), epsi06 ) * zswi(ji,jj) ! mean ice thickness 232 END DO 233 END DO 234 WHERE( z2d(:,:) < 0.2 .AND. z2d(:,:) > 0. ) ; z2da = wfx_bog 235 ELSEWHERE ; z2da = 0._wp 236 END WHERE 237 CALL iom_put( "vfxthin", ( wfx_opw + z2da ) * ztmp ) 238 ENDIF 239 235 240 !-------------------------------- 236 241 ! Output values for each category
Note: See TracChangeset
for help on using the changeset viewer.