Changeset 1238 for trunk/LMDZ.MARS/libf


Ignore:
Timestamp:
May 7, 2014, 12:35:26 PM (11 years ago)
Author:
emillour
Message:

Mars GCM and common dynamics:

Common dynamics:

  • correction in inidissip (only matters in Martian case)
  • added correction in addfi on theta to account for surface pressure change.

Mars GCM:
Some fixes and updates:

  • addfi (dyn3d): Add correction on theta when surface pressure changes
  • vdif_cd (phymars): Correction for coefficients in stable nighttime case
  • jthermcalc (aeronomars): Fix for some pathological cases (further investigations on the origin of these is needed)

EM

Location:
trunk/LMDZ.MARS/libf
Files:
3 edited

Legend:

Unmodified
Added
Removed

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

    r1119 r1238  
    158158      limup=1.e26
    159159
     160      ! Ehouarn: sanity check
     161      ! test that auxcolinp in monotonously increasing from 1 to nlayermx
     162      do j=1,nlayermx-1
     163        if (auxcolinp(j).gt.auxcolinp(j+1)) then
     164          !there is a problem
     165          write(*,*) "jthermcalc error: "
     166          write(*,*) "auxcolinp() not increasing with altitude index!"
     167          write(*,*) "j=",j," auxcolinp(j)=",auxcolinp(j)
     168          write(*,*) "                auxcolinp(j+1)=",auxcolinp(j+1)
     169          ! Quick fix:
     170          if (j==1) then
     171            auxcolinp(j)=auxcolinp(j+1)/2.
     172          else
     173            ! compute it as a geometric mean from encompassing values
     174            auxcolinp(j)=sqrt(auxcolinp(j-1)*auxcolinp(j+1))
     175          endif
     176          write(*,*) " Quick fixed to auxcolinp(j)=",auxcolinp(j)
     177        endif
     178      enddo
    160179
    161180c     Interpolations

  • trunk/LMDZ.MARS/libf/dyn3d/addfi.F

    r38 r1238  
    104104       ENDDO
    105105      ENDDO
    106 c
     106!***********************
     107! Correction on teta due to surface pressure changes
     108      DO k = 1,llm
     109        DO j = 1,ip1jmp1
     110           pteta(j,k)= pteta(j,k)*(1+pdpfi(j)*pdt/pps(j))**kappa
     111        ENDDO
     112      ENDDO
     113!***********************
    107114
    108115! increment covariant zonal wind

  • trunk/LMDZ.MARS/libf/phymars/vdif_cd.F

    r1236 r1238  
    195195         else
    196196c For Ri>Ric, we consider Ri->Infinity => no turbulent mixing at surface
    197             fm(ig)=0.
    198             fh(ig)=0.
     197!            fm(ig)=0.
     198!            fh(ig)=0.
     199            fm(ig)=1.
     200            fh(ig)=1.
    199201         endif
    200202c Unstable case :
Note: See TracChangeset for help on using the changeset viewer.