Changeset 3686 for trunk/LMDZ.PLUTO

Timestamp:

Mar 14, 2025, 5:59:49 PM (5 months ago)

Author:

debatzbr

Message:

Minor corrections for the opacity calculation
BBT

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• optcv.F90 (modified) (3 diffs)
• optcv_pluto_mod.F90 (modified) (1 diff)

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

@@ -132 +132 @@
   dpr(:)        = 0.0
   lkcoef(:,:)   = 0.0
+  DTAUKV(:,:,:) = 0.0 
 
   do K=2,L_LEVELS
@@ -365 +366 @@
 
   ! Total extinction optical depths
-
-  DO NG=1,L_NGAUSS       ! full gauss loop
+  DO NG=1,L_NGAUSS ! full gauss loop
      DO NW=1,L_NSPECTV       
+        TAUV(1,NW,NG)=0.0D0
         TAUCUMV(1,NW,NG)=0.0D0
+        
         DO K=2,L_LEVELS
            TAUCUMV(K,NW,NG)=TAUCUMV(K-1,NW,NG)+DTAUKV(K,NW,NG)
@@ -378 +380 @@
         TAUV(L,NW,NG)=TAUCUMV(2*L_NLAYRAD+1,NW,NG)
      END DO            
-  END DO                 ! end full gauss loop
-
-
-
+  END DO ! end full gauss loop
 
 end subroutine optcv

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

@@ -260 +260 @@
 
         ! Last level
-
         L              = L_NLAYRAD
         K              = 2*L+1
-        DTAUV(L,nw,ng) = DTAUKV(K,NW,NG)
+             DTAUV(L,nw,ng) = DTAUKV(K,NW,NG)
 
         WBARV(L,NW,NG) = ctemp(L,NW) / DTAUV(L,NW,NG)
-        !print*, 'TB22 : WBARV(L)=',WBARV(L,NW,NG),NW,NG
-        !print*, 'TB22 : ctemp(L)=',ctemp(L,NW),NW
-        !print*, 'TB22 : dtauv(L)=',DTAUV(L,NW,NG),NW,NG
-     END DO                 ! NW spectral loop
-  END DO                    ! NG Gauss loop
+     END DO ! NW spectral loop
+  END DO ! NG Gauss loop
 
   ! Total extinction optical depths
-
-  DO NG=1,L_NGAUSS       ! full gauss loop
+  DO NG=1,L_NGAUSS ! full gauss loop
      DO NW=1,L_NSPECTV
         TAUV(1,NW,NG)=0.0D0
+        TAUCUMV(1,NW,NG)=0.0D0
+        
         DO L=1,L_NLAYRAD
            TAUV(L+1,NW,NG)=TAUV(L,NW,NG)+DTAUV(L,NW,NG)
         END DO
 
-        TAUCUMV(1,NW,NG)=0.0D0
         DO K=2,L_LEVELS
            TAUCUMV(K,NW,NG)=TAUCUMV(K-1,NW,NG)+DTAUKV(K,NW,NG)
         END DO
      END DO
-  END DO                 ! end full gauss loop
-
+  END DO ! end full gauss loop
 
   return
 
-
 end subroutine optcv_pluto