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

LMDZ5/trunk/libf/phylmd/conf_phys_m.F90

@@ -110 +110 @@
     integer,SAVE        :: iflag_radia_omp
     integer,SAVE        :: iflag_rrtm_omp
+    integer,SAVE        :: iflag_albedo_omp !albedo SB
+    logical,save        :: ok_chlorophyll_omp ! albedo SB  
     integer,SAVE        :: NSW_omp
     integer,SAVE        :: iflag_cldth_omp, ip_ebil_phy_omp
@@ -889 +891 @@
     NSW_omp = 6
     call getin('NSW',NSW_omp)
+!albedo SB >>>
+    iflag_albedo_omp = 0
+    call getin('iflag_albedo',iflag_albedo_omp)
+
+    ok_chlorophyll_omp=.false.
+    call getin('ok_chlorophyll',ok_chlorophyll_omp)
+!albedo SB <<<
 
     !
@@ -2135 +2144 @@
     write(lunout,*)' iflag_rrtm = ', iflag_rrtm
     write(lunout,*)' NSW = ', NSW
+    write(lunout,*)' iflag_albedo = ', iflag_albedo !albedo SB
+    write(lunout,*)' ok_chlorophyll =',ok_chlorophyll ! albedo SB
     write(lunout,*)' iflag_ratqs = ', iflag_ratqs
     write(lunout,*)' seuil_inversion = ', seuil_inversion