Changeset 3707 for LMDZ6/branches


Ignore:
Timestamp:
Jun 11, 2020, 11:09:39 AM (4 years ago)
Author:
adurocher
Message:

New generic implementation for compression

Location:
LMDZ6/branches/Optimisation_LMDZ/libf/phylmd
Files:
1 added
2 edited
1 moved

Legend:

Unmodified
Added
Removed
  • LMDZ6/branches/Optimisation_LMDZ/libf/phylmd/cv3a_compress_old.F90

    r3706 r3707  
    1 SUBROUTINE cv3a_compress(len, nloc, ncum, nd, ntra, compress, &
     1SUBROUTINE cv3a_compress_old(len, nloc, ncum, nd, ntra, compress, &
    22                         iflag1, nk1, icb1, icbs1, &
    33                         plcl1, tnk1, qnk1, gznk1, hnk1, unk1, vnk1, &
     
    247247
    248248  RETURN
    249 END SUBROUTINE cv3a_compress
     249END SUBROUTINE cv3a_compress_old
  • LMDZ6/branches/Optimisation_LMDZ/libf/phylmd/cv3a_uncompress.F90

    r3496 r3707  
    1 SUBROUTINE cv3a_uncompress(nloc, len, ncum, nd, ntra, idcum, compress, &
     1SUBROUTINE cv3a_uncompress_old(nloc, len, ncum, nd, ntra, idcum, compress, &
    22                           iflag, kbas, ktop, &
    33                           precip, cbmf, plcl, plfc, wbeff, sig, w0, ptop2, &
     
    332332
    333333  RETURN
    334 END SUBROUTINE cv3a_uncompress
    335 
     334END SUBROUTINE cv3a_uncompress_old
     335
  • LMDZ6/branches/Optimisation_LMDZ/libf/phylmd/cva_driver.F90

    r3496 r3707  
    4141  USE print_control_mod, ONLY: prt_level, lunout
    4242  USE add_phys_tend_mod, ONLY: fl_cor_ebil
     43  USE cv3a_compress_mod
    4344  IMPLICIT NONE
    4445
     
    542543!$OMP THREADPRIVATE(igout)
    543544
     545  type(compress_data_t) :: compress_data
     546  type(array_list) :: cv3a_compress_list, cv3a_uncompress_list
    544547
    545548! print *, 't1, t1_wake ',(k,t1(1,k),t1_wake(1,k),k=1,nd)
     
    843846!   If the fraction of convective points is larger than comp_threshold, then compression
    844847!   is assumed useless.
    845 !
    846   compress = ncum .lt. len*comp_threshold
    847 !
    848   IF (.not. compress) THEN
    849     DO i = 1,len
    850       idcum(i) = i
    851     ENDDO
    852   ENDIF
    853 !
    854 !>jyg
    855         if (prt_level >= 9) &
    856              PRINT *, 'cva_driver -> cv3a_compress'
    857       CALL cv3a_compress(len, nloc, ncum, nd, ntra, compress, &
    858                          iflag1, nk1, icb1, icbs1, &
    859                          plcl1, tnk1, qnk1, gznk1, hnk1, unk1, vnk1, &
    860                          wghti1, pbase1, buoybase1, &
    861                          t1, q1, qs1, t1_wake, q1_wake, qs1_wake, s1_wake, &
    862                          u1, v1, gz1, th1, th1_wake, &
    863                          tra1, &
    864                          h1, lv1, lf1, cpn1, p1, ph1, tv1, tp1, tvp1, clw1, &
    865                          h1_wake, lv1_wake, lf1_wake, cpn1_wake, tv1_wake, &
    866                          sig1, w01, ptop21, &
    867                          Ale1, Alp1, omega1, &
    868                          iflag, nk, icb, icbs, &
    869                          plcl, tnk, qnk, gznk, hnk, unk, vnk, &
    870                          wghti, pbase, buoybase, &
    871                          t, q, qs, t_wake, q_wake, qs_wake, s_wake, &
    872                          u, v, gz, th, th_wake, &
    873                          tra, &
    874                          h, lv, lf, cpn, p, ph, tv, tp, tvp, clw, &
    875                          h_wake, lv_wake, lf_wake, cpn_wake, tv_wake, &
    876                          sig, w0, ptop2, &
    877                          Ale, Alp, omega)
     848  compress = ncum <= len*comp_threshold
     849
     850  if( compress ) then
     851    compress_mode = COMPRESS_MODE_COMPRESS
     852  else
     853    compress_mode = COMPRESS_MODE_COPY
     854  endif
     855
     856  if (prt_level >= 9) &
     857        PRINT *, 'cva_driver -> cv3a_compress'
     858
     859  call add_array_i1(cv3a_compress_list, iflag1, iflag)
     860  call add_array_i1(cv3a_compress_list, nk1, nk)
     861  call add_array_i1(cv3a_compress_list, icb1, icb)
     862  call add_array_i1(cv3a_compress_list, icbs1, icbs)
     863  call add_array_r1(cv3a_compress_list, plcl1, plcl)
     864  call add_array_r1(cv3a_compress_list, tnk1, tnk)
     865  call add_array_r1(cv3a_compress_list, qnk1, qnk)
     866  call add_array_r1(cv3a_compress_list, gznk1, gznk)
     867  call add_array_r1(cv3a_compress_list, hnk1, hnk)
     868  call add_array_r1(cv3a_compress_list, unk1, unk)
     869  call add_array_r1(cv3a_compress_list, vnk1, vnk)
     870  call add_array_r2(cv3a_compress_list, wghti1, wghti)
     871  call add_array_r1(cv3a_compress_list, pbase1, pbase)
     872  call add_array_r1(cv3a_compress_list, buoybase1, buoybase)
     873  call add_array_r2(cv3a_compress_list, th1, th)
     874  call add_array_r2(cv3a_compress_list, t1, t)
     875  call add_array_r2(cv3a_compress_list, q1, q)
     876  call add_array_r2(cv3a_compress_list, qs1, qs)
     877  call add_array_r2(cv3a_compress_list, t1_wake, t_wake)
     878  call add_array_r2(cv3a_compress_list, q1_wake, q_wake)
     879  call add_array_r2(cv3a_compress_list, qs1_wake, qs_wake)
     880  call add_array_r1(cv3a_compress_list, s1_wake, s_wake)
     881  call add_array_r2(cv3a_compress_list, u1, u)
     882  call add_array_r2(cv3a_compress_list, v1, v)
     883  call add_array_r2(cv3a_compress_list, gz1, gz)
     884  call add_array_r2(cv3a_compress_list, h1, h)
     885  call add_array_r2(cv3a_compress_list, th1_wake, th_wake)
     886  !call add_array_r3(cv3a_compress_list, tra1, tra) !
     887  call add_array_r2(cv3a_compress_list, lv1, lv)
     888  call add_array_r2(cv3a_compress_list, lf1, lf)
     889  call add_array_r2(cv3a_compress_list, cpn1, cpn)
     890  call add_array_r2(cv3a_compress_list, p1, p)
     891  call add_array_r2(cv3a_compress_list, ph1, ph)
     892  call add_array_r2(cv3a_compress_list, tv1, tv)
     893  call add_array_r2(cv3a_compress_list, tp1, tp)
     894  call add_array_r2(cv3a_compress_list, tvp1, tvp)
     895  call add_array_r2(cv3a_compress_list, clw1, clw)
     896  call add_array_r2(cv3a_compress_list, h1_wake, h_wake)
     897  call add_array_r2(cv3a_compress_list, lv1_wake, lv_wake)
     898  call add_array_r2(cv3a_compress_list, lf1_wake, lf_wake)
     899  call add_array_r2(cv3a_compress_list, cpn1_wake, cpn_wake)
     900  call add_array_r2(cv3a_compress_list, tv1_wake, tv_wake)
     901  call add_array_r2(cv3a_compress_list, sig1, sig)
     902  call add_array_r1(cv3a_compress_list, sig1(:,nd), sig(:,nd))
     903  call add_array_r2(cv3a_compress_list, w01, w0)
     904  call add_array_r1(cv3a_compress_list, ptop21, ptop2)
     905  call add_array_r1(cv3a_compress_list, Ale1, Ale)
     906  call add_array_r1(cv3a_compress_list, Alp1, Alp)
     907  call add_array_r2(cv3a_compress_list, omega1, omega)
     908 
     909  call cv3a_compress(len, (iflag1 == 0), cv3a_compress_list, compress_data)
     910  ncum = compress_data%ncum
     911      ! IF (.not. compress) THEN
     912      !   DO i = 1,len
     913      !     idcum(i) = i
     914      !   ENDDO
     915      ! ENDIF
     916      ! CALL cv3a_compress_old(len, nloc, ncum, nd, ntra, compress, &
     917      !                    iflag1, nk1, icb1, icbs1, &
     918      !                    plcl1, tnk1, qnk1, gznk1, hnk1, unk1, vnk1, &
     919      !                    wghti1, pbase1, buoybase1, &
     920      !                    t1, q1, qs1, t1_wake, q1_wake, qs1_wake, s1_wake, &
     921      !                    u1, v1, gz1, th1, th1_wake, &
     922      !                    tra1, &
     923      !                    h1, lv1, lf1, cpn1, p1, ph1, tv1, tp1, tvp1, clw1, &
     924      !                    h1_wake, lv1_wake, lf1_wake, cpn1_wake, tv1_wake, &
     925      !                    sig1, w01, ptop21, &
     926      !                    Ale1, Alp1, omega1, &
     927      !                    iflag, nk, icb, icbs, &
     928      !                    plcl, tnk, qnk, gznk, hnk, unk, vnk, &
     929      !                    wghti, pbase, buoybase, &
     930      !                    t, q, qs, t_wake, q_wake, qs_wake, s_wake, &
     931      !                    u, v, gz, th, th_wake, &
     932      !                    tra, &
     933      !                    h, lv, lf, cpn, p, ph, tv, tp, tvp, clw, &
     934      !                    h_wake, lv_wake, lf_wake, cpn_wake, tv_wake, &
     935      !                    sig, w0, ptop2, &
     936      !                    Ale, Alp, omega)
    878937
    879938! print*,'tv ',tv
     
    885944        if (prt_level >= 9) &
    886945             PRINT *, 'cva_driver -> cv_compress'
     946      ! TODO : new compress interface could be used here too
    887947      CALL cv_compress(len, nloc, ncum, nd, &
    888948                       iflag1, nk1, icb1, &
     
    11921252
    11931253    IF (iflag_con==3) THEN
    1194         if (prt_level >= 9) &
    1195              PRINT *, 'cva_driver -> cv3a_uncompress'
    1196       CALL cv3a_uncompress(nloc, len, ncum, nd, ntra, idcum, compress, &
    1197                            iflag, icb, inb, &
    1198                            precip, cbmf, plcl, plfc, wbeff, sig, w0, ptop2, &
    1199                            ft, fq, fu, fv, ftra, &
    1200                            sigd, ma, mip, vprecip, vprecipi, upwd, dnwd, dnwd0, &
    1201                            qcondc, wd, cape, cin, &
    1202                            tvp, &
    1203                            ftd, fqd, &
    1204                            Plim1, plim2, asupmax, supmax0, &
    1205                            asupmaxmin, &
    1206                            da, phi, mp, phi2, d1a, dam, sigij, &         ! RomP
    1207                            qta, clw, elij, evap, ep, epmlmMm, eplaMm, &  ! RomP
    1208                            wdtrainA, wdtrainS, wdtrainM, &                         ! RomP
    1209                            qtc, sigt, epmax_diag, & ! epmax_cape
    1210                            iflag1, kbas1, ktop1, &
    1211                            precip1, cbmf1, plcl1, plfc1, wbeff1, sig1, w01, ptop21, &
    1212                            ft1, fq1, fu1, fv1, ftra1, &
    1213                            sigd1, ma1, mip1, vprecip1, vprecipi1, upwd1, dnwd1, dnwd01, &
    1214                            qcondc1, wd1, cape1, cin1, &
    1215                            tvp1, &
    1216                            ftd1, fqd1, &
    1217                            Plim11, plim21, asupmax1, supmax01, &
    1218                            asupmaxmin1, &
    1219                            da1, phi1, mp1, phi21, d1a1, dam1, sigij1,  &       ! RomP
    1220                            qta1, clw1, elij1, evap1, ep1, epmlmMm1, eplaMm1, & ! RomP
    1221                            wdtrainA1, wdtrainS1, wdtrainM1,                  & ! RomP
    1222                            qtc1, sigt1, epmax_diag1) ! epmax_cape
     1254      if (prt_level >= 9) &
     1255            PRINT *, 'cva_driver -> cv3a_uncompress'
     1256      call add_array_i1(cv3a_uncompress_list, iflag, iflag1)
     1257      call add_array_i1(cv3a_uncompress_list, icb, icb1)
     1258      call add_array_i1(cv3a_uncompress_list, inb, ktop1)
     1259      call add_array_r1(cv3a_uncompress_list, precip, precip1)
     1260      call add_array_r1(cv3a_uncompress_list, cbmf, cbmf1)
     1261      call add_array_r1(cv3a_uncompress_list, plcl, plcl1)
     1262      call add_array_r1(cv3a_uncompress_list, plfc, plfc1)
     1263      call add_array_r1(cv3a_uncompress_list, wbeff, wbeff1)
     1264      call add_array_r2(cv3a_uncompress_list, sig, sig1)
     1265      call add_array_r2(cv3a_uncompress_list, w0, w01)
     1266      call add_array_r1(cv3a_uncompress_list, ptop2, ptop21)
     1267      call add_array_r2(cv3a_uncompress_list, ft, ft1)
     1268      call add_array_r2(cv3a_uncompress_list, fq, fq1)
     1269      call add_array_r2(cv3a_uncompress_list, fu, fu1)
     1270      call add_array_r2(cv3a_uncompress_list, fv, fv1)
     1271      !call add_array_r3(cv3a_uncompress_list, ftra, ftra1)
     1272      call add_array_r1(cv3a_uncompress_list, sigd, sigd1)
     1273      call add_array_r2(cv3a_uncompress_list, ma, ma1)
     1274      call add_array_r2(cv3a_uncompress_list, mip, mip1)
     1275      call add_array_r2(cv3a_uncompress_list, vprecip, vprecip1)
     1276      call add_array_r2(cv3a_uncompress_list, vprecipi, vprecipi1)
     1277      call add_array_r2(cv3a_uncompress_list, upwd, upwd1)
     1278      call add_array_r2(cv3a_uncompress_list, dnwd, dnwd1)
     1279      call add_array_r2(cv3a_uncompress_list, dnwd0, dnwd01)
     1280      call add_array_r2(cv3a_uncompress_list, qcondc, qcondc1)
     1281      call add_array_r1(cv3a_uncompress_list, wd, wd1)
     1282      call add_array_r1(cv3a_uncompress_list, cape, cape1)
     1283      call add_array_r1(cv3a_uncompress_list, cin, cin1)
     1284      call add_array_r2(cv3a_uncompress_list, tvp, tvp1)
     1285      call add_array_r2(cv3a_uncompress_list, ftd, ftd1)
     1286      call add_array_r2(cv3a_uncompress_list, fqd, fqd1)
     1287      call add_array_r1(cv3a_uncompress_list, Plim1, Plim11)
     1288      call add_array_r1(cv3a_uncompress_list, plim2, plim21)
     1289      call add_array_r2(cv3a_uncompress_list, asupmax, asupmax1)
     1290      call add_array_r1(cv3a_uncompress_list, supmax0, supmax01)
     1291      call add_array_r1(cv3a_uncompress_list, asupmaxmin, asupmaxmin1)
     1292      call add_array_r2(cv3a_uncompress_list, da, da1)
     1293      call add_array_r3(cv3a_uncompress_list, phi, phi1)
     1294      call add_array_r2(cv3a_uncompress_list, mp, mp1)
     1295      call add_array_r3(cv3a_uncompress_list, phi2, phi21)
     1296      call add_array_r2(cv3a_uncompress_list, d1a, d1a1)
     1297      call add_array_r2(cv3a_uncompress_list, dam, dam1)
     1298      call add_array_r3(cv3a_uncompress_list, sigij, sigij1)
     1299      call add_array_r2(cv3a_uncompress_list, qta, qta1)
     1300      call add_array_r2(cv3a_uncompress_list, clw, clw1)
     1301      call add_array_r3(cv3a_uncompress_list, elij, elij1)
     1302      call add_array_r2(cv3a_uncompress_list, evap, evap1)
     1303      call add_array_r2(cv3a_uncompress_list, ep, ep1)
     1304      call add_array_r3(cv3a_uncompress_list, epmlmMm, epmlmMm1)
     1305      call add_array_r2(cv3a_uncompress_list, eplaMm, eplaMm1)
     1306      call add_array_r2(cv3a_uncompress_list, wdtrainA, wdtrainA1)
     1307      call add_array_r2(cv3a_uncompress_list, wdtrainS, wdtrainS1)
     1308      call add_array_r2(cv3a_uncompress_list, wdtrainM, wdtrainM1)
     1309      call add_array_r2(cv3a_uncompress_list, qtc, qtc1)
     1310      call add_array_r2(cv3a_uncompress_list, sigt, sigt1)
     1311      call add_array_r1(cv3a_uncompress_list, epmax_diag, epmax_diag1)
     1312      call add_array_r1(cv3a_uncompress_list, sig(:,nd), sig1(:,nd))
     1313      call cv3a_uncompress(len, compress_data, cv3a_uncompress_list)
     1314      ! CALL cv3a_uncompress_old(nloc, len, ncum, nd, ntra, idcum, compress, &
     1315      !                      iflag, icb, inb, &
     1316      !                      precip, cbmf, plcl, plfc, wbeff, sig, w0, ptop2, &
     1317      !                      ft, fq, fu, fv, ftra, &
     1318      !                      sigd, ma, mip, vprecip, vprecipi, upwd, dnwd, dnwd0, &
     1319      !                      qcondc, wd, cape, cin, &
     1320      !                      tvp, &
     1321      !                      ftd, fqd, &
     1322      !                      Plim1, plim2, asupmax, supmax0, &
     1323      !                      asupmaxmin, &
     1324      !                      da, phi, mp, phi2, d1a, dam, sigij, &         ! RomP
     1325      !                      qta, clw, elij, evap, ep, epmlmMm, eplaMm, &  ! RomP
     1326      !                      wdtrainA, wdtrainS, wdtrainM, &                         ! RomP
     1327      !                      qtc, sigt, epmax_diag, & ! epmax_cape
     1328      !                      iflag1, kbas1, ktop1, &
     1329      !                      precip1, cbmf1, plcl1, plfc1, wbeff1, sig1, w01, ptop21, &
     1330      !                      ft1, fq1, fu1, fv1, ftra1, &
     1331      !                      sigd1, ma1, mip1, vprecip1, vprecipi1, upwd1, dnwd1, dnwd01, &
     1332      !                      qcondc1, wd1, cape1, cin1, &
     1333      !                      tvp1, &
     1334      !                      ftd1, fqd1, &
     1335      !                      Plim11, plim21, asupmax1, supmax01, &
     1336      !                      asupmaxmin1, &
     1337      !                      da1, phi1, mp1, phi21, d1a1, dam1, sigij1,  &       ! RomP
     1338      !                      qta1, clw1, elij1, evap1, ep1, epmlmMm1, eplaMm1, & ! RomP
     1339      !                      wdtrainA1, wdtrainS1, wdtrainM1,                  & ! RomP
     1340      !                      qtc1, sigt1, epmax_diag1) ! epmax_cape
    12231341!   
    12241342      IF (prt_level >= 10) THEN
Note: See TracChangeset for help on using the changeset viewer.