Ignore:
Timestamp:
Apr 8, 2009, 12:47:23 PM (16 years ago)
Author:
jghattas
Message:

Ajoute des variables q2m et t2m comme argument d'entrees dans ORCHIDEE. Pour cela, ajoute de la calcul de ces variables avec stdlevvar avant calcul sur sous-surface.

Location:
LMDZ4/branches/LMDZ4-dev/libf/phylmd
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/branches/LMDZ4-dev/libf/phylmd/pbl_surface_mod.F90

    r1142 r1144  
    765765       r_co2_ppm(:) = co2_ppm
    766766
     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
    767789!****************************************************************************************
    768790!
    … …  
    787809               AcoefU, AcoefV, BcoefU, BcoefV, &
    788810               ypsref, yu1, yv1, yrugoro, pctsrf, &
     811               ylwdown, yq2m, yt2m, &
    789812               ysnow, yqsol, yagesno, ytsoil, &
    790813               yz0_new, yalb1_new, yalb2_new, yevap, yfluxsens, yfluxlat, &
    791814               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               
    794817     
    795818       CASE(is_lic)

  • LMDZ4/branches/LMDZ4-dev/libf/phylmd/surf_land_mod.F90

    r1132 r1144  
    1515       AcoefU, AcoefV, BcoefU, BcoefV, &
    1616       pref, u1, v1, rugoro, pctsrf, &
     17       lwdown_m, q2m, t2m, &
    1718       snow, qsol, agesno, tsoil, &
    1819       z0_new, alb1_new, alb2_new, evap, fluxsens, fluxlat, &
    1920       qsurf, tsurf_new, dflux_s, dflux_l, &
    20        flux_u1, flux_v1, &
    21        lwdown_m)
     21       flux_u1, flux_v1 )
    2222
    2323    USE dimphy
    … …  
    5858    REAL, DIMENSION(klon), INTENT(IN)       :: rugoro
    5959    REAL, DIMENSION(klon,nbsrf), INTENT(IN) :: pctsrf
    60 
    6160    REAL, DIMENSION(klon), INTENT(IN)       :: lwdown_m  ! downwelling longwave radiation at mean surface
    6261                                                         ! corresponds to previous sollwdown
     62    REAL, DIMENSION(klon), INTENT(IN)       :: q2m, t2m
    6363
    6464! In/Output variables
    … …  
    129129            cdragh, AcoefH, AcoefQ, BcoefH, BcoefQ, &
    130130            precip_rain, precip_snow, lwdown_m, swnet, swdown, &
    131             pref_tmp, &
     131            pref_tmp, q2m, t2m, &
    132132            evap, fluxsens, fluxlat, &             
    133133            tsol_rad, tsurf_new, alb1_new, alb2_new, &

  • LMDZ4/branches/LMDZ4-dev/libf/phylmd/surf_land_orchidee_mod.F90

    r1132 r1144  
    3636       tq_cdrag, petAcoef, peqAcoef, petBcoef, peqBcoef, &
    3737       precip_rain, precip_snow, lwdown, swnet, swdown, &
    38        ps, &
     38       ps, q2m, t2m, &
    3939       evap, fluxsens, fluxlat, &             
    4040       tsol_rad, tsurf_new, alb1_new, alb2_new, &
    … …  
    120120    REAL, DIMENSION(klon), INTENT(IN)         :: precip_rain, precip_snow
    121121    REAL, DIMENSION(klon), INTENT(IN)         :: lwdown, swnet, swdown, ps
    122     REAL, DIMENSION(klon)                     :: swdown_vrai
     122    REAL, DIMENSION(klon), INTENT(IN)         :: q2m, t2m
    123123
    124124! Parametres de sortie
    … …  
    133133    INTEGER                                   :: ij, jj, igrid, ireal, index
    134134    INTEGER                                   :: error
     135    REAL, DIMENSION(klon)                     :: swdown_vrai
    135136    CHARACTER (len = 20)                      :: modname = 'surf_land_orchidee'
    136137    CHARACTER (len = 80)                      :: abort_message
    … …  
    391392               evap, fluxsens, fluxlat, coastalflow, riverflow, &
    392393               tsol_rad, tsurf_new, qsurf, albedo_out, emis_new, z0_new, &
    393                lon_scat, lat_scat)
     394               lon_scat, lat_scat, q2m, t2m)
    394395#endif         
    395396       ENDIF
    … …  
    415416            evap(1:knon), fluxsens(1:knon), fluxlat(1:knon), coastalflow(1:knon), riverflow(1:knon), &
    416417            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)
    418419#endif       
    419420    ENDIF

  • LMDZ4/branches/LMDZ4-dev/libf/phylmd/surf_land_orchidee_noopenmp_mod.F90

    r1138 r1144  
    4040       tq_cdrag, petAcoef, peqAcoef, petBcoef, peqBcoef, &
    4141       precip_rain, precip_snow, lwdown, swnet, swdown, &
    42        ps, &
     42       ps, q2m, t2m, &
    4343       evap, fluxsens, fluxlat, &             
    4444       tsol_rad, tsurf_new, alb1_new, alb2_new, &
    … …  
    121121    REAL, DIMENSION(klon), INTENT(IN)         :: precip_rain, precip_snow
    122122    REAL, DIMENSION(klon), INTENT(IN)         :: lwdown, swnet, swdown, ps
    123     REAL, DIMENSION(klon)                     :: swdown_vrai
     123    REAL, DIMENSION(klon), INTENT(IN)         :: q2m, t2m
    124124
    125125! Parametres de sortie
    … …  
    134134    INTEGER                                   :: ij, jj, igrid, ireal, index
    135135    INTEGER                                   :: error
     136    REAL, DIMENSION(klon)                     :: swdown_vrai
    136137    CHARACTER (len = 20)                      :: modname = 'surf_land_orchidee'
    137138    CHARACTER (len = 80)                      :: abort_message
    … …  
    390391               evap, fluxsens, fluxlat, coastalflow, riverflow, &
    391392               tsol_rad, tsurf_new, qsurf, albedo_out, emis_new, z0_new, &
    392                lon_scat, lat_scat)
     393               lon_scat, lat_scat, q2m, t2m)
    393394
    394395#else         
    … …  
    402403               evap(1:knon), fluxsens(1:knon), fluxlat(1:knon), coastalflow(1:knon), riverflow(1:knon), &
    403404               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)
    405406#endif
    406407         
    … …  
    426427            evap, fluxsens, fluxlat, coastalflow, riverflow, &
    427428            tsol_rad, tsurf_new, qsurf, albedo_out, emis_new, z0_new, &
    428             lon_scat, lat_scat)
     429            lon_scat, lat_scat, q2m, t2m)
    429430       
    430431#else
    … …  
    438439            evap(1:knon), fluxsens(1:knon), fluxlat(1:knon), coastalflow(1:knon), riverflow(1:knon), &
    439440            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)
    441442#endif
    442443       
Note: See TracChangeset for help on using the changeset viewer.