Changeset 3695 for trunk


Ignore:
Timestamp:
Mar 19, 2025, 2:27:25 PM (4 months ago)
Author:
debatzbr
Message:

Calculate the visible opacity everywhere
BBT

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

Legend:

Unmodified
Added
Removed
  • TabularUnified trunk/LMDZ.PLUTO/libf/phypluto/callcorrk.F90

    r3687 r3695  
    927927!-----------------------------------------------------------------------
    928928
     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
    929934         if(fract(ig) .ge. 1.0e-4) then ! Only during daylight.
    930935            if((ngrid.eq.1).and.(global1d))then
     
    938943            endif
    939944
    940             call optcv(dtauv,tauv,taucumv,plevrad,                 &
    941                  qxvaer,qsvaer,gvaer,wbarv,cosbv,tauray,tauaero,   &
    942                  tmid,pmid,taugsurf,qvar,muvarrad)
    943 
    944945            call sfluxv(dtauv,tauv,taucumv,albv,dwnv,wbarv,cosbv,   &
    945946                 acosz,stel_fract,nfluxtopv,fluxtopvdn,nfluxoutv_nu,&
     
    953954            nfluxoutv_nu(:) = 0.0d0
    954955            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
    961960         end if
    962961
     
    10681067            do nw=1,L_NSPECTV
    10691068               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)
    10711070               enddo
    10721071            enddo
    10731072            do nw=1,L_NSPECTI
    10741073               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)
    10761075               enddo
    10771076            enddo
  • TabularUnified trunk/LMDZ.PLUTO/libf/phypluto/optci.F90

    r3184 r3695  
    133133  !     spectral interval, NW, and each Gauss point, NG.
    134134
     135  dtaui(:,:,:) = 0.0
     136  taucumi(:,:,:) = 0.0
     137
    135138  taugsurf(:,:) = 0.0
    136139  dpr(:)        = 0.0
  • TabularUnified trunk/LMDZ.PLUTO/libf/phypluto/optcv.F90

    r3691 r3695  
    128128  !     Calculate the continuum opacities, i.e., those that do not depend on
    129129  !     NG, the Gauss index.
     130
     131  dtauv(:,:,:)   = 0.0
     132  tauv(:,:,:)    = 0.0
     133  taucumv(:,:,:) = 0.0
    130134
    131135  taugsurf(:,:) = 0.0
  • TabularUnified trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90

    r3694 r3695  
    22012201      !   call write_output('dtaui'//str2,'Layer optical thickness attenuation in IR band '//str2,'',int_dtaui(:,nlayer:1:-1,nw))
    22022202      !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     
    22042208      ! Temporary inclusions for heating diagnostics.
    22052209      if (.not.fast) then
Note: See TracChangeset for help on using the changeset viewer.