Changeset 2475 for branches/nemo_v3_3_beta/NEMOGCM/NEMO/LIM_SRC_2/ice_2.F90
- Timestamp:
- 2010-12-17T07:45:43+01:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/nemo_v3_3_beta/NEMOGCM/NEMO/LIM_SRC_2/ice_2.F90
r2370 r2475 15 15 IMPLICIT NONE 16 16 PRIVATE 17 18 INTEGER , PUBLIC :: numit !: ice iteration index 19 REAL(wp), PUBLIC :: rdt_ice !: ice time step 17 20 18 ! !* Share parameters namelist (namicerun read in iceini)*21 ! !!* namicerun read in iceini * 19 22 CHARACTER(len=32) , PUBLIC :: cn_icerst_in = "restart_ice_in" !: suffix of ice restart name (input) 20 23 CHARACTER(len=32) , PUBLIC :: cn_icerst_out = "restart_ice" !: suffix of ice restart name (output) … … 22 25 LOGICAL , PUBLIC :: ln_limdmp = .FALSE. !: Ice damping 23 26 LOGICAL , PUBLIC :: ln_nicep = .TRUE. !: flag grid points output (T) or not (F) 24 REAL(wp) , PUBLIC :: hsndif = 0.e0 !: computation of temp. in snow (0) or not (9999) 25 REAL(wp) , PUBLIC :: hicdif = 0.e0 !: computation of temp. in ice (0) or not (9999) 26 REAL(wp), DIMENSION(2), PUBLIC :: acrit = (/ 1.e-06 , 1.e-06 /) !: minimum fraction for leads in 27 ! !: north and south hemisphere 28 !!* ice-dynamic namelist (namicedyn) * 29 INTEGER , PUBLIC :: nbiter = 1 !: number of sub-time steps for relaxation 30 INTEGER , PUBLIC :: nbitdr = 250 !: maximum number of iterations for relaxation 31 INTEGER , PUBLIC :: nevp = 360 !: number of EVP subcycling iterations 32 INTEGER , PUBLIC :: telast = 3600 !: timescale for EVP elastic waves 33 REAL(wp), PUBLIC :: rdt_ice !: ice time step 34 REAL(wp), PUBLIC :: epsd = 1.0e-20 !: tolerance parameter for dynamic 35 REAL(wp), PUBLIC :: alpha = 0.5 !: coefficient for semi-implicit coriolis 36 REAL(wp), PUBLIC :: dm = 0.6e+03 !: diffusion constant for dynamics 37 REAL(wp), PUBLIC :: om = 0.5 !: relaxation constant 38 REAL(wp), PUBLIC :: resl = 5.0e-05 !: maximum value for the residual of relaxation 39 REAL(wp), PUBLIC :: cw = 5.0e-03 !: drag coefficient for oceanic stress 40 REAL(wp), PUBLIC :: angvg = 0.e0 !: turning angle for oceanic stress 41 REAL(wp), PUBLIC :: pstar = 1.0e+04 !: first bulk-rheology parameter 42 REAL(wp), PUBLIC :: c_rhg = 20.e0 !: second bulk-rhelogy parameter 43 REAL(wp), PUBLIC :: etamn = 0.e+07 !: minimun value for viscosity 44 REAL(wp), PUBLIC :: creepl = 2.e-08 !: creep limit 45 REAL(wp), PUBLIC :: ecc = 2.e0 !: eccentricity of the elliptical yield curve 46 REAL(wp), PUBLIC :: ahi0 = 350.e0 !: sea-ice hor. eddy diffusivity coeff. (m2/s) 47 REAL(wp), PUBLIC :: alphaevp = 1.e0 !: coefficient for the solution of EVP int. stresses 27 REAL(wp) , PUBLIC :: hsndif = 0._wp !: snow temp. computation (0) or not (9999) 28 REAL(wp) , PUBLIC :: hicdif = 0._wp !: ice temp. computation (0) or not (9999) 29 REAL(wp), DIMENSION(2), PUBLIC :: acrit = (/ 1.e-6_wp , 1.e-6_wp /) !: minimum lead fraction in the 2 hemisphere 30 31 ! !!* ice-dynamic namelist (namicedyn) * 32 INTEGER , PUBLIC :: nbiter = 1 !: number of sub-time steps for relaxation 33 INTEGER , PUBLIC :: nbitdr = 250 !: maximum number of iterations for relaxation 34 INTEGER , PUBLIC :: nevp = 360 !: number of EVP subcycling iterations 35 INTEGER , PUBLIC :: telast = 3600 !: timescale for EVP elastic waves 36 REAL(wp), PUBLIC :: epsd = 1.0e-20_wp !: tolerance parameter for dynamic 37 REAL(wp), PUBLIC :: alpha = 0.5_wp !: coefficient for semi-implicit coriolis 38 REAL(wp), PUBLIC :: dm = 0.6e+03_wp !: diffusion constant for dynamics 39 REAL(wp), PUBLIC :: om = 0.5_wp !: relaxation constant 40 REAL(wp), PUBLIC :: resl = 5.0e-05_wp !: maximum value for the residual of relaxation 41 REAL(wp), PUBLIC :: cw = 5.0e-03_wp !: drag coefficient for oceanic stress 42 REAL(wp), PUBLIC :: angvg = 0._wp !: turning angle for oceanic stress 43 REAL(wp), PUBLIC :: pstar = 1.0e+04_wp !: first bulk-rheology parameter 44 REAL(wp), PUBLIC :: c_rhg = 20._wp !: second bulk-rhelogy parameter 45 REAL(wp), PUBLIC :: etamn = 0._wp !: minimun value for viscosity 46 REAL(wp), PUBLIC :: creepl = 2.e-08_wp !: creep limit 47 REAL(wp), PUBLIC :: ecc = 2._wp !: eccentricity of the elliptical yield curve 48 REAL(wp), PUBLIC :: ahi0 = 350._wp !: sea-ice hor. eddy diffusivity coeff. (m2/s) 49 REAL(wp), PUBLIC :: alphaevp = 1._wp !: coefficient for the solution of EVP int. stresses 48 50 49 REAL(wp), PUBLIC :: usecc2 !: = 1.0 / ( ecc * ecc )50 REAL(wp), PUBLIC :: rhoco !: = rau0 * cw51 REAL(wp), PUBLIC :: sangvg, cangvg !: sin and cos of the turning angle for ocean stress52 REAL(wp), PUBLIC :: pstarh !: pstar / 2.051 REAL(wp), PUBLIC :: usecc2 !: = 1.0 / ( ecc * ecc ) 52 REAL(wp), PUBLIC :: rhoco !: = rau0 * cw 53 REAL(wp), PUBLIC :: sangvg, cangvg !: sin and cos of the turning angle for ocean stress 54 REAL(wp), PUBLIC :: pstarh !: pstar / 2.0 53 55 54 56 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: ahiu , ahiv !: hor. diffusivity coeff. at ocean U- and V-points (m2/s) … … 58 60 !!* Ice Rheology 59 61 # if defined key_lim2_vp 60 ! 62 ! !!* VP rheology * 61 63 LOGICAL , PUBLIC :: lk_lim2_vp = .TRUE. !: Visco-Plactic reology flag 62 64 ! … … 64 66 ! 65 67 # else 66 ! 68 ! !!* EVP rheology * 67 69 LOGICAL , PUBLIC:: lk_lim2_vp = .FALSE. !: Visco-Plactic reology flag 68 70 ! … … 96 98 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: rdmicif !: Variation of ice mass 97 99 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: qldif !: heat balance of the lead (or of the open ocean) 98 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: qcmif !: Energy needed to bring the ocean surface layer until its freezing100 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: qcmif !: Energy needed to freeze the ocean surface layer 99 101 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: fdtcn !: net downward heat flux from the ice to the ocean 100 102 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: qdtcn !: energy from the ice to the ocean point (at a factor 2) 101 103 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: thcm !: part of the solar energy used in the lead heat budget 102 104 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: fstric !: Solar flux transmitted trough the ice 103 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: ffltbif !: Arraylinked with the max heat contained in brine pockets (?)105 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: ffltbif !: linked with the max heat contained in brine pockets (?) 104 106 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: fscmbq !: Linked with the solar flux below the ice (?) 105 107 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: fsbbq !: Also linked with the solar flux below the ice (?) 106 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: qfvbq !: Arrayused to store energy in case of toral lateral ablation (?)108 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: qfvbq !: used to store energy in case of toral lateral ablation (?) 107 109 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: dmgwi !: Variation of the mass of snow ice 108 109 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: u_ice, v_ice !: two components of the ice velocity at I-point (m/s) 110 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: u_oce, v_oce !: two components of the ocean velocity at I-point (m/s) 110 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: u_ice, v_ice !: two components of the ice velocity at I-point (m/s) 111 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: u_oce, v_oce !: two components of the ocean velocity at I-point (m/s) 111 112 112 113 REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jplayersp1) :: tbif !: Temperature inside the ice/snow layer
Note: See TracChangeset
for help on using the changeset viewer.