Changeset 1960 for LMDZ5/trunk/libf/phy1d/lmdz1d.F

Timestamp:

Feb 7, 2014, 5:20:25 PM (11 years ago)

Author:

fhourdin

Message:

Nettoyage du code 1D pour limiter les warning en mode debug.
Essentiellement des declarations inutiles.

Cleaning of the 1D code in order to limit warnings in debug mode.
Essentially supressing unused declared variables.

File:

• LMDZ5/trunk/libf/phy1d/lmdz1d.F (modified) (11 diffs)

LMDZ5/trunk/libf/phy1d/lmdz1d.F

@@ -81 +81 @@
       integer :: an
  
-!
-      real :: paire    = 1.     ! aire de la maille
-!**      common /flux_arp/fsens,flat,ok_flux_surf
-
 !---------------------------------------------------------------------
 !  Declarations related to forcing and initial profiles 
@@ -92 +88 @@
         integer llm700,nq1,nq2
         INTEGER, PARAMETER :: nlev_max=1000, nqmx=1000
-        real timestep, frac, timeit
+        real timestep, frac
         real height(nlev_max),tttprof(nlev_max),qtprof(nlev_max),
      .              uprof(nlev_max),vprof(nlev_max),e12prof(nlev_max),
@@ -100 +96 @@
      .              qprof(nlev_max,nqmx)
 
-        real    :: fff
 c        integer :: forcing_type
         logical :: forcing_les     = .false.
@@ -138 +133 @@
       real :: pzero=1.e5
       real :: play (llm),zlay (llm),sig_s(llm),plev(llm+1)
-      real :: playd(llm),zlayd(llm),ap_amma(llm+1),bp_amma(llm+1),poub
 
 !---------------------------------------------------------------------
@@ -144 +138 @@
 !---------------------------------------------------------------------
 
-      integer :: iq
       real :: phi(llm)
       real :: teta(llm),tetal(llm),temp(llm),u(llm),v(llm),w(llm)
@@ -152 +145 @@
       real :: sfdt, cfdt
       real :: du_phys(llm),dv_phys(llm),dt_phys(llm)
-      real :: du_dyn(llm),dv_dyn(llm),dt_dyn(llm)
-      real :: dt_cooling(llm),d_t_cool(llm),d_th_adv(llm)
-      real :: dq_cooling(llm),d_q_cool(llm)
-      real :: tmpvar(llm)
+      real :: dt_dyn(llm)
+      real :: dt_cooling(llm),d_th_adv(llm)
       real :: alpha
+      real :: ttt
 
       REAL, ALLOCATABLE, DIMENSION(:,:):: q
@@ -205 +197 @@
 !  Fichiers et d'autres variables
 !---------------------------------------------------------------------
-      real ttt,bow,q1
-      integer :: ierr,k,l,i,it=1,mxcalc
+      integer :: k,l,i,it=1,mxcalc
       integer jjmp1
       parameter (jjmp1=jjm+1-1/jjm)
@@ -231 +222 @@
 !---------------------------------------------------------------------
 cAl1
-        call conf_unicol(99)
+        call conf_unicol
 cAl1 moves this gcssold var from common fcg_gcssold to 
         Turb_fcg_gcssold = xTurb_fcg_gcssold
@@ -375 +366 @@
       itau_phy = 0
       call ymds2ju(annee_ref,mois,day_ref,heure,day)
-      day_ini = day
+      day_ini = int(day)
       day_end = day_ini + fnday
 
@@ -859 +850 @@
 !
         du_age(1:mxcalc)= -2.*sfdt/timestep*
-     :          (sfdt*(u(1:mxcalc)-ug(1:mxcalc)) -
-     :           cfdt*(v(1:mxcalc)-vg(1:mxcalc))  )
+     s          (sfdt*(u(1:mxcalc)-ug(1:mxcalc)) -
+     s           cfdt*(v(1:mxcalc)-vg(1:mxcalc))  )
 !!     : fcoriolis*(v(1:mxcalc)-vg(1:mxcalc))
 !
        dv_age(1:mxcalc)= -2.*sfdt/timestep*
-     :          (cfdt*(u(1:mxcalc)-ug(1:mxcalc)) +
-     :           sfdt*(v(1:mxcalc)-vg(1:mxcalc))  )
+     s          (cfdt*(u(1:mxcalc)-ug(1:mxcalc)) +
+     s           sfdt*(v(1:mxcalc)-vg(1:mxcalc))  )
 !!     : -fcoriolis*(u(1:mxcalc)-ug(1:mxcalc))
 !
@@ -879 +870 @@
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!
         u(1:mxcalc)=u(1:mxcalc) + timestep*(
-     :              du_phys(1:mxcalc)
-     :             +du_age(1:mxcalc) )            
+     s              du_phys(1:mxcalc)
+     s             +du_age(1:mxcalc) )            
         v(1:mxcalc)=v(1:mxcalc) + timestep*(
-     :               dv_phys(1:mxcalc)
-     :             +dv_age(1:mxcalc) )
+     s              dv_phys(1:mxcalc)
+     s             +dv_age(1:mxcalc) )
         q(1:mxcalc,:)=q(1:mxcalc,:)+timestep*(
-     :                dq(1:mxcalc,:)
-     :               +d_q_adv(1:mxcalc,:) )
+     s                dq(1:mxcalc,:)
+     s               +d_q_adv(1:mxcalc,:) )
 
         if (prt_level.ge.1) then