Changeset 2138


Ignore:
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
Files:
8 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.GENERIC/README

    r2133 r2138  
    14821482 outputs dtaui/v(altitude) in diagfi.nc for every narrowband (could be done with one var
    14831483 but would require to be able to have writediag in 5D)
     1484 EDIT (22/05/19) : To have correct calculations, output is now exp(-tau), you need to postproc it to have tau.

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

    r2134 r2138  
    934934            int_dtauv(ig,l,nw) = 0.0d0
    935935             DO k=1,L_NGAUSS
    936               int_dtauv(ig,l,nw)= int_dtauv(ig,l,nw) + dtauv(l,nw,k)*gweight(k)
     936              ! Output exp(-tau) because gweight ponderates exp and not tau itself
     937              int_dtauv(ig,l,nw)= int_dtauv(ig,l,nw) + exp(-dtauv(l,nw,k))*gweight(k)
    937938             ENDDO
    938939          enddo
    … …  
    940941           int_dtaui(ig,l,nw) = 0.0d0
    941942             DO k=1,L_NGAUSS
    942               int_dtaui(ig,l,nw)= int_dtaui(ig,l,nw) + dtaui(l,nw,k)*gweight(k)
     943              ! Output exp(-tau) because gweight ponderates exp and not tau itself
     944              int_dtaui(ig,l,nw)= int_dtaui(ig,l,nw) + exp(-dtaui(l,nw,k))*gweight(k)
    943945             ENDDO
    944946          enddo

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

    r2133 r2138  
    327327     write(*,*)" specOLR = ",specOLR
    328328
    329      write(*,*)"Output diagnostic optical thickness?"
     329     write(*,*)"Output diagnostic optical thickness attenuation exp(-tau)?"
    330330     diagdtau=.false.
    331331     call getin_p("diagdtau",diagdtau)

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

    r2133 r2138  
    21942194
    21952195        ! Diagnostics of optical thickness
     2196        ! Warning this is exp(-tau), I let you postproc with -log to have tau itself - JVO 19
    21962197        if (diagdtau) then               
    21972198          do nw=1,L_NSPECTV
    21982199            write(str2,'(i2.2)') nw
    2199             call writediagfi(ngrid,'dtauv'//str2,'Layer optical thickness in VI band '//str2,'',1,int_dtauv(:,nlayer:1:-1,nw))
     2200            call writediagfi(ngrid,'dtauv'//str2,'Layer optical thickness attenuation in VI band '//str2,'',1,int_dtauv(:,nlayer:1:-1,nw))
    22002201          enddo
    22012202          do nw=1,L_NSPECTI
    22022203            write(str2,'(i2.2)') nw
    2203             call writediagfi(ngrid,'dtaui'//str2,'Layer optical thickness in IR band '//str2,'',1,int_dtaui(:,nlayer:1:-1,nw))
     2204            call writediagfi(ngrid,'dtaui'//str2,'Layer optical thickness attenuation in IR band '//str2,'',1,int_dtaui(:,nlayer:1:-1,nw))
    22042205          enddo
    22052206        endif

  • trunk/LMDZ.TITAN/README

    r2133 r2138  
    15621562 outputs dtaui/v(altitude) in diagfi.nc for every narrowband (could be done with one var
    15631563 but would require to be able to have writediag in 5D)
     1564 EDIT (22/05/19) : To have correct calculations, output is now exp(-tau), you need to postproc it to have tau.

  • trunk/LMDZ.TITAN/libf/phytitan/callcorrk.F90

    r2133 r2138  
    543543            int_dtauv(ig,l,nw) = 0.0d0
    544544             DO k=1,L_NGAUSS
    545               int_dtauv(ig,l,nw)= int_dtauv(ig,l,nw) + dtauv(l,nw,k)*gweight(k)
     545              ! Output exp(-tau) because gweight ponderates exp and not tau itself
     546              int_dtauv(ig,l,nw)= int_dtauv(ig,l,nw) + exp(-dtauv(l,nw,k))*gweight(k)
    546547             ENDDO
    547548          enddo
    … …  
    549550           int_dtaui(ig,l,nw) = 0.0d0
    550551             DO k=1,L_NGAUSS
    551               int_dtaui(ig,l,nw)= int_dtaui(ig,l,nw) + dtaui(l,nw,k)*gweight(k)
     552              ! Output exp(-tau) because gweight ponderates exp and not tau itself
     553              int_dtaui(ig,l,nw)= int_dtaui(ig,l,nw) + exp(-dtaui(l,nw,k))*gweight(k)
    552554             ENDDO
    553555          enddo

  • trunk/LMDZ.TITAN/libf/phytitan/inifis_mod.F90

    r2133 r2138  
    337337     write(*,*)" specOLR = ",specOLR
    338338
    339      write(*,*)"Output diagnostic optical thickness?"
     339     write(*,*)"Output diagnostic optical thickness attenuation exp(-tau)?"
    340340     diagdtau=.false.
    341341     call getin_p("diagdtau",diagdtau)

  • trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90

    r2133 r2138  
    15191519
    15201520        ! Diagnostics of optical thickness
     1521        ! Warning this is exp(-tau), I let you postproc with -log to have tau itself - JVO 19
    15211522        if (diagdtau) then               
    15221523          do nw=1,L_NSPECTV
    15231524            write(str2,'(i2.2)') nw
    1524             call writediagfi(ngrid,'dtauv'//str2,'Layer optical thickness in VI band '//str2,'',1,int_dtauv(:,nlayer:1:-1,nw))
     1525            call writediagfi(ngrid,'dtauv'//str2,'Layer optical thickness attenuation in VI band '//str2,'',1,int_dtauv(:,nlayer:1:-1,nw))
    15251526          enddo
    15261527          do nw=1,L_NSPECTI
    15271528            write(str2,'(i2.2)') nw
    1528             call writediagfi(ngrid,'dtaui'//str2,'Layer optical thickness in IR band '//str2,'',1,int_dtaui(:,nlayer:1:-1,nw))
     1529            call writediagfi(ngrid,'dtaui'//str2,'Layer optical thickness attenuation in IR band '//str2,'',1,int_dtaui(:,nlayer:1:-1,nw))
    15291530          enddo
    15301531        endif
Note: See TracChangeset for help on using the changeset viewer.