Newer
Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
!-----------------------------------------------------------------------
&namberg ! iceberg parameters (default: OFF)
!-----------------------------------------------------------------------
ln_icebergs = .false. ! activate iceberg floats (force =F with "key_agrif")
!
! ! diagnostics:
ln_bergdia = .true. ! Calculate budgets
nn_verbose_level = 0 ! Turn on more verbose output if level > 0
nn_verbose_write = 15 ! Timesteps between verbose messages
nn_sample_rate = 1 ! Timesteps between sampling for trajectory storage
!
! ! iceberg setting:
! ! Initial mass required for an iceberg of each class
rn_initial_mass = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11
! ! Proportion of calving mass to apportion to each class
rn_distribution = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02
! ! Ratio between effective and real iceberg mass (non-dim)
! ! i.e. number of icebergs represented at a point
rn_mass_scaling = 2000., 200., 50., 20., 10., 5., 2., 1., 1., 1.
! thickness of newly calved bergs (m)
rn_initial_thickness = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250.
!
rn_rho_bergs = 850. ! Density of icebergs
rn_LoW_ratio = 1.5 ! Initial ratio L/W for newly calved icebergs
ln_operator_splitting = .true. ! Use first order operator splitting for thermodynamics
rn_bits_erosion_fraction = 0. ! Fraction of erosion melt flux to divert to bergy bits
rn_sicn_shift = 0. ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1)
ln_passive_mode = .false. ! iceberg - ocean decoupling
nn_test_icebergs = 10 ! Create test icebergs of this class (-1 = no)
! ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2)
rn_test_box = 108.0, 116.0, -66.0, -58.0
ln_use_calving = .false. ! Use calving data even when nn_test_icebergs > 0
rn_speed_limit = 0. ! CFL speed limit for a berg
ln_M2016 = .false. ! use Merino et al. (2016) modification (use of 3d ocean data instead of only sea surface data)
ln_icb_grd = .false. ! ground icb when icb bottom level hit oce bottom level (need ln_M2016 to be activated)
cn_dir = './' ! root directory for the calving data location
!___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
sn_icb = 'calving' , -1. ,'calvingmask', .true. , .true. , 'yearly' , '' , '' , ''
/