- Timestamp:
- 2010-05-03T13:59:46+02:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/DEV_r1837_mass_heat_salt_fluxes/NEMO/LIM_SRC_2/limdyn_2.F90
r1855 r1857 17 17 !!---------------------------------------------------------------------- 18 18 USE dom_oce ! ocean space and time domain 19 USE sbc_oce ! surface boundary condition : ocean20 USE phycst ! physical constants21 USE ice_2 ! LIM-2 sea-ice variables22 USE dom_ice_2 ! LIM-2 domain23 USE limistate_2 ! LIM-2 initial state24 USE limrhg_2 ! LIM-2rheology19 USE sbc_oce ! 20 USE phycst ! 21 USE ice_2 ! 22 USE dom_ice_2 ! 23 USE limistate_2 ! 24 USE limrhg_2 ! ice rheology 25 25 26 26 USE lbclnk ! … … 34 34 PUBLIC lim_dyn_2 ! routine called by sbc_ice_lim 35 35 36 !! * Module variables 36 37 REAL(wp) :: rone = 1.e0 ! constant value 37 38 38 39 # include "vectopt_loop_substitute.h90" 39 40 !!---------------------------------------------------------------------- 40 !! NEMO/LIM 3.3, UCL-LOCEAN-IPSL (2010)41 !! LIM 2.0, UCL-LOCEAN-IPSL (2006) 41 42 !! $Id$ 42 43 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) … … 58 59 !! - treatment of the case if no ice dynamic 59 60 !!--------------------------------------------------------------------- 60 INTEGER, INTENT(in) :: kt ! number of iteration61 !! 62 INTEGER :: ji, jj ! dummy loop indices63 INTEGER :: i_j1, i_jpj ! Starting/ending j-indices for rheology64 REAL(wp) :: zcoef ! temporary scalar61 INTEGER, INTENT(in) :: kt ! number of iteration 62 !! 63 INTEGER :: ji, jj ! dummy loop indices 64 INTEGER :: i_j1, i_jpj ! Starting/ending j-indices for rheology 65 REAL(wp) :: zcoef ! temporary scalar 65 66 REAL(wp), DIMENSION(jpj) :: zind ! i-averaged indicator of sea-ice 66 67 REAL(wp), DIMENSION(jpj) :: zmsk ! i-averaged of tmask … … 143 144 ENDIF 144 145 145 IF(ln_ctl) CALL prt_ctl( tab2d_1=u_ice, clinfo1=' lim_dyn : u_ice :', tab2d_2=v_ice, clinfo2=' v_ice :')146 IF(ln_ctl) CALL prt_ctl(tab2d_1=u_ice , clinfo1=' lim_dyn : u_ice :', tab2d_2=v_ice , clinfo2=' v_ice :') 146 147 147 148 ! computation of friction velocity … … 178 179 CALL lbc_lnk( ust2s, 'T', 1. ) ! T-point 179 180 ! 180 IF(ln_ctl) CALL prt_ctl( tab2d_1=ust2s, clinfo1=' lim_dyn : ust2s :')181 IF(ln_ctl) CALL prt_ctl(tab2d_1=ust2s , clinfo1=' lim_dyn : ust2s :') 181 182 182 183 END SUBROUTINE lim_dyn_2
Note: See TracChangeset
for help on using the changeset viewer.