Changeset 1668 for trunk/LMDZ.TITAN/libf/phytitan/initracer.F

Timestamp:

Feb 17, 2017, 3:44:15 PM (8 years ago)

Author:

jvatant

Message:

Another round of cleaning of dust and dummy tracers
JVO

File:

• trunk/LMDZ.TITAN/libf/phytitan/initracer.F (modified) (5 diffs)

trunk/LMDZ.TITAN/libf/phytitan/initracer.F

@@ -8 +8 @@
 c   --------
 c   Initialization related to tracer 
-c   (transported dust, water, chemical species, ice...)
+c   (chemical species, ice...)
 c
 c   Name of the tracer
@@ -38 +38 @@
 c  alpha_lift(nq)  ! saltation vertical flux/horiz flux ratio (m-1)
 c  alpha_devil(nq) ! lifting coeeficient by dust devil
-c  doubleq           ! if method with mass (iq=1) and number(iq=2) mixing ratio
-c  varian            ! Characteristic variance of log-normal distribution
 c-----------------------------------------------------------------------
 
@@ -53 +51 @@
        ALLOCATE(alpha_devil(nq))
        ALLOCATE(qextrhor(nq))
-       ALLOCATE(igcm_dustbin(nq))
        !! initialization
        alpha_lift(:)=0.
@@ -72 +69 @@
       ! 0. initialize tracer indexes to zero:
       ! NB: igcm_* indexes are commons in 'tracer.h'
-      igcm_co=0
-      igcm_o=0
-      igcm_o1d=0
-      igcm_o2=0
-      igcm_o3=0
-      igcm_h=0
-      igcm_h2=0
-      igcm_oh=0
-      igcm_ho2=0
-      igcm_h2o2=0
       igcm_n2=0
-      igcm_ar=0
-      igcm_ar_n2=0
 
       write(*,*) 'initracer: noms() ', noms
 
 
-      !print*,'Setting dustbin = 0 in initracer.F'
-      !dustbin=0
+      ! 1. find chemistry tracers
+      count = 0.    
 
-      ! 1. find dust tracers
-      count=0
-!      if (dustbin.gt.0) then
-!        do iq=1,nq
-!          txt=" "
-!          write(txt,'(a4,i2.2)')'dust',count+1   
-!          if (noms(iq).eq.txt) then
-!            count=count+1
-!            igcm_dustbin(count)=iq
-!            mmol(iq)=100.
-!          endif
-!        enddo !do iq=1,nq
-!      endif ! of if (dustbin.gt.0)
-
-
-!      if (doubleq) then
-!        do iq=1,nq
-!          if (noms(iq).eq."dust_mass") then
-!            igcm_dust_mass=iq
-!            count=count+1
-!          endif
-!          if (noms(iq).eq."dust_number") then
-!            igcm_dust_number=iq
-!            count=count+1
-!          endif
-!        enddo
-!      endif ! of if (doubleq)
-      ! 2. find chemistry and water tracers
-     
+ 
       ! check that we identified all tracers:
       if (count.ne.nq) then
@@ -150 +107 @@
       write(*,*) 'alpha_devil = ', alpha_devil
       write(*,*) 'radius  = ', radius
-!      if(doubleq) then
-!        write(*,*) 'reff_lift (um) =  ', reff_lift
-!        write(*,*) 'size distribution variance  = ', varian
-!        write(*,*) 'r3n_q , ref_r0 : ', r3n_q , ref_r0
-!      end if
+      
       write(*,*) 'Qext  = ', qext 
       write(*,*)