Changeset 2133 for trunk/LMDZ.GENERIC/libf/phystd/optcv.F90

Timestamp:

Apr 30, 2019, 10:34:11 AM (6 years ago)

Author:

jvatant

Message:

Fic a bug in r2131 the writediagfi cannot be called
in RT routines because of the iradia it goes crazy with outputs freq ...
--JVO

File:

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

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

@@ -10 +10 @@
 
   use radinc_h, only: L_NLAYRAD, L_NLEVRAD, L_LEVELS, L_NSPECTV, L_NGAUSS, L_REFVAR, NAERKIND
-  use radcommon_h, only: gasv, tlimit, wrefVAR, Cmk, tgasref, pfgasref,wnov,scalep,indv,glat_ig,gweight
+  use radcommon_h, only: gasv, tlimit, wrefVAR, Cmk, tgasref, pfgasref,wnov,scalep,indv,glat_ig
   use gases_h, only: gfrac, ngasmx, igas_H2, igas_H2O, igas_He, igas_N2
   use comcstfi_mod, only: g, r, mugaz
-  use callkeys_mod, only: kastprof,continuum,graybody,H2Ocont_simple,callgasvis,diagdtau
+  use callkeys_mod, only: kastprof,continuum,graybody,H2Ocont_simple,callgasvis
 
   implicit none
@@ -92 +92 @@
   real*8 dz(L_LEVELS)
 
-  ! variables for diags of optical depth
-  real*8 INT_DTAU(L_NLAYRAD,L_NSPECTV)
-  character*2 str2
 
   integer igas, jgas
@@ -379 +376 @@
 
 
-  !  Optical thickness for 1D diagnostics (added by JVO)
-  if (diagdtau) then ! diagtau can be true only if 1D
-    do l=1,L_NLAYRAD
-       do nw=1,L_NSPECTV
-        INT_DTAU(L,NW) = 0.0d+0
-          DO NG=1,L_NGAUSS
-             INT_DTAU(L,NW)= INT_DTAU(L,NW) + dtauv(L,nw,ng)*gweight(NG)
-          enddo
-       enddo
-    enddo
-    do nw=1,L_NSPECTV
-      write(str2,'(i2.2)') nw
-      call writediagfi(1,'dtauv'//str2,'Layer optical thickness in VI band '//str2,'',1,int_dtau(L_NLAYRAD:1:-1,nw))
-    enddo
-  endif