Changeset 1331

Timestamp:

Aug 26, 2014, 2:46:46 PM (10 years ago)

Author:

emillour

Message:

Mars GCM:

• Bug fix: Enforce recomputation of reaction rates at every call of routine prodsandlosses (in paramfoto_compact.F)

FGG

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/aeronomars/paramfoto_compact.F (modified) (4 diffs)

trunk/LMDZ.MARS/README

@@ -2137 +2137 @@
 the physics. This is to be consistent with what is done in LMDZ.COMMON.
 
+== 26/08/2014 == FGG
+- Bug fix: Enforce recomputation of reaction rates at every call of
+  routine prodsandlosses (in paramfoto_compact.F)

trunk/LMDZ.MARS/libf/aeronomars/paramfoto_compact.F

@@ -1283 +1283 @@
          jion(2,i,1)=jion(2,i,1) +
      $        jfotsout(inter,2,i) * fluxtop(inter) * (1.-efdiso2(inter))
-!         if(ig.eq.1.and.i.eq.39)write(*,*)'paramfoto_compact/1355',
-!     $        jion(2,i,1),jfotsout(inter,2,i),fluxtop(inter)
          
 c     O3P
@@ -2971 +2969 @@
 ccccccccccccccc CODE STARTS 
 
-      !Initialization
-      if (firstcall) then
-        firstcall= .false.
+      !Initialization'
+!      if (firstcall) then
+!        firstcall= .false.
         do j=1,nreact
            Lco2(i,j) = 0.d0
@@ -3089 +3087 @@
         Pelecttot(i) = 0.0d0
         Lelecttot(i) = 0.0d0
-      endif
-
-
+!      endif
+
+      
 
       !!! Productions and losses reaction by reaction
@@ -3567 +3565 @@
          Pnoplus(i,62)     = ch62 *  co2plusxinput * noxinput
          Pco2(i,62)        = Pnoplus(i,62) !ch62 *  co2plusxinput * noxinput
-
+         
 
 c     R63 :   CO2+  +  N    -->  NO  + CO+