Changeset 1144
- Timestamp:
- Apr 8, 2009, 12:47:23 PM (15 years ago)
- Location:
- LMDZ4/branches/LMDZ4-dev/libf/phylmd
- Files:
-
- 4 edited
-
pbl_surface_mod.F90 (modified) (2 diffs)
-
surf_land_mod.F90 (modified) (3 diffs)
-
surf_land_orchidee_mod.F90 (modified) (5 diffs)
-
surf_land_orchidee_noopenmp_mod.F90 (modified) (7 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/LMDZ4-dev/libf/phylmd/pbl_surface_mod.F90
r1142 r1144 765 765 r_co2_ppm(:) = co2_ppm 766 766 767 768 !**************************************************************************************** 769 ! 770 ! Calulate t2m and q2m for the case of calculation at land grid points 771 ! t2m and q2m are needed as input to ORCHIDEE 772 ! 773 !**************************************************************************************** 774 IF (nsrf == is_ter) THEN 775 776 DO i = 1, knon 777 zgeo1(i) = RD * yt(i,1) / (0.5*(ypaprs(i,1)+ypplay(i,1))) & 778 * (ypaprs(i,1)-ypplay(i,1)) 779 END DO 780 781 ! Calculate the temperature et relative humidity at 2m and the wind at 10m 782 CALL stdlevvar(klon, knon, is_ter, zxli, & 783 yu(:,1), yv(:,1), yt(:,1), yq(:,1), zgeo1, & 784 yts, yqsurf, yrugos, ypaprs(:,1), ypplay(:,1), & 785 yt2m, yq2m, yt10m, yq10m, yu10m, yustar) 786 787 END IF 788 767 789 !**************************************************************************************** 768 790 ! … … 787 809 AcoefU, AcoefV, BcoefU, BcoefV, & 788 810 ypsref, yu1, yv1, yrugoro, pctsrf, & 811 ylwdown, yq2m, yt2m, & 789 812 ysnow, yqsol, yagesno, ytsoil, & 790 813 yz0_new, yalb1_new, yalb2_new, yevap, yfluxsens, yfluxlat, & 791 814 yqsurf, ytsurf_new, y_dflux_t, y_dflux_q, & 792 y_flux_u1, y_flux_v1 , &793 ylwdown)815 y_flux_u1, y_flux_v1 ) 816 794 817 795 818 CASE(is_lic) -
LMDZ4/branches/LMDZ4-dev/libf/phylmd/surf_land_mod.F90
r1132 r1144 15 15 AcoefU, AcoefV, BcoefU, BcoefV, & 16 16 pref, u1, v1, rugoro, pctsrf, & 17 lwdown_m, q2m, t2m, & 17 18 snow, qsol, agesno, tsoil, & 18 19 z0_new, alb1_new, alb2_new, evap, fluxsens, fluxlat, & 19 20 qsurf, tsurf_new, dflux_s, dflux_l, & 20 flux_u1, flux_v1, & 21 lwdown_m) 21 flux_u1, flux_v1 ) 22 22 23 23 USE dimphy … … 58 58 REAL, DIMENSION(klon), INTENT(IN) :: rugoro 59 59 REAL, DIMENSION(klon,nbsrf), INTENT(IN) :: pctsrf 60 61 60 REAL, DIMENSION(klon), INTENT(IN) :: lwdown_m ! downwelling longwave radiation at mean surface 62 61 ! corresponds to previous sollwdown 62 REAL, DIMENSION(klon), INTENT(IN) :: q2m, t2m 63 63 64 64 ! In/Output variables … … 129 129 cdragh, AcoefH, AcoefQ, BcoefH, BcoefQ, & 130 130 precip_rain, precip_snow, lwdown_m, swnet, swdown, & 131 pref_tmp, &131 pref_tmp, q2m, t2m, & 132 132 evap, fluxsens, fluxlat, & 133 133 tsol_rad, tsurf_new, alb1_new, alb2_new, & -
LMDZ4/branches/LMDZ4-dev/libf/phylmd/surf_land_orchidee_mod.F90
r1132 r1144 36 36 tq_cdrag, petAcoef, peqAcoef, petBcoef, peqBcoef, & 37 37 precip_rain, precip_snow, lwdown, swnet, swdown, & 38 ps, &38 ps, q2m, t2m, & 39 39 evap, fluxsens, fluxlat, & 40 40 tsol_rad, tsurf_new, alb1_new, alb2_new, & … … 120 120 REAL, DIMENSION(klon), INTENT(IN) :: precip_rain, precip_snow 121 121 REAL, DIMENSION(klon), INTENT(IN) :: lwdown, swnet, swdown, ps 122 REAL, DIMENSION(klon) :: swdown_vrai122 REAL, DIMENSION(klon), INTENT(IN) :: q2m, t2m 123 123 124 124 ! Parametres de sortie … … 133 133 INTEGER :: ij, jj, igrid, ireal, index 134 134 INTEGER :: error 135 REAL, DIMENSION(klon) :: swdown_vrai 135 136 CHARACTER (len = 20) :: modname = 'surf_land_orchidee' 136 137 CHARACTER (len = 80) :: abort_message … … 391 392 evap, fluxsens, fluxlat, coastalflow, riverflow, & 392 393 tsol_rad, tsurf_new, qsurf, albedo_out, emis_new, z0_new, & 393 lon_scat, lat_scat )394 lon_scat, lat_scat, q2m, t2m) 394 395 #endif 395 396 ENDIF … … 415 416 evap(1:knon), fluxsens(1:knon), fluxlat(1:knon), coastalflow(1:knon), riverflow(1:knon), & 416 417 tsol_rad(1:knon), tsurf_new(1:knon), qsurf(1:knon), albedo_out(1:knon,:), emis_new(1:knon), z0_new(1:knon), & 417 lon_scat, lat_scat )418 lon_scat, lat_scat, q2m, t2m) 418 419 #endif 419 420 ENDIF -
LMDZ4/branches/LMDZ4-dev/libf/phylmd/surf_land_orchidee_noopenmp_mod.F90
r1138 r1144 40 40 tq_cdrag, petAcoef, peqAcoef, petBcoef, peqBcoef, & 41 41 precip_rain, precip_snow, lwdown, swnet, swdown, & 42 ps, &42 ps, q2m, t2m, & 43 43 evap, fluxsens, fluxlat, & 44 44 tsol_rad, tsurf_new, alb1_new, alb2_new, & … … 121 121 REAL, DIMENSION(klon), INTENT(IN) :: precip_rain, precip_snow 122 122 REAL, DIMENSION(klon), INTENT(IN) :: lwdown, swnet, swdown, ps 123 REAL, DIMENSION(klon) :: swdown_vrai123 REAL, DIMENSION(klon), INTENT(IN) :: q2m, t2m 124 124 125 125 ! Parametres de sortie … … 134 134 INTEGER :: ij, jj, igrid, ireal, index 135 135 INTEGER :: error 136 REAL, DIMENSION(klon) :: swdown_vrai 136 137 CHARACTER (len = 20) :: modname = 'surf_land_orchidee' 137 138 CHARACTER (len = 80) :: abort_message … … 390 391 evap, fluxsens, fluxlat, coastalflow, riverflow, & 391 392 tsol_rad, tsurf_new, qsurf, albedo_out, emis_new, z0_new, & 392 lon_scat, lat_scat )393 lon_scat, lat_scat, q2m, t2m) 393 394 394 395 #else … … 402 403 evap(1:knon), fluxsens(1:knon), fluxlat(1:knon), coastalflow(1:knon), riverflow(1:knon), & 403 404 tsol_rad(1:knon), tsurf_new(1:knon), qsurf(1:knon), albedo_out(1:knon,:), emis_new(1:knon), z0_new(1:knon), & 404 lon_scat, lat_scat )405 lon_scat, lat_scat, q2m, t2m) 405 406 #endif 406 407 … … 426 427 evap, fluxsens, fluxlat, coastalflow, riverflow, & 427 428 tsol_rad, tsurf_new, qsurf, albedo_out, emis_new, z0_new, & 428 lon_scat, lat_scat )429 lon_scat, lat_scat, q2m, t2m) 429 430 430 431 #else … … 438 439 evap(1:knon), fluxsens(1:knon), fluxlat(1:knon), coastalflow(1:knon), riverflow(1:knon), & 439 440 tsol_rad(1:knon), tsurf_new(1:knon), qsurf(1:knon), albedo_out(1:knon,:), emis_new(1:knon), z0_new(1:knon), & 440 lon_scat, lat_scat )441 lon_scat, lat_scat, q2m, t2m) 441 442 #endif 442 443
Note: See TracChangeset
for help on using the changeset viewer.
