Changeset 5791 for LMDZ6/branches/contrails/libf/phylmd/cv3a_compress.f90

Timestamp:

Jul 28, 2025, 7:23:15 PM (7 days ago)

Author:

aborella

Message:

Merge with trunk r5789

Location:

LMDZ6/branches/contrails

Files:

• . (modified) (1 prop)
• libf/phylmd/cv3a_compress.f90 (modified) (4 diffs)

LMDZ6/branches/contrails/libf/phylmd/cv3a_compress.f90

@@ +1 @@
+MODULE cv3a_compress_mod
+
+CONTAINS
+
 SUBROUTINE cv3a_compress(len, nloc, ncum, nd, ntra, compress, &
                          iflag1, nk1, icb1, icbs1, &
@@ -92 +96 @@
   INTEGER i, k, nn, j
 
-  CHARACTER (LEN=20) :: modname = 'cv3a_compress'
+  CHARACTER (LEN=20),PARAMETER :: modname = 'cv3a_compress'
   CHARACTER (LEN=80) :: abort_message
 
@@ -191 +195 @@
   ELSE  !(compress)
 !
-      ncum = len
-!
-      wghti(:,1:nl+1) = wghti1(:,1:nl+1)
-      t(:,1:nl+1) = t1(:,1:nl+1)
-      q(:,1:nl+1) = q1(:,1:nl+1)
-      qs(:,1:nl+1) = qs1(:,1:nl+1)
-      t_wake(:,1:nl+1) = t1_wake(:,1:nl+1)
-      q_wake(:,1:nl+1) = q1_wake(:,1:nl+1)
-      qs_wake(:,1:nl+1) = qs1_wake(:,1:nl+1)
-      u(:,1:nl+1) = u1(:,1:nl+1)
-      v(:,1:nl+1) = v1(:,1:nl+1)
-      gz(:,1:nl+1) = gz1(:,1:nl+1)
-      th(:,1:nl+1) = th1(:,1:nl+1)
-      th_wake(:,1:nl+1) = th1_wake(:,1:nl+1)
-      h(:,1:nl+1) = h1(:,1:nl+1)
-      lv(:,1:nl+1) = lv1(:,1:nl+1)
-      lf(:,1:nl+1) = lf1(:,1:nl+1)
-      cpn(:,1:nl+1) = cpn1(:,1:nl+1)
-      p(:,1:nl+1) = p1(:,1:nl+1)
-      ph(:,1:nl+1) = ph1(:,1:nl+1)
-      tv(:,1:nl+1) = tv1(:,1:nl+1)
-      tp(:,1:nl+1) = tp1(:,1:nl+1)
-      tvp(:,1:nl+1) = tvp1(:,1:nl+1)
-      clw(:,1:nl+1) = clw1(:,1:nl+1)
-      h_wake(:,1:nl+1) = h1_wake(:,1:nl+1)
-      lv_wake(:,1:nl+1) = lv1_wake(:,1:nl+1)
-      lf_wake(:,1:nl+1) = lf1_wake(:,1:nl+1)
-      cpn_wake(:,1:nl+1) = cpn1_wake(:,1:nl+1)
-      tv_wake(:,1:nl+1) = tv1_wake(:,1:nl+1)
-      sig(:,1:nl+1) = sig1(:,1:nl+1)
-      w0(:,1:nl+1) = w01(:,1:nl+1)
-      omega(:,1:nl+1) = omega1(:,1:nl+1)
-!
-      s_wake(:) = s1_wake(:)
-      iflag(:) = iflag1(:)
-      nk(:) = nk1(:)
-      icb(:) = icb1(:)
-      icbs(:) = icbs1(:)
-      plcl(:) = plcl1(:)
-      tnk(:) = tnk1(:)
-      qnk(:) = qnk1(:)
-      gznk(:) = gznk1(:)
-      hnk(:) = hnk1(:)
-      unk(:) = unk1(:)
-      vnk(:) = vnk1(:)
-      pbase(:) = pbase1(:)
-      buoybase(:) = buoybase1(:)
-      sig(:, nd) = sig1(:, nd)
-      ptop2(:) = ptop2(:)
-      Ale(:) = Ale1(:)
-      Alp(:) = Alp1(:)
+    wghti(:,1:nl+1) = wghti1(:,1:nl+1)
+    t(:,1:nl+1) = t1(:,1:nl+1)
+    q(:,1:nl+1) = q1(:,1:nl+1)
+    qs(:,1:nl+1) = qs1(:,1:nl+1)
+    t_wake(:,1:nl+1) = t1_wake(:,1:nl+1)
+    q_wake(:,1:nl+1) = q1_wake(:,1:nl+1)
+    qs_wake(:,1:nl+1) = qs1_wake(:,1:nl+1)
+    u(:,1:nl+1) = u1(:,1:nl+1)
+    v(:,1:nl+1) = v1(:,1:nl+1)
+    gz(:,1:nl+1) = gz1(:,1:nl+1)
+    th(:,1:nl+1) = th1(:,1:nl+1)
+    th_wake(:,1:nl+1) = th1_wake(:,1:nl+1)
+    h(:,1:nl+1) = h1(:,1:nl+1)
+    lv(:,1:nl+1) = lv1(:,1:nl+1)
+    lf(:,1:nl+1) = lf1(:,1:nl+1)
+    cpn(:,1:nl+1) = cpn1(:,1:nl+1)
+    p(:,1:nl+1) = p1(:,1:nl+1)
+    ph(:,1:nl+1) = ph1(:,1:nl+1)
+    tv(:,1:nl+1) = tv1(:,1:nl+1)
+    tp(:,1:nl+1) = tp1(:,1:nl+1)
+    tvp(:,1:nl+1) = tvp1(:,1:nl+1)
+    clw(:,1:nl+1) = clw1(:,1:nl+1)
+    h_wake(:,1:nl+1) = h1_wake(:,1:nl+1)
+    lv_wake(:,1:nl+1) = lv1_wake(:,1:nl+1)
+    lf_wake(:,1:nl+1) = lf1_wake(:,1:nl+1)
+    cpn_wake(:,1:nl+1) = cpn1_wake(:,1:nl+1)
+    tv_wake(:,1:nl+1) = tv1_wake(:,1:nl+1)
+    sig(:,1:nl+1) = sig1(:,1:nl+1)
+    w0(:,1:nl+1) = w01(:,1:nl+1)
+    omega(:,1:nl+1) = omega1(:,1:nl+1)
+
+    s_wake(:) = s1_wake(:)
+    iflag(:) = iflag1(:)
+    nk(:) = nk1(:)
+    icb(:) = icb1(:)
+    icbs(:) = icbs1(:)
+    plcl(:) = plcl1(:)
+    tnk(:) = tnk1(:)
+    qnk(:) = qnk1(:)
+    gznk(:) = gznk1(:)
+    hnk(:) = hnk1(:)
+    unk(:) = unk1(:)
+    vnk(:) = vnk1(:)
+    pbase(:) = pbase1(:)
+    buoybase(:) = buoybase1(:)
+    sig(:, nd) = sig1(:, nd)
+    ptop2(:) = ptop2(:)
+    Ale(:) = Ale1(:)
+    Alp(:) = Alp1(:)
 !
   ENDIF !(compress)
@@ -248 +250 @@
   RETURN
 END SUBROUTINE cv3a_compress
+
+END MODULE cv3a_compress_mod