- Timestamp:
- 2020-10-01T13:33:30+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r13383_HPC-02_Daley_Tiling/src/ICE/icestp.F90
r13216 r13553 201 201 IF( lrst_ice ) CALL ice_rst_write( kt ) ! -- Ice restart file 202 202 ! 203 IF( ln_icectl ) CALL ice_ctl( kt ) ! -- alerts in case of model crash203 IF( ln_icectl ) CALL ice_ctl( kt ) ! -- Control checks 204 204 ! 205 205 ENDIF ! End sea-ice time step only … … 224 224 INTEGER, INTENT(in) :: Kbb, Kmm, Kaa 225 225 ! 226 INTEGER :: ji, jj,ierr226 INTEGER :: ierr 227 227 !!---------------------------------------------------------------------- 228 228 IF(lwp) WRITE(numout,*) … … 252 252 IF( ierr /= 0 ) CALL ctl_stop('STOP', 'ice_init : unable to allocate ice arrays') 253 253 ! 254 CALL ice_itd_init ! ice thickness distribution initialization255 !256 CALL ice_thd_init ! set ice thermodynics parameters (clem: important to call it first for melt ponds)257 !258 ! ! Initial sea-ice state259 IF( .NOT. ln_rstart ) THEN ! start from rest: sea-ice deduced from sst260 CALL ice_istate_init261 CALL ice_istate( nit000, Kbb, Kmm, Kaa )262 ELSE ! start from a restart file263 CALL ice_rst_read( Kbb, Kmm, Kaa )264 ENDIF265 CALL ice_var_glo2eqv266 CALL ice_var_agg(1)267 !268 CALL ice_sbc_init ! set ice-ocean and ice-atm. coupling parameters269 !270 CALL ice_dyn_init ! set ice dynamics parameters271 !272 CALL ice_update_init ! ice surface boundary condition273 !274 CALL ice_alb_init ! ice surface albedo275 !276 CALL ice_dia_init ! initialization for diags277 !278 fr_i (:,:) = at_i(:,:) ! initialisation of sea-ice fraction279 tn_ice(:,:,:) = t_su(:,:,:) ! initialisation of surface temp for coupled simu280 !281 254 ! ! set max concentration in both hemispheres 282 255 WHERE( gphit(:,:) > 0._wp ) ; rn_amax_2d(:,:) = rn_amax_n ! NH 283 256 ELSEWHERE ; rn_amax_2d(:,:) = rn_amax_s ! SH 284 257 END WHERE 285 258 ! 259 CALL diag_set0 ! set diag of mass, heat and salt fluxes to 0: needed for Agrif child grids 260 ! 261 CALL ice_itd_init ! ice thickness distribution initialization 262 ! 263 CALL ice_thd_init ! set ice thermodynics parameters (clem: important to call it first for melt ponds) 264 ! 265 ! ! Initial sea-ice state 266 CALL ice_istate_init 267 IF ( ln_rstart .OR. nn_iceini_file == 2 ) THEN 268 CALL ice_rst_read( Kbb, Kmm, Kaa ) ! start from a restart file 269 ELSE 270 CALL ice_istate( nit000, Kbb, Kmm, Kaa ) ! start from rest or read a file 271 ENDIF 272 CALL ice_var_glo2eqv 273 CALL ice_var_agg(1) 274 ! 275 CALL ice_sbc_init ! set ice-ocean and ice-atm. coupling parameters 276 ! 277 CALL ice_dyn_init ! set ice dynamics parameters 278 ! 279 CALL ice_update_init ! ice surface boundary condition 280 ! 281 CALL ice_alb_init ! ice surface albedo 282 ! 283 CALL ice_dia_init ! initialization for diags 284 ! 285 fr_i (:,:) = at_i(:,:) ! initialisation of sea-ice fraction 286 tn_ice(:,:,:) = t_su(:,:,:) ! initialisation of surface temp for coupled simu 287 ! 286 288 IF( ln_rstart ) CALL iom_close( numrir ) ! close input ice restart file 287 289 ! … … 366 368 v_s_b (:,:,:) = v_s (:,:,:) ! snow volume 367 369 sv_i_b(:,:,:) = sv_i(:,:,:) ! salt content 368 oa_i_b(:,:,:) = oa_i(:,:,:) ! areal age content369 370 e_s_b (:,:,:,:) = e_s (:,:,:,:) ! snow thermal energy 370 371 e_i_b (:,:,:,:) = e_i (:,:,:,:) ! ice thermal energy … … 375 376 h_i_b(:,:,:) = 0._wp 376 377 h_s_b(:,:,:) = 0._wp 377 END WHERE378 379 WHERE( a_ip(:,:,:) >= epsi20 )380 h_ip_b(:,:,:) = v_ip(:,:,:) / a_ip(:,:,:) ! ice pond thickness381 ELSEWHERE382 h_ip_b(:,:,:) = 0._wp383 378 END WHERE 384 379 ! … … 424 419 hfx_res(:,:) = 0._wp ; hfx_sub(:,:) = 0._wp 425 420 hfx_spr(:,:) = 0._wp ; hfx_dif(:,:) = 0._wp 426 hfx_err_rem(:,:) = 0._wp427 421 hfx_err_dif(:,:) = 0._wp 428 422 wfx_err_sub(:,:) = 0._wp … … 445 439 diag_trp_ei(:,:) = 0._wp ; diag_trp_es(:,:) = 0._wp 446 440 diag_trp_sv(:,:) = 0._wp 447 441 448 442 END SUBROUTINE diag_set0 449 443
Note: See TracChangeset
for help on using the changeset viewer.