Changeset 2543 for trunk/LMDZ.GENERIC/libf/phystd/optci.F90

Timestamp:

Jul 5, 2021, 4:04:28 PM (3 years ago)

Author:

aslmd

Message:

Generic GCM:

Adding k-coefficients mixing on the fly
Working with MordernTrac?

JVO + YJ

File:

• trunk/LMDZ.GENERIC/libf/phystd/optci.F90 (modified) (3 diffs)

trunk/LMDZ.GENERIC/libf/phystd/optci.F90

@@ -15 +15 @@
   use comcstfi_mod, only: g, r, mugaz
   use callkeys_mod, only: kastprof,continuum,graybody
+  use recombin_corrk_mod, only: corrk_recombin, gasi_recomb
   implicit none
 
@@ -247 +248 @@
               ! transfer on the tested simulations !
 
-              tmpk = GASI(MT(K):MT(K)+1,MP(K):MP(K)+1,1,NW,NG)
+              IF (corrk_recombin) THEN ! added by JVO
+                tmpk = GASI_RECOMB(MT(K):MT(K)+1,MP(K):MP(K)+1,1,NW,NG) ! contains the mix of recombined species
+              ELSE
+                tmpk = GASI(MT(K):MT(K)+1,MP(K):MP(K)+1,1,NW,NG)
+              ENDIF
 
               KCOEF(1) = tmpk(1,1) ! KCOEF(1) = GASI(MT(K),MP(K),1,NW,NG)
@@ -256 +261 @@
            else
 
-              tmpkvar = GASI(MT(K):MT(K)+1,MP(K):MP(K)+1,NVAR(K):NVAR(K)+1,NW,NG)
+              IF (corrk_recombin) THEN ! added by JVO
+                tmpkvar = GASI_RECOMB(MT(K):MT(K)+1,MP(K):MP(K)+1,NVAR(K):NVAR(K)+1,NW,NG)
+              ELSE
+                tmpkvar = GASI(MT(K):MT(K)+1,MP(K):MP(K)+1,NVAR(K):NVAR(K)+1,NW,NG)
+              ENDIF
 
               KCOEF(1) = tmpkvar(1,1,1) + WRATIO(K) *  &