- Timestamp:
- 2017-09-18T16:54:04+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icealb.F90
r8531 r8534 2 2 !!====================================================================== 3 3 !! *** MODULE icealb *** 4 !! Ocean forcing: bulk thermohaline forcing of theice4 !! Atmospheric forcing: Albedo over sea ice 5 5 !!===================================================================== 6 6 !! History : 4.0 ! 2017-07 (C. Rousset) Split ice and ocean albedos … … 8 8 #if defined key_lim3 9 9 !!---------------------------------------------------------------------- 10 !! 'key_lim3' LIM 3.0sea-ice model11 !!---------------------------------------------------------------------- 12 !! ice_alb : albedo for ice (clear and overcast skies)13 !! alb_init: initialisation of albedo computation10 !! 'key_lim3' ESIM sea-ice model 11 !!---------------------------------------------------------------------- 12 !! ice_alb : albedo for ice (clear and overcast skies) 13 !! ice_alb_init : initialisation of albedo computation 14 14 !!---------------------------------------------------------------------- 15 15 USE ice, ONLY: jpl ! sea-ice: number of categories … … 19 19 USE lib_mpp ! MPP library 20 20 USE lib_fortran ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined) 21 USE timing ! Timing 21 22 22 23 IMPLICIT NONE … … 109 110 REAL(wp), DIMENSION(jpi,jpj,jpl) :: zafrac_pnd ! relative ice fraction (effective) 110 111 !!--------------------------------------------------------------------- 111 112 ! 113 IF( nn_timing == 1 ) CALL timing_start('icealb') 114 ! 112 115 !----------------------------------------------------- 113 116 ! Snow-free albedo (no ice thickness dependence yet) … … 147 150 ENDIF 148 151 ! 149 !!gm: optimization ( rn_alb_smlt - rn_alb_imlt ) * r1_c2 can be computed one for all150 !!gm before the DO-loop below151 152 DO jl = 1, jpl 152 153 DO jj = 1, jpj … … 301 302 END SELECT 302 303 ! 304 IF( nn_timing == 1 ) CALL timing_stop('icealb') 305 ! 303 306 END SUBROUTINE ice_alb 304 307 … … 343 346 #else 344 347 !!---------------------------------------------------------------------- 345 !! Default option Dummy module NO LIM 3.0sea-ice model348 !! Default option Dummy module NO ESIM sea-ice model 346 349 !!---------------------------------------------------------------------- 347 350 #endif
Note: See TracChangeset
for help on using the changeset viewer.