Changeset 2952 for LMDZ5/trunk/libf/phylmd/surf_land_orchidee_mod.F90

Timestamp:

Jul 14, 2017, 10:23:40 PM (7 years ago)

Author:

Laurent Fairhead

Message:

Parametrization of drag by copses
Need version 4465 of ORCHIDEE at least

6. Cheruy

Location:

LMDZ5/trunk

Files:

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

LMDZ5/trunk/libf/phylmd/surf_land_orchidee_mod.F90

@@ -3 +3 @@
 #ifndef ORCHIDEE_NOOPENMP
 #ifndef ORCHIDEE_NOZ0H
+#ifndef ORCHIDEE_NOFREIN
 !
 ! This module controles the interface towards the model ORCHIDEE.
 !
 ! Compatibility with ORCHIDIEE :
-! The current version can be used with ORCHIDEE/trunk from revision 3525.
-! This interface is used if none of the cpp keys ORCHIDEE_NOOPENMP or ORCHIDEE_NOZ0H is set.
+! The current version can be used with ORCHIDEE/trunk from revision 4465.
+! This interface is used if none of the cpp keys ORCHIDEE_NOOPENMP, 
+! ORCHIDEE_NOZ0H or ORCHIDEE_NOFREIN is set.
 !
 ! Subroutines in this module : surf_land_orchidee
@@ -43 +45 @@
        evap, fluxsens, fluxlat, &              
        tsol_rad, tsurf_new, alb1_new, alb2_new, &
-       emis_new, z0m_new, z0h_new, qsurf)
+       emis_new, z0m_new, z0h_new, qsurf, &
+       veget, lai, height )
+
 
     USE mod_surf_para
@@ -108 +112 @@
 !
     INCLUDE "YOMCST.h"
+    INCLUDE "dimpft.h"
+
+
   
 !
@@ -138 +145 @@
     REAL, DIMENSION(klon), INTENT(OUT)        :: alb1_new, alb2_new
     REAL, DIMENSION(klon), INTENT(OUT)        :: emis_new, z0m_new, z0h_new
+    REAL, DIMENSION(klon,nvm_lmdz), INTENT(OUT) :: veget
+    REAL, DIMENSION(klon,nvm_lmdz), INTENT(OUT) :: lai
+    REAL, DIMENSION(klon,nvm_lmdz), INTENT(OUT) :: height
+
 
 ! Local
@@ -409 +420 @@
                precip_rain, precip_snow, lwdown, swnet, swdown, ps, &
                evap, fluxsens, fluxlat, coastalflow, riverflow, &
-               tsol_rad, tsurf_new, qsurf, albedo_out, emis_new, z0m_new, &
-               lon_scat, lat_scat, q2m, t2m, z0h_new)
+               tsol_rad, tsurf_new, qsurf, albedo_out, emis_new, z0m_new, &    
+               lon_scat, lat_scat, q2m, t2m, z0h_new, nvm_orch)
 #endif         
        ENDIF
@@ -426 +437 @@
     IF (knon > 0) THEN
 #ifdef CPP_VEGET    
+       IF (nvm_orch .NE. nvm_lmdz ) THEN
+          abort_message='Pb de dimensiosn PFT: nvm_orch et nvm_lmdz differents.'
+          CALL abort_physic(modname,abort_message,1)
+       ENDIF
+
        CALL intersurf_main_gathered (itime+itau_phy, nbp_lon, nbp_lat, knon, ktindex, dtime,  &
             lrestart_read, lrestart_write, lalo, &
@@ -434 +450 @@
             evap(1:knon), fluxsens(1:knon), fluxlat(1:knon), coastalflow(1:knon), riverflow(1:knon), &
             tsol_rad(1:knon), tsurf_new(1:knon), qsurf(1:knon), albedo_out(1:knon,:), emis_new(1:knon), z0m_new(1:knon), &
-            lon_scat, lat_scat, q2m, t2m, z0h_new(1:knon), coszang=yrmu0(1:knon))
+            lon_scat, lat_scat, q2m, t2m, z0h_new(1:knon),&
+            veget(1:knon,:),lai(1:knon,:),height(1:knon,:),&
+            coszang=yrmu0(1:knon))
 #endif       
     ENDIF
@@ -664 +682 @@
 #endif
 #endif
+#endif
 END MODULE surf_land_orchidee_mod