Changeset 1457

Timestamp:

Nov 22, 2010, 4:19:25 PM (14 years ago)

Author:

musat

Message:

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

Location:

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

Files:

• phyetat0.F (modified) (3 diffs)
• phyredem.F (modified) (3 diffs)
• physiq.F (modified) (4 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
+      USE carbon_cycle_mod,ONLY :
+     &     carbon_cycle_tr, carbon_cycle_cpl, co2_send
 
       IMPLICIT none
@@ -133 +134 @@
 
        
-
-         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
-
+      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 )
 
 c
@@ -1078 +1054 @@
 
          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
 

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

@@ -11 +11 @@
       USE pbl_surface_mod,  ONLY : pbl_surface_final
       USE phys_state_var_mod
-      USE iostart  
+      USE iostart
       USE traclmdz_mod, ONLY : traclmdz_to_restart
       USE infotrac
       USE control_mod
+      USE carbon_cycle_mod, ONLY : carbon_cycle_cpl, co2_send
 
       IMPLICIT none
@@ -97 +98 @@
       CALL put_var("controle","Parametres de controle",tab_cntrl)
 c
+
       CALL put_field("longitude",
      .               "Longitudes de la grille physique",rlon)
@@ -335 +337 @@
             CALL put_field("trs_"//tname(iiq),"",trs(:,it))
          END DO
+         IF (carbon_cycle_cpl) THEN
+            IF (.NOT. ALLOCATED(co2_send)) THEN
+               ! This is the case of create_etat0_limit, ce0l
+               ALLOCATE(co2_send(klon))
+               co2_send(:) = co2_ppm0
+            END IF
+            CALL put_field("co2_send","co2_ppm for coupling",co2_send)
+         END IF
       END IF
 

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

@@ -1250 +1250 @@
 cym Attention pbase pas initialise dans concvl !!!!
           pbase=0
-          paire_ter(:)=0.    
 cIM 180608
 c         pmflxr=0.
@@ -1849 +1848 @@
 !!!   jyg 17 Sep 2010 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
       call ymds2ju(year_cur, mth_eq, day_eq,0., jD_eq)
-      day_since_equinox = (jD_cur + jH_cur) - jD_eq
+      day_since_equinox = (jD_cur + jH_cur) - jD_eq 
 !
-!   choix entre calcul de la longitude solaire vraie ou valeur fixee a
+!   choix entre calcul de la longitude solaire vraie ou valeur fixee a 
 !   solarlong0
       if (solarlong0<-999.) then
@@ -1857 +1856 @@
 ! calcul selon la routine utilisee pour les planetes
         call solarlong(day_since_equinox, zlongi, dist)
-       else
+       else    
 ! calcul selon la routine utilisee pour l'AR4
         CALL orbite(REAL(days_elapsed+1),zlongi,dist)
-       endif
-      else
+       endif   
+      else    
            zlongi=solarlong0  ! longitude solaire vraie
-           dist=1.            ! distance au soleil / moyenne
-      endif
-!!!
-
+           dist=1.            ! distance au soleil / moyenne 
+      endif   
       if(prt_level.ge.1)                                                &
      &    write(lunout,*)'Longitude solaire ',zlongi,solarlong0,dist
@@ -3379 +3376 @@
      I                   cdragh,coefh,u1,v1,ftsol,pctsrf,
      I                   frac_impa, frac_nucl,
-     I                   pphis,airephy,dtime,itap)
+     I                   pphis,airephy,dtime,itap,
+     I                   rlon,rlat,qx(:,:,ivap),da,phi,mp,upwd,dnwd)