Changeset 2213 for trunk/LMDZ.MARS/libf/aeronomars/jthermcalc.F

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

File:

• trunk/LMDZ.MARS/libf/aeronomars/jthermcalc.F (modified) (2 diffs)

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