Changeset 2213

Timestamp:

Jan 8, 2020, 12:01:49 PM (5 years ago)

Author:

emillour

Message:

Mars GCM:
Bug fix in calchim (only had an impact when computations include the
ionosphere), parameters should not be reset at every call.
EM

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/aeronomars/calchim_mod.F90 (modified) (2 diffs)
• libf/aeronomars/jthermcalc.F (modified) (2 diffs)

trunk/LMDZ.MARS/README

@@ -2814 +2814 @@
 Manual User Update: Old stuff removed from part4: Running the model
 
-
+== 08/01/2020 == EM
+Bug fix in calchim (only had an impact when computations include the
+ionosphere), parameters should not be reset at every call.
+
+

trunk/LMDZ.MARS/libf/aeronomars/calchim_mod.F90

@@ -191 +191 @@
 !     for output:
 
-      logical :: output             ! to issue calls to writediagfi and stats
+      logical,save :: output        ! to issue calls to writediagfi and stats
       real :: jo3_3d(ngrid,nlayer)  ! Photodissociation rate O3->O1D (s-1)
       real :: jh2o_3d(ngrid,nlayer)  ! Photodissociation rate H2O->H+OH (s-1)
@@ -200 +200 @@
 
 !=======================================================================
+!     initialization of the chemistry (first call only)
+!=======================================================================
+
+      if (firstcall) then
+
+!=======================================================================
 !     main dashboard for the chemistry
 !=======================================================================
 
-      unichim = .true.     ! true : unified chemistry ! false : separate models in lower and upper atmosphere
-      jonline = .true.     ! true : on-line calculation of photodissociation rates ! false : lookup table
-      ionchem = .false.    ! switch for ionospheric chemistry
-      depos   = .false.    ! switch for dry deposition
-      output  = .false.    ! true : triggers writing of specific outputs (photolysis rates, emission rates, etc)
-
-!=======================================================================
-!     initialization of the chemistry (first call only)
-!=======================================================================
-
-      if (firstcall) then
+        unichim = .true.     ! true : unified chemistry ! false : separate models in lower and upper atmosphere
+        jonline = .true.     ! true : on-line calculation of photodissociation rates ! false : lookup table
+        ionchem = .false.    ! switch for ionospheric chemistry
+        depos   = .false.    ! switch for dry deposition
+        output  = .false.    ! true : triggers writing of specific outputs (photolysis rates, emission rates, etc)
 
          if (photochem) then

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

@@ -1377 +1377 @@
 C escout(nlayer) on pressure grid p(nlayer).
 C
-      real*8 wm(nlayer),wp(nlayer),p(nlayer)
-      integer nm(nlayer)
-      real*8 pin(nl)
-      real*8 limup,limdown
-      integer nl,nlayer,n1,n,np,nini
+      real*8,intent(out) :: wm(nlayer),wp(nlayer) ! interpolation weights
+      integer,intent(out) :: nm(nlayer) ! index of nearest point
+      real*8,intent(in) :: pin(nl),p(nlayer)
+      real*8,intent(in) :: limup,limdown
+      integer,intent(in) :: nl,nlayer
+      integer :: n1,n,np,nini
       nini=1
       do n1=1,nlayer
@@ -1400 +1401 @@
          endif
       enddo
-      return
+
       end