Changeset 2808 for trunk/LMDZ.MARS/libf

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/libf

Files:

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

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)