Changeset 2131 for trunk/LMDZ.GENERIC/libf

Timestamp:

Apr 29, 2019, 2:52:42 PM (6 years ago)

Author:

jvatant

Message:

+ Add diagnostics of optical thickness, in 1D, if 'diagdtau' key is activated, it

outputs dtaui/v(altitude) in diagfi.nc for every narrowband (could be done with one var
but would require to be able to have writediag in 5D

--JVO

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• callkeys_mod.F90 (modified) (1 diff)
• inifis_mod.F90 (modified) (1 diff)
• optci.F90 (modified) (3 diffs)
• optcv.F90 (modified) (3 diffs)

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

@@ -20 +20 @@
       logical,save :: specOLR
       logical,save :: kastprof
-!$OMP THREADPRIVATE(enertest,nonideal,meanOLR,kastprof)
+      logical,save :: diagdtau
+!$OMP THREADPRIVATE(enertest,nonideal,meanOLR,kastprof,diagdtau)
       logical,save :: newtonian
       logical,save :: check_cpp_match

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

@@ -327 +327 @@
      write(*,*)" specOLR = ",specOLR
 
+     write(*,*)"Output diagnostic optical thickness?"
+     diagdtau=.false.
+     call getin_p("diagdtau",diagdtau)
+     write(*,*)" diagdtau = ",diagdtau
+     ! sanity check
+     if (diagdtau.and.(ngrid.ne.1)) then
+       print*,"Diagnostic optical thickness can be output in 1D only !"
+       print*,"Start again with diagdtau=.false."
+       stop
+     endif
+
      write(*,*)"Operate in kastprof mode?"
      kastprof=.false.

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

@@ -11 +11 @@
   use radinc_h, only: L_LEVELS, L_NLAYRAD, L_NSPECTI, L_NGAUSS, &
                       L_NLEVRAD, L_REFVAR, naerkind
-  use radcommon_h, only: gasi,tlimit,wrefVAR,Cmk,tgasref,pfgasref,wnoi,scalep,indi,glat_ig
+  use radcommon_h, only: gasi,tlimit,wrefVAR,Cmk,tgasref,pfgasref,wnoi,scalep,indi,glat_ig,gweight
   use gases_h, only: gfrac, ngasmx, igas_N2, igas_He, igas_H2O, igas_H2
   use comcstfi_mod, only: g, r, mugaz
-  use callkeys_mod, only: kastprof,continuum,graybody,H2Ocont_simple
+  use callkeys_mod, only: kastprof,continuum,graybody,H2Ocont_simple,diagdtau
   implicit none
 
@@ -84 +84 @@
   !real*8 rho !! see test below
 
+  ! variables for diags of optical depth
+  real *8 INT_DTAU(L_NLAYRAD,L_NSPECTI)
+  character*2 str2
+
   integer igas, jgas
 
@@ -449 +453 @@
 !  print*,DTAUI
 !  call abort
+
+
+  !  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_NSPECTI
+        INT_DTAU(L,NW) = 0.0d+0
+          DO NG=1,L_NGAUSS
+             INT_DTAU(L,NW)= INT_DTAU(L,NW) + dtaui(L,nw,ng)*gweight(NG)
+          enddo
+       enddo
+    enddo
+    do nw=1,L_NSPECTI
+      write(str2,'(i2.2)') nw
+      call writediagfi(1,'dtaui'//str2,'Layer optical thickness in IR band '//str2,'',1,int_dtau(L_NLAYRAD:1:-1,nw))
+    enddo
+  endif
+  
   
 

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
+  use radcommon_h, only: gasv, tlimit, wrefVAR, Cmk, tgasref, pfgasref,wnov,scalep,indv,glat_ig,gweight
   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
+  use callkeys_mod, only: kastprof,continuum,graybody,H2Ocont_simple,callgasvis,diagdtau
 
   implicit none
@@ -91 +91 @@
   ! variables for k in units m^-1
   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
@@ -375 +379 @@
 
 
+  !  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
+
+
 end subroutine optcv