Changeset 1788 for trunk/LMDZ.TITAN/libf/phytitan/optcv.F90

Timestamp:

Sep 26, 2017, 12:43:42 PM (7 years ago)

Author:

jvatant

Message:

Get rid of all the old-generic dummy aerosol scheme ( just left scatterers_h for compilation )
as the new microphysics for Titan will be plugged in
-> even removed sedimentation ( will be done in the microphysical model )
--JVO

File:

• trunk/LMDZ.TITAN/libf/phytitan/optcv.F90 (modified) (8 diffs)

trunk/LMDZ.TITAN/libf/phytitan/optcv.F90

@@ -1 +1 @@
 SUBROUTINE OPTCV(DTAUV,TAUV,TAUCUMV,PLEV,  &
-     QXVAER,QSVAER,GVAER,WBARV,COSBV,       &
-     TAURAY,TAUAERO,TMID,PMID,TAUGSURF)
+     WBARV,COSBV,TAURAY,TMID,PMID,TAUGSURF)
 
   use radinc_h
@@ -46 +45 @@
   real*8 COSBV(L_NLAYRAD,L_NSPECTV,L_NGAUSS)
   real*8 WBARV(L_NLAYRAD,L_NSPECTV,L_NGAUSS)
-
-  ! for aerosols
-  real*8  QXVAER(L_LEVELS,L_NSPECTV,NAERKIND)
-  real*8  QSVAER(L_LEVELS,L_NSPECTV,NAERKIND)
-  real*8  GVAER(L_LEVELS,L_NSPECTV,NAERKIND)
-  real*8  TAUAERO(L_LEVELS,NAERKIND)
-  real*8  TAUAEROLK(L_LEVELS,L_NSPECTV,NAERKIND)
-  real*8  TAEROS(L_LEVELS,L_NSPECTV,NAERKIND)
 
   ! Titan customisation
@@ -76 +67 @@
 
   real*8 taugsurf(L_NSPECTV,L_NGAUSS-1)
-  real*8 DCONT,DAERO
+  real*8 DCONT
   real*8 DRAYAER
   double precision wn_cont, p_cont, p_air, T_cont, dtemp, dtempc
@@ -129 +120 @@
   end do                    ! levels
 
-  ! Spectral dependance of aerosol absorption
-  do iaer=1,naerkind
-     do NW=1,L_NSPECTV
-        do K=2,L_LEVELS
-           TAEROS(K,NW,IAER) = TAUAERO(K,IAER) * QXVAER(K,NW,IAER)
-        end do                    ! levels
-     end do
-  end do
-  
   ! Rayleigh scattering
   do NW=1,L_NSPECTV
@@ -158 +140 @@
         DRAYAER = TRAY(K,NW)
         !     DRAYAER is Tau RAYleigh scattering, plus AERosol opacity
-        do iaer=1,naerkind
-           DRAYAER = DRAYAER + TAEROS(K,NW,IAER)
-        end do
-
         DRAYAER = DRAYAER + DHAZE_T(K,NW) ! Titan's aerosol
 
@@ -275 +253 @@
   !     we need to calculate the scattering albedo and asymmetry factors
 
-  do iaer=1,naerkind
-    DO NW=1,L_NSPECTV
-      DO K=2,L_LEVELS
-           TAUAEROLK(K,NW,IAER) = TAUAERO(K,IAER) * QSVAER(K,NW,IAER) ! effect of scattering albedo
-      ENDDO
-    ENDDO
-  end do
-  
   ! Haze scattering
   DO NW=1,L_NSPECTV
@@ -294 +264 @@
      DO L=1,L_NLAYRAD-1
         K              = 2*L+1
-        atemp(L,NW) = SUM(GVAER(K,NW,1:naerkind) * TAUAEROLK(K,NW,1:naerkind))+SUM(GVAER(K+1,NW,1:naerkind) * TAUAEROLK(K+1,NW,1:naerkind)) &
-                    + ASF_T(K,NW)*DHAZES_T(K,NW) + ASF_T(K+1,NW)*DHAZES_T(K+1,NW)
-        btemp(L,NW) = SUM(TAUAEROLK(K,NW,1:naerkind)) + SUM(TAUAEROLK(K+1,NW,1:naerkind)) &
-                    + DHAZES_T(K,NW) + DHAZES_T(K+1,NW)
+        atemp(L,NW) = ASF_T(K,NW)*DHAZES_T(K,NW) + ASF_T(K+1,NW)*DHAZES_T(K+1,NW)
+        btemp(L,NW) = DHAZES_T(K,NW) + DHAZES_T(K+1,NW)
         ctemp(L,NW) = btemp(L,NW) + 0.9999*(TRAY(K,NW) + TRAY(K+1,NW)) ! JVO 2017 : does this 0.999 is really meaningful ?
         btemp(L,NW) = btemp(L,NW) + TRAY(K,NW) + TRAY(K+1,NW)
@@ -306 +274 @@
      L           = L_NLAYRAD
      K           = 2*L+1
-     atemp(L,NW) = SUM(GVAER(K,NW,1:naerkind) * TAUAEROLK(K,NW,1:naerkind)) &
-                 + ASF_T(K,NW)*DHAZES_T(K,NW)
-     btemp(L,NW) = SUM(TAUAEROLK(K,NW,1:naerkind)) &
-                 + DHAZES_T(K,NW)
+     atemp(L,NW) = ASF_T(K,NW)*DHAZES_T(K,NW)
+     btemp(L,NW) = DHAZES_T(K,NW)
      ctemp(L,NW) = btemp(L,NW) + 0.9999*TRAY(K,NW) ! JVO 2017 : does this 0.999 is really meaningful ?
      btemp(L,NW) = btemp(L,NW) + TRAY(K,NW)