Changeset 601 for trunk/LMDZ.GENERIC

Timestamp:

Mar 27, 2012, 6:27:23 PM (13 years ago)

Author:

aslmd

Message:

GENERIC: minor modifications to jupiter def. plus output profile at each timestep for possible restarts and tests. plus allowed the model to keep on integrations if out of allowed continuum temperature, a lot of warnings are given to the user.

Location:

trunk/LMDZ.GENERIC

Files:

• deftank/jupiter1d.callphys.def (modified) (1 diff)
• deftank/jupiter1d.rcm1d.def (modified) (1 diff)
• libf/phystd/interpolateH2H2.F90 (modified) (3 diffs)
• libf/phystd/rcm1d.F (modified) (2 diffs)

trunk/LMDZ.GENERIC/deftank/jupiter1d.callphys.def

@@ -141 +141 @@
 ## CO2 options 
 ## ~~~~~~~~~~~
-# gas is non-ideal CO2 ?
-nonideal      = .false.
 # call CO2 condensation ?
 co2cond       = .false.

trunk/LMDZ.GENERIC/deftank/jupiter1d.rcm1d.def

@@ -48 +48 @@
 # Simulate global averaged conditions ? 
 global1d       = .true.
-# Latitude (deg) [only used if force_cpp = F]
+# Latitude (deg) [only used if global1d = F]
 latitude       = 0.0
-# Solar Zenith angle (deg) [only used if force_cpp = T]
+# Solar Zenith angle (deg) [only used if global1d = T]
 szangle        = 45. 
 # Force specific heat capacity and molecular mass values

trunk/LMDZ.GENERIC/libf/phystd/interpolateH2H2.F90

@@ -180 +180 @@
          write(*,*) 'Warning from bilinearH2H2:'
          write(*,*) 'Outside continuum temperature range!'
+         write(*,*) y,y_arr(1),y_arr(nY)
          if(y.lt.y_arr(1))then
             y=y_arr(1)+0.01
@@ -186 +187 @@
             y=y_arr(nY)-0.01
          endif
-      else
+      endif
+      !else
 
 !     in the y (temperature) direction 2nd
@@ -200 +202 @@
             goto 20
          endif
-      endif
+      !endif
       
       f11=f2d_arr(a,b)

trunk/LMDZ.GENERIC/libf/phystd/rcm1d.F

@@ -117 +117 @@
 
       saveprofile=.false.
+      saveprofile=.true.
 
 c ----------------------------------------
@@ -806 +807 @@
         END DO
 
-      ENDDO   ! fin de la boucle temporelle
-
 c    ========================================================
 c    GESTION DES SORTIE
 c    ========================================================
-
-! save temperature profile
       if(saveprofile)then
          OPEN(12,file='profile.out',form='formatted')
+         write(12,*) tsurf
          DO ilayer=1,nlayermx
-            write(12,*) temp(ilayer), play(ilayer)
+            write(12,*) temp(ilayer) !, play(ilayer) !AS12 only temp so that iprofile=8 can be used
          ENDDO
          CLOSE(12)
       endif
 
+
+      ENDDO   ! fin de la boucle temporelle
 
 c    ========================================================