Changeset 8407
- Timestamp:
- 2017-08-07T16:45:49+02:00 (8 years ago)
- Location:
- branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO
- Files:
-
- 3 added
- 3 deleted
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/ice.F90
r8360 r8407 74 74 !! in icestp.F90, the routines that compute the changes in the ice | 75 75 !! state variables are called | 76 !! - lim_dyn: ice dynamics |76 !! - ice_rhg : ice dynamics | 77 77 !! - lim_trp : ice transport | 78 78 !! - lim_itd_me : mechanical redistribution (ridging and rafting) | … … 208 208 REAL(wp), PUBLIC :: rn_crhg !: determines changes in ice strength 209 209 LOGICAL , PUBLIC :: ln_icestr_bvf !: use brine volume to diminish ice strength 210 ! -- limdyn & limrhg -- !210 ! -- limdyn & icerhg -- ! 211 211 REAL(wp), PUBLIC :: rn_ishlat !: lateral boundary condition for sea-ice 212 212 REAL(wp), PUBLIC :: rn_cio !: drag coefficient for oceanic stress -
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icestp.F90
r8404 r8407 32 32 USE phycst ! Define parameters for the routines 33 33 USE eosbn2 ! equation of state 34 USE limdyn! Ice dynamics34 USE icerhg ! Ice dynamics 35 35 USE limtrp ! Ice transport 36 36 USE limthd ! Ice thermodynamics … … 60 60 61 61 USE bdy_oce , ONLY: ln_bdy 62 USE bdyice _lim ! unstructured open boundary data (bdy_ice_lim routine)62 USE bdyice ! unstructured open boundary data 63 63 # if defined key_agrif 64 64 USE agrif_ice … … 106 106 ! 3 = bulk formulation, 107 107 ! 4 = Pure Coupled formulation) 108 !!109 108 INTEGER :: jl ! dummy loop index 110 REAL(wp), DIMENSION(jpi,jpj,jpl) :: zalb_os, zalb_cs ! ice albedo under overcast/clear sky111 REAL(wp), DIMENSION(jpi,jpj) :: zutau_ice, zvtau_ice112 109 !!---------------------------------------------------------------------- 113 110 … … 158 155 ! 159 156 IF( nn_limdyn /= 0 ) THEN ! -- Ice dynamics 160 CALL lim_dyn( kt ) ! rheology157 CALL ice_rhg( kt ) ! rheology 161 158 ELSE 162 159 u_ice(:,:) = rn_uice * umask(:,:,1) ! or prescribed velocity … … 176 173 IF( .NOT. Agrif_Root() ) CALL agrif_interp_lim3('T') 177 174 #endif 178 IF( ln_limthd .AND. ln_bdy ) CALL bdy_ice _lim( kt )! -- bdy ice thermo175 IF( ln_limthd .AND. ln_bdy ) CALL bdy_ice( kt ) ! -- bdy ice thermo 179 176 ! previous lead fraction and ice volume for flux calculations 180 177 CALL lim_var_glo2eqv ! ht_i and ht_s for ice albedo calculation … … 278 275 IF( ierr /= 0 ) CALL ctl_stop('STOP', 'ice_init : unable to allocate ice arrays') 279 276 ! 280 CALL lim_dyn_init ! set ice dynamics parameters277 CALL ice_rhg_init ! set ice dynamics parameters 281 278 ! 282 279 CALL ice_itd_init ! ice thickness distribution initialization
Note: See TracChangeset
for help on using the changeset viewer.