Changeset 803 for LMDZ4/trunk
- Timestamp:
- Jul 3, 2007, 5:18:23 PM (17 years ago)
- Location:
- LMDZ4/trunk/libf/phylmd
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/trunk/libf/phylmd/ocean_cpl_mod.F90
r793 r803 12 12 USE calcul_fluxs_mod, ONLY : calcul_fluxs 13 13 USE climb_wind_mod, ONLY : calcul_wind_flux 14 USE mod_clvent, ONLY : calcul_flux_vent15 USE surface_data, ONLY : newwind16 14 17 15 IMPLICIT NONE … … 160 158 tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l) 161 159 162 !jg : temporary test, only calcul_flux_wind should be saved in the future163 160 ! Calcultate the flux of u and v at surface 164 IF (newwind) THEN 165 ! calculate flux with new methode 166 CALL calcul_wind_flux(knon, dtime, taux, tauy) 167 ELSE 168 ! use results from old clvent. jg: TO BE REMOVED 169 CALL calcul_flux_vent(klon, taux, tauy) 170 ENDIF 171 161 CALL calcul_wind_flux(knon, dtime, taux, tauy) 162 172 163 173 164 !**************************************************************************************** … … 312 303 313 304 ! Calcultate the flux of u and v at surface 314 IF (newwind) THEN 315 ! calculate flux with new methode 316 CALL calcul_wind_flux(knon, dtime, taux, tauy) 317 ELSE 318 ! use results from old clvent. jg: TO BE REMOVED 319 CALL calcul_flux_vent(klon, taux, tauy) 320 END IF 321 305 CALL calcul_wind_flux(knon, dtime, taux, tauy) 306 322 307 !**************************************************************************************** 323 308 ! Flux ocean-atmosphere useful for "slab" ocean but here calculated only for printing -
LMDZ4/trunk/libf/phylmd/pbl_surface_mod.F90
r793 r803 13 13 USE mod_phys_lmdz_para, ONLY : mpi_size 14 14 USE ioipsl 15 USE surface_data, ONLY : ocean, ok_veget , debug_surf, newwind15 USE surface_data, ONLY : ocean, ok_veget 16 16 USE surf_land_mod, ONLY : surf_land 17 17 USE surf_landice_mod, ONLY : surf_landice … … 22 22 USE climb_wind_mod, ONLY : climb_wind_down, climb_wind_up 23 23 USE coef_diff_turb_mod, ONLY : coef_diff_turb 24 25 !jg+ temporary26 USE mod_clvent, ONLY : clvent, save_flux27 !jg-28 24 29 25 IMPLICIT NONE … … 455 451 456 452 457 !jg+ temporary test458 REAL, DIMENSION(klon,klev) :: y_flux_u_old, y_flux_v_old459 REAL, DIMENSION(klon,klev) :: y_d_u_old, y_d_v_old460 !jg-461 462 453 !**************************************************************************************** 463 454 ! End of declarations … … 698 689 ycoefm, ycoefh) 699 690 700 !jg+701 !****************************************************************************************702 ! => Old method703 ! Calculer la diffusion des vitesses "u" et "v"704 ! Output can be used : y_d_u_old, y_flux_u_old, y_d_v_old, y_flux_v_old705 !706 !****************************************************************************************707 CALL clvent(knon,dtime,yu1,yv1,ycoefm,yt,yu,ypaprs,ypplay,ydelp, &708 y_d_u_old, y_flux_u_old)709 CALL clvent(knon,dtime,yu1,yv1,ycoefm,yt,yv,ypaprs,ypplay,ydelp, &710 y_d_v_old, y_flux_v_old)711 ! save_flux est utile pour pouvoir utilise calcul_flux_vent plus tard712 CALL save_flux(klon, y_flux_u_old(:,1), y_flux_v_old(:,1))713 !jg-714 715 691 !**************************************************************************************** 716 692 ! … … 848 824 y_flux_u, y_flux_v, y_d_u, y_d_v) 849 825 850 !jg+ temporary for testing851 ! Use the results from old method852 IF (.NOT. newwind) THEN853 y_flux_u(:,:) = y_flux_u_old(:,:)854 y_flux_v(:,:) = y_flux_v_old(:,:)855 y_d_u(:,:) = y_d_u_old(:,:)856 y_d_v(:,:) = y_d_v_old(:,:)857 ENDIF858 !jg-859 860 826 DO j = 1, knon 861 827 y_dflux_t(j) = y_dflux_t(j) * ypct(j) -
LMDZ4/trunk/libf/phylmd/surface_data.F90
r781 r803 10 10 REAL, PARAMETER :: calsno=1./(2.3867e+06*.15) 11 11 12 LOGICAL, SAVE :: debug_surf=.FALSE.13 !$OMP THREADPRIVATE(debug_surf)14 LOGICAL, SAVE :: newwind=.FALSE.15 !$OMP THREADPRIVATE(newwind)16 17 12 LOGICAL, SAVE :: ok_veget ! true for use of vegetation model ORCHIDEE 18 13 !$OMP THREADPRIVATE(ok_veget)
Note: See TracChangeset
for help on using the changeset viewer.