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

Timestamp:

Nov 28, 2018, 12:30:55 PM (6 years ago)

Author:

jvatant

Message:

Add a simple seasonal+latitudinal variation model (season_haze.F90) for the vertical mean profile
of haze. Based on Karkoscha,2016 . Activable through "seashaze" key (default=T).
--JVO

File:

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

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

@@ -1 +1 @@
 SUBROUTINE OPTCV(PQO,NLAY,PLEV,TMID,PMID,  &
-     DTAUV,TAUV,TAUCUMV,WBARV,COSBV,TAURAY,TAUGSURF)
+     DTAUV,TAUV,TAUCUMV,WBARV,COSBV,TAURAY,TAUGSURF,SEASHAZEFACT)
 
   use radinc_h
@@ -6 +6 @@
   use gases_h
   use comcstfi_mod, only: r
-  use callkeys_mod, only: continuum,graybody,callgasvis,callclouds,callmufi,uncoupl_optic_haze
+  use callkeys_mod, only: continuum,graybody,callgasvis,callclouds,callmufi,seashaze,uncoupl_optic_haze
   use tracer_h, only: nmicro,nice
   use MMP_OPTICS
@@ -48 +48 @@
   REAL*8, INTENT(IN)  :: TMID(L_LEVELS), PMID(L_LEVELS)
   REAL*8, INTENT(IN)  :: TAURAY(L_NSPECTV)
+  REAL*8, INTENT(IN)  :: SEASHAZEFACT(L_LEVELS)
   
   REAL*8, INTENT(OUT) :: DTAUV(L_NLAYRAD,L_NSPECTV,L_NGAUSS)
@@ -184 +185 @@
           ! Call fixed vertical haze profile of extinction - same for all columns
           call disr_haze(dz(k),plev(k),wnov(nw),dhaze_T(k,nw),SSA_T(k,nw),ASF_T(k,nw))
+          if (seashaze) dhaze_T(k,nw) = dhaze_T(k,nw)*seashazefact(k)
         ENDIF