Changeset 1945 for trunk/LMDZ.VENUS/libf

Timestamp:

Jun 8, 2018, 10:40:22 AM (7 years ago)

Author:

slebonnois

Message:

Correction of a bug in radlwsw Venus

File:

• trunk/LMDZ.VENUS/libf/phyvenus/radlwsw.F (modified) (2 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/radlwsw.F

@@ -381 +381 @@
        
 c -------- Probleme aux bords
-       if ((ips0.eq.0).and.(psurfve(1,lat).gt.paprs(j,1))) then
-              ips0  = 1
-              print*,'Extrapolation! ig=',j,'  ips0=',ips0
-              factp = (paprs(j,1)         -psurfve(ips0,lat))
-     .               /(psurfve(ips0+1,lat)-psurfve(ips0,lat))
-       endif
-       if ((ips0.eq.0).and.(psurfve(nbpsve(lat),lat).le.paprs(j,1)))
+c surf press lower than the lowest surf pres in matrices
+       if ((ips0.eq.0).and.(psurfve(nbpsve(lat),lat).gt.paprs(j,1)))
      . then
               ips0  = nbpsve(lat)-1
@@ -394 +389 @@
      .               /(psurfve(ips0+1,lat)-psurfve(ips0,lat))
        endif
+c surf press higher than the highest surf pres in matrices
+       if ((ips0.eq.0).and.(psurfve(1,lat).le.paprs(j,1))) then
+              ips0  = 1
+              print*,'Extrapolation! ig=',j,'  ips0=',ips0
+              factp = (paprs(j,1)         -psurfve(ips0,lat))
+     .               /(psurfve(ips0+1,lat)-psurfve(ips0,lat))
+       endif
+
+c this has to be done for ips0A and ips0B also...
+       if (latA.eq.lat) then
+         ips0A = ips0
+       else
+         if (latA.gt.0) then
+          if ((ips0A.eq.0).and.
+     .        (psurfve(nbpsve(latA),latA).gt.paprs(j,1))) then
+              ips0A  = nbpsve(latA)-1
+          endif
+          if ((ips0A.eq.0).and.(psurfve(1,latA).le.paprs(j,1))) then
+              ips0A  = 1
+          endif
+         endif
+       endif
+       if (latB.eq.lat) then
+         ips0B = ips0
+       else
+         if (latB.gt.0) then
+          if ((ips0B.eq.0).and.
+     .        (psurfve(nbpsve(latB),latB).gt.paprs(j,1))) then
+              ips0B  = nbpsve(latB)-1
+          endif
+          if ((ips0B.eq.0).and.(psurfve(1,latB).le.paprs(j,1))) then
+              ips0B  = 1
+          endif
+         endif
+       endif       
+
 c ---------