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

Timestamp:

Jan 11, 2017, 3:33:51 PM (8 years ago)

Author:

jvatant

Message:

+ Major clean of the new LMDZ.TITAN from too-generic options and routines (water, co2, ocean, surface type ...)
+ From this revision LMDZ.TITAN begins to be really separated from LMDZ.GENERIC
+ Partial desactivation of aerosols, only the dummy case is still enabled to keep the code running ( new aerosol routines to come in followings commits )

JVO

File:

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

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

@@ -7 +7 @@
   use gases_h
   use comcstfi_mod, only: g, r, mugaz
-  use callkeys_mod, only: kastprof,continuum,graybody,H2Ocont_simple,callgasvis
+  use callkeys_mod, only: continuum,graybody,callgasvis
 
   implicit none
@@ -106 +106 @@
 
      ! if we have continuum opacities, we need dz
-     if(kastprof)then
-        dz(k) = dpr(k)*(1000.0d0*8.3145d0/muvar(k))*TMID(K)/(g*PMID(K))
-        U(k)  = Cmk*DPR(k)*mugaz/muvar(k) 
-     else
-        dz(k) = dpr(k)*R*TMID(K)/(glat_ig*PMID(K))*mugaz/muvar(k)
-        U(k)  = Cmk*DPR(k)*mugaz/muvar(k)     ! only Cmk line in optci.F  
-            !JL13 the mugaz/muvar factor takes into account water meanmolecular weight if water is present
-     endif
+
+      dz(k) = dpr(k)*R*TMID(K)/(glat_ig*PMID(K))*mugaz/muvar(k)
+      U(k)  = Cmk*DPR(k)*mugaz/muvar(k)     ! only Cmk line in optci.F  
+          !JL13 the mugaz/muvar factor takes into account water meanmolecular weight if water is present
+     
 
      call tpindex(PMID(K),TMID(K),QVAR(K),pfgasref,tgasref,WREFVAR, &
@@ -186 +183 @@
                        indv(nw,igas,jgas) = interm
                        ! should be irrelevant in the visible
-                    elseif(jgas.eq.igas_He)then 
-                       interm = indv(nw,igas,jgas)
-                       call interpolateH2He(wn_cont,T_cont,p_cross,p_cont,dtempc,.false.,interm)
-                       indv(nw,igas,jgas) = interm
                     endif
                     dtemp = dtemp + dtempc
                  enddo
-
-              elseif(igas.eq.igas_H2O.and.T_cont.gt.200.0)then
-
-                 p_air = dble(PMID(k)*scalep) - p_cont ! note assumes background is air!
-                 if(H2Ocont_simple)then
-                    call interpolateH2Ocont_PPC(wn_cont,T_cont,p_cont,p_air,dtemp,.false.)
-                 else
-                    interm = indv(nw,igas,igas)
-                    call interpolateH2Ocont_CKD(wn_cont,T_cont,p_cont,p_air,dtemp,.false.,interm)
-                    indv(nw,igas,igas) = interm
-                 endif
 
               endif