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

Timestamp:

Oct 5, 2017, 12:15:17 PM (7 years ago)

Author:

jvatant

Message:

+ Found a bug for the limit cases in disr_haze leading to uninit. value.

Fixed by init to zero but should be relooked at.

+ Also added omp instructions in disr_haze
--JVO

File:

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

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

@@ -94 +94 @@
       indv = -9999 ! this initial value means "to be calculated"
   ENDIF
+  
+  ! Some initialisation beacause there's a pb with disr_haze at the limits (nw=1)
+  ! I should check this - For now we set vars to zero : better than nans - JVO 2017
+  
+  dhaze_t(:,:) = 0.
+  ssa_t(:,:) = 0.
+  asf_t(:,:) = 0.
+
 
   !=======================================================================
@@ -135 +143 @@
 
         !================= Titan customisation ========================================
-        call disr_haze(dz(k),plev(k),wnov(nw),dhaze_T(k,nw),SSA_T(k,nw),ASF_T(k,nw))
+        call disr_haze(dz(k),plev(k),wnov(nw),dhaze_T(k,nw),SSA_T(k,nw),ASF_T(k,nw))
         ! =============================================================================
-
+         
         DRAYAER = TRAY(K,NW)
         !     DRAYAER is Tau RAYleigh scattering, plus AERosol opacity