Changeset 1301 for trunk/LMDZ.VENUS/libf/phyvenus/phytrac_emiss.F

Timestamp:

Jun 26, 2014, 12:25:31 PM (11 years ago)

Author:

slebonnois

Message:

SL: many bug corrections in phyvenus, some cleaning, and a new ksi matrix format for Venus IR

File:

• trunk/LMDZ.VENUS/libf/phyvenus/phytrac_emiss.F (modified) (4 diffs)

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

@@ -77 +77 @@
       integer,parameter :: nbsrc=2,nblat=5,nblon=4
       integer,parameter :: Nemiss=1   ! duree emission (Ed)
-      real,save :: source_volcan(nbsrc)
+      integer,save :: Nemiss(nbsrc)      ! duration emission (Ed)
+      real,save :: source_volcan(nbsrc)  ! flux emission (kg/s)
       real,save :: lat_volcan(nblat),lon_volcan(nblon)
       real,save :: area_emiss(nblat,nblon)
@@ -131 +132 @@
          source_volcan(1) = 1.
          source_volcan(2) = 1000.
+c duration in Ed
+         Nemiss(1) = 1
+         Nemiss(2) = 10
 c localisation volcan
          lat_volcan(1) =  70.
@@ -173 +177 @@
            do ilat  = 1,nblat
             do ilon  = 1,nblon
-             it=min(iemiss*ilat*ilon,nqtot)
+             it=(iemiss-1)*nblat*nblon+(ilat-1)*nblon+ilon
+             it=min(it,nqtot)
+             deltatr(i,1,it) = 0.
 
              if (i .eq. ig_volcan(ilat,ilon)) then
 
 c source appliquee pendant Nemiss Ed
-               if (timesimu .lt. 86400.*Nemiss) then
+               if (timesimu .lt. 86400.*Nemiss(iemiss)) then
 
 c source en kg/kg/s
@@ -185 +191 @@
            tr_seri(i,1,it) = tr_seri(i,1,it) + deltatr(i,1,it)*pdtphys
 
-               else
-           deltatr(i,1,it) = 0.
                end if  ! duree emission
-
              end if ! i localisation
             end do