Changeset 921 for trunk/NEMO/LIM_SRC_3/ice.F90
- Timestamp:
- 2008-05-13T10:28:52+02:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/LIM_SRC_3/ice.F90
r904 r921 217 217 s_i_min = 0.1 , & !: minimum ice salinity (ppt) 218 218 s_i_0 = 3.5 , & !: 1st sal. value for the computation of sal .prof. 219 !: (ppt)219 !: (ppt) 220 220 s_i_1 = 4.5 , & !: 2nd sal. value for the computation of sal .prof. 221 !: (ppt)221 !: (ppt) 222 222 sal_G = 5.00 , & !: restoring salinity for gravity drainage 223 !: (ppt)223 !: (ppt) 224 224 sal_F = 2.50 , & !: restoring salinity for flushing 225 !: (ppt)225 !: (ppt) 226 226 time_G = 1.728e+06,&!: restoring time constant for gravity drainage 227 !: (= 20 days, in s)227 !: (= 20 days, in s) 228 228 time_F = 8.640e+05,&!: restoring time constant for gravity drainage 229 !: (= 10 days, in s)229 !: (= 10 days, in s) 230 230 bulk_sal = 4.0 !: bulk salinity (ppt) in case of constant salinity 231 231 232 232 INTEGER , PUBLIC :: & !!: ** ice-salinity namelist (namicesal) ** 233 233 num_sal = 1 , & !: salinity configuration used in the model 234 !: 1 - s constant in space and time235 !: 2 - prognostic salinity (s(z,t))236 !: 3 - salinity profile, constant in time237 !: 4 - salinity variations affect only ice238 ! thermodynamics234 !: 1 - s constant in space and time 235 !: 2 - prognostic salinity (s(z,t)) 236 !: 3 - salinity profile, constant in time 237 !: 4 - salinity variations affect only ice 238 ! thermodynamics 239 239 sal_prof = 1 , & !: salinity profile or not 240 240 thcon_i_swi = 1 !: thermal conductivity of Untersteiner (1964) (1) or 241 241 !: Pringle et al (2007) (2) 242 242 243 243 REAL(wp), PUBLIC :: & !!: ** ice-mechanical redistribution namelist (namiceitdme) … … 249 249 astar = 0.05 , & !!: equivalent of G* for an exponential participation function 250 250 Hstar = 100.0 , & !!: thickness that determines the maximal thickness of ridged 251 !!: ice251 !!: ice 252 252 hparmeter = 0.75, & !!: threshold thickness (m) for rafting / ridging 253 253 Craft = 5.0 , & !!: coefficient for smoothness of the hyperbolic tangent in rafting … … 256 256 betas = 1.0 , & !:: coef. for partitioning of snowfall between leads and sea ice 257 257 kappa_i = 1.0 , & !!: coefficient for the extinction of radiation 258 !!: Grenfell et al. (2006) (m-1)258 !!: Grenfell et al. (2006) (m-1) 259 259 nconv_i_thd = 50 , & !!: maximal number of iterations for heat diffusion 260 260 maxer_i_thd = 1.0e-4 !!: maximal tolerated error (C) for heat diffusion … … 264 264 raftswi = 1, & !!: rafting of ice or not 265 265 partfun_swi = 1, & !!: participation function Thorndike et al. JGR75 (0) 266 !!: or Lipscomb et al. JGR07 (1)266 !!: or Lipscomb et al. JGR07 (1) 267 267 transfun_swi = 0, & !!: transfer function of Hibler, MWR80 (0) 268 !!: or Lipscomb et al., 2007 (1)268 !!: or Lipscomb et al., 2007 (1) 269 269 brinstren_swi = 0 !!: use brine volume to diminish ice strength 270 270 … … 301 301 t_bo , & !: Sea-Ice bottom temperature (Kelvin) 302 302 hicifp , & !: Ice production/melting 303 !obsolete... can be removed303 !obsolete... can be removed 304 304 frld , & !: Leads fraction = 1-a/totalarea REFERS TO LEAD FRACTION everywhere 305 !: except in the OUTPUTS!!!!305 !: except in the OUTPUTS!!!! 306 306 pfrld , & !: Leads fraction at previous time 307 307 phicif , & !: Old ice thickness … … 328 328 fheat_res, & !: Residual heat flux due to correction of ice thickness 329 329 fhmec !: Heat flux due to snow loss during compression 330 330 331 331 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: & !: 332 332 albege , & !: Albedo of the snow or ice (only for outputs) … … 334 334 tauc !: Cloud optical depth 335 335 336 ! temporary arrays for dummy version of the code336 ! temporary arrays for dummy version of the code 337 337 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: & !: 338 338 dh_i_surf2D, dh_i_bott2D, fstbif, fsup2D, focea2D, q_s … … 354 354 sm_i , & !: Sea-Ice Bulk salinity (ppt) 355 355 smv_i , & !: Sea-Ice Bulk salinity times volume per area (ppt.m) 356 !: this is an extensive variable that has to be transported356 !: this is an extensive variable that has to be transported 357 357 o_i , & !: Sea-Ice Age (days) 358 358 ov_i , & !: Sea-Ice Age times volume per area (days.m) … … 401 401 !!-------------------------------------------------------------------------- 402 402 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: & !: 403 403 sxopw, syopw, sxxopw, syyopw, sxyopw !: open water in sea ice 404 404 405 405 REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpl) :: & !: 406 407 408 409 410 411 406 sxice, syice, sxxice, syyice, sxyice, & !: ice thickness moments for advection 407 sxsn, sysn, sxxsn, syysn, sxysn, & !: snow thickness 408 sxa, sya, sxxa, syya, sxya, & !: lead fraction 409 sxc0, syc0, sxxc0, syyc0, sxyc0, & !: snow thermal content 410 sxsal, sysal, sxxsal, syysal, sxysal, & !: ice salinity 411 sxage, syage, sxxage, syyage, sxyage !: ice age 412 412 413 413 REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jkmax,jpl) :: & !: 414 414 sxe , sye , sxxe , syye , sxye !: ice layers heat content 415 415 416 416 !!-------------------------------------------------------------------------- … … 446 446 REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jkmax,jpl) :: & !: 447 447 d_e_i_thd, d_e_i_trp 448 448 449 449 REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: & !: ice velocity 450 450 d_u_ice_dyn, d_v_ice_dyn … … 459 459 INTEGER, PUBLIC, DIMENSION(jpm,2) :: & !: 460 460 ice_cat_bounds !: Matrix containing the integer upper and 461 461 !: lower boundaries of ice thickness categories 462 462 463 463 ! REMOVE … … 474 474 REAL(wp), PUBLIC, DIMENSION(0:jpl,jpm) :: & !: 475 475 hi_max_typ !: Boundary of ice thickness categories 476 !:in thickness space (same but specific for each ice type) 476 !:in thickness space (same but specific for each ice type) 477 478 !!-------------------------------------------------------------------------- 479 !! * Ice Run 480 !!-------------------------------------------------------------------------- 481 !! Namelist namicerun read in iceini 482 LOGICAL , PUBLIC :: & !!! ** init namelist (namicerun) ** 483 ln_limdyn = .TRUE., & !: flag for ice dynamics (T) or not (F) 484 ln_nicep = .TRUE. !: flag for sea-ice points output (T) or not (F) 485 REAL(wp), PUBLIC :: & !: 486 hsndif = 0.e0 , & !: computation of temp. in snow (0) or not (9999) 487 hicdif = 0.e0 , & !: computation of temp. in ice (0) or not (9999) 488 cai = 1.40e-3 , & !: atmospheric drag over sea ice 489 cao = 1.00e-3 !: atmospheric drag over ocean 490 REAL(wp), PUBLIC, DIMENSION(2) :: & !: 491 acrit = (/ 1.e-06 , 1.e-06 /) !: minimum fraction for leads in 492 ! ! north and south hemisphere 477 493 478 494 !!--------------------------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.