Changeset 2227 for LMDZ5/trunk/libf/phylmd/change_srf_frac_mod.F90

Timestamp:

Mar 12, 2015, 12:07:43 PM (9 years ago)

Author:

Laurent Fairhead

Message:

New ocean albedo.

To activate the new scheme, put iflag_albedo=1 in physiq.def

To activate chlorophyll concentration effect on albedo,
put ok_chlorophyll=y in def file

and download file named chlorophyll.nc
chlorophyll.nc has the same dimension as the model grid with 12 months data,
(i=lon, j=lat, L=1:12) and can be degraded from the original file of dimension
i=1:4320 , j=1:2160 , L=1:12
ada:/workgpfs/rech/psl/rpsl949/clima/chlor_seasonal_clim_seawifs.nc

For 96X96 resolution, chlorophyll.nc file is in
ada:/workgpfs/rech/psl/rpsl949/clima/chlorophyll.nc

19. Baek

File:

• LMDZ5/trunk/libf/phylmd/change_srf_frac_mod.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/change_srf_frac_mod.F90

@@ -1 +1 @@
 !
-! $Header$
+! $Id$
 !
 MODULE change_srf_frac_mod
@@ -12 +12 @@
 
   SUBROUTINE change_srf_frac(itime, dtime, jour, &
-       pctsrf, alb1, alb2, tsurf, ustar, u10m, v10m, pbl_tke)
+!albedo SB >>>
+!       pctsrf, alb1, alb2, tsurf, ustar, u10m, v10m, pbl_tke)
+        pctsrf, alb_dir, alb_dif, tsurf, ustar, u10m, v10m, pbl_tke)
+!albedo SB <<<
+   
+
+
 !
 ! This subroutine is called from physiq.F at each timestep. 
@@ -32 +38 @@
     INCLUDE "iniprint.h"
     INCLUDE "YOMCST.h"
+!albedo SB >>>
+    include "clesphys.h"
+!albedo SB <<<
+
+
 
 ! Input arguments
@@ -43 +54 @@
    
     REAL, DIMENSION(klon,nbsrf), INTENT(INOUT) :: pctsrf ! sub-surface fraction
-    REAL, DIMENSION(klon,nbsrf), INTENT(INOUT) :: alb1   ! albedo first interval in SW spektrum
-    REAL, DIMENSION(klon,nbsrf), INTENT(INOUT) :: alb2   ! albedo second interval in SW spektrum
+!albedo SB >>>
+!   REAL, DIMENSION(klon,nbsrf), INTENT(INOUT) :: alb1   ! albedo first interval in SW spektrum
+!   REAL, DIMENSION(klon,nbsrf), INTENT(INOUT) :: alb2   ! albedo second interval in SW spektrum
+    REAL, DIMENSION(klon,nsw,nbsrf), INTENT(INOUT) :: alb_dir,alb_dif
+!albedo SB <<<
+
     REAL, DIMENSION(klon,nbsrf), INTENT(INOUT) :: tsurf
     REAL, DIMENSION(klon,nbsrf), INTENT(INOUT) :: ustar
@@ -160 +175 @@
 ! 
 !****************************************************************************************
-       CALL pbl_surface_newfrac(itime, pctsrf, pctsrf_old, tsurf, alb1, alb2, ustar, u10m, v10m, pbl_tke)
+
+!albedo SB >>>
+! CALL pbl_surface_newfrac(itime, pctsrf, pctsrf_old, tsurf, alb1, alb2, ustar,
+! u10m, v10m, pbl_tke)
+       CALL pbl_surface_newfrac(itime, pctsrf, pctsrf_old, tsurf, alb_dir,alb_dif, ustar, u10m, v10m, pbl_tke)
+!albedo SB <<<
+
+
 
     ELSE