- Timestamp:
- 2017-09-06T19:02:43+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icealb.F90
r8498 r8505 13 13 !! alb_init : initialisation of albedo computation 14 14 !!---------------------------------------------------------------------- 15 USE ice , ONLY : jpl ! number of ice category15 USE ice, ONLY: jpl ! sea-ice: number of categories 16 16 USE phycst ! physical constants 17 17 ! … … 23 23 PRIVATE 24 24 25 PUBLIC ice_alb ! routine called in iceforcing.F90 and iceupdate.F90 25 PUBLIC ice_alb_init ! called in icestp 26 PUBLIC ice_alb ! called in iceforcing.F90 and iceupdate.F90 26 27 27 28 REAL(wp), PUBLIC, PARAMETER :: rn_alb_oce = 0.066 !: ocean or lead albedo (Pegau and Paulson, Ann. Glac. 2001) 28 29 29 INTEGER :: albd_init = 0 ! control flag for initialization30 30 REAL(wp) , PARAMETER :: rc1 = 0.05 ! snow thickness (only for nn_ice_alb=0) 31 31 REAL(wp) , PARAMETER :: rc2 = 0.10 ! " " … … 33 33 REAL(wp) , PARAMETER :: r1_c1 = 1. / rc1 34 34 REAL(wp) , PARAMETER :: r1_c2 = 1. / rc2 35 36 ! !!* namelist namsbc_alb * 37 INTEGER :: nn_ice_alb 38 REAL(wp) :: rn_alb_sdry, rn_alb_smlt, rn_alb_idry, rn_alb_imlt, rn_alb_dpnd 35 ! 36 ! ** albedo namelist (namicealb) 37 INTEGER :: nn_ice_alb ! type of albedo scheme: 0: Shine & Henderson-Sellers (JGR 1985) 38 ! ! 1: "home made" based on Brandt et al. (JClim 2005) 39 ! ! and Grenfell & Perovich (JGR 2004) 40 ! ! 2: Same as 1 but with melt ponds 41 REAL(wp) :: rn_alb_sdry ! dry snow albedo 42 REAL(wp) :: rn_alb_smlt ! melting snow albedo 43 REAL(wp) :: rn_alb_idry ! dry ice albedo 44 REAL(wp) :: rn_alb_imlt ! bare puddled ice albedo 45 REAL(wp) :: rn_alb_dpnd ! ponded ice albedo 39 46 40 47 !!---------------------------------------------------------------------- … … 102 109 REAL(wp), DIMENSION(jpi,jpj,jpl) :: zafrac_pnd ! relative ice fraction (effective) 103 110 !!--------------------------------------------------------------------- 104 105 IF( albd_init == 0 ) CALL alb_init ! initialization106 111 107 112 !----------------------------------------------------- … … 299 304 300 305 301 SUBROUTINE alb_init306 SUBROUTINE ice_alb_init 302 307 !!---------------------------------------------------------------------- 303 308 !! *** ROUTINE alb_init *** … … 311 316 NAMELIST/namicealb/ nn_ice_alb, rn_alb_sdry, rn_alb_smlt, rn_alb_idry, rn_alb_imlt, rn_alb_dpnd 312 317 !!---------------------------------------------------------------------- 313 !314 albd_init = 1 ! indicate that the initialization has been done315 318 ! 316 319 REWIND( numnam_ice_ref ) ! Namelist namicealb in reference namelist : Albedo parameters … … 336 339 ENDIF 337 340 ! 338 END SUBROUTINE alb_init341 END SUBROUTINE ice_alb_init 339 342 340 343 #else
Note: See TracChangeset
for help on using the changeset viewer.