Changeset 2040 for trunk

Timestamp:

Nov 13, 2018, 6:03:56 PM (6 years ago)

Author:

jvatant

Message:

Fix a bug on a haze opt properties in case of vertical extension.
--JVO

File:

• trunk/LMDZ.TITAN/libf/phytitan/disr_haze.F90 (modified) (3 diffs)

trunk/LMDZ.TITAN/libf/phytitan/disr_haze.F90

@@ -106 +106 @@
 else if ((ip.ne.nblev_PL) .and. (iw.eq.nbwl_PL)) then
 
-taeros =  ext_PL(ip,iw)**(1.-factp)  * ext_PL(ip+1,iw)**factp
-        
+taeros =  ext_PL(ip,iw)**(1.-factp) * ext_PL(ip+1,iw)**factp
 ssa    =  ssa_PL(ip,iw)**(1.-factp) * ssa_PL(ip+1,iw)**factp
-        
 cbar   =  asf_PL(ip,iw)**(1.-factp) * asf_PL(ip+1,iw)**factp
 
@@ -132 +130 @@
            / ( ext_PL(ip-5,iw)*(1.-factw) + ext_PL(ip-5,iw+1)  *factw ) )
            
-     taeros = ext_PL(ip,iw)*(1.-factw) + ext_PL(ip,iw+1)  *factw
+     taeros =  ext_PL(ip,iw)*(1.-factw) + ext_PL(ip,iw+1)*factw
+     ssa    =  ssa_PL(ip,iw)*(1.-factw) + ssa_PL(ip,iw+1)*factw
+     cbar   =  asf_PL(ip,iw)*(1.-factw) + asf_PL(ip,iw+1)*factw
 
   else if (iw.eq.nbwl_PL) then
@@ -138 +138 @@
      fact_t = log10( ext_PL(ip,iw) / ext_PL(ip-5,iw) )
     
-     taeros = ext_PL(ip,iw)
+     taeros =  ext_PL(ip,iw)
+     ssa    =  ssa_PL(ip,iw)
+     cbar   =  asf_PL(ip,iw)
 
   endif