Changeset 815 for LMDZ4/trunk/libf/phytherm/pbl_surface_mod.F90

Timestamp:

Jul 14, 2007, 10:31:47 PM (17 years ago)

Author:

Laurent Fairhead

Message:

Mise a jour des routines de phytherm vers la version V3 de LMDZ4
LF

File:

• LMDZ4/trunk/libf/phytherm/pbl_surface_mod.F90 (modified) (5 diffs)

LMDZ4/trunk/libf/phytherm/pbl_surface_mod.F90

@@ -13 +13 @@
   USE mod_phys_lmdz_para,  ONLY : mpi_size
   USE ioipsl
-  USE surface_data,        ONLY : ocean, ok_veget, debug_surf, newwind
+  USE surface_data,        ONLY : ocean, ok_veget
   USE surf_land_mod,       ONLY : surf_land
   USE surf_landice_mod,    ONLY : surf_landice
@@ -22 +22 @@
   USE climb_wind_mod,      ONLY : climb_wind_down, climb_wind_up
   USE coef_diff_turb_mod,  ONLY : coef_diff_turb
-
-!jg+ temporary
-  USE mod_clvent, ONLY : clvent, save_flux
-!jg-
 
   IMPLICIT NONE
@@ -455 +451 @@
 
 
-!jg+ temporary test
-    REAL, DIMENSION(klon,klev)         :: y_flux_u_old, y_flux_v_old
-    REAL, DIMENSION(klon,klev)         :: y_d_u_old, y_d_v_old
-!jg-
-    
 !****************************************************************************************
 ! End of declarations
@@ -698 +689 @@
             ycoefm, ycoefh)
        
-!jg+
-!****************************************************************************************
-!    => Old method
-! Calculer la diffusion des vitesses "u" et "v"   
-! Output can be used : y_d_u_old, y_flux_u_old, y_d_v_old, y_flux_v_old
-!
-!****************************************************************************************
-       CALL clvent(knon,dtime,yu1,yv1,ycoefm,yt,yu,ypaprs,ypplay,ydelp, &
-            y_d_u_old, y_flux_u_old)
-       CALL clvent(knon,dtime,yu1,yv1,ycoefm,yt,yv,ypaprs,ypplay,ydelp, &
-            y_d_v_old, y_flux_v_old)
-       ! save_flux est utile pour pouvoir utilise calcul_flux_vent plus tard
-       CALL save_flux(klon, y_flux_u_old(:,1), y_flux_v_old(:,1))
-!jg-
-
 !****************************************************************************************
 ! 
@@ -848 +824 @@
             y_flux_u, y_flux_v, y_d_u, y_d_v)
 
-!jg+ temporary for testing 
-! Use the results from old method
-       IF (.NOT. newwind) THEN
-          y_flux_u(:,:) = y_flux_u_old(:,:)
-          y_flux_v(:,:) = y_flux_v_old(:,:)
-          y_d_u(:,:)    = y_d_u_old(:,:)
-          y_d_v(:,:)    = y_d_v_old(:,:)
-       ENDIF
-!jg-
-   
        DO j = 1, knon
           y_dflux_t(j) = y_dflux_t(j) * ypct(j)