Changeset 3823
- Timestamp:
- Feb 5, 2021, 10:00:13 AM (4 years ago)
- Location:
- LMDZ6/branches/IPSL-CM6A-MR/libf/phylmd
- Files:
-
- 4 added
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/IPSL-CM6A-MR/libf/phylmd/pbl_surface_mod.F90
r3435 r3823 23 23 USE climb_wind_mod, ONLY : climb_wind_down, climb_wind_up 24 24 USE coef_diff_turb_mod, ONLY : coef_diff_turb 25 USE ioipsl_getin_p_mod, ONLY : getin_p 26 USE cdrag_mod 27 USE stdlevvar_mod 25 28 USE wx_pbl_mod, ONLY : wx_pbl_init, wx_pbl_final, & 26 29 !! wx_pbl_fuse_no_dts, wx_pbl_split_no_dts, & … … 42 45 INTEGER, SAVE :: iflag_pbl_surface_t2m_bug 43 46 !$OMP THREADPRIVATE(iflag_pbl_surface_t2m_bug) 47 integer, save :: iflag_new_t2mq2m 48 !$OMP THREADPRIVATE(iflag_new_t2mq2m) 44 49 !FC 45 50 ! integer, save :: iflag_frein … … 177 182 alb_dir_m, alb_dif_m, zxsens, zxevap, & 178 183 alb3_lic, runoff, snowhgt, qsnow, to_ice, sissnow, & 179 zxtsol, zxfluxlat, zt2m, qsat2m, 184 zxtsol, zxfluxlat, zt2m, qsat2m, zn2mout, & 180 185 d_t, d_q, d_u, d_v, d_t_diss, & 181 186 !!! nrlmd+jyg le 02/05/2011 et le 20/02/2012 … … 394 399 REAL, DIMENSION(klon), INTENT(OUT) :: zxfluxlat ! latent flux, mean for each grid point 395 400 REAL, DIMENSION(klon), INTENT(OUT) :: zt2m ! temperature at 2m, mean for each grid point 401 INTEGER, DIMENSION(klon, 6), INTENT(OUT) :: zn2mout ! nb of times t2m is out of the [tsol,temp] 396 402 REAL, DIMENSION(klon), INTENT(OUT) :: qsat2m 397 403 REAL, DIMENSION(klon, klev), INTENT(OUT) :: d_t ! change in temperature … … 537 543 REAL, DIMENSION(klon) :: y_flux_u1, y_flux_v1 538 544 REAL, DIMENSION(klon) :: yt2m, yq2m, yu10m 545 INTEGER, DIMENSION(klon, nbsrf, 6) :: yn2mout, yn2mout_x, yn2mout_w 546 INTEGER, DIMENSION(klon, nbsrf, 6) :: n2mout, n2mout_x, n2mout_w 539 547 REAL, DIMENSION(klon) :: yustar 540 548 REAL, DIMENSION(klon) :: ywstar … … 838 846 839 847 IF (first_call) THEN 848 849 iflag_new_t2mq2m=0 850 CALL getin_p('iflag_new_t2mq2m',iflag_new_t2mq2m) 851 print*,'pbl_iflag_new_t2mq2m=',iflag_new_t2mq2m 852 840 853 print*,'PBL SURFACE AVEC GUSTINESS' 841 854 first_call=.FALSE. … … 903 916 zxfluxlat(:)=0. 904 917 zt2m(:)=0. ; zq2m(:)=0. ; qsat2m(:)=0. ; rh2m(:)=0. 918 zn2mout(:,:)=0 ; 905 919 d_t(:,:)=0. ; d_t_diss(:,:)=0. ; d_q(:,:)=0. ; d_u(:,:)=0. ; d_v(:,:)=0. 906 920 zcoefh(:,:,:)=0. ; zcoefm(:,:,:)=0. … … 1821 1835 1822 1836 ! Calculate the temperature et relative humidity at 2m and the wind at 10m 1837 IF (iflag_new_t2mq2m==1) THEN 1838 CALL stdlevvarn(klon, knon, is_ter, zxli, & 1839 yu(:,1), yv(:,1), yt(:,1), yq(:,1), zgeo1, & 1840 yts, yqsurf, yz0m, yz0h, ypaprs(:,1), ypplay(:,1), & 1841 yt2m, yq2m, yt10m, yq10m, yu10m, yn2mout(:, nsrf, :)) 1842 ELSE 1823 1843 CALL stdlevvar(klon, knon, is_ter, zxli, & 1824 1844 yu(:,1), yv(:,1), yt(:,1), yq(:,1), zgeo1, & 1825 1845 yts, yqsurf, yz0m, yz0h, ypaprs(:,1), ypplay(:,1), & 1826 1846 yt2m, yq2m, yt10m, yq10m, yu10m, yustar) 1827 1847 ENDIF 1828 1848 ENDIF 1829 1849 … … 2614 2634 !!! jyg le 07/02/2012 2615 2635 IF (iflag_split .eq.0) THEN 2636 IF (iflag_new_t2mq2m==1) THEN 2637 CALL stdlevvarn(klon, knon, nsrf, zxli, & 2638 uzon, vmer, tair1, qair1, zgeo1, & 2639 tairsol, qairsol, yz0m, yz0h_oupas, psfce, patm, & 2640 yt2m, yq2m, yt10m, yq10m, yu10m, yn2mout(:, nsrf, :)) 2641 ELSE 2616 2642 CALL stdlevvar(klon, knon, nsrf, zxli, & 2617 2643 uzon, vmer, tair1, qair1, zgeo1, & 2618 2644 tairsol, qairsol, yz0m, yz0h_oupas, psfce, patm, & 2619 2645 yt2m, yq2m, yt10m, yq10m, yu10m, yustar) 2646 ENDIF 2620 2647 ELSE !(iflag_split .eq.0) 2648 IF (iflag_new_t2mq2m==1) THEN 2649 CALL stdlevvarn(klon, knon, nsrf, zxli, & 2650 uzon_x, vmer_x, tair1_x, qair1_x, zgeo1_x, & 2651 tairsol_x, qairsol, yz0m, yz0h_oupas, psfce, patm, & 2652 yt2m_x, yq2m_x, yt10m_x, yq10m_x, yu10m_x, yn2mout_x(:, nsrf, :)) 2653 CALL stdlevvarn(klon, knon, nsrf, zxli, & 2654 uzon_w, vmer_w, tair1_w, qair1_w, zgeo1_w, & 2655 tairsol_w, qairsol, yz0m, yz0h_oupas, psfce, patm, & 2656 yt2m_w, yq2m_w, yt10m_w, yq10m_w, yu10m_w, yn2mout_w(:, nsrf, :)) 2657 ELSE 2621 2658 CALL stdlevvar(klon, knon, nsrf, zxli, & 2622 2659 uzon_x, vmer_x, tair1_x, qair1_x, zgeo1_x, & … … 2627 2664 tairsol_w, qairsol, yz0m, yz0h_oupas, psfce, patm, & 2628 2665 yt2m_w, yq2m_w, yt10m_w, yq10m_w, yu10m_w, yustar_w) 2666 ENDIF 2629 2667 !!! 2630 2668 ENDIF ! (iflag_split .eq.0) … … 2640 2678 u10m(i,nsrf)=(yu10m(j) * uzon(j))/SQRT(uzon(j)**2+vmer(j)**2) 2641 2679 v10m(i,nsrf)=(yu10m(j) * vmer(j))/SQRT(uzon(j)**2+vmer(j)**2) 2680 ! 2681 DO k = 1, 6 2682 n2mout(i,nsrf,k) = yn2mout(j,nsrf,k) 2683 END DO 2684 ! 2642 2685 ENDDO 2643 2686 ELSE !(iflag_split .eq.0) … … 2650 2693 u10m_x(i,nsrf)=(yu10m_x(j) * uzon_x(j))/SQRT(uzon_x(j)**2+vmer_x(j)**2) 2651 2694 v10m_x(i,nsrf)=(yu10m_x(j) * vmer_x(j))/SQRT(uzon_x(j)**2+vmer_x(j)**2) 2695 ! 2696 DO k = 1, 6 2697 n2mout_x(i,nsrf,k) = yn2mout_x(j,nsrf,k) 2698 END DO 2699 ! 2652 2700 ENDDO 2653 2701 DO j=1, knon … … 2663 2711 u10m(i,nsrf) = u10m_x(i,nsrf) + wake_s(i)*(u10m_w(i,nsrf)-u10m_x(i,nsrf)) 2664 2712 v10m(i,nsrf) = v10m_x(i,nsrf) + wake_s(i)*(v10m_w(i,nsrf)-v10m_x(i,nsrf)) 2713 ! 2714 DO k = 1, 6 2715 n2mout_w(i,nsrf,k) = yn2mout_w(j,nsrf,k) 2716 END DO 2717 ! 2665 2718 ENDDO 2666 2719 !!! … … 2919 2972 ! 2920 2973 zxtsol(:) = 0.0 ; zxfluxlat(:) = 0.0 2921 zt2m(:) = 0.0 ; zq2m(:) = 0.0 2974 zt2m(:) = 0.0 ; zq2m(:) = 0.0 ; zn2mout(:,:) = 0 2922 2975 zustar(:)=0.0 ; zu10m(:) = 0.0 ; zv10m(:) = 0.0 2923 2976 s_pblh(:) = 0.0 ; s_plcl(:) = 0.0 … … 2972 3025 zt2m(i) = zt2m(i) + t2m(i,nsrf) * pctsrf(i,nsrf) 2973 3026 zq2m(i) = zq2m(i) + q2m(i,nsrf) * pctsrf(i,nsrf) 3027 ! 3028 DO k = 1, 6 3029 zn2mout(i,k) = zn2mout(i,k) + n2mout(i,nsrf,k) * pctsrf(i,nsrf) 3030 ENDDO 3031 ! 2974 3032 zustar(i) = zustar(i) + ustar(i,nsrf) * pctsrf(i,nsrf) 2975 3033 wstar(i,is_ave)=wstar(i,is_ave)+wstar(i,nsrf)*pctsrf(i,nsrf) … … 3003 3061 zt2m(i) = zt2m(i) + (t2m_x(i,nsrf)+wake_s(i)*(t2m_w(i,nsrf)-t2m_x(i,nsrf))) * pctsrf(i,nsrf) 3004 3062 zq2m(i) = zq2m(i) + q2m_x(i,nsrf) * pctsrf(i,nsrf) 3063 ! 3064 DO k = 1, 6 3065 zn2mout(i,k) = zn2mout(i,k) + n2mout_x(i,nsrf,k) * pctsrf(i,nsrf) 3066 ENDDO 3067 ! 3005 3068 zustar(i) = zustar(i) + ustar_x(i,nsrf) * pctsrf(i,nsrf) 3006 3069 wstar(i,is_ave)=wstar(i,is_ave)+wstar_x(i,nsrf)*pctsrf(i,nsrf) … … 3081 3144 DO nsrf = 1, nbsrf 3082 3145 DO i = 1, klon 3083 zxqsurf(i) = zxqsurf(i) + qsurf(i,nsrf) * pctsrf(i,nsrf)3146 zxqsurf(i) = zxqsurf(i) + MAX(qsurf(i,nsrf),0.0) * pctsrf(i,nsrf) 3084 3147 zxsnow(i) = zxsnow(i) + snow(i,nsrf) * pctsrf(i,nsrf) 3085 3148 ENDDO -
LMDZ6/branches/IPSL-CM6A-MR/libf/phylmd/phys_local_var_mod.F90
r3662 r3823 380 380 REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: t2m_min_mon, t2m_max_mon 381 381 !$OMP THREADPRIVATE(t2m_min_mon, t2m_max_mon) 382 REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zq2m_cor, zt2m_cor383 !$OMP THREADPRIVATE(zq2m_cor, zt2m_cor)384 REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zu10m_cor, zv10m_cor385 !$OMP THREADPRIVATE(zu10m_cor, zv10m_cor)386 REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zrh2m_cor, zqsat2m_cor387 !$OMP THREADPRIVATE(zrh2m_cor, zqsat2m_cor)388 382 REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: weak_inversion 389 383 !$OMP THREADPRIVATE(weak_inversion) … … 475 469 REAL, ALLOCATABLE, SAVE, DIMENSION(:) :: p_tropopause, z_tropopause, t_tropopause 476 470 !$OMP THREADPRIVATE(p_tropopause, z_tropopause, t_tropopause) 471 472 INTEGER,ALLOCATABLE,SAVE,DIMENSION(:,:) :: zn2mout 473 !$OMP THREADPRIVATE(zn2mout) 477 474 478 475 #ifdef CPP_StratAer … … 760 757 ALLOCATE(zt2m_min_mon(klon), zt2m_max_mon(klon)) 761 758 ALLOCATE(t2m_min_mon(klon), t2m_max_mon(klon)) 762 ALLOCATE(zq2m_cor(klon), zt2m_cor(klon), zu10m_cor(klon), zv10m_cor(klon))763 ALLOCATE(zrh2m_cor(klon), zqsat2m_cor(klon))764 759 ALLOCATE(sens(klon), flwp(klon), fiwp(klon)) 765 760 ALLOCATE(alp_bl_conv(klon), alp_bl_det(klon)) … … 831 826 ALLOCATE (z_tropopause(klon)) 832 827 ALLOCATE (t_tropopause(klon)) 828 829 ALLOCATE(zn2mout(klon,6)) 833 830 834 831 #ifdef CPP_StratAer … … 1057 1054 DEALLOCATE(zt2m_min_mon, zt2m_max_mon) 1058 1055 DEALLOCATE(t2m_min_mon, t2m_max_mon) 1059 DEALLOCATE(zq2m_cor, zt2m_cor, zu10m_cor, zv10m_cor)1060 DEALLOCATE(zrh2m_cor, zqsat2m_cor)1061 1056 DEALLOCATE(sens, flwp, fiwp) 1062 1057 DEALLOCATE(alp_bl_conv,alp_bl_det) … … 1119 1114 DEALLOCATE (z_tropopause) 1120 1115 DEALLOCATE (t_tropopause) 1116 1117 DEALLOCATE(zn2mout) 1121 1118 1122 1119 #ifdef CPP_StratAer -
LMDZ6/branches/IPSL-CM6A-MR/libf/phylmd/phys_output_ctrlout_mod.F90
r3702 r3823 271 271 ctrl_out((/ 10, 6, 10, 10, 10, 10, 11, 11, 11, 11/), & 272 272 't2m_sic', "Temp 2m "//clnsurf(4), "K", (/ ('', i=1, 10) /)) /) 273 274 TYPE(ctrl_out), SAVE :: o_nt2mout = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11, 11/), & 275 'nt2mout', 'Nbt2m out of range complete computation', '-', (/ ('', i=1, 10) /)) 276 TYPE(ctrl_out), SAVE :: o_nq2mout = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11, 11/), & 277 'nq2mout', 'Nbq2m out of range complete computation', '-', (/ ('', i=1, 10) /)) 278 TYPE(ctrl_out), SAVE :: o_nu2mout = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11, 11/), & 279 'nu2mout', 'Nbu2m out of range complete computation', '-', (/ ('', i=1, 10) /)) 280 281 TYPE(ctrl_out), SAVE :: o_nt2moutfg = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11, 11/), & 282 'nt2moutfg', 'Nbt2m out of range complete/fgRi1 computation', '-', (/ ('', i=1, 10) /)) 283 TYPE(ctrl_out), SAVE :: o_nq2moutfg = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11, 11/), & 284 'nq2moutfg', 'Nbq2m out of range complete/fgRi1 computation', '-', (/ ('', i=1, 10) /)) 285 TYPE(ctrl_out), SAVE :: o_nu2moutfg = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11, 11/), & 286 'nu2moutfg', 'Nbu2m out of range complete/fgRi1 computation', '-', (/ ('', i=1, 10) /)) 273 287 274 288 TYPE(ctrl_out), SAVE :: o_gusts = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11, 11/), & -
LMDZ6/branches/IPSL-CM6A-MR/libf/phylmd/phys_output_write_mod.F90
r3702 r3823 38 38 o_t2m, o_t2m_min, o_t2m_max, & 39 39 o_t2m_min_mon, o_t2m_max_mon, & 40 o_nt2mout, o_nt2moutfg, & 41 o_nq2mout, o_nq2moutfg, & 42 o_nu2mout, o_nu2moutfg, & 40 43 o_q2m, o_ustar, o_u10m, o_v10m, & 41 44 o_wind10m, o_wind10max, o_wind100m, o_gusts, o_sicf, & … … 249 252 250 253 USE phys_local_var_mod, ONLY: zxfluxlat, slp, ptstar, pt0, zxtsol, zt2m, & 251 zt2m_cor,zq2m_cor,zu10m_cor,zv10m_cor, zrh2m_cor, zqsat2m_cor, & 252 t2m_min_mon, t2m_max_mon, evap, & 254 zn2mout, t2m_min_mon, t2m_max_mon, evap, & 253 255 l_mixmin,l_mix, & 254 256 zu10m, zv10m, zq2m, zustar, zxqsurf, & … … 676 678 CALL histwrite_phy(o_slp, slp) 677 679 CALL histwrite_phy(o_tsol, zxtsol) 678 CALL histwrite_phy(o_t2m, zt2m _cor)679 CALL histwrite_phy(o_t2m_min, zt2m _cor)680 CALL histwrite_phy(o_t2m_max, zt2m _cor)680 CALL histwrite_phy(o_t2m, zt2m) 681 CALL histwrite_phy(o_t2m_min, zt2m) 682 CALL histwrite_phy(o_t2m_max, zt2m) 681 683 CALL histwrite_phy(o_t2m_max_mon, t2m_max_mon) 682 684 CALL histwrite_phy(o_t2m_min_mon, t2m_min_mon) … … 684 686 IF (vars_defined) THEN 685 687 DO i=1, klon 686 zx_tmp_fi2d(i)=SQRT(zu10m_cor(i)*zu10m_cor(i)+zv10m_cor(i)*zv10m_cor(i)) 688 zx_tmp_fi2d(i)=real(zn2mout(i,1)) 689 ENDDO 690 ENDIF 691 CALL histwrite_phy(o_nt2mout, zx_tmp_fi2d) 692 693 IF (vars_defined) THEN 694 DO i=1, klon 695 zx_tmp_fi2d(i)=real(zn2mout(i,2)) 696 ENDDO 697 ENDIF 698 CALL histwrite_phy(o_nt2moutfg, zx_tmp_fi2d) 699 700 IF (vars_defined) THEN 701 DO i=1, klon 702 zx_tmp_fi2d(i)=real(zn2mout(i,3)) 703 ENDDO 704 ENDIF 705 CALL histwrite_phy(o_nq2mout, zx_tmp_fi2d) 706 707 IF (vars_defined) THEN 708 DO i=1, klon 709 zx_tmp_fi2d(i)=real(zn2mout(i,4)) 710 ENDDO 711 ENDIF 712 CALL histwrite_phy(o_nq2moutfg, zx_tmp_fi2d) 713 714 IF (vars_defined) THEN 715 DO i=1, klon 716 zx_tmp_fi2d(i)=real(zn2mout(i,5)) 717 ENDDO 718 ENDIF 719 CALL histwrite_phy(o_nu2mout, zx_tmp_fi2d) 720 721 IF (vars_defined) THEN 722 DO i=1, klon 723 zx_tmp_fi2d(i)=real(zn2mout(i,6)) 724 ENDDO 725 ENDIF 726 CALL histwrite_phy(o_nu2moutfg, zx_tmp_fi2d) 727 728 IF (vars_defined) THEN 729 DO i=1, klon 730 zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i)) 687 731 ENDDO 688 732 ENDIF … … 691 735 IF (vars_defined) THEN 692 736 DO i=1, klon 693 zx_tmp_fi2d(i)=SQRT(zu10m _cor(i)*zu10m_cor(i)+zv10m_cor(i)*zv10m_cor(i))737 zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i)) 694 738 ENDDO 695 739 ENDIF … … 771 815 ENDIF 772 816 CALL histwrite_phy(o_sicf, zx_tmp_fi2d) 773 CALL histwrite_phy(o_q2m, zq2m_cor) 774 CALL histwrite_phy(o_ustar, zustar) 775 CALL histwrite_phy(o_u10m, zu10m_cor) 776 CALL histwrite_phy(o_v10m, zv10m_cor) 817 CALL histwrite_phy(o_q2m, zq2m) 818 IF (vars_defined) zx_tmp_fi2d = zustar 819 CALL histwrite_phy(o_ustar, zx_tmp_fi2d) 820 CALL histwrite_phy(o_u10m, zu10m) 821 CALL histwrite_phy(o_v10m, zv10m) 777 822 778 823 IF (vars_defined) THEN … … 1412 1457 ! CALL histwrite_phy(o_rh2m_max, zx_tmp_fi2d) 1413 1458 1414 CALL histwrite_phy(o_qsat2m, zqsat2m_cor)1459 CALL histwrite_phy(o_qsat2m, qsat2m) 1415 1460 CALL histwrite_phy(o_tpot, tpot) 1416 1461 CALL histwrite_phy(o_tpote, tpote) -
LMDZ6/branches/IPSL-CM6A-MR/libf/phylmd/physiq_mod.F90
r3666 r3823 119 119 cdragm, cdragh, & 120 120 zustar, zu10m, zv10m, rh2m, qsat2m, & 121 zq2m, zt2m, weak_inversion, & 122 zq2m_cor,zt2m_cor,zu10m_cor,zv10m_cor, & ! pour corriger d'un bug 123 zrh2m_cor,zqsat2m_cor, & 121 zq2m, zt2m, zn2mout, weak_inversion, & 124 122 zt2m_min_mon, zt2m_max_mon, & ! pour calcul_divers.h 125 123 t2m_min_mon, t2m_max_mon, & ! pour calcul_divers.h … … 603 601 !$OMP THREADPRIVATE(iflag_alp_wk_cond) 604 602 605 INTEGER, SAVE :: iflag_bug_t2m_ipslcm61=1 !606 !$OMP THREADPRIVATE(iflag_bug_t2m_ipslcm61)607 INTEGER, SAVE :: iflag_bug_t2m_stab_ipslcm61=-1 !608 !$OMP THREADPRIVATE(iflag_bug_t2m_stab_ipslcm61)609 610 603 REAL t_w(klon,klev),q_w(klon,klev) ! temperature and moisture profiles in the wake region 611 604 REAL t_x(klon,klev),q_x(klon,klev) ! temperature and moisture profiles in the off-wake region … … 1314 1307 tau_gl=86400.*tau_gl 1315 1308 WRITE(lunout,*) 'debut physiq_mod tau_gl=',tau_gl 1316 1317 iflag_bug_t2m_ipslcm61 = 11318 CALL getin_p('iflag_bug_t2m_ipslcm61', iflag_bug_t2m_ipslcm61)1319 iflag_bug_t2m_stab_ipslcm61 = -11320 CALL getin_p('iflag_bug_t2m_stab_ipslcm61', iflag_bug_t2m_stab_ipslcm61)1321 1309 1322 1310 CALL getin_p('iflag_alp_wk_cond', iflag_alp_wk_cond) … … 2529 2517 !albedo SB <<< 2530 2518 albsol3_lic,runoff, snowhgt, qsnow, to_ice, sissnow, & 2531 zxtsol, zxfluxlat, zt2m, qsat2m, &2519 zxtsol, zxfluxlat, zt2m, qsat2m, zn2mout, & 2532 2520 d_t_vdf, d_q_vdf, d_u_vdf, d_v_vdf, d_t_diss, & 2533 2521 !nrlmd< … … 2576 2564 !>jyg 2577 2565 ENDIF 2578 2579 !add limitation for t,q at and wind at 10m2580 if ( iflag_bug_t2m_ipslcm61 == 0 ) THEN2581 CALL borne_var_surf( klon,klev,nbsrf, &2582 iflag_bug_t2m_stab_ipslcm61, &2583 t_seri(:,1),q_seri(:,1),u_seri(:,1),v_seri(:,1), &2584 ftsol,zxqsurf,pctsrf,paprs, &2585 t2m, q2m, u10m, v10m, &2586 zt2m_cor, zq2m_cor, zu10m_cor, zv10m_cor, &2587 zrh2m_cor, zqsat2m_cor)2588 ELSE2589 zt2m_cor(:)=zt2m(:)2590 zq2m_cor(:)=zq2m(:)2591 zu10m_cor(:)=zu10m(:)2592 zv10m_cor(:)=zv10m(:)2593 zqsat2m_cor=999.9992594 ENDIF2595 2566 2596 2567 !---------------------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.