Changeset 1160 for trunk/LMDZ.VENUS/libf/phyvenus/physiq.F

Timestamp:

Jan 16, 2014, 9:49:42 AM (11 years ago)

Author:

slebonnois

Message:

SL: update for tracers management in phytrac, Venus

File:

• trunk/LMDZ.VENUS/libf/phyvenus/physiq.F (modified) (1 diff)

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

@@ -669 +669 @@
 
       if (iflag_trac.eq.1) then
-      call phytrac (
-     I                   itap, gmtime,
-     I                   debut,lafin,
-     I                   nqmax,
-     I                   nlon,nlev,dtime,
-     I                   u,v,t,paprs,pplay,
-     I                   rlatd,
-     I                   rlond,presnivs,pphis,pphi,
-     I                   falbe,
+
+       if (tr_scheme.eq.1) then
+! Case 1: pseudo-chemistry with relaxation toward fixed profile
+         call phytrac_relax (debut,lafin,nqmax,
+     I                   nlon,nlev,dtime,pplay,
      O                   tr_seri)
+
+       elseif (tr_scheme.eq.2) then
+! Case 2: surface emission
+! For the moment, inspired from Mars version
+! However, the variable 'source' could be used in physiq
+! so the call to phytrac_emiss could be to initialise it.
+         call phytrac_emiss ( (rjourvrai+gmtime)*RDAY,
+     I                   debut,lafin,nqmax,
+     I                   nlon,nlev,dtime,paprs,
+     I                   rlatd,rlond,
+     O                   tr_seri)
+       elseif (tr_scheme.eq.3) then
+! Case 3: Full chemistry
+!        call phytrac_chem ( ?? )
+         print*,"Chemistry not yet implemented..." 
+         print*,"See Aurelien Stolzenbach" 
+       endif
       endif