Changeset 2138 for trunk/LMDZ.GENERIC/libf

Timestamp:

May 22, 2019, 4:15:51 PM (6 years ago)

Author:

jvatant

Message:

Correct outptuts of dtaui/v with correct ponderations
of weights, but now it is attenuation exp(-tau)
--JVO

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• callcorrk.F90 (modified) (2 diffs)
• inifis_mod.F90 (modified) (1 diff)
• physiq_mod.F90 (modified) (1 diff)

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

@@ -934 +934 @@
             int_dtauv(ig,l,nw) = 0.0d0
              DO k=1,L_NGAUSS
-              int_dtauv(ig,l,nw)= int_dtauv(ig,l,nw) + dtauv(l,nw,k)*gweight(k)
+              ! Output exp(-tau) because gweight ponderates exp and not tau itself
+              int_dtauv(ig,l,nw)= int_dtauv(ig,l,nw) + exp(-dtauv(l,nw,k))*gweight(k)
              ENDDO
           enddo
@@ -940 +941 @@
            int_dtaui(ig,l,nw) = 0.0d0
              DO k=1,L_NGAUSS
-              int_dtaui(ig,l,nw)= int_dtaui(ig,l,nw) + dtaui(l,nw,k)*gweight(k)
+              ! Output exp(-tau) because gweight ponderates exp and not tau itself
+              int_dtaui(ig,l,nw)= int_dtaui(ig,l,nw) + exp(-dtaui(l,nw,k))*gweight(k)
              ENDDO
           enddo

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

@@ -327 +327 @@
      write(*,*)" specOLR = ",specOLR
 
-     write(*,*)"Output diagnostic optical thickness?"
+     write(*,*)"Output diagnostic optical thickness attenuation exp(-tau)?"
      diagdtau=.false.
      call getin_p("diagdtau",diagdtau)

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

@@ -2194 +2194 @@
 
         ! Diagnostics of optical thickness
+        ! Warning this is exp(-tau), I let you postproc with -log to have tau itself - JVO 19
         if (diagdtau) then                
           do nw=1,L_NSPECTV
             write(str2,'(i2.2)') nw
-            call writediagfi(ngrid,'dtauv'//str2,'Layer optical thickness in VI band '//str2,'',1,int_dtauv(:,nlayer:1:-1,nw))
+            call writediagfi(ngrid,'dtauv'//str2,'Layer optical thickness attenuation in VI band '//str2,'',1,int_dtauv(:,nlayer:1:-1,nw))
           enddo
           do nw=1,L_NSPECTI
             write(str2,'(i2.2)') nw
-            call writediagfi(ngrid,'dtaui'//str2,'Layer optical thickness in IR band '//str2,'',1,int_dtaui(:,nlayer:1:-1,nw))
+            call writediagfi(ngrid,'dtaui'//str2,'Layer optical thickness attenuation in IR band '//str2,'',1,int_dtaui(:,nlayer:1:-1,nw))
           enddo
         endif