Changeset 2808

Timestamp:

Oct 20, 2022, 10:09:46 AM (2 years ago)

Author:

emillour

Message:

Mars GCM:
Bug fix in jthermcalc_e107 (wrong indexes used in some computations involving
NO2 and H2).
AM

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/aeronomars/jthermcalc_e107.F (modified) (2 diffs)
• libf/phymars/physiq_mod.F (modified) (1 diff)

trunk/LMDZ.MARS/README

@@ -3753 +3753 @@
 Add reference def files used to generate MCDv6.1 in deftank
 (as a side note MCDv6.1 was generated using r2740 of the GCM)
+
+== 20/10/2022 == AM
+Bug fix in jthermcalc_e107 (wrong indexes used in some computations involving
+NO2 and H2).

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

@@ -238 +238 @@
      $        (wm(i)*auxjo3p(ind+1) + wp(i)*auxjo3p(ind))
          !H2 interpolated coefficient
-         jfotsout(indexint,5,auxi) = jfotsout(indexint,5,auxi) *
+         jfotsout(indexint,5,auxi) = jfotsout(indexint,5,nlayer) *
      $        (wm(i)*auxjh2(ind+1) + wp(i)*auxjh2(ind))
       enddo
@@ -790 +790 @@
                !NO2 interpolated coefficient
                jfotsout(indexint,13,auxi)=
-     $              jfotsout(indexint,13,auxi) *
+     $              jfotsout(indexint,13,nlayer) *
      $              (wm(i)*auxjno2(ind+1)+wp(i)*auxjno2(ind)) * 
      $              cortemp(i)

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -2999 +2999 @@
 c        call WRITEDIAGFI(ngrid,"pplay","Pressure","Pa",3,zplay)
 c        call WRITEDIAGFI(ngrid,"pplev","Pressure","Pa",3,zplev)
+         call writediagfi(ngrid,"pphi","Geopotential","m2s-2",3,
+     &                    pphi)
+         call writediagfi(ngrid,"phisfi","Surface geopotential",
+     &                    "m2s-2",2,phisfi)
          call WRITEDIAGFI(ngrid,"tsurf","Surface temperature","K",2,
      &                  tsurf)