Changeset 1456 for LMDZ5/branches/LMDZ5V2.0-dev/libf/phylmd/phyetat0.F

Timestamp:

Nov 22, 2010, 2:16:44 PM (14 years ago)

Author:

musat

Message:

phyetat0, phyredem: correction dimension verticale pbl_tke: pbl_tke(:,1:klev+1,:)
physiq: pour pouvoir fixer la longitude solaire avec la nouvelle orbite
JYG/IM

File:

• LMDZ5/branches/LMDZ5V2.0-dev/libf/phylmd/phyetat0.F (modified) (5 diffs)

LMDZ5/branches/LMDZ5V2.0-dev/libf/phylmd/phyetat0.F

@@ -21 +21 @@
       USE infotrac
       USE traclmdz_mod,    ONLY : traclmdz_from_restart
-      USE carbon_cycle_mod,ONLY :
-     &     carbon_cycle_tr, carbon_cycle_cpl, co2_send
+      USE carbon_cycle_mod,ONLY : carbon_cycle_tr, carbon_cycle_cpl
 
       IMPLICIT none
@@ -134 +133 @@
 
        
-      clesphy0(1)=tab_cntrl( 5 )
-      clesphy0(2)=tab_cntrl( 6 )
-      clesphy0(3)=tab_cntrl( 7 )
-      clesphy0(4)=tab_cntrl( 8 )
-      clesphy0(5)=tab_cntrl( 9 )
-      clesphy0(6)=tab_cntrl( 10 )
-      clesphy0(7)=tab_cntrl( 11 )
-      clesphy0(8)=tab_cntrl( 12 )
+
+         IF( clesphy0(1).NE.tab_cntrl( 5 ) )  THEN
+             clesphy0(1)=tab_cntrl( 5 )
+         ENDIF
+
+         IF( clesphy0(2).NE.tab_cntrl( 6 ) )  THEN
+             clesphy0(2)=tab_cntrl( 6 )
+         ENDIF
+
+         IF( clesphy0(3).NE.tab_cntrl( 7 ) )  THEN
+             clesphy0(3)=tab_cntrl( 7 )
+         ENDIF
+
+         IF( clesphy0(4).NE.tab_cntrl( 8 ) )  THEN
+             clesphy0(4)=tab_cntrl( 8 )
+         ENDIF
+
+         IF( clesphy0(5).NE.tab_cntrl( 9 ) )  THEN
+             clesphy0(5)=tab_cntrl( 9 )
+         ENDIF
+
+         IF( clesphy0(6).NE.tab_cntrl( 10 ) )  THEN
+             clesphy0(6)=tab_cntrl( 10 )
+         ENDIF
+
+         IF( clesphy0(7).NE.tab_cntrl( 11 ) )  THEN
+             clesphy0(7)=tab_cntrl( 11 )
+         ENDIF
+
+         IF( clesphy0(8).NE.tab_cntrl( 12 ) )  THEN
+             clesphy0(8)=tab_cntrl( 12 )
+         ENDIF
+
 
 c
@@ -817 +841 @@
           ENDIF
           WRITE(str2,'(i2.2)') nsrf
-          CALL get_field("TKE"//str2,pbl_tke(:,1:klev,nsrf),found)
+          CALL get_field("TKE"//str2,pbl_tke(:,1:klev+1,nsrf),found)
           IF (.NOT. found) THEN
             PRINT*, "phyetat0: <TKE"//str2//"> est absent"
@@ -824 +848 @@
           xmin = 1.0E+20
           xmax = -1.0E+20
-          DO k = 1, klev
+          DO k = 1, klev+1
             DO i = 1, klon
               xmin = MIN(pbl_tke(i,k,nsrf),xmin)
@@ -1054 +1078 @@
 
          END DO
+         
          CALL traclmdz_from_restart(trs)
-
-         IF (carbon_cycle_cpl) THEN
-            ALLOCATE(co2_send(klon), stat=ierr)
-            IF (ierr /= 0) CALL abort_gcm
-     &           ('phyetat0','pb allocation co2_send',1)
-            CALL get_field("co2_send",co2_send,found)
-            IF (.NOT. found) THEN
-               PRINT*,"phyetat0: Le champ <co2_send> est absent"
-               PRINT*,"Initialisation uniforme a co2_ppm=",co2_ppm
-               co2_send(:) = co2_ppm
-            END IF
-         END IF
       END IF