Changeset 3685
- Timestamp:
- Mar 14, 2025, 3:23:23 PM (4 months ago)
- Location:
- trunk/LMDZ.PLUTO/libf/phypluto
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.PLUTO/libf/phypluto/callcorrk.F90
r3684 r3685 104 104 REAL,INTENT(OUT) :: OSR_nu(ngrid,L_NSPECTV) ! Outgoing SW radiation in each band (Normalized to the band width (W/m2/cm-1). 105 105 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 (). 106 108 REAL,INTENT(OUT) :: tau_col(ngrid) ! Diagnostic from aeropacity. 107 109 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 113 110 114 111 -
trunk/LMDZ.PLUTO/libf/phypluto/callcorrk_pluto_mod.F90
r3683 r3685 5 5 CONTAINS 6 6 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) 13 14 14 15 use radinc_h … … 102 103 REAL fluxtop_sw(ngrid) ! outgoing LW flux to space (W/m2) 103 104 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 (). 104 107 105 108 !----------------------------------------------------------------------- … … 745 748 end do 746 749 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 747 767 ! ********************************************************** 748 768 ! NON NLTE correction in Pluto atmosphere -
trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90
r3684 r3685 1022 1022 zzlay,zzlev,tsurf,fract,dist_star,dtau_aer, & 1023 1023 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) 1026 1027 albedo_equivalent(1:ngrid)=albedo(1:ngrid,1) 1027 1028 fluxrad_sky(1:ngrid)=emis(1:ngrid)*fluxsurf_lw(1:ngrid)+ &
Note: See TracChangeset
for help on using the changeset viewer.