Changeset 2227 for LMDZ5/trunk/libf/phylmd/phys_state_var_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/phys_state_var_mod.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/phys_state_var_mod.F90

@@ -30 +30 @@
       REAL, ALLOCATABLE, SAVE :: falb1(:,:), falb2(:,:)
 !$OMP THREADPRIVATE(falb1, falb2)
+
+!albedo SB >>>
+      REAL, ALLOCATABLE, SAVE :: falb_dif(:,:,:), falb_dir(:,:,:)
+      real, allocatable, save :: chl_con(:)
+!$OMP THREADPRIVATE(falb_dir,falb_dif,chl_con)
+!albedo SB <<<
+
+
       REAL, ALLOCATABLE, SAVE :: rain_fall(:), snow_fall(:)
 !$OMP THREADPRIVATE( rain_fall, snow_fall)
@@ -261 +269 @@
 !$OMP THREADPRIVATE(albsol1,albsol2)
 
+!albedo SB >>>
+      REAL,ALLOCATABLE,SAVE :: albsol_dif(:,:),albsol_dir(:,:)
+!$OMP THREADPRIVATE(albsol_dif,albsol_dir)
+!albedo SB <<<
+
+
       REAL, ALLOCATABLE, SAVE:: wo(:, :, :)
       ! column-density of ozone in a layer, in kilo-Dobsons
@@ -404 +418 @@
       ALLOCATE(falb1(klon,nbsrf))
       ALLOCATE(falb2(klon,nbsrf))
+!albedo SB >>> 
+      ALLOCATE(falb_dir(klon,nsw,nbsrf),falb_dif(klon,nsw,nbsrf))
+      ALLOCATE(chl_con(klon))
+!albedo SB <<<
       ALLOCATE(rain_fall(klon))
       ALLOCATE(snow_fall(klon))
@@ -501 +519 @@
       ALLOCATE(paire_ter(klon))
       ALLOCATE(albsol1(klon), albsol2(klon))
+!albedo SB >>>
+      ALLOCATE(albsol_dir(klon,nsw),albsol_dif(klon,nsw))
+!albedo SB <<<
 
       if (read_climoz <= 1) then
@@ -634 +655 @@
       deallocate(paire_ter)
       deallocate(albsol1, albsol2)
+!albedo SB >>>
+      deallocate(albsol_dir,albsol_dif,falb_dir,falb_dif,chl_con)
+!albedo SB <<<
       deallocate(wo)
       deallocate(clwcon0,rnebcon0)