Changeset 3695 for trunk/LMDZ.PLUTO/libf/phypluto/callcorrk.F90

Timestamp:

Mar 19, 2025, 2:27:25 PM (4 months ago)

Author:

debatzbr

Message:

Calculate the visible opacity everywhere
BBT

File:

• trunk/LMDZ.PLUTO/libf/phypluto/callcorrk.F90 (modified) (4 diffs)

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

@@ -927 +927 @@
 !-----------------------------------------------------------------------
 
+         ! Call everywhere for diagnostics.
+         call optcv(dtauv,tauv,taucumv,plevrad,                 &
+                 qxvaer,qsvaer,gvaer,wbarv,cosbv,tauray,tauaero,   &
+                 tmid,pmid,taugsurf,qvar,muvarrad)
+
          if(fract(ig) .ge. 1.0e-4) then ! Only during daylight.
             if((ngrid.eq.1).and.(global1d))then
@@ -938 +943 @@
             endif
 
-            call optcv(dtauv,tauv,taucumv,plevrad,                 &
-                 qxvaer,qsvaer,gvaer,wbarv,cosbv,tauray,tauaero,   &
-                 tmid,pmid,taugsurf,qvar,muvarrad)
-
             call sfluxv(dtauv,tauv,taucumv,albv,dwnv,wbarv,cosbv,   &
                  acosz,stel_fract,nfluxtopv,fluxtopvdn,nfluxoutv_nu,&
@@ -953 +954 @@
             nfluxoutv_nu(:) = 0.0d0
             nfluxgndv_nu(:) = 0.0d0
-            fmnetv_nu(:,:)=0.0
-            do l=1,L_NLAYRAD
-               fmnetv(l)=0.0d0
-               fluxupv(l)=0.0d0
-               fluxdnv(l)=0.0d0
-            end do
+            fmnetv_nu(:,:)  = 0.0d0
+            fmnetv(:)       = 0.0d0
+            fluxupv(:)      = 0.0d0
+            fluxdnv(:)      = 0.0d0
          end if
 
@@ -1068 +1067 @@
             do nw=1,L_NSPECTV
                do k=1,L_NGAUSS
-                  int_dtauv(ig,l,nw)= int_dtauv(ig,l,nw) + exp(-dtauv(l,nw,k))*gweight(k)
+                  int_dtauv(ig,l,nw) = int_dtauv(ig,l,nw) + exp(-dtauv(l,nw,k))*gweight(k)
                enddo
             enddo
             do nw=1,L_NSPECTI
                do k=1,L_NGAUSS
-                  int_dtaui(ig,l,nw)= int_dtaui(ig,l,nw) + exp(-dtaui(l,nw,k))*gweight(k)
+                  int_dtaui(ig,l,nw) = int_dtaui(ig,l,nw) + exp(-dtaui(l,nw,k))*gweight(k)
                enddo
             enddo