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

Timestamp:

Feb 8, 2019, 5:45:36 PM (6 years ago)

Author:

jvatant

Message:

Update radiative transfer with some of the latest updates in the generic model
Cf r2056 by AB and r1987-1991 by JL
--JVO

File:

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

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

@@ -154 +154 @@
   ! Rayleigh scattering
   do NW=1,L_NSPECTV
-     do K=2,L_LEVELS
+     TRAY(1:4,NW)   = 1d-30
+     do K=5,L_LEVELS
         TRAY(K,NW)   = TAURAY(NW) * DPR(K)
      end do                    ! levels
@@ -189 +190 @@
           if (seashaze) dhaze_T(k,nw) = dhaze_T(k,nw)*seashazefact(k)
         ENDIF
+        
+        !JL18 It seems to be good to have aerosols in the first "radiative layer" of the gcm in the IR
+        !   but visible does not handle very well diffusion in first layer.
+        !   The tauaero and tauray are thus set to 0 (a small value for rayleigh because the code crashes otherwise)
+        !   in the 4 first semilayers in optcv, but not optci.
+        !   This solves random variations of the sw heating at the model top. 
+        if (k<5)  dhaze_T(K,:) = 0.0
          
         DRAYAER = TRAY(K,NW)
@@ -311 +319 @@
 
   ! Haze scattering
+            !JL18 It seems to be good to have aerosols in the first "radiative layer" of the gcm in the IR
+            !   but not in the visible
+            !   The dhaze_s is thus set to 0 in the 4 first semilayers in optcv, but not optci.
+            !   This solves random variations of the sw heating at the model top. 
   DO NW=1,L_NSPECTV
-    DO K=2,L_LEVELS
+      DHAZES_T(1:4,NW) = 0.d0
+    DO K=5,L_LEVELS
       DHAZES_T(K,NW) = DHAZE_T(K,NW) * SSA_T(K,NW) ! effect of scattering albedo on haze
     ENDDO
@@ -365 +378 @@
   DO NG=1,L_NGAUSS       ! full gauss loop
      DO NW=1,L_NSPECTV       
-        TAUV(1,NW,NG)=0.0D0
-        DO L=1,L_NLAYRAD
-           TAUV(L+1,NW,NG)=TAUV(L,NW,NG)+DTAUV(L,NW,NG)
-        END DO
-
         TAUCUMV(1,NW,NG)=0.0D0
         DO K=2,L_LEVELS
            TAUCUMV(K,NW,NG)=TAUCUMV(K-1,NW,NG)+DTAUKV(K,NW,NG)
         END DO
+
+        DO L=1,L_NLAYRAD
+           TAUV(L,NW,NG)=TAUCUMV(2*L,NW,NG)
+        END DO
+        TAUV(L,NW,NG)=TAUCUMV(2*L_NLAYRAD+1,NW,NG)
      END DO            
   END DO                 ! end full gauss loop