Changeset 2413 for LMDZ5/trunk/libf/phylmd/surf_ocean_mod.F90

Timestamp:

Dec 18, 2015, 8:27:39 PM (8 years ago)

Author:

oboucher

Message:

limiting iflag_rrtm=0 to NSW=2 case (NSW=4 and 6 were possible with slightly different averaging procedures)
limiting iflag_rrtm choices to 0 and 1
limiting iflag_albedo choices to 0 and 1
cleaning up surface ocean albedo parametrisations
fmagic and pmagic corrections are now done in surf_ocean_albedo.F90
fmagic and pmagic corrections also possible for iflag_rrtm=1 now

File:

• LMDZ5/trunk/libf/phylmd/surf_ocean_mod.F90 (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/surf_ocean_mod.F90

@@ -175 +175 @@
 
 !******************************************************************************
-! Calculate albedo
+! Calculate ocean surface albedo
 !******************************************************************************
 !albedo SB >>>
-  if(iflag_albedo==1)then
-    call ocean_albedo(knon,rmu0,knindex,windsp,SFRWL,alb_dir_new,alb_dif_new)
-  else
+IF (iflag_albedo==0) THEN
+!--old parametrizations of ocean surface albedo
+!
     IF (cycle_diurne) THEN
+!
        CALL alboc_cd(rmu0,alb_eau)
+!
+!--ad-hoc correction for model radiative balance tuning
+!--now outside alboc_cd routine
+       alb_eau(:) = fmagic*alb_eau(:) + pmagic
+       alb_eau=MIN(MAX(alb_eau,0.0),1.0)
+!
     ELSE
+!
        CALL alboc(REAL(jour),rlat,alb_eau)
+!--ad-hoc correction for model radiative balance tuning
+!--now outside alboc routine
+       alb_eau(:) = fmagic*alb_eau(:) + pmagic
+       alb_eau=MIN(MAX(alb_eau(i),0.04),0.60)
+!
     ENDIF
-
+!
     DO i =1, knon
-      do  k=1,nsw
+      DO  k=1,nsw
        alb_dir_new(i,k) = alb_eau(knindex(i))
-      enddo
+      ENDDO
     ENDDO
-!IM 09122015 beg
-!IM next line corresponds to Sunghye Baek change that she made when she added new ocean albedo
-!    alb_dif_new=0.05 !alb_dir_new
-!IM next line corresponds to the old values of LMDZ5A/IPSLCM5A versions 
+!IM 09122015 next line corresponds to the old way of doing in LMDZ5A/IPSLCM5A versions 
+!albedo for diffuse radiation is taken the same as for direct radiation
      alb_dif_new=alb_dir_new
 !IM 09122015 end
-endif
-
+!
+ELSE IF (iflag_albedo==1) THEN
+!--new parametrization of ocean surface albedo by Sunghye Baek
+!--albedo for direct and diffuse radiation are different
+!
+    CALL ocean_albedo(knon,rmu0,knindex,windsp,SFRWL,alb_dir_new,alb_dif_new)
+!
+!--ad-hoc correction for model radiative balance tuning
+    alb_dir_new(:,:) = fmagic*alb_dir_new(:,:) + pmagic
+    alb_dir_new=MIN(MAX(alb_dir_new,0.0),1.0)
+    alb_dif_new=MIN(MAX(alb_dif_new,0.0),1.0)
+!
+ENDIF
 !albedo SB <<<