Changeset 1498 for trunk/LMDZ.GENERIC

Timestamp:

Nov 13, 2015, 7:21:01 PM (9 years ago)

Author:

mturbet

Message:

spectral snow albedo added

Location:

trunk/LMDZ.GENERIC

Files:

• README (modified) (1 diff)
• deftank/callphys.GJ581d (modified) (2 diffs)
• deftank/callphys.earlymars (modified) (2 diffs)
• deftank/callphys.earth (modified) (2 diffs)
• deftank/earth1d_Dry.callphys.def (modified) (2 diffs)
• libf/phystd/albedo_calcv.F90 (added)
• libf/phystd/albedo_snow_calc.F90 (added)
• libf/phystd/callcorrk.F90 (modified) (1 diff)
• libf/phystd/callkeys_mod.F90 (modified) (1 diff)
• libf/phystd/inifis.F (modified) (1 diff)
• libf/phystd/physiq.F90 (modified) (3 diffs)

trunk/LMDZ.GENERIC/README

@@ -1106 +1106 @@
 == 05/11/2015 == MT
 - Minor Albedo correction of bugs in rcm1d.
+
+== 13/11/2015 == MT
+- Spectral Snow Albebdo implemented in the code.

trunk/LMDZ.GENERIC/deftank/callphys.GJ581d

@@ -126 +126 @@
 # Include hydrology ?
 hydrology     = .true.
+# Spectral Dependant Albedo ?
+albedo_spectral_mode=.false.
+# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+# If albedo_spectral_mode=.true., albedosnow becomes the 0.5 micron snow albedo.
+#
+# albedosnow = 0.95  (0.73 Sun-integrated) for fresh snow.
+#            = 0.50  (0.39 Sun-integrated) for dirty snow.
+#            = 0.645 (0.50 Sun-integrated) for 'realistic' snow.
+# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 # H2O snow (and ice) albedo ?
-albedosnow    = 0.5
+albedosnow    = 0.50
 # Maximum sea ice thickness ?
 maxicethick   = 0.05
@@ -135 +144 @@
 ## CO2 options 
 ## ~~~~~~~~~~~
+# Co2 ice albedo ?
+albedoco2ice   = 0.5
 # gas is non-ideal CO2 ?
 nonideal      = .false.

trunk/LMDZ.GENERIC/deftank/callphys.earlymars

@@ -133 +133 @@
 # Include hydrology ?
 hydrology     = .true.
+# Spectral Dependant Albedo ?
+albedo_spectral_mode=.false.
+# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+# If albedo_spectral_mode=.true., albedosnow becomes the 0.5 micron snow albedo.
+#
+# albedosnow = 0.95  (0.73 Sun-integrated) for fresh snow.
+#            = 0.50  (0.39 Sun-integrated) for dirty snow.
+#            = 0.645 (0.50 Sun-integrated) for 'realistic' snow.
+# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 # H2O snow (and ice) albedo ?
-albedosnow    = 0.5
+albedosnow    = 0.50
 # Maximum sea ice thickness ?
 maxicethick   = 0.05
@@ -146 +155 @@
 ## CO2 options 
 ## ~~~~~~~~~~~
+# Co2 ice albedo ?
+albedoco2ice   = 0.5
 # gas is non-ideal CO2 ?
 nonideal      = .false.

trunk/LMDZ.GENERIC/deftank/callphys.earth

@@ -132 +132 @@
 # WATER: hydrology ?
 hydrology     = .true.
+# Spectral Dependant Albedo ?
+albedo_spectral_mode=.false.
+# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+# If albedo_spectral_mode=.true., albedosnow becomes the 0.5 micron snow albedo.
+#
+# albedosnow = 0.95  (0.73 Sun-integrated) for fresh snow.
+#            = 0.50  (0.39 Sun-integrated) for dirty snow.
+#            = 0.78  (0.60 Sun-integrated) for 'realistic' snow.
+# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 # H2O snow (and ice) albedo ?
-albedosnow    = 0.6
+albedosnow    = 0.60
 # Maximum sea ice thickness ?
 maxicethick   = 10.
@@ -143 +152 @@
 ## CO2 options 
 ## ~~~~~~~~~~~
+# Co2 ice albedo ?
+albedoco2ice   = 0.5
 # gas is non-ideal CO2 ?
 nonideal      = .false.

trunk/LMDZ.GENERIC/deftank/earth1d_Dry.callphys.def

@@ -167 +167 @@
 # Include hydrology ?
 hydrology     = .false.
+# Spectral Dependant Albedo ?
+albedo_spectral_mode=.false.
+# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+# If albedo_spectral_mode=.true., albedosnow becomes the 0.5 micron snow albedo.
+#
+# albedosnow = 0.95  (0.73 Sun-integrated) for fresh snow.
+#            = 0.50  (0.39 Sun-integrated) for dirty snow.
+#            = 0.645 (0.50 Sun-integrated) for 'realistic' snow.
+# ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 # H2O snow (and ice) albedo ?
-albedosnow    = 0.5
+albedosnow    = 0.50
 # Maximum sea ice thickness ?
 maxicethick   = 0.5
@@ -178 +187 @@
 ## CO2 options 
 ## ~~~~~~~~~~~
+# Co2 ice albedo ?
+albedoco2ice   = 0.5
 # call CO2 condensation ?
 co2cond       = .false.

trunk/LMDZ.GENERIC/libf/phystd/callcorrk.F90

@@ -231 +231 @@
          call sugas_corrk       ! Set up gaseous absorption properties.
          call suaer_corrk       ! Set up aerosol optical properties.
-
+        
 
          if((igcm_h2o_vap.eq.0) .and. varactive)then

trunk/LMDZ.GENERIC/libf/phystd/callkeys_mod.F90

@@ -39 +39 @@
       logical ok_slab_sic
       logical ok_slab_heat_transp
+      logical albedo_spectral_mode
 
       integer iddist

trunk/LMDZ.GENERIC/libf/phystd/inifis.F

@@ -614 +614 @@
          call getin_p("icetstep",icetstep)
          write(*,*) " icetstep = ",icetstep
+         
+         write(*,*) "Spectral Dependant albedo ?"
+         albedo_spectral_mode=.false. ! default value
+         call getin_p("albedo_spectral_mode",albedo_spectral_mode)
+         write(*,*) " albedo_spectral_mode = ",albedo_spectral_mode
 
          write(*,*) "Snow albedo ?"
+         write(*,*) "If albedo_spectral_mode=.true., then this "
+         write(*,*) "corresponds to the 0.5 microns snow albedo."
          albedosnow=0.5         ! default value
          call getin_p("albedosnow",albedosnow)

trunk/LMDZ.GENERIC/libf/phystd/physiq.F90

@@ -539 +539 @@
          call surfini(ngrid,nq,qsurf,albedo,albedo_bareground,albedo_snow_SPECTV,albedo_co2_ice_SPECTV) 
          
-         
 !        Initialize orbital calculation. 
 !        ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
@@ -928 +927 @@
                ! Net atmospheric radiative heating rate (K.s-1)
                dtrad(1:ngrid,1:nlayer)=zdtsw(1:ngrid,1:nlayer)+zdtlw(1:ngrid,1:nlayer)
+               
+               ! Late initialization of the Ice Spectral Albedo. We needed the visible bands to do that ! 
+               if (firstcall .and. albedo_spectral_mode) then
+                  call spectral_albedo_calc(albedo_snow_SPECTV)
+               endif
 
             elseif(newtonian)then
@@ -951 +955 @@
                print*,'------------WARNING---WARNING------------' ! by MT2015.
                print*,'You are in corrk=false mode, '
-               print*,'and the surface albedo is taken equal to its value at 0.5 micrometers ...'
-               if ( (10000./BWNV(1) .le. 0.5) .or. (10000./BWNV(L_NSPECTV) .ge. 0.5) ) then
-                  print*,'0.5 microns band doesnt match your visible bands ! Abort !'
-                  call abort
-               endif
-               DO nw=1, L_NSPECTV-1
-                  if ( (10000./BWNV(nw) .ge. 0.5) .and. (10000./BWNV(nw+1) .lt. 0.5) ) then
-                     fluxsurfabs_sw(1:ngrid) = fluxtop_dn(1:ngrid)*(1.-albedo(1:ngrid,nw))
-                     fluxrad_sky(1:ngrid)    = fluxsurfabs_sw(1:ngrid)
-                  endif
-               ENDDO
+               print*,'and the surface albedo is taken equal to the first visible spectral value'               
+               
+               fluxsurfabs_sw(1:ngrid) = fluxtop_dn(1:ngrid)*(1.-albedo(1:ngrid,1))
+               fluxrad_sky(1:ngrid)    = fluxsurfabs_sw(1:ngrid)
                fluxtop_lw(1:ngrid)  = emis(1:ngrid)*sigma*tsurf(1:ngrid)**4