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

Timestamp:

Oct 17, 2023, 9:40:14 AM (13 months ago)

Author:

slebonnois

Message:

BdeBatz? : Cleans microphysics and makes few corrections for physics

File:

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

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

@@ -25 +25 @@
   !     -------
   !     Adapted from the NASA Ames code by R. Wordsworth (2009)
-  !     Clean and adaptation to Titan by J. Vatant d'Ollone (2016-17)
-  !     Clean and correction to Titan by B. de Batz de Trenquelléon (2022)
-  !     New optics added for Titan's clouds by B. de Batz de Trenquelléon (2023)
-  !     The whole routine needs to be redone...
+  !     
+  !     Modified
+  !     --------
+  !     J. Vatant d'Ollone (2016-17)
+  !         --> Clean and adaptation to Titan
+  !     B. de Batz de Trenquelléon (2022-2023)
+  !         --> Clean and correction to Titan
+  !         --> New optics added for Titan's clouds
   !     
   !==================================================================
@@ -116 +120 @@
   ! Variables for new optics
   integer iq, iw, FTYPE, CTYPE
-  real*8 m0as,m0af,m0ccn,m3as,m3af,m3ccn,m3ices,m3cld
+  real*8 m0as,m0af,m0ccn,m3as,m3af,m3ccn,m3cld
   real*8 dtauaer_s,dtauaer_f,dtau_ccn,dtau_cld
   real*8,save :: rhoaer_s(L_NSPECTV),ssa_s(L_NSPECTV),asf_s(L_NSPECTV)
@@ -202 +206 @@
                m3as = pqmo(ilay,2) / 2.0
                m3af = pqmo(ilay,4) / 2.0
-               ! Cut-off
-               IF (ilay .lt. 19) THEN
-                  m3as = pqmo(19,2) / 2.0 * (zlev(ilay+1)-zlev(ilay)) / (zlev(20)-zlev(19))
-                  m3af = pqmo(19,4) / 2.0 * (zlev(ilay+1)-zlev(ilay)) / (zlev(20)-zlev(19))
+               ! Cut-off (here for p = 2.7e3Pa / alt = 70km)
+               IF (ilay .lt. 23) THEN
+                  m3as = pqmo(23,2) / 2.0 * (zlev(ilay+1)-zlev(ilay)) / (zlev(24)-zlev(23))
+                  m3af = pqmo(23,4) / 2.0 * (zlev(ilay+1)-zlev(ilay)) / (zlev(24)-zlev(23))
                ENDIF
                
@@ -226 +230 @@
                m0as  = pqmo(ilay,1) / 2.0
                m3as  = pqmo(ilay,2) / 2.0
-               ! If not callclouds : must have a cut-off
+               ! If not callclouds : must have a cut-off (here for p = 2.7e3Pa / alt = 70km)
                IF (.NOT. callclouds) THEN
-                  IF (ilay .lt. 19) THEN
-                     m0as = pqmo(19,1) / 2.0 * (zlev(ilay+1)-zlev(ilay)) / (zlev(20)-zlev(19))
-                     m3as = pqmo(19,2) / 2.0 * (zlev(ilay+1)-zlev(ilay)) / (zlev(20)-zlev(19))
+                  IF (ilay .lt. 23) THEN
+                     m0as = pqmo(23,1) / 2.0 * (zlev(ilay+1)-zlev(ilay)) / (zlev(24)-zlev(23))
+                     m3as = pqmo(23,2) / 2.0 * (zlev(ilay+1)-zlev(ilay)) / (zlev(24)-zlev(23))
                   ENDIF
                ENDIF
@@ -240 +244 @@
                m0af  = pqmo(ilay,3) / 2.0
                m3af  = pqmo(ilay,4) / 2.0
-               ! If not callclouds : must have a cut-off
+               ! If not callclouds : must have a cut-off (here for p = 2.7e3Pa / alt = 70km)
                IF (.NOT. callclouds) THEN
-                  IF (ilay .lt. 19) THEN
-                     m0af = pqmo(19,3) / 2.0 * (zlev(ilay+1)-zlev(ilay)) / (zlev(20)-zlev(19))
-                     m3af = pqmo(19,4) / 2.0 * (zlev(ilay+1)-zlev(ilay)) / (zlev(20)-zlev(19))
+                  IF (ilay .lt. 23) THEN
+                     m0af = pqmo(23,3) / 2.0 * (zlev(ilay+1)-zlev(ilay)) / (zlev(24)-zlev(23))
+                     m3af = pqmo(23,4) / 2.0 * (zlev(ilay+1)-zlev(ilay)) / (zlev(24)-zlev(23))
                   ENDIF
                ENDIF
@@ -273 +277 @@
                m0ccn = pqmo(ilay,5) / 2.0
                m3ccn = pqmo(ilay,6) / 2.0
-               m3ices = 0.0d0
-               m3cld  = m3ccn
+               m3cld = 0.0d0
                
                ! Clear / Dark column method :
@@ -285 +288 @@
                IF (CDCOLUMN == 0) THEN
                   DO iq = 2, nice
-                     m3ices = m3ices + (pqmo(ilay,ices_indx(iq)) / 2.0)
                      m3cld  = m3cld + (pqmo(ilay,ices_indx(iq)) / 2.0)
                   ENDDO
@@ -293 +295 @@
                ELSEIF (CDCOLUMN == 1) THEN
                   DO iq = 1, nice
-                     m3ices = m3ices + (pqmo(ilay,ices_indx(iq)) / 2.0)
                      m3cld  = m3cld + (pqmo(ilay,ices_indx(iq)) / 2.0)
                   ENDDO
@@ -302 +303 @@
                   WRITE(*,*) 'WE USE DARK COLUMN ...'
                   DO iq = 1, nice
-                     m3ices = m3ices + (pqmo(ilay,ices_indx(iq)) / 2.0)
                      m3cld  = m3cld + (pqmo(ilay,ices_indx(iq)) / 2.0)
                   ENDDO
@@ -309 +309 @@
 
                ! For small dropplets, opacity of nucleus dominates
-               !ssa_cld(nw) = (ssa_ccn(nw)*m3ccn + ssa_cld(nw)*m3ices) / (m3ccn + m3ices)
                ssa_cld(nw) = (ssa_ccn(nw)*m3ccn + ssa_cld(nw)*m3cld) / (m3ccn + m3cld)
                ssa_cld(nw) = ssa_cld(nw) * Fssa