Changeset 601 for trunk/LMDZ.GENERIC/libf/phystd

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/libf/phystd

Files:

• interpolateH2H2.F90 (modified) (3 diffs)
• rcm1d.F (modified) (2 diffs)

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    ========================================================