Changeset 5099 for LMDZ6/branches/Amaury_dev/libf/phylmd/surf_ocean_mod.F90
- Timestamp:
- Jul 22, 2024, 9:29:09 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/surf_ocean_mod.F90
r5093 r5099 1 ! 1 2 2 ! $Id$ 3 ! 3 4 4 MODULE surf_ocean_mod 5 5 … … 7 7 8 8 CONTAINS 9 ! 9 10 10 !****************************************************************************** 11 ! 11 12 12 SUBROUTINE surf_ocean(rlon, rlat, swnet, lwnet, alb1, & 13 13 windsp, rmu0, fder, tsurf_in, & … … 45 45 USE limit_read_mod 46 46 USE config_ocean_skin_m, ONLY: activate_ocean_skin 47 ! 47 48 48 ! This subroutine will make a call to ocean_XXX_noice according to the ocean mode (force, 49 49 ! slab or couple). The calculations of albedo and rugosity for the ocean surface are … … 196 196 !****************************************************************************** 197 197 ! Calculate total net radiance at surface 198 ! 198 199 199 !****************************************************************************** 200 200 radsol(1:klon) = 0.0 ! initialisation a priori inutile … … 294 294 IF (iflag_albedo==0) THEN 295 295 !--old parametrizations of ocean surface albedo 296 ! 296 297 297 IF (iflag_cycle_diurne>=1) THEN 298 ! 298 299 299 CALL alboc_cd(rmu0,alb_eau) 300 ! 300 301 301 !--ad-hoc correction for model radiative balance tuning 302 302 !--now outside alboc_cd routine 303 303 alb_eau(1:klon) = fmagic*alb_eau(1:klon) + pmagic 304 304 alb_eau(1:klon)=MIN(MAX(alb_eau(1:klon),0.0),1.0) 305 ! 305 306 306 ELSE 307 ! 307 308 308 CALL alboc(REAL(jour),rlat,alb_eau) 309 309 !--ad-hoc correction for model radiative balance tuning … … 311 311 alb_eau(1:klon) = fmagic*alb_eau(1:klon) + pmagic 312 312 alb_eau(1:klon)=MIN(MAX(alb_eau(1:klon),0.04),0.60) 313 ! 313 314 314 ENDIF 315 ! 315 316 316 DO i =1, knon 317 317 DO k=1,nsw … … 323 323 alb_dif_new(1:knon,:)=alb_dir_new(1:knon,:) 324 324 !IM 09122015 end 325 ! 325 326 326 ELSE IF (iflag_albedo==1) THEN 327 327 !--new parametrization of ocean surface albedo by Sunghye Baek 328 328 !--albedo for direct and diffuse radiation are different 329 ! 329 330 330 CALL ocean_albedo(knon,rmu0,knindex,windsp,SFRWL,alb_dir_new,alb_dif_new) 331 ! 331 332 332 !--ad-hoc correction for model radiative balance tuning 333 333 alb_dir_new(1:knon,:) = fmagic*alb_dir_new(1:knon,:) + pmagic … … 335 335 alb_dir_new(1:knon,:)=MIN(MAX(alb_dir_new(1:knon,:),0.0),1.0) 336 336 alb_dif_new(1:knon,:)=MIN(MAX(alb_dif_new(1:knon,:),0.0),1.0) 337 ! 337 338 338 ELSE IF (iflag_albedo==2) THEN 339 339 ! F. Codron albedo read from limit.nc … … 399 399 END SUBROUTINE surf_ocean 400 400 !**************************************************************************** 401 ! 401 402 402 END MODULE surf_ocean_mod
Note: See TracChangeset
for help on using the changeset viewer.