Changeset 4569


Ignore:
Timestamp:
Jun 13, 2023, 5:14:32 PM (18 months ago)
Author:
jyg
Message:

replacing (:) by (1:knon) in operations on compressed variables

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90

    r4545 r4569  
    13881388!   Provisional : set ybeta to standard values
    13891389       IF (nsrf .NE. is_ter) THEN
    1390            ybeta(:) = 1.
     1390           ybeta(1:knon) = 1.
    13911391       ELSE
    13921392           IF (iflag_split .EQ. 0) THEN
    1393               ybeta(:) = 1.
     1393              ybeta(1:knon) = 1.
    13941394           ELSE
    13951395             DO j = 1, knon
    … …  
    16131613     ENDIF
    16141614
    1615         IF (prt_level >=10) print *,'clcdrag -> ycdragh ', ycdragh
     1615        IF (prt_level >=10) print *,'cdrag -> ycdragh ', ycdragh(1:knon)
    16161616       ELSE  !(iflag_split .eq.0)
    16171617
    … …  
    16611661            ycdragm_w, ycdragh_w, zri1_w, pref_w )
    16621662!
    1663         zgeo1(:) = wake_s(:)*zgeo1_w(:) + (1.-wake_s(:))*zgeo1_x(:)
     1663!!!bug !!        zgeo1(:) = wake_s(:)*zgeo1_w(:) + (1.-wake_s(:))*zgeo1_x(:)
     1664        zgeo1(1:knon) = wake_s(1:knon)*zgeo1_w(1:knon) + (1.-wake_s(1:knon))*zgeo1_x(1:knon)
    16641665
    16651666! --- special Dice. JYG+MPL 25112013 puis BOMEX
    … …  
    16871688!!! nrlmd & jyg les 02/05/2011, 13/06/2011, 05/02/2012
    16881689      IF (prt_level >=10) THEN
    1689       print *,' args coef_diff_turb: yu ',  yu 
    1690       print *,' args coef_diff_turb: yv ',  yv 
    1691       print *,' args coef_diff_turb: yq ',  yq 
    1692       print *,' args coef_diff_turb: yt ',  yt 
    1693       print *,' args coef_diff_turb: yts ', yts 
    1694       print *,' args coef_diff_turb: yz0m ', yz0m 
    1695       print *,' args coef_diff_turb: yqsurf ', yqsurf 
    1696       print *,' args coef_diff_turb: ycdragm ', ycdragm
    1697       print *,' args coef_diff_turb: ycdragh ', ycdragh
    1698       print *,' args coef_diff_turb: ytke ', ytke
     1690      print *,' args coef_diff_turb: yu ',  yu(1:knon,:) 
     1691      print *,' args coef_diff_turb: yv ',  yv(1:knon,:)   
     1692      print *,' args coef_diff_turb: yq ',  yq(1:knon,:)   
     1693      print *,' args coef_diff_turb: yt ',  yt(1:knon,:)   
     1694      print *,' args coef_diff_turb: yts ', yts(1:knon)
     1695      print *,' args coef_diff_turb: yz0m ', yz0m(1:knon)
     1696      print *,' args coef_diff_turb: yqsurf ', yqsurf(1:knon) 
     1697      print *,' args coef_diff_turb: ycdragm ', ycdragm(1:knon)
     1698      print *,' args coef_diff_turb: ycdragh ', ycdragh(1:knon)
     1699      print *,' args coef_diff_turb: ytke ', ytke(1:knon,:)   
    16991700       ENDIF
    17001701
    … …  
    17241725       ENDIF ! iflag_pbl >= 50
    17251726
    1726         IF (prt_level >=10) print *,'coef_diff_turb -> ycoefh ',ycoefh
     1727        IF (prt_level >=10) print *,'coef_diff_turb -> ycoefh ',ycoefh(1:knon,:)
    17271728
    17281729
    … …  
    17311732     
    17321733      IF (prt_level >=10) THEN
    1733       print *,' args coef_diff_turb: yu_x ',  yu_x 
    1734       print *,' args coef_diff_turb: yv_x ',  yv_x 
    1735       print *,' args coef_diff_turb: yq_x ',  yq_x 
    1736       print *,' args coef_diff_turb: yt_x ',  yt_x 
    1737       print *,' args coef_diff_turb: yts_x ', yts_x 
    1738       print *,' args coef_diff_turb: yqsurf ', yqsurf 
    1739       print *,' args coef_diff_turb: ycdragm_x ', ycdragm_x
    1740       print *,' args coef_diff_turb: ycdragh_x ', ycdragh_x
    1741       print *,' args coef_diff_turb: ytke_x ', ytke_x
     1734      print *,' args coef_diff_turb: yu_x ',  yu_x(1:knon,:)     
     1735      print *,' args coef_diff_turb: yv_x ',  yv_x(1:knon,:)     
     1736      print *,' args coef_diff_turb: yq_x ',  yq_x(1:knon,:)     
     1737      print *,' args coef_diff_turb: yt_x ',  yt_x(1:knon,:)     
     1738      print *,' args coef_diff_turb: yts_x ', yts_x(1:knon)
     1739      print *,' args coef_diff_turb: yqsurf ', yqsurf(1:knon) 
     1740      print *,' args coef_diff_turb: ycdragm_x ', ycdragm_x(1:knon)
     1741      print *,' args coef_diff_turb: ycdragh_x ', ycdragh_x(1:knon)
     1742      print *,' args coef_diff_turb: ytke_x ', ytke_x(1:knon,:)   
    17421743      ENDIF
    17431744
    … …  
    17681769        ENDIF ! iflag_pbl >= 50
    17691770
    1770         IF (prt_level >=10) print *,'coef_diff_turb -> ycoefh_x ',ycoefh_x
     1771        IF (prt_level >=10) print *,'coef_diff_turb -> ycoefh_x ',ycoefh_x(1:knon,:)
    17711772!
    17721773      IF (prt_level >=10) THEN
    1773       print *,' args coef_diff_turb: yu_w ',  yu_w 
    1774       print *,' args coef_diff_turb: yv_w ',  yv_w 
    1775       print *,' args coef_diff_turb: yq_w ',  yq_w 
    1776       print *,' args coef_diff_turb: yt_w ',  yt_w 
    1777       print *,' args coef_diff_turb: yts_w ', yts_w 
    1778       print *,' args coef_diff_turb: yqsurf ', yqsurf 
    1779       print *,' args coef_diff_turb: ycdragm_w ', ycdragm_w
    1780       print *,' args coef_diff_turb: ycdragh_w ', ycdragh_w
    1781       print *,' args coef_diff_turb: ytke_w ', ytke_w
     1774      print *,' args coef_diff_turb: yu_w ',  yu_w(1:knon,:)
     1775      print *,' args coef_diff_turb: yv_w ',  yv_w(1:knon,:) 
     1776      print *,' args coef_diff_turb: yq_w ',  yq_w(1:knon,:) 
     1777      print *,' args coef_diff_turb: yt_w ',  yt_w(1:knon,:) 
     1778      print *,' args coef_diff_turb: yts_w ', yts_w(1:knon)
     1779      print *,' args coef_diff_turb: yqsurf ', yqsurf(1:knon) 
     1780      print *,' args coef_diff_turb: ycdragm_w ', ycdragm_w(1:knon)
     1781      print *,' args coef_diff_turb: ycdragh_w ', ycdragh_w(1:knon)
     1782      print *,' args coef_diff_turb: ytke_w ', ytke_w(1:knon,:)
    17821783      ENDIF
    17831784     
    … …  
    18071808
    18081809
    1809         IF (prt_level >=10) print *,'coef_diff_turb -> ycoefh_w ',ycoefh_w
     1810        IF (prt_level >=10) print *,'coef_diff_turb -> ycoefh_w ',ycoefh_w(1:knon,:)
    18101811
    18111812!!!jyg le 10/04/2013
    … …  
    20432044!!!
    20442045       IF (prt_level >=10) THEN
    2045          PRINT *,'pbl_surface (merge->): yt(1,:) ',yt(1,:)
    2046          PRINT *,'pbl_surface (merge->): yq(1,:) ',yq(1,:)
    2047          PRINT *,'pbl_surface (merge->): yu(1,:) ',yu(1,:)
    2048          PRINT *,'pbl_surface (merge->): yv(1,:) ',yv(1,:)
    2049          PRINT *,'pbl_surface (merge->): AcoefH(1), AcoefQ(1), AcoefU(1), AcoefV(1) ', &
    2050                                          AcoefH(1), AcoefQ(1), AcoefU(1), AcoefV(1)
    2051          PRINT *,'pbl_surface (merge->): BcoefH(1), BcoefQ(1), BcoefU(1), BcoefV(1) ', &
    2052                                          BcoefH(1), BcoefQ(1), BcoefU(1), BcoefV(1)
     2046         DO i = 1, min(1,knon)
     2047           PRINT *,'pbl_surface (merge->): yt(1,:) ',yt(i,:)
     2048           PRINT *,'pbl_surface (merge->): yq(1,:) ',yq(i,:)
     2049           PRINT *,'pbl_surface (merge->): yu(1,:) ',yu(i,:)
     2050           PRINT *,'pbl_surface (merge->): yv(1,:) ',yv(i,:)
     2051           PRINT *,'pbl_surface (merge->): AcoefH(1), AcoefQ(1), AcoefU(1), AcoefV(1) ', &
     2052                                           AcoefH(i), AcoefQ(i), AcoefU(i), AcoefV(i)
     2053           PRINT *,'pbl_surface (merge->): BcoefH(1), BcoefQ(1), BcoefU(1), BcoefV(1) ', &
     2054                                           BcoefH(i), BcoefQ(i), BcoefU(i), BcoefV(i)
     2055         ENDDO
    20532056
    20542057       ENDIF
    … …  
    22022205               ydt_ds(:knon), ytkt(:knon), ytks(:knon), ytaur(:knon), ysss)
    22032206      IF (prt_level >=10) THEN
    2204           print *,'arg de surf_ocean: ycdragh ',ycdragh
    2205           print *,'arg de surf_ocean: ycdragm ',ycdragm
    2206           print *,'arg de surf_ocean: yt ', yt
    2207           print *,'arg de surf_ocean: yq ', yq
    2208           print *,'arg de surf_ocean: yts ', yts
    2209           print *,'arg de surf_ocean: AcoefH ',AcoefH
    2210           print *,'arg de surf_ocean: AcoefQ ',AcoefQ
    2211           print *,'arg de surf_ocean: BcoefH ',BcoefH
    2212           print *,'arg de surf_ocean: BcoefQ ',BcoefQ
    2213           print *,'arg de surf_ocean: yevap ',yevap
    2214           print *,'arg de surf_ocean: yfluxsens ',yfluxsens
    2215           print *,'arg de surf_ocean: yfluxlat ',yfluxlat
    2216           print *,'arg de surf_ocean: ytsurf_new ',ytsurf_new
     2207          print *,'arg de surf_ocean: ycdragh ',ycdragh(1:knon)
     2208          print *,'arg de surf_ocean: ycdragm ',ycdragm(1:knon)
     2209          print *,'arg de surf_ocean: yt ', yt(1:knon,:)
     2210          print *,'arg de surf_ocean: yq ', yq(1:knon,:)
     2211          print *,'arg de surf_ocean: yts ', yts(1:knon)
     2212          print *,'arg de surf_ocean: AcoefH ',AcoefH(1:knon)
     2213          print *,'arg de surf_ocean: AcoefQ ',AcoefQ(1:knon)
     2214          print *,'arg de surf_ocean: BcoefH ',BcoefH(1:knon)
     2215          print *,'arg de surf_ocean: BcoefQ ',BcoefQ(1:knon)
     2216          print *,'arg de surf_ocean: yevap ',yevap(1:knon)
     2217          print *,'arg de surf_ocean: yfluxsens ',yfluxsens(1:knon)
     2218          print *,'arg de surf_ocean: yfluxlat ',yfluxlat(1:knon)
     2219          print *,'arg de surf_ocean: ytsurf_new ',ytsurf_new(1:knon)
    22172220       ENDIF
    22182221! Special DICE MPL 05082013 puis BOMEX MPL 20150410
    … …  
    22642267
    22652268       IF (evap0>=0.) THEN
    2266           yevap(:)=evap0
    2267           yevap(:)=RLVTT*evap0
     2269          yevap(1:knon)=evap0
     2270          yevap(1:knon)=RLVTT*evap0
    22682271       ENDIF
    22692272
    … …  
    24722475       IF (prt_level >= 10) THEN
    24732476         print *,'pbl_surface, ybeta , yevap, yevap_pot ', &
    2474                                ybeta , yevap, yevap_pot
     2477                               ybeta(1:knon) , yevap(1:knon), yevap_pot(1:knon)
    24752478       ENDIF  ! (prt_level >= 10)
    24762479!
    … …  
    29572960       IF (prt_level >=10) THEN
    29582961         print *, 'pbl_surface tendencies for w: d_t_w, d_t_x, d_t ', &
    2959           d_t_w(:,1), d_t_x(:,1), d_t(:,1)
     2962          d_t_w(1:knon,1), d_t_x(1:knon,1), d_t(1:knon,1)
    29602963       ENDIF
    29612964
    … …  
    32033206       IF (prt_level >=10) THEN
    32043207         print *, 'T2m, q2m, RH2m ', &
    3205           t2m, q2m, rh2m
     3208          t2m(1:knon,:), q2m(1:knon,:), rh2m(1:knon)
    32063209       ENDIF
    32073210
    … …  
    32163219            ytherm,ytrmb1,ytrmb2,ytrmb3,ylcl)
    32173220          IF (prt_level >=10) THEN
    3218        print *,' Arg. de HBTM: yt2m ',yt2m
    3219        print *,' Arg. de HBTM: yt10m ',yt10m
    3220        print *,' Arg. de HBTM: yq2m ',yq2m
    3221        print *,' Arg. de HBTM: yq10m ',yq10m
    3222        print *,' Arg. de HBTM: yustar ',yustar
    3223        print *,' Arg. de HBTM: y_flux_t ',y_flux_t
    3224        print *,' Arg. de HBTM: y_flux_q ',y_flux_q
    3225        print *,' Arg. de HBTM: yu ',yu
    3226        print *,' Arg. de HBTM: yv ',yv
    3227        print *,' Arg. de HBTM: yt ',yt
    3228        print *,' Arg. de HBTM: yq ',yq
     3221       print *,' Arg. de HBTM: yt2m ',yt2m(1:knon)
     3222       print *,' Arg. de HBTM: yt10m ',yt10m(1:knon)
     3223       print *,' Arg. de HBTM: yq2m ',yq2m(1:knon)
     3224       print *,' Arg. de HBTM: yq10m ',yq10m(1:knon)
     3225       print *,' Arg. de HBTM: yustar ',yustar(1:knon)
     3226       print *,' Arg. de HBTM: y_flux_t ',y_flux_t(1:knon,:)
     3227       print *,' Arg. de HBTM: y_flux_q ',y_flux_q(1:knon,:)
     3228       print *,' Arg. de HBTM: yu ',yu(1:knon,:)
     3229       print *,' Arg. de HBTM: yv ',yv(1:knon,:)
     3230       print *,' Arg. de HBTM: yt ',yt(1:knon,:)
     3231       print *,' Arg. de HBTM: yq ',yq(1:knon,:)
    32293232          ENDIF
    32303233       ELSE  ! (iflag_split .eq.0)
    … …  
    32353238            ytherm_x,ytrmb1_x,ytrmb2_x,ytrmb3_x,ylcl_x)
    32363239          IF (prt_level >=10) THEN
    3237        print *,' Arg. de HBTM: yt2m_x ',yt2m_x
    3238        print *,' Arg. de HBTM: yt10m_x ',yt10m_x
    3239        print *,' Arg. de HBTM: yq2m_x ',yq2m_x
    3240        print *,' Arg. de HBTM: yq10m_x ',yq10m_x
    3241        print *,' Arg. de HBTM: yustar_x ',yustar_x
    3242        print *,' Arg. de HBTM: y_flux_t_x ',y_flux_t_x
    3243        print *,' Arg. de HBTM: y_flux_q_x ',y_flux_q_x
    3244        print *,' Arg. de HBTM: yu_x ',yu_x
    3245        print *,' Arg. de HBTM: yv_x ',yv_x
    3246        print *,' Arg. de HBTM: yt_x ',yt_x
    3247        print *,' Arg. de HBTM: yq_x ',yq_x
     3240       print *,' Arg. de HBTM: yt2m_x ',yt2m_x(1:knon)
     3241       print *,' Arg. de HBTM: yt10m_x ',yt10m_x(1:knon)
     3242       print *,' Arg. de HBTM: yq2m_x ',yq2m_x(1:knon)
     3243       print *,' Arg. de HBTM: yq10m_x ',yq10m_x(1:knon)
     3244       print *,' Arg. de HBTM: yustar_x ',yustar_x(1:knon)
     3245       print *,' Arg. de HBTM: y_flux_t_x ',y_flux_t_x(1:knon,:)
     3246       print *,' Arg. de HBTM: y_flux_q_x ',y_flux_q_x(1:knon,:)
     3247       print *,' Arg. de HBTM: yu_x ',yu_x(1:knon,:)
     3248       print *,' Arg. de HBTM: yv_x ',yv_x(1:knon,:)
     3249       print *,' Arg. de HBTM: yt_x ',yt_x(1:knon,:)
     3250       print *,' Arg. de HBTM: yq_x ',yq_x(1:knon,:)
    32483251          ENDIF
    32493252        CALL HBTM(knon, ypaprs, ypplay, &
    … …  
    32743277        ENDDO
    32753278        IF (prt_level >=10) THEN
    3276           print *, 'After HBTM: pblh ', pblh
    3277           print *, 'After HBTM: plcl ', plcl
    3278           print *, 'After HBTM: cteiCL ', cteiCL
     3279          print *, 'After HBTM: pblh ', pblh(1:knon,:)
     3280          print *, 'After HBTM: plcl ', plcl(1:knon,:)
     3281          print *, 'After HBTM: cteiCL ', cteiCL(1:knon,:)
    32793282        ENDIF
    32803283       ELSE  !(iflag_split .eq.0)
    … …  
    32943297        ENDDO
    32953298        IF (prt_level >=10) THEN
    3296           print *, 'After HBTM: pblh_x ', pblh_x
    3297           print *, 'After HBTM: plcl_x ', plcl_x
    3298           print *, 'After HBTM: cteiCL_x ', cteiCL_x
     3299          print *, 'After HBTM: pblh_x ', pblh_x(1:knon,:)
     3300          print *, 'After HBTM: plcl_x ', plcl_x(1:knon,:)
     3301          print *, 'After HBTM: cteiCL_x ', cteiCL_x(1:knon,:)
    32993302        ENDIF
    33003303        DO j=1, knon
    … …  
    33133316        ENDDO
    33143317        IF (prt_level >=10) THEN
    3315           print *, 'After HBTM: pblh_w ', pblh_w
    3316           print *, 'After HBTM: plcl_w ', plcl_w
    3317           print *, 'After HBTM: cteiCL_w ', cteiCL_w
     3318          print *, 'After HBTM: pblh_w ', pblh_w(1:knon,:)
     3319          print *, 'After HBTM: plcl_w ', plcl_w(1:knon,:)
     3320          print *, 'After HBTM: cteiCL_w ', cteiCL_w(1:knon,:)
    33183321        ENDIF
    33193322!!!
Note: See TracChangeset for help on using the changeset viewer.