Changeset 5167 for trunk/NEMOGCM/NEMO/OPA_SRC/SBC
- Timestamp:
- 2015-03-24T18:35:00+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim.F90
r5146 r5167 184 184 numit = numit + nn_fsbc ! Ice model time step 185 185 ! 186 CALL sbc_lim_ update! Store previous ice values186 CALL sbc_lim_bef ! Store previous ice values 187 187 188 188 CALL sbc_lim_diag0 ! set diag of mass, heat and salt fluxes to 0 … … 202 202 203 203 #if defined key_bdy 204 CALL lim_var_glo2eqv205 204 CALL bdy_ice_lim( kt ) ! bdy ice thermo 206 CALL lim_var_zapsmall207 CALL lim_var_agg(1)208 205 IF( ln_icectl ) CALL lim_prt( kt, iiceprt, jiceprt, 1, ' - ice thermo bdy - ' ) 209 206 #endif … … 212 209 ENDIF 213 210 214 CALL sbc_lim_ update! Store previous ice values211 CALL sbc_lim_bef ! Store previous ice values 215 212 216 213 ! ---------------------------------------------- 217 214 ! ice thermodynamics 218 215 ! ---------------------------------------------- 219 CALL lim_var_glo2eqv220 216 CALL lim_var_agg(1) 221 217 … … 248 244 ! 249 245 IF( lrst_ice ) CALL lim_rst_write( kt ) ! Ice restart file 250 CALL lim_var_glo2eqv ! ??? 251 ! 252 IF( ln_icectl ) CALL lim_ctl( kt ) ! alerts in case of model crash 246 ! 247 IF( ln_icectl ) CALL lim_ctl( kt ) ! alerts in case of model crash 253 248 ! 254 249 CALL wrk_dealloc( jpi,jpj,jpl, zalb_os, zalb_cs, zalb_ice ) … … 389 384 r1_nlay_s = 1._wp / REAL( nlay_s, wp ) 390 385 ! 386 #if defined key_bdy 387 IF( lwp .AND. ln_limdiahsb ) CALL ctl_warn('online conservation check activated but it does not work with BDY') 388 #endif 389 ! 391 390 END SUBROUTINE ice_run 392 391 … … 555 554 END SUBROUTINE ice_lim_flx 556 555 557 SUBROUTINE sbc_lim_ update558 !!---------------------------------------------------------------------- 559 !! *** ROUTINE sbc_lim_ update***556 SUBROUTINE sbc_lim_bef 557 !!---------------------------------------------------------------------- 558 !! *** ROUTINE sbc_lim_bef *** 560 559 !! 561 560 !! ** purpose : store ice variables at "before" time step … … 571 570 v_ice_b(:,:) = v_ice(:,:) 572 571 573 END SUBROUTINE sbc_lim_ update572 END SUBROUTINE sbc_lim_bef 574 573 575 574 SUBROUTINE sbc_lim_diag0 … … 602 601 hfx_spr(:,:) = 0._wp ; hfx_dif(:,:) = 0._wp 603 602 hfx_err(:,:) = 0._wp ; hfx_err_rem(:,:) = 0._wp 604 hfx_err_dif(:,:) = 0._wp 603 hfx_err_dif(:,:) = 0._wp ; 605 604 606 605 afx_tot(:,:) = 0._wp ; 607 606 afx_dyn(:,:) = 0._wp ; afx_thd(:,:) = 0._wp 608 607 609 diag_heat_dhc(:,:) = 0._wp ; 608 diag_heat(:,:) = 0._wp ; diag_smvi(:,:) = 0._wp ; 609 diag_vice(:,:) = 0._wp ; diag_vsnw(:,:) = 0._wp ; 610 610 611 611 END SUBROUTINE sbc_lim_diag0 … … 636 636 637 637 fice_ice_ave (:,:) = 0.0_wp 638 WHERE ( at_i (:,:) .GT.0.0_wp ) fice_ice_ave (:,:) = fice_cell_ave ( ptab (:,:,:)) / at_i (:,:)638 WHERE ( at_i (:,:) > 0.0_wp ) fice_ice_ave (:,:) = fice_cell_ave ( ptab (:,:,:)) / at_i (:,:) 639 639 640 640 END FUNCTION fice_ice_ave
Note: See TracChangeset
for help on using the changeset viewer.