Changeset 919 for trunk/LMDZ.GENERIC

Timestamp:

Mar 28, 2013, 5:29:00 PM (12 years ago)

Author:

jleconte

Message:

28/03/2013 == JL

• small bug in previous commit found in debug mode

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

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

trunk/LMDZ.GENERIC/libf/phystd/bilinearbig.F90

@@ -30 +30 @@
       if ((x.lt.x_arr(2)).or.(x.gt.x_arr(nX-2))) then
          ind=-1
-         return
       else
         i=1

trunk/LMDZ.GENERIC/libf/phystd/optcv.F90

@@ -289 +289 @@
 
   DO NW=1,L_NSPECTV
-     DO L=1,L_NLAYRAD
+     DO L=1,L_NLAYRAD-1
         K              = 2*L+1
         atemp(L,NW) = SUM(GVAER(K,NW,1:naerkind) * TAUAEROLK(K,NW,1:naerkind))+SUM(GVAER(K+1,NW,1:naerkind) * TAUAEROLK(K+1,NW,1:naerkind))
@@ -297 +297 @@
         COSBV(L,NW,1:L_NGAUSS) = atemp(L,NW)/btemp(L,NW)
      END DO ! L vertical loop
+     
+     !last level
+     L              = L_NLAYRAD
+     K              = 2*L+1
+     atemp(L,NW)    = SUM(GVAER(K,NW,1:naerkind) * TAUAEROLK(K,NW,1:naerkind))
+     btemp(L,NW) = SUM(TAUAEROLK(K,NW,1:naerkind))
+     ctemp(L,NW) = btemp(L,NW) + 0.9999*TRAY(K,NW)
+     btemp(L,NW) = btemp(L,NW) + TRAY(K,NW)
+     COSBV(L,NW,1:L_NGAUSS) = atemp(L,NW)/btemp(L,NW)
+     
+     
   END DO                    ! NW spectral loop
 
@@ -313 +324 @@
         L              = L_NLAYRAD
         K              = 2*L+1
-        DTAUV(L,nw,ng) = DTAUKV(K,NW,NG)
-        WBARV(L,nw,ng) = ctemp(L,NW) / DTAUV(L,nw,ng)
+        DTAUV(L,nw,ng) = DTAUKV(K,NW,NG)
+
+        WBARV(L,NW,NG) = ctemp(L,NW) / DTAUV(L,NW,NG)
      END DO                 ! NW spectral loop
   END DO                    ! NG Gauss loop