- Timestamp:
- Mar 19, 2025, 2:27:25 PM (4 months ago)
- Location:
- trunk/LMDZ.PLUTO/libf/phypluto
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified trunk/LMDZ.PLUTO/libf/phypluto/callcorrk.F90 ¶
r3687 r3695 927 927 !----------------------------------------------------------------------- 928 928 929 ! Call everywhere for diagnostics. 930 call optcv(dtauv,tauv,taucumv,plevrad, & 931 qxvaer,qsvaer,gvaer,wbarv,cosbv,tauray,tauaero, & 932 tmid,pmid,taugsurf,qvar,muvarrad) 933 929 934 if(fract(ig) .ge. 1.0e-4) then ! Only during daylight. 930 935 if((ngrid.eq.1).and.(global1d))then … … 938 943 endif 939 944 940 call optcv(dtauv,tauv,taucumv,plevrad, &941 qxvaer,qsvaer,gvaer,wbarv,cosbv,tauray,tauaero, &942 tmid,pmid,taugsurf,qvar,muvarrad)943 944 945 call sfluxv(dtauv,tauv,taucumv,albv,dwnv,wbarv,cosbv, & 945 946 acosz,stel_fract,nfluxtopv,fluxtopvdn,nfluxoutv_nu,& … … 953 954 nfluxoutv_nu(:) = 0.0d0 954 955 nfluxgndv_nu(:) = 0.0d0 955 fmnetv_nu(:,:)=0.0 956 do l=1,L_NLAYRAD 957 fmnetv(l)=0.0d0 958 fluxupv(l)=0.0d0 959 fluxdnv(l)=0.0d0 960 end do 956 fmnetv_nu(:,:) = 0.0d0 957 fmnetv(:) = 0.0d0 958 fluxupv(:) = 0.0d0 959 fluxdnv(:) = 0.0d0 961 960 end if 962 961 … … 1068 1067 do nw=1,L_NSPECTV 1069 1068 do k=1,L_NGAUSS 1070 int_dtauv(ig,l,nw) = int_dtauv(ig,l,nw) + exp(-dtauv(l,nw,k))*gweight(k)1069 int_dtauv(ig,l,nw) = int_dtauv(ig,l,nw) + exp(-dtauv(l,nw,k))*gweight(k) 1071 1070 enddo 1072 1071 enddo 1073 1072 do nw=1,L_NSPECTI 1074 1073 do k=1,L_NGAUSS 1075 int_dtaui(ig,l,nw) = int_dtaui(ig,l,nw) + exp(-dtaui(l,nw,k))*gweight(k)1074 int_dtaui(ig,l,nw) = int_dtaui(ig,l,nw) + exp(-dtaui(l,nw,k))*gweight(k) 1076 1075 enddo 1077 1076 enddo -
TabularUnified trunk/LMDZ.PLUTO/libf/phypluto/optci.F90 ¶
r3184 r3695 133 133 ! spectral interval, NW, and each Gauss point, NG. 134 134 135 dtaui(:,:,:) = 0.0 136 taucumi(:,:,:) = 0.0 137 135 138 taugsurf(:,:) = 0.0 136 139 dpr(:) = 0.0 -
TabularUnified trunk/LMDZ.PLUTO/libf/phypluto/optcv.F90 ¶
r3691 r3695 128 128 ! Calculate the continuum opacities, i.e., those that do not depend on 129 129 ! NG, the Gauss index. 130 131 dtauv(:,:,:) = 0.0 132 tauv(:,:,:) = 0.0 133 taucumv(:,:,:) = 0.0 130 134 131 135 taugsurf(:,:) = 0.0 -
TabularUnified trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90 ¶
r3694 r3695 2201 2201 ! call write_output('dtaui'//str2,'Layer optical thickness attenuation in IR band '//str2,'',int_dtaui(:,nlayer:1:-1,nw)) 2202 2202 !enddo 2203 2203 call write_output('dtauv_01','Layer optical thickness attenuation in VI band','',int_dtauv(:,nlayer:1:-1,1)) 2204 call write_output('dtauv_23','Layer optical thickness attenuation in VI band','',int_dtauv(:,nlayer:1:-1,23)) 2205 call write_output('dtaui_01','Layer optical thickness attenuation in IR band','',int_dtaui(:,nlayer:1:-1,1)) 2206 call write_output('dtaui_17','Layer optical thickness attenuation in IR band','',int_dtaui(:,nlayer:1:-1,17)) 2207 2204 2208 ! Temporary inclusions for heating diagnostics. 2205 2209 if (.not.fast) then
Note: See TracChangeset
for help on using the changeset viewer.