Changeset 7596 for branches/2016/v3_6_CMIP6_ice_diagnostics/NEMOGCM
- Timestamp:
- 2017-01-20T18:50:18+01:00 (7 years ago)
- Location:
- branches/2016/v3_6_CMIP6_ice_diagnostics/NEMOGCM
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/v3_6_CMIP6_ice_diagnostics/NEMOGCM/CONFIG/SHARED/namelist_ice_lim3_ref
r7537 r7596 98 98 ! 0: k = k0 + beta.S/T (Untersteiner, 1964) 99 99 ! 1: k = k0 + beta1.S/T - beta2.T (Pringle et al., 2007) 100 rn_ rcdsn= 0.31 ! thermal conductivity of the snow (0.31 W/m/K, Maykut and Untersteiner, 1971). Obs: 0.1-0.5 (Lecomte et al, JAMES 2013)100 rn_cdsn = 0.31 ! thermal conductivity of the snow (0.31 W/m/K, Maykut and Untersteiner, 1971). Obs: 0.1-0.5 (Lecomte et al, JAMES 2013) 101 101 nn_monocat = 0 ! virtual ITD mono-category parameterizations (1, jpl = 1 only) or not (0) 102 102 ! 2: simple piling instead of ridging --- temporary option -
branches/2016/v3_6_CMIP6_ice_diagnostics/NEMOGCM/CONFIG/SHARED/namelist_ref
r7537 r7596 504 504 ! 1: "home made" based on Brandt et al. (JClim 2005) and Grenfell & Perovich (JGR 2004), 505 505 ! giving cloud-sky albedo 506 rn_alb_ drsno = 0.80! dry snow albedo : 0.80 (nn_ice_alb = 0); 0.85 (nn_ice_alb = 1); obs 0.85-0.87 (cloud-sky)507 rn_alb_ mlsno = 0.65 ! melting snow albedo : 0.65 ( '' ) ; 0.75 ( '' ) ; obs 0.72-0.82 ( '' )508 rn_alb_ drice = 0.72! dry ice albedo : 0.72 ( '' ) ; 0.60 ( '' ) ; obs 0.54-0.65 ( '' )509 rn_alb_ mlice = 0.53! bare puddled ice albedo : 0.53 ( '' ) ; 0.50 ( '' ) ; obs 0.49-0.58 ( '' )506 rn_alb_sdry = 0.85 ! dry snow albedo : 0.80 (nn_ice_alb = 0); 0.85 (nn_ice_alb = 1); obs 0.85-0.87 (cloud-sky) 507 rn_alb_smlt = 0.75 ! melting snow albedo : 0.65 ( '' ) ; 0.75 ( '' ) ; obs 0.72-0.82 ( '' ) 508 rn_alb_idry = 0.60 ! dry ice albedo : 0.72 ( '' ) ; 0.60 ( '' ) ; obs 0.54-0.65 ( '' ) 509 rn_alb_imlt = 0.50 ! bare puddled ice albedo : 0.53 ( '' ) ; 0.50 ( '' ) ; obs 0.49-0.58 ( '' ) 510 510 / 511 511 !----------------------------------------------------------------------- -
branches/2016/v3_6_CMIP6_ice_diagnostics/NEMOGCM/NEMO/LIM_SRC_3/ice.F90
r7537 r7596 212 212 REAL(wp), PUBLIC :: rn_betas !: coef. for partitioning of snowfall between leads and sea ice 213 213 REAL(wp), PUBLIC :: rn_kappa_i !: coef. for the extinction of radiation Grenfell et al. (2006) [1/m] 214 REAL(wp), PUBLIC :: rn_ rcdsn!: thermal conductivity of the snow [W/m/K]214 REAL(wp), PUBLIC :: rn_cdsn !: thermal conductivity of the snow [W/m/K] 215 215 REAL(wp), PUBLIC :: nn_conv_dif !: maximal number of iterations for heat diffusion 216 216 REAL(wp), PUBLIC :: rn_terr_dif !: maximal tolerated error (C) for heat diffusion -
branches/2016/v3_6_CMIP6_ice_diagnostics/NEMOGCM/NEMO/LIM_SRC_3/limthd.F90
r7537 r7596 645 645 NAMELIST/namicethd/ rn_hnewice, ln_frazil, rn_maxfrazb, rn_vfrazb, rn_Cfrazb, & 646 646 & rn_himin, rn_betas, rn_kappa_i, nn_conv_dif, rn_terr_dif, nn_ice_thcon, & 647 & rn_ rcdsn, nn_monocat, ln_it_qnsice647 & rn_cdsn, nn_monocat, ln_it_qnsice 648 648 !!------------------------------------------------------------------- 649 649 ! … … 684 684 WRITE(numout,*)' maximal err. on T for heat diffusion computation rn_terr_dif = ', rn_terr_dif 685 685 WRITE(numout,*)' switch for comp. of thermal conductivity in the ice nn_ice_thcon = ', nn_ice_thcon 686 WRITE(numout,*)' thermal conductivity of the snow rn_ rcdsn = ', rn_rcdsn686 WRITE(numout,*)' thermal conductivity of the snow rn_cdsn = ', rn_cdsn 687 687 WRITE(numout,*)' check heat conservation in the ice/snow con_i = ', con_i 688 688 WRITE(numout,*)' virtual ITD mono-category parameterizations (1) or not nn_monocat = ', nn_monocat -
branches/2016/v3_6_CMIP6_ice_diagnostics/NEMOGCM/NEMO/LIM_SRC_3/limthd_dif.F90
r7537 r7596 376 376 377 377 ! Effective thickness he (zhe) 378 zfac = 1._wp / ( rn_ rcdsn + zkimean )379 zratio_s = rn_ rcdsn * zfac378 zfac = 1._wp / ( rn_cdsn + zkimean ) 379 zratio_s = rn_cdsn * zfac 380 380 zratio_i = zkimean * zfac 381 381 zhe = zratio_s * ht_i_1d(ji) + zratio_i * ht_s_1d(ji) … … 400 400 DO ji = kideb, kiut 401 401 zfac = 1. / MAX( epsi10 , zh_s(ji) ) 402 zkappa_s(ji,0) = zghe(ji) * rn_ rcdsn * zfac403 zkappa_s(ji,nlay_s) = zghe(ji) * rn_ rcdsn * zfac402 zkappa_s(ji,0) = zghe(ji) * rn_cdsn * zfac 403 zkappa_s(ji,nlay_s) = zghe(ji) * rn_cdsn * zfac 404 404 END DO 405 405 406 406 DO jk = 1, nlay_s-1 407 407 DO ji = kideb , kiut 408 zkappa_s(ji,jk) = zghe(ji) * 2.0 * rn_ rcdsn / MAX( epsi10, 2.0 * zh_s(ji) )408 zkappa_s(ji,jk) = zghe(ji) * 2.0 * rn_cdsn / MAX( epsi10, 2.0 * zh_s(ji) ) 409 409 END DO 410 410 END DO … … 422 422 zkappa_i(ji,0) = zghe(ji) * ztcond_i(ji,0) * zfac 423 423 zkappa_i(ji,nlay_i) = zghe(ji) * ztcond_i(ji,nlay_i) * zfac 424 zkappa_s(ji,nlay_s) = zghe(ji) * zghe(ji) * 2.0 * rn_ rcdsn * ztcond_i(ji,0) / &425 & MAX( epsi10, ( zghe(ji) * ztcond_i(ji,0) * zh_s(ji) + zghe(ji) * rn_ rcdsn * zh_i(ji) ) )424 zkappa_s(ji,nlay_s) = zghe(ji) * zghe(ji) * 2.0 * rn_cdsn * ztcond_i(ji,0) / & 425 & MAX( epsi10, ( zghe(ji) * ztcond_i(ji,0) * zh_s(ji) + zghe(ji) * rn_cdsn * zh_i(ji) ) ) 426 426 zkappa_i(ji,0) = zkappa_s(ji,nlay_s) * isnow(ji) + zkappa_i(ji,0) * ( 1._wp - isnow(ji) ) 427 427 END DO … … 738 738 !--- Snow-ice interfacial temperature (diagnostic SIMIP) 739 739 DO ji = kideb, kiut 740 zfac = 1. / MAX( epsi10 , rn_ rcdsn * zh_i(ji) + ztcond_i(ji,1) * zh_s(ji) )741 t_si_1d(ji) = ( rn_ rcdsn * zh_i(ji) * t_s_1d(ji,1) + &740 zfac = 1. / MAX( epsi10 , rn_cdsn * zh_i(ji) + ztcond_i(ji,1) * zh_s(ji) ) 741 t_si_1d(ji) = ( rn_cdsn * zh_i(ji) * t_s_1d(ji,1) + & 742 742 & ztcond_i(ji,1) * zh_s(ji) * t_i_1d(ji,1) ) * zfac 743 743 END DO -
branches/2016/v3_6_CMIP6_ice_diagnostics/NEMOGCM/NEMO/OPA_SRC/SBC/albedo.F90
r7537 r7596 39 39 ! !!* namelist namsbc_alb 40 40 INTEGER :: nn_ice_alb 41 REAL(wp) :: rn_alb_ drsno, rn_alb_mlsno, rn_alb_drice, rn_alb_mlice41 REAL(wp) :: rn_alb_sdry, rn_alb_smlt, rn_alb_idry, rn_alb_imlt 42 42 43 43 !!---------------------------------------------------------------------- … … 101 101 IF( albd_init == 0 ) CALL albedo_init ! initialization 102 102 103 ralb_sf = rn_alb_ drsno! dry snow104 ralb_sm = rn_alb_ mlsno! melting snow105 ralb_if = rn_alb_ drice! bare frozen ice106 ralb_im = rn_alb_ mlice! bare puddled ice103 ralb_sf = rn_alb_sdry ! dry snow 104 ralb_sm = rn_alb_smlt ! melting snow 105 ralb_if = rn_alb_idry ! bare frozen ice 106 ralb_im = rn_alb_imlt ! bare puddled ice 107 107 108 108 SELECT CASE ( nn_ice_alb ) … … 252 252 !!---------------------------------------------------------------------- 253 253 INTEGER :: ios ! Local integer output status for namelist read 254 NAMELIST/namsbc_alb/ nn_ice_alb, rn_alb_ drsno, rn_alb_mlsno, rn_alb_drice , rn_alb_mlice254 NAMELIST/namsbc_alb/ nn_ice_alb, rn_alb_sdry, rn_alb_smlt, rn_alb_idry , rn_alb_imlt 255 255 !!---------------------------------------------------------------------- 256 256 ! … … 272 272 WRITE(numout,*) ' Namelist namsbc_alb : albedo ' 273 273 WRITE(numout,*) ' choose the albedo parameterization nn_ice_alb = ', nn_ice_alb 274 WRITE(numout,*) ' albedo of dry snow rn_alb_ drsno = ', rn_alb_drsno275 WRITE(numout,*) ' albedo of melting snow rn_alb_ mlsno = ', rn_alb_mlsno276 WRITE(numout,*) ' albedo of dry ice rn_alb_ drice = ', rn_alb_drice277 WRITE(numout,*) ' albedo of bare puddled ice rn_alb_ mlice = ', rn_alb_mlice274 WRITE(numout,*) ' albedo of dry snow rn_alb_sdry = ', rn_alb_sdry 275 WRITE(numout,*) ' albedo of melting snow rn_alb_smlt = ', rn_alb_smlt 276 WRITE(numout,*) ' albedo of dry ice rn_alb_idry = ', rn_alb_idry 277 WRITE(numout,*) ' albedo of bare puddled ice rn_alb_imlt = ', rn_alb_imlt 278 278 ENDIF 279 279 !
Note: See TracChangeset
for help on using the changeset viewer.