Changeset 13284 for NEMO/releases/r4.0/r4.0-HEAD/src/ICE/icesbc.F90
- Timestamp:
- 2020-07-09T17:12:23+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/releases/r4.0/r4.0-HEAD/src/ICE/icesbc.F90
r11575 r13284 116 116 INTEGER :: ji, jj, jl ! dummy loop index 117 117 REAL(wp) :: zmiss_val ! missing value retrieved from xios 118 REAL(wp), DIMENSION(jpi,jpj,jpl) :: zalb_os, zalb_cs ! ice albedo under overcast/clear sky 119 REAL(wp), DIMENSION(:,:) , ALLOCATABLE :: zalb, zmsk00 ! 2D workspace 118 REAL(wp), DIMENSION(:,:), ALLOCATABLE :: zalb, zmsk00 ! 2D workspace 120 119 !!-------------------------------------------------------------------- 121 120 ! … … 131 130 CALL iom_miss_val( "icetemp", zmiss_val ) 132 131 133 ! --- cloud-sky and overcast-sky ice albedos --- ! 134 CALL ice_alb( t_su, h_i, h_s, ln_pnd_alb, a_ip_frac, h_ip, zalb_cs, zalb_os ) 135 136 ! albedo depends on cloud fraction because of non-linear spectral effects 137 !!gm cldf_ice is a real, DOCTOR naming rule: start with cd means CHARACTER passed in argument ! 138 alb_ice(:,:,:) = ( 1. - cldf_ice ) * zalb_cs(:,:,:) + cldf_ice * zalb_os(:,:,:) 139 ! 132 ! --- ice albedo --- ! 133 CALL ice_alb( t_su, h_i, h_s, ln_pnd_alb, a_ip_eff, h_ip, cloud_fra, alb_ice ) 134 140 135 ! 141 136 SELECT CASE( ksbc ) !== fluxes over sea ice ==! … … 281 276 INTEGER :: ios, ioptio ! Local integer 282 277 !! 283 NAMELIST/namsbc/ rn_cio, rn_blow_s, nn_flxdist, ln_cndflx, ln_cndemulate278 NAMELIST/namsbc/ rn_cio, nn_snwfra, rn_snwblow, nn_flxdist, ln_cndflx, ln_cndemulate, nn_qtrice 284 279 !!------------------------------------------------------------------- 285 280 ! … … 297 292 WRITE(numout,*) '~~~~~~~~~~~~~~~~' 298 293 WRITE(numout,*) ' Namelist namsbc:' 299 WRITE(numout,*) ' drag coefficient for oceanic stress rn_cio = ', rn_cio 300 WRITE(numout,*) ' coefficient for ice-lead partition of snowfall rn_blow_s = ', rn_blow_s 301 WRITE(numout,*) ' Multicategory heat flux formulation nn_flxdist = ', nn_flxdist 302 WRITE(numout,*) ' Use conduction flux as surface condition ln_cndflx = ', ln_cndflx 303 WRITE(numout,*) ' emulate conduction flux ln_cndemulate = ', ln_cndemulate 294 WRITE(numout,*) ' drag coefficient for oceanic stress rn_cio = ', rn_cio 295 WRITE(numout,*) ' fraction of ice covered by snow (options 0,1,2) nn_snwfra = ', nn_snwfra 296 WRITE(numout,*) ' coefficient for ice-lead partition of snowfall rn_snwblow = ', rn_snwblow 297 WRITE(numout,*) ' Multicategory heat flux formulation nn_flxdist = ', nn_flxdist 298 WRITE(numout,*) ' Use conduction flux as surface condition ln_cndflx = ', ln_cndflx 299 WRITE(numout,*) ' emulate conduction flux ln_cndemulate = ', ln_cndemulate 300 WRITE(numout,*) ' solar flux transmitted thru the surface scattering layer nn_qtrice = ', nn_qtrice 301 WRITE(numout,*) ' = 0 Grenfell and Maykut 1977' 302 WRITE(numout,*) ' = 1 Lebrun 2019' 304 303 ENDIF 305 304 !
Note: See TracChangeset
for help on using the changeset viewer.