Changeset 3685


Ignore:
Timestamp:
Mar 14, 2025, 3:23:23 PM (4 months ago)
Author:
debatzbr
Message:

callcorrk_pluto_mod: new diagnostics of optical thickness
BBT

Location:
trunk/LMDZ.PLUTO/libf/phypluto
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.PLUTO/libf/phypluto/callcorrk.F90

    r3684 r3685  
    104104      REAL,INTENT(OUT) :: OSR_nu(ngrid,L_NSPECTV)         ! Outgoing SW radiation in each band (Normalized to the band width (W/m2/cm-1).
    105105      REAL,INTENT(OUT) :: GSR_nu(ngrid,L_NSPECTV)         ! Surface SW radiation in each band (Normalized to the band width (W/m2/cm-1).
     106      REAL,INTENT(OUT) :: int_dtaui(ngrid,nlayer,L_NSPECTI) ! VI optical thickness of layers within narrowbands for diags ().
     107      REAL,INTENT(OUT) :: int_dtauv(ngrid,nlayer,L_NSPECTV) ! IR optical thickness of layers within narrowbands for diags ().
    106108      REAL,INTENT(OUT) :: tau_col(ngrid)                  ! Diagnostic from aeropacity.
    107109      REAL,INTENT(OUT) :: albedo_equivalent(ngrid)        ! Spectrally Integrated Albedo. For Diagnostic. By MT2015
    108       REAL,INTENT(OUT) :: int_dtaui(ngrid,nlayer,L_NSPECTI) ! VI optical thickness of layers within narrowbands for diags ().
    109       REAL,INTENT(OUT) :: int_dtauv(ngrid,nlayer,L_NSPECTV) ! IR optical thickness of layers within narrowbands for diags ().
    110 
    111 
    112 
    113110
    114111

  • trunk/LMDZ.PLUTO/libf/phypluto/callcorrk_pluto_mod.F90

    r3683 r3685  
    55CONTAINS
    66
    7     subroutine callcorrk_pluto(icount,ngrid,nlayer,pq,nq,qsurf,   &
    8           albedo,emis,mu0,pplev,pplay,pt,   &
    9           zzlay,zzlev,tsurf,fract,dist_star,dtau_aer,    &
    10           dtlw,dtsw,fluxsurf_lw,    &
    11           fluxsurf_sw,fluxtop_lw,fluxtop_sw,fluxtop_dn, &
    12           reffrad,tau_col,ptime,pday,firstcall,lastcall)
     7    subroutine callcorrk_pluto(icount,ngrid,nlayer,pq,nq,qsurf,&
     8          albedo,emis,mu0,pplev,pplay,pt,                      &
     9          zzlay,zzlev,tsurf,fract,dist_star,dtau_aer,          &
     10          dtlw,dtsw,fluxsurf_lw,                               &
     11          fluxsurf_sw,fluxtop_lw,fluxtop_sw,fluxtop_dn,        &
     12          reffrad,int_dtaui,int_dtauv,tau_col,                 &
     13          ptime,pday,firstcall,lastcall)
    1314
    1415      use radinc_h
    … …  
    102103      REAL fluxtop_sw(ngrid)    ! outgoing LW flux to space (W/m2)
    103104      REAL fluxtop_dn(ngrid)    ! incident top of atmosphere SW flux (W/m2)
     105      REAL int_dtaui(ngrid,nlayer,L_NSPECTI) ! VI optical thickness of layers within narrowbands for diags ().
     106      REAL int_dtauv(ngrid,nlayer,L_NSPECTV) ! IR optical thickness of layers within narrowbands for diags ().
    104107
    105108!-----------------------------------------------------------------------
    … …  
    745748         end do
    746749
     750         ! Optical thickness diagnostics
     751         ! Output exp(-tau) because gweight ponderates exp and not tau itself
     752         int_dtauv(ig,:,:) = 0.0d0
     753         int_dtaui(ig,:,:) = 0.0d0
     754         do l=1,L_NLAYRAD
     755            do nw=1,L_NSPECTV
     756               do k=1,L_NGAUSS
     757                  int_dtauv(ig,l,nw)= int_dtauv(ig,l,nw) + exp(-dtauv(l,nw,k))*gweight(k)
     758               enddo
     759            enddo
     760            do nw=1,L_NSPECTI
     761               do k=1,L_NGAUSS
     762                  int_dtaui(ig,l,nw)= int_dtaui(ig,l,nw) + exp(-dtaui(l,nw,k))*gweight(k)
     763               enddo
     764            enddo
     765         enddo
     766
    747767! **********************************************************
    748768!     NON NLTE correction in Pluto atmosphere

  • trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90

    r3684 r3685  
    10221022                               zzlay,zzlev,tsurf,fract,dist_star,dtau_aer,       &
    10231023                               zdtlw,zdtsw,fluxsurf_lw,fluxsurf_sw,fluxtop_lw,   &
    1024                                fluxtop_sw,fluxtop_dn,reffrad,tau_col,ptime,pday, &
    1025                                firstcall,lastcall)
     1024                               fluxtop_sw,fluxtop_dn,reffrad,                    &
     1025                               int_dtaui,int_dtauv,tau_col,                      &
     1026                               ptime,pday,firstcall,lastcall)
    10261027                  albedo_equivalent(1:ngrid)=albedo(1:ngrid,1)
    10271028                  fluxrad_sky(1:ngrid)=emis(1:ngrid)*fluxsurf_lw(1:ngrid)+       &
Note: See TracChangeset for help on using the changeset viewer.