| 76 | {{{ |
| 77 | !----------------------------------------------------------------------- |
| 78 | &namsbc ! Surface Boundary Condition (surface module) |
| 79 | !----------------------------------------------------------------------- |
| 80 | nn_isf = 1 ! 0=no isf / 1 = presence of ISF / 2 = bg03 parametrisation / 3 = rnf file for isf / 4 = prescribed melt rate |
| 81 | / |
| 82 | !----------------------------------------------------------------------- |
| 83 | &namsbc_isf ! Top boundary layer (ISF) |
| 84 | !----------------------------------------------------------------------- |
| 85 | ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! |
| 86 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! |
| 87 | ! nn_isf == 4 |
| 88 | sn_qisf = 'rnfisf' , -12 ,'sohflisf', .false. , .true. , 'yearly' , '' , '' |
| 89 | sn_fwfisf = 'rnfisf' , -12 ,'sowflisf', .false. , .true. , 'yearly' , '' , '' |
| 90 | ! nn_isf == 3 |
| 91 | sn_rnfisf = 'runoffs' , -12 ,'sofwfisf', .false. , .true. , 'yearly' , '' , '' |
| 92 | ! nn_isf == 2 and 3 |
| 93 | sn_depmax_isf = 'runoffs' , -12 ,'sozisfmax' , .false. , .true. , 'yearly' , '' , '' |
| 94 | sn_depmin_isf = 'runoffs' , -12 ,'sozisfmin' , .false. , .true. , 'yearly' , '' , '' |
| 95 | ! nn_isf == 2 |
| 96 | sn_Leff_isf = 'rnfisf' , 0 ,'Leff' , .false. , .true. , 'yearly' , '' , '' |
| 97 | ! for all case |
| 98 | ln_divisf = .true. ! apply isf melting as a mass flux or in the salinity trend. (maybe I should remove this option as for runoff?) |
| 99 | ! only for nn_isf = 1 or 2 |
| 100 | rn_gammat0 = 1.0e-4 ! gammat coefficient used in blk formula |
| 101 | rn_gammas0 = 1.0e-4 ! gammas coefficient used in blk formula |
| 102 | ! only for nn_isf = 1 |
| 103 | nn_isfblk = 1 ! 1 ISOMIP ; 2 conservative (3 equation formulation, Jenkins et al. 1991 ??) |
| 104 | rn_hisf_tbl = 30. ! thickness of the top boundary layer (Losh et al. 2008) |
| 105 | ! 0 => thickness of the tbl = thickness of the first wet cell |
| 106 | ln_conserve = .true. ! conservative case (take into account meltwater advection) |
| 107 | nn_gammablk = 0 ! 0 = cst Gammat (= gammat/s) |
| 108 | ! 1 = velocity dependend Gamma (u* * gammat/s) (Jenkins et al. 2010) |
| 109 | ! if you want to keep the cd as in global config, adjust rn_gammat0 to compensate |
| 110 | ! 2 = velocity and stability dependent Gamma Holland et al. 1999 |
| 111 | / |
| 112 | !----------------------------------------------------------------------- |
| 113 | &nambfr ! bottom/top friction |
| 114 | !----------------------------------------------------------------------- |
| 115 | rn_tfri1 = 4.e-4 ! top drag coefficient (linear case) |
| 116 | rn_tfri2 = 2.5e-3 ! top drag coefficient (non linear case). Minimum coeft if ln_loglayer=T |
| 117 | rn_tfri2_max = 1.e-1 ! max. top drag coefficient (non linear case and ln_loglayer=T) |
| 118 | rn_tfeb2 = 0.0 ! top turbulent kinetic energy background (m2/s2) |
| 119 | rn_tfrz0 = 3.e-3 ! top roughness [m] if ln_loglayer=T |
| 120 | ln_tfr2d = .false. ! horizontal variation of the top friction coef (read a 2D mask file ) |
| 121 | rn_tfrien = 50. ! local multiplying factor of tfr (ln_tfr2d=T) |
| 122 | / |
| 123 | }}} |