Changeset 1279 for LMDZ4/trunk/libf/phylmd/surf_land_orchidee_mod.F90

Timestamp:

Dec 10, 2009, 10:02:56 AM (14 years ago)

Author:

Laurent Fairhead

Message:

Merged LMDZ4-dev branch changes r1241:1278 into the trunk
Running trunk and LMDZ4-dev in LMDZOR configuration on local
machine (sequential) and SX8 (4-proc) yields identical results
(restart and restartphy are identical binarily)
Log history from r1241 to r1278 is available by switching to
source:LMDZ4/branches/LMDZ4-dev-20091210

Location:

LMDZ4/trunk

Files:

• . (modified) (1 prop)
• libf/phylmd/surf_land_orchidee_mod.F90 (modified) (5 diffs)

LMDZ4/trunk/libf/phylmd/surf_land_orchidee_mod.F90

@@ -40 +40 @@
        tsol_rad, tsurf_new, alb1_new, alb2_new, &
        emis_new, z0_new, qsurf)
-   USE mod_surf_para
-   USE mod_synchro_omp
+
+    USE mod_surf_para
+    USE mod_synchro_omp
    
+USE carbon_cycle_mod, ONLY : carbon_cycle_cpl, fco2_land_inst, fco2_lu_inst
+
 !    
 ! Cette routine sert d'interface entre le modele atmospherique et le 
@@ -67 +70 @@
 !   spechum      humidite specifique 1ere couche
 !   epot_air     temp pot de l'air
-!   ccanopy      concentration CO2 canopee
+!   ccanopy      concentration CO2 canopee, correspond au co2_send de 
+!                carbon_cycle_mod ou valeur constant co2_ppm
 !   tq_cdrag     cdrag
 !   petAcoef     coeff. A de la resolution de la CL pour t
@@ -134 +138 @@
     INTEGER                                   :: error
     REAL, DIMENSION(klon)                     :: swdown_vrai
+    REAL, DIMENSION(klon)                     :: fco2_land_comp  ! sur grille compresse
+    REAL, DIMENSION(klon)                     :: fco2_lu_comp    ! sur grille compresse
     CHARACTER (len = 20)                      :: modname = 'surf_land_orchidee'
     CHARACTER (len = 80)                      :: abort_message
@@ -341 +347 @@
           CALL abort_gcm(modname,abort_message,1)
        ENDIF
-
+!
+! Allocate variables needed for carbon_cycle_mod
+!
+       IF (carbon_cycle_cpl) THEN
+          IF (.NOT. ALLOCATED(fco2_land_inst)) THEN
+             ALLOCATE(fco2_land_inst(klon),stat=error)
+             IF (error /= 0)  CALL abort_gcm(modname,'Pb in allocation fco2_land_inst',1)
+             
+             ALLOCATE(fco2_lu_inst(klon),stat=error)
+             IF(error /=0) CALL abort_gcm(modname,'Pb in allocation fco2_lu_inst',1)
+          END IF
+       END IF
+       
     ENDIF                          ! (fin debut) 
+ 
 
 ! 
@@ -443 +462 @@
     
     IF (debut) CALL Finalize_surf_para
-    
+
+    
+! JG : TEMPORAIRE!!!! Les variables fco2_land_comp et fco2_lu_comp seront plus tard en sortie d'ORCHIDEE
+!      ici mis a valeur quelquonque pour test. Ces variables sont sur la grille compresse avec uniquement des points de terres
+
+    fco2_land_comp(:) = 1.
+    fco2_lu_comp(:)   = 10.
+
+! Decompress variables for the module carbon_cycle_mod
+    IF (carbon_cycle_cpl) THEN
+       fco2_land_inst(:)=0.
+       fco2_lu_inst(:)=0.
+       
+       DO igrid = 1, knon
+          ireal = knindex(igrid)
+          fco2_land_inst(ireal) = fco2_land_comp(igrid)
+          fco2_lu_inst(ireal)   = fco2_lu_comp(igrid)
+       END DO
+    END IF
+
   END SUBROUTINE surf_land_orchidee
 !