Changeset 13334 for NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/ICE/iceistate.F90
- Timestamp:
- 2020-07-22T16:20:32+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/ICE/iceistate.F90
r13295 r13334 33 33 USE fldread ! read input fields 34 34 35 # if defined key_agrif 36 USE agrif_oce 37 USE agrif_ice 38 USE agrif_ice_interp 39 # endif 35 USE agrif_oce ! initial state interpolation 36 USE agrif_ice_interp 40 37 41 38 IMPLICIT NONE … … 177 174 178 175 IF( ln_iceini ) THEN 179 ! !---------------! 180 181 IF( Agrif_Root() ) THEN 182 176 177 #if defined key_agrif 178 IF ( ( Agrif_Root() ).OR.(.NOT.ln_init_chfrpar ) ) THEN 179 #endif 180 ! !---------------! 183 181 IF( ln_iceini_file )THEN ! Read a file ! 184 182 ! !---------------! … … 376 374 t1_ice(:,:,:) = t_i (:,:,1,:) 377 375 ! 378 379 #if defined key_agrif 376 #if defined key_agrif 380 377 ELSE 381 382 Agrif_SpecialValue = -9999. 383 Agrif_UseSpecialValue = .TRUE. 384 CALL Agrif_init_variable(tra_iceini_id,procname=interp_tra_ice) 385 use_sign_north = .TRUE. 386 sign_north = -1. 387 CALL Agrif_init_variable(u_iceini_id ,procname=interp_u_ice) 388 CALL Agrif_init_variable(v_iceini_id ,procname=interp_v_ice) 389 Agrif_SpecialValue = 0._wp 390 use_sign_north = .FALSE. 391 Agrif_UseSpecialValue = .FALSE. 392 ! lbc ???? 393 ! Here we know : a_i, v_i, v_s, sv_i, oa_i, a_ip, v_ip, t_su, e_s, e_i 394 CALL ice_var_glo2eqv 395 CALL ice_var_zapsmall 396 CALL ice_var_agg(2) 397 398 ! Melt ponds 399 WHERE( a_i > epsi10 ) 400 a_ip_frac(:,:,:) = a_ip(:,:,:) / a_i(:,:,:) 401 ELSEWHERE 402 a_ip_frac(:,:,:) = 0._wp 403 END WHERE 404 WHERE( a_ip > 0._wp ) ! ??????? 405 h_ip(:,:,:) = v_ip(:,:,:) / a_ip(:,:,:) 406 ELSEWHERE 407 h_ip(:,:,:) = 0._wp 408 END WHERE 409 410 tn_ice(:,:,:) = t_su(:,:,:) 411 t1_ice(:,:,:) = t_i (:,:,1,:) 378 CALL agrif_istate_ice 379 ENDIF 412 380 #endif 413 ENDIF ! Agrif_Root414 381 ENDIF ! ln_iceini 415 382 !
Note: See TracChangeset
for help on using the changeset viewer.