Ignore:
Timestamp:
May 12, 2012, 8:10:08 AM (13 years ago)
Author:
emillour
Message:

Mars GCM:

  • updated high atmosphere photochemistry (jthermcalc.F, param_v4.h, iono.h, paramfoto_compact.F, param_read.F , thermosphere.F).
  • minor change in calchim.F90 (to not use maxloc(zycol, dim = 2) function which seems to be a problem for g95) .
  • minor bug fix in perosat.F; set tendency on pdqscloud for h2o2 to zero if none is computed.
  • in "moldiff.F", changed "tridag" to "tridag_sp", "LUBKSB" to "LUBKSB_SP" and "LUDCMP" to "LUDCMP_SP" to avoid possible conflicts with same routines defined in "moldiff_red.F". Added use of automatic-sized array in "tridag" and "tridag_sp" local array "gam" (to avoid using an a priori oversized local array).

FGG+JYC+EM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.MARS/libf/aeronomars/moldiff_red.F90

    r645 r658  
    726726
    727727      subroutine tridag(a,b,c,r,u,n)
    728       parameter (nmax=4000)
     728!      parameter (nmax=4000)
    729729!      dimension gam(nmax),a(n),b(n),c(n),r(n),u(n)
    730       real*8 gam(nmax),a(n),b(n),c(n),r(n),u(n)
     730      real*8 gam(n),a(n),b(n),c(n),r(n),u(n)
    731731      if(b(1).eq.0.)then
    732732        stop 'tridag: error: b(1)=0 !!! '
Note: See TracChangeset for help on using the changeset viewer.