Changeset 1238 for trunk/LMDZ.MARS/libf/aeronomars

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

File:

• trunk/LMDZ.MARS/libf/aeronomars/jthermcalc.F (modified) (1 diff)

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

@@ -158 +158 @@
       limup=1.e26
 
+      ! Ehouarn: sanity check
+      ! test that auxcolinp in monotonously increasing from 1 to nlayermx
+      do j=1,nlayermx-1
+        if (auxcolinp(j).gt.auxcolinp(j+1)) then
+          !there is a problem
+          write(*,*) "jthermcalc error: "
+          write(*,*) "auxcolinp() not increasing with altitude index!"
+          write(*,*) "j=",j," auxcolinp(j)=",auxcolinp(j)
+          write(*,*) "                auxcolinp(j+1)=",auxcolinp(j+1)
+          ! Quick fix:
+          if (j==1) then
+            auxcolinp(j)=auxcolinp(j+1)/2.
+          else
+            ! compute it as a geometric mean from encompassing values
+            auxcolinp(j)=sqrt(auxcolinp(j-1)*auxcolinp(j+1))
+          endif
+          write(*,*) " Quick fixed to auxcolinp(j)=",auxcolinp(j)
+        endif
+      enddo
 
 c     Interpolations