Changeset 695 for trunk/LMDZ.MARS/libf/phymars/nlte_calc.F

Timestamp:

Jun 5, 2012, 1:41:39 PM (12 years ago)

Author:

acolaitis

Message:

GCM: modifications to allow for nesting compilation. Transparent to GCM user. MESOSCALE: handling of functions and modules for nesting compilation. Full nested configuration with new physics is now fully compiling.

File:

• trunk/LMDZ.MARS/libf/phymars/nlte_calc.F (modified) (9 diffs)

trunk/LMDZ.MARS/libf/phymars/nlte_calc.F

@@ -155 +155 @@
                                                             
 c vibr. temp of the bending mode :              
-      if (isot.eq.1) call interdp ( tvtbs,zyd,nzy, v626t1,zld,nl, 1 )  
-      if (isot.eq.2) call interdp ( tvtbs,zyd,nzy, v628t1,zld,nl, 1 )  
-      if (isot.eq.3) call interdp ( tvtbs,zyd,nzy, v636t1,zld,nl, 1 )  
-      if (isot.eq.4) call interdp ( tvtbs,zyd,nzy, v627t1,zld,nl, 1 )  
+      if (isot.eq.1) call interdp (tvtbs,zyd,nzy,v626t1,zld,nl,1)  
+      if (isot.eq.2) call interdp (tvtbs,zyd,nzy,v628t1,zld,nl,1)  
+      if (isot.eq.3) call interdp (tvtbs,zyd,nzy,v636t1,zld,nl,1)  
+      if (isot.eq.4) call interdp (tvtbs,zyd,nzy,v627t1,zld,nl,1)  
       
 c 2nd: correccion a la n10(i) (cantidad de absorbente en el lower state)
@@ -683 +683 @@
 c***********************************************************************
                                                             
-      call mzescape (ig,taustar21,tauinf210,tauii210,1,2,irw_mztf,imu) 
-      call mzescape (ig,taustar31,tauinf310,tauii310,1,3,irw_mztf,imu)
-      call mzescape (ig,taustar41,tauinf410,tauii410,1,4,irw_mztf,imu)
+      call mzescape (ig,taustar21,tauinf210,tauii210,1,2
+     & ,irw_mztf,imu) 
+      call mzescape (ig,taustar31,tauinf310,tauii310,1,3
+     & ,irw_mztf,imu)
+      call mzescape (ig,taustar41,tauinf410,tauii410,1,4
+     & ,irw_mztf,imu)
                                                             
       istyle = 2                                     
@@ -915 +918 @@
       enddo      
 c vibr. temp of the bending mode :  
-      if (isot.eq.1) call interdp ( tvtbs,zyd,nzy, v626t1,zld,nl, 1 )  
-      if (isot.eq.2) call interdp ( tvtbs,zyd,nzy, v628t1,zld,nl, 1 )  
-      if (isot.eq.3) call interdp ( tvtbs,zyd,nzy, v636t1,zld,nl, 1 )  
-      if (isot.eq.4) call interdp ( tvtbs,zyd,nzy, v627t1,zld,nl, 1 )  
+      if (isot.eq.1) call interdp(tvtbs,zyd,nzy, v626t1,zld,nl,1)  
+      if (isot.eq.2) call interdp(tvtbs,zyd,nzy, v628t1,zld,nl,1)  
+      if (isot.eq.3) call interdp(tvtbs,zyd,nzy, v636t1,zld,nl,1)  
+      if (isot.eq.4) call interdp(tvtbs,zyd,nzy, v627t1,zld,nl,1)  
         !if (isot.eq.5) call interdp ( tvtbs,zxd,nz, vcot1,zld,nl, 1 )  
         
@@ -2037 +2040 @@
             
 c vibr. temp of the bending mode :  
-      if (isot.eq.1) call interdp ( tvtbs,zyd,nzy, v626t1,zld,nl, 1 )  
-      if (isot.eq.2) call interdp ( tvtbs,zyd,nzy, v628t1,zld,nl, 1 )  
-      if (isot.eq.3) call interdp ( tvtbs,zyd,nzy, v636t1,zld,nl, 1 )  
-      if (isot.eq.4) call interdp ( tvtbs,zyd,nzy, v627t1,zld,nl, 1 )  
+      if (isot.eq.1) call interdp(tvtbs,zyd,nzy, v626t1,zld,nl,1)  
+      if (isot.eq.2) call interdp(tvtbs,zyd,nzy, v628t1,zld,nl,1)  
+      if (isot.eq.3) call interdp(tvtbs,zyd,nzy, v636t1,zld,nl,1)  
+      if (isot.eq.4) call interdp(tvtbs,zyd,nzy, v627t1,zld,nl,1)  
         !if (isot.eq.5) call interdp ( tvtbs,zxd,nz, vcot1,zld,nl, 1 )  
             
@@ -2533 +2536 @@
                                                 
 ! tvtbs is the bending mode of the molecule. used in xqv.   
-      if (isot.eq.1) call interdp (tvtbs,zyd,nzy, v626t1,zld,nl, 1 )  
-      if (isot.eq.2) call interdp (tvtbs,zyd,nzy, v628t1,zld,nl, 1 )  
-      if (isot.eq.3) call interdp (tvtbs,zyd,nzy, v636t1,zld,nl, 1 )  
-      if (isot.eq.4) call interdp (tvtbs,zyd,nzy, v627t1,zld,nl, 1 )  
-      if (isot.eq.5) call interdp (tvtbs,zyd,nzy, vcot1,zld,nl, 1 )   
+      if (isot.eq.1) call interdp (tvtbs,zyd,nzy,v626t1,zld,nl,1)  
+      if (isot.eq.2) call interdp (tvtbs,zyd,nzy,v628t1,zld,nl,1)  
+      if (isot.eq.3) call interdp (tvtbs,zyd,nzy,v636t1,zld,nl,1)  
+      if (isot.eq.4) call interdp (tvtbs,zyd,nzy,v627t1,zld,nl,1)  
+      if (isot.eq.5) call interdp (tvtbs,zyd,nzy,vcot1,zld,nl,1)   
       
 ! tvt0 is the lower level of the transition. used in xlower.            
       if (ib.eq.2 .or. ib.eq.3 .or. ib.eq.4 .or. ib.eq.15) then  
-         if (isot.eq.1) call interdp (tvt0,zyd,nzy, v626t1,zld,nl, 1 ) 
-         if (isot.eq.2) call interdp (tvt0,zyd,nzy, v628t1,zld,nl, 1 ) 
-         if (isot.eq.3) call interdp (tvt0,zyd,nzy, v636t1,zld,nl, 1 ) 
-         if (isot.eq.4) call interdp (tvt0,zyd,nzy, v627t1,zld,nl, 1 ) 
+         if (isot.eq.1) call interdp(tvt0,zyd,nzy,v626t1,zld,nl,1) 
+         if (isot.eq.2) call interdp(tvt0,zyd,nzy,v628t1,zld,nl,1) 
+         if (isot.eq.3) call interdp(tvt0,zyd,nzy,v636t1,zld,nl,1) 
+         if (isot.eq.4) call interdp(tvt0,zyd,nzy,v627t1,zld,nl,1) 
       elseif (ib.eq.6 .or. ib.eq.8 .or. ib.eq.10     
      @        .or. ib.eq.13 .or. ib.eq.14                  
      @        .or. ib.eq.17 .or. ib.eq.19 .or. ib.eq.20) then          
-         if (isot.eq.1) call interdp ( tvt0,zyd,nzy, v626t2,zld,nl, 1 )
-         if (isot.eq.2) call interdp ( tvt0,zyd,nzy, v628t2,zld,nl, 1 )
-         if (isot.eq.3) call interdp ( tvt0,zyd,nzy, v636t2,zld,nl, 1 )
+         if (isot.eq.1) call interdp(tvt0,zyd,nzy,v626t2,zld,nl,1)
+         if (isot.eq.2) call interdp(tvt0,zyd,nzy,v628t2,zld,nl,1)
+         if (isot.eq.3) call interdp(tvt0,zyd,nzy,v636t2,zld,nl,1)
          if (isot.eq.4) then  
             call interdp ( tvt0,zyd,nzy, v627t2,zld,nl, 1 )
@@ -2566 +2569 @@
 c tvt1 is the upper level of the transition.    
       if (ib.eq.13 .or. ib.eq.14) then 
-         if (isot.eq.1) call interdp ( tvt1,zyd,nzy, v626t4,zld,nl, 1 )
-         if (isot.eq.2) call interdp ( tvt1,zyd,nzy, v628t4,zld,nl, 1 ) 
-         if (isot.eq.3) call interdp ( tvt1,zyd,nzy, v636t4,zld,nl, 1 )
-         if (isot.eq.4) call interdp ( tvt1,zyd,nzy, v627t4,zld,nl, 1 )
+         if (isot.eq.1) call interdp(tvt1,zyd,nzy,v626t4,zld,nl,1)
+         if (isot.eq.2) call interdp(tvt1,zyd,nzy,v628t4,zld,nl,1) 
+         if (isot.eq.3) call interdp(tvt1,zyd,nzy,v636t4,zld,nl,1)
+         if (isot.eq.4) call interdp(tvt1,zyd,nzy,v627t4,zld,nl,1)
       end if
       
@@ -2586 +2589 @@
 
             !!! vt of the 3 levels in (020)  (see pag. 36a-sn1 for this)       
-            xalfa = 1.d0/2.d0* exp( dble(-ee*(nu12_1000-nu(1,2))/ty(i)) )      
-            xbeta = 1.d0/2.d0* exp( dble(-ee*(nu12_0200-nu(1,2))/ty(i)) )      
+            xalfa = 1.d0/2.d0*exp(dble(-ee*(nu12_1000-nu(1,2))/ty(i)))      
+            xbeta = 1.d0/2.d0*exp(dble(-ee*(nu12_0200-nu(1,2))/ty(i)))      
             xtv0200 = dble( - ee * nu12_0200 ) /       
      @           ( log( xbeta/(1.d0+xalfa+xbeta) ) - 
@@ -2603 +2606 @@
      @           ( log(2./(1.d0+xalfa+xbeta)) - ee*nu(1,2)/tvt0(i) ) )     
             xtv1000 = dble( - ee * nu12_1000 /       
-     @           ( log(4.*xalfa/(1.d0+xalfa+xbeta))-ee*nu(1,2)/tvt0(i)))
+     @           (log(4.*xalfa/(1.d0+xalfa+xbeta))-ee*nu(1,2)/tvt0(i)))
 
 !            if ( icurt_pop.eq.30 ) then 
@@ -4013 +4016 @@
             !write (*,*)  ' cm15um_hb/11 '                                  
             !write (*,*)  ' ib, ist, irw, imu =', ib, ist, irw_mztf, imu     
-               call mztf ( ig,cax1,cax2,cax3,v1,v2, ib,ist,irw_mztf,imu,      
+               call mztf(ig,cax1,cax2,cax3,v1,v2,ib,ist,irw_mztf,imu,
      @              itauout,icfout,itableout)                    
 !         else