Changeset 9987 for branches/UKMO/dev_r5518_obs_oper_update_icethick/NEMOGCM/NEMO/LIM_SRC_3/limwri.F90
- Timestamp:
- 2018-07-23T11:33:03+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_obs_oper_update_icethick/NEMOGCM/NEMO/LIM_SRC_3/limwri.F90
r7960 r9987 157 157 ENDIF 158 158 159 IF ( iom_use( "icecolf" ) ) THEN 160 DO jj = 1, jpj 161 DO ji = 1, jpi 162 rswitch = MAX( 0._wp , SIGN( 1._wp , at_i(ji,jj) ) ) 163 z2d(ji,jj) = hicol(ji,jj) * rswitch 164 END DO 165 END DO 166 CALL iom_put( "icecolf" , z2d ) ! frazil ice collection thickness 167 ENDIF 159 IF ( iom_use( "icecolf" ) ) CALL iom_put( "icecolf", hicol ) ! frazil ice collection thickness 168 160 169 161 CALL iom_put( "isst" , sst_m ) ! sea surface temperature … … 190 182 CALL iom_put( "destrp" , diag_trp_es ) ! advected snw enthalpy (W/m2) 191 183 192 CALL iom_put( "sfxbog" , sfx_bog * rday ) ! salt flux from b rines193 CALL iom_put( "sfxbom" , sfx_bom * rday ) ! salt flux from b rines194 CALL iom_put( "sfxsum" , sfx_sum * rday ) ! salt flux from brines195 CALL iom_put( "sfxsni" , sfx_sni * rday ) ! salt flux from brines196 CALL iom_put( "sfxopw" , sfx_opw * rday ) ! salt flux from brines184 CALL iom_put( "sfxbog" , sfx_bog * rday ) ! salt flux from bottom growth 185 CALL iom_put( "sfxbom" , sfx_bom * rday ) ! salt flux from bottom melt 186 CALL iom_put( "sfxsum" , sfx_sum * rday ) ! salt flux from surface melt 187 CALL iom_put( "sfxsni" , sfx_sni * rday ) ! salt flux from snow ice formation 188 CALL iom_put( "sfxopw" , sfx_opw * rday ) ! salt flux from open water formation 197 189 CALL iom_put( "sfxdyn" , sfx_dyn * rday ) ! salt flux from ridging rafting 198 CALL iom_put( "sfxres" , sfx_res * rday ) ! salt flux from limupdate (resultant)190 CALL iom_put( "sfxres" , sfx_res * rday ) ! salt flux from residual 199 191 CALL iom_put( "sfxbri" , sfx_bri * rday ) ! salt flux from brines 192 CALL iom_put( "sfxsub" , sfx_sub * rday ) ! salt flux from sublimation 200 193 CALL iom_put( "sfx" , sfx * rday ) ! total salt flux 201 194 … … 235 228 CALL iom_put ('hfxdhc' , diag_heat(:,:) ) ! Heat content variation in snow and ice 236 229 CALL iom_put ('hfxspr' , hfx_spr(:,:) ) ! Heat content of snow precip 230 231 232 IF ( iom_use( "vfxthin" ) ) THEN ! ice production for open water + thin ice (<20cm) => comparable to observations 233 DO jj = 1, jpj 234 DO ji = 1, jpi 235 z2d(ji,jj) = vt_i(ji,jj) / MAX( at_i(ji,jj), epsi06 ) * zswi(ji,jj) ! mean ice thickness 236 END DO 237 END DO 238 WHERE( z2d(:,:) < 0.2 .AND. z2d(:,:) > 0. ) ; z2da = wfx_bog 239 ELSEWHERE ; z2da = 0._wp 240 END WHERE 241 CALL iom_put( "vfxthin", ( wfx_opw + z2da ) * ztmp ) 242 ENDIF 237 243 238 244 !-------------------------------- … … 311 317 !! 312 318 !! History : 313 !! 4. 1! 2013-06 (C. Rousset)319 !! 4.0 ! 2013-06 (C. Rousset) 314 320 !!---------------------------------------------------------------------- 315 321 INTEGER, INTENT( in ) :: kt ! ocean time-step index)
Note: See TracChangeset
for help on using the changeset viewer.