Changeset 1498 for trunk/LMDZ.GENERIC/libf/phystd

Timestamp:

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

Author:

mturbet

Message:

spectral snow albedo added

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• albedo_calcv.F90 (added)
• albedo_snow_calc.F90 (added)
• callcorrk.F90 (modified) (1 diff)
• callkeys_mod.F90 (modified) (1 diff)
• inifis.F (modified) (1 diff)
• physiq.F90 (modified) (3 diffs)

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