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

Timestamp:

Nov 30, 2018, 12:55:49 AM (6 years ago)

Author:

jvatant

Message:

Major radiative transfer contribution : Add the 'corrk_recombin' option that allows to use
correlated-k for single species instead of pre-mix and enables more flexiblity for variable species.
-> Algorithm inspired from Lacis and Oinas 1991 and Amundsen et al 2016

+ Added 'recombin_corrk.F90' - Important improvements in 'sugas_corrk.90' and 'callcork.F90'

• Must have the desired variable species as tracers -> TBD : Enable to force composition even if no tracers

• To have decent CPU time recombining is not done on all gridpoints and wavelenghts but we calculate a gasi/v_recomb variable on the reference corrk-k T,P grid (only for T,P values who match the atmospheric conditions ) which is then processed as a standard pre-mix in optci/v routines, but updated every time tracers on the ref P grid have varied > 1%.

READ CAREFULY :

• In case of 'corrk_recombin', the variable L_NREFVAR doesn't have the same meaning as before and doesn't stand for the different mixing ratios but the different species.
• Input corr-k should be found in corrkdir within 'corrk_gcm_IR/VI_XXX.dat' and can contain a 'fixed' specie ( compulsory if you include self-broadening ) that MUST have been created with correct mixing ratios, or a variable specie for which mixing ratio MUST have been set to 1 ( no self-broadening then, assume it's a trace sepecie ) -> You can't neither have CIA of variable species included upstream in the corr-k

File:

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

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

@@ -1 @@
-SUBROUTINE OPTCV(PQO,NLAY,PLEV,TMID,PMID,  &
+SUBROUTINE OPTCV(PQMO,NLAY,PLEV,TMID,PMID,  &
      DTAUV,TAUV,TAUCUMV,WBARV,COSBV,TAURAY,TAUGSURF,SEASHAZEFACT)
 
   use radinc_h
-  use radcommon_h, only: gasv,tlimit,Cmk,gzlat_ig,tgasref,pfgasref,wnov,scalep,indv,gweight
+  use radcommon_h, only: gasv,gasv_recomb,tlimit,Cmk,gzlat_ig, &
+                         tgasref,pfgasref,wnov,scalep,indv,gweight
   use gases_h
   use comcstfi_mod, only: r
-  use callkeys_mod, only: continuum,graybody,callgasvis,callclouds,callmufi,seashaze,uncoupl_optic_haze
+  use callkeys_mod, only: continuum,graybody,callgasvis,corrk_recombin,  &
+                          callclouds,callmufi,seashaze,uncoupl_optic_haze
   use tracer_h, only: nmicro,nice
   use MMP_OPTICS
@@ -43 +45 @@
   ! Input/Output
   !==========================================================
-  REAL*8, INTENT(IN)  :: PQO(nlay,nmicro) ! Tracers (X/m2).
-  INTEGER, INTENT(IN) :: NLAY             ! Number of pressure layers (for pqo)
+  REAL*8, INTENT(IN)  :: PQMO(nlay,nmicro)  ! Tracers for microphysics optics (X/m2).
+  INTEGER, INTENT(IN) :: NLAY               ! Number of pressure layers (for pqmo)
   REAL*8, INTENT(IN)  :: PLEV(L_LEVELS)
   REAL*8, INTENT(IN)  :: TMID(L_LEVELS), PMID(L_LEVELS)
@@ -167 +169 @@
         ! as long as the microphysics only isn't fully debugged -- JVO 01/18
         IF (callmufi .AND. (.NOT. uncoupl_optic_haze)) THEN
-          m0as = pqo(ilay,1)
-          m3as = pqo(ilay,2)
-          m0af = pqo(ilay,3)
-          m3af = pqo(ilay,4)
+          m0as = pqmo(ilay,1)
+          m3as = pqmo(ilay,2)
+          m0af = pqmo(ilay,3)
+          m3af = pqmo(ilay,4)
 
           IF (.NOT.mmp_sph_optics_vis(m0as,m3as,nw,ext_s,sca_s,ssa_s,asf_s)) &
@@ -267 +269 @@
            ! transfer on the tested simulations !
 
-           tmpk = GASV(MT(K):MT(K)+1,MP(K):MP(K)+1,1,NW,NG)
+           if (corrk_recombin) then 
+             tmpk = GASV_RECOMB(MT(K):MT(K)+1,MP(K):MP(K)+1,NW,NG)
+           else
+             tmpk = GASV(MT(K):MT(K)+1,MP(K):MP(K)+1,1,NW,NG)
+           endif
               
            KCOEF(1) = tmpk(1,1) ! KCOEF(1) = GASV(MT(K),MP(K),1,NW,NG)
@@ -302 +308 @@
   !     Now the full treatment for the layers, where besides the opacity
   !     we need to calculate the scattering albedo and asymmetry factors
+  ! ======================================================================
 
   ! Haze scattering