Changeset 11896
- Timestamp:
- 2019-11-13T12:40:52+01:00 (5 years ago)
- Location:
- NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/ISOMIP+
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/ISOMIP+/EXPREF/namelist_cfg
r11889 r11896 62 62 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 63 63 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! 64 sn_tem = 'TS_init', -12. , 'Tinit', .true. , .true. , 'yearly' , '' , '' , '' 65 sn_sal = 'TS_init', -12. , 'Sinit', .true. , .true. , 'yearly' , '' , '' , '' 64 sn_tem = 'nemo_base_WARM.nc' , -12. , 'Tinit', .false. , .true. , 'yearly' , '' , '' , '' 65 sn_sal = 'nemo_base_WARM.nc' , -12. , 'Sinit', .false. , .true. , 'yearly' , '' , '' , '' 66 67 sn_dmpt = 'nemo_base_WARM.nc' , -12. , 'Tinit', .false. , .true. , 'yearly' , '' , '' , '' 68 sn_dmps = 'nemo_base_WARM.nc' , -12. , 'Sinit', .false. , .true. , 'yearly' , '' , '' , '' 66 69 / 67 70 !----------------------------------------------------------------------- … … 111 114 112 115 ln_usr = .true. ! user defined formulation (T => check usrdef_sbc) 116 nn_fwb = 1 113 117 / 114 118 !----------------------------------------------------------------------- … … 181 185 ! ! ad15 = velocity dependend Gamma (u* * gammat/s) (Jenkins et al. 2010) 182 186 ! ! hj99 = velocity and stability dependent Gamma (Holland et al. 1999) 183 rn_gammat0 = 0.0 45! gammat coefficient used in blk formula184 rn_gammas0 = 1.29e-3! gammas coefficient used in blk formula187 rn_gammat0 = 0.0215 ! gammat coefficient used in blk formula 188 rn_gammas0 = 0.614e-3 ! gammas coefficient used in blk formula 185 189 rn_vtide = 0.01 ! tidal velocity [m/s] 186 190 ! … … 227 231 ! ! free slip ! partial slip ! no slip ! strong slip 228 232 rn_shlat = 2. ! shlat = 0 ! 0 < shlat < 2 ! shlat = 2 ! 2 < shlat 229 ln_vorlat = .false. ! consistency of vorticity boundary condition with analytical Eqs.230 233 / 231 234 !----------------------------------------------------------------------- … … 270 273 !----------------------------------------------------------------------- 271 274 rn_Cd0 = 2.5e-3 ! drag coefficient [-] 275 rn_ke0 = 0.0e-3 ! background kinetic energy [m2/s2] (non-linear cases) 272 276 / 273 277 !----------------------------------------------------------------------- … … 275 279 !----------------------------------------------------------------------- 276 280 rn_Cd0 = 2.5e-3 ! drag coefficient [-] 281 rn_ke0 = 0.0e-3 ! background kinetic energy [m2/s2] (non-linear cases) 277 282 / 278 283 !----------------------------------------------------------------------- … … 300 305 &nameos ! ocean Equation Of Seawater (default: NO selection) 301 306 !----------------------------------------------------------------------- 302 ln_eos80 = .true. ! = Use EOS80 307 ln_teos10 = .false. ! = Use TEOS-10 308 ln_eos80 = .false. ! = Use EOS80 309 ln_leos = .true. ! = Use S-EOS (simplified Eq.) 310 ! 311 ! ! S-EOS coefficients (ln_seos=T): 312 ! ! rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 313 ! ! L-EOS coefficients (ln_seos=T): 314 ! ! rd(T,S,Z)*rau0 = rau0*(-a0*dT+b0*dS) 315 rn_a0 = 3.7330e-5 ! thermal expension coefficient 316 rn_b0 = 7.8430e-4 ! saline expension coefficient 303 317 / 304 318 !----------------------------------------------------------------------- -
NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/ISOMIP+/MY_SRC/sbcfwb.F90
r11889 r11896 105 105 ! 106 106 IF( MOD( kt-1, kn_fsbc ) == 0 ) THEN 107 y_fwfnow(1) = local_sum( e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf_cav(:,:) + fwfisf_par(:,:) - snwice_fmass(:,:) ) ) 108 CALL mpp_delay_sum( 'sbcfwb', 'fwb', y_fwfnow(:), z_fwfprv(:), kt == nitend - nn_fsbc + 1 ) 109 z_fwfprv(1) = z_fwfprv(1) / area 110 zcoef = z_fwfprv(1) * rcp 111 emp(:,:) = emp(:,:) - z_fwfprv(1) * tmask(:,:,1) 112 qns(:,:) = qns(:,:) + zcoef * sst_m(:,:) * tmask(:,:,1) ! account for change to the heat budget due to fw correction 113 ENDIF 114 ! 115 CASE ( 4 ) !== global mean fwf set to zero (ISOMIP case) ==! 116 ! 117 IF( MOD( kt-1, kn_fsbc ) == 0 ) THEN 107 118 z_fwf = glob_sum( 'sbcfwb', e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf_cav(:,:) + fwfisf_par(:,:) - snwice_fmass(:,:) ) ) 108 119 z_fwf = z_fwf / area 109 120 zcoef = z_fwf * rcp 110 emp(:,:) = emp(:,:) - z_fwf * tmask(:,:,1)111 qns(:,:) = qns(:,:) + zcoef * sst_m(:,:) * tmask(:,:,1) ! account for change to the heat budget due to fw correction112 sfx(:,:) = sfx(:,:) + z_fwf * sss_m(:,:) * tmask(:,:,1)121 emp(:,:) = emp(:,:) - z_fwf * tmask(:,:,1) ! (Eq. 34 AD2015) 122 qns(:,:) = qns(:,:) + zcoef * ( -1.9 ) * tmask(:,:,1) ! (Eq. 35 AD2015) ! could be sst_m if we don't want any bouyancy fluxes 123 sfx(:,:) = sfx(:,:) + z_fwf * ( 33.8 ) * tmask(:,:,1) ! (Eq. 36 AD2015) ! could be sss_m if we don't want any bouyancy fluxes 113 124 ENDIF 114 125 !
Note: See TracChangeset
for help on using the changeset viewer.