Changeset 5113 for LMDZ6/branches/Amaury_dev/libf/dyn3dmem/caldyn_loc.f90

Timestamp:

Jul 24, 2024, 1:17:08 PM (2 months ago)

Author:

abarral

Message:

Rename modules in misc from *_mod > lmdz_*
Put cbrt.f90, ch*.f90, pch*.f90 in new lmdz_libmath_pch.f90

File:

• LMDZ6/branches/Amaury_dev/libf/dyn3dmem/caldyn_loc.f90 (modified) (5 diffs)

LMDZ6/branches/Amaury_dev/libf/dyn3dmem/caldyn_loc.f90

@@ -66 +66 @@
   !--------------------------------
 
-  ! ! compute contravariant winds ucont() and vcont
+  ! compute contravariant winds ucont() and vcont
   CALL covcont_loc  (llm, ucov, vcov, ucont, vcont)
-  ! ! compute pressure p()
+  ! compute pressure p()
   CALL pression_loc (ip1jmp1, ap, bp, ps, p)
   !ym      CALL psextbar (   ps   , psexbarxy                          )
   !$OMP BARRIER
-  ! ! compute mass in each atmospheric mesh: masse()
+  ! compute mass in each atmospheric mesh: masse()
   CALL massdair_loc (p, masse)
-  ! ! compute X and Y-averages of mass, massebx() and masseby()
+  ! compute X and Y-averages of mass, massebx() and masseby()
   CALL massbar_loc  (masse, massebx, masseby)
-  ! ! compute XY-average of mass, massebxy()
+  ! compute XY-average of mass, massebxy()
   CALL massbarxy_loc(masse, massebxy)
-  ! ! compute mass fluxes pbaru() and pbarv()
+  ! compute mass fluxes pbaru() and pbarv()
   CALL flumass_loc  (massebx, masseby, vcont, ucont, pbaru, pbarv)
-  ! ! compute dteta() , horizontal converging flux of theta
+  ! compute dteta() , horizontal converging flux of theta
   CALL dteta1_loc   (teta, pbaru, pbarv, dteta)
-  ! ! compute convm(), horizontal converging flux of mass
+  ! compute convm(), horizontal converging flux of mass
   CALL convmas1_loc  (pbaru, pbarv, convm)
   !$OMP BARRIER
@@ -105 +105 @@
   ijb = ij_begin
   ije = ij_end
-  ! ! compute pressure variation due to mass convergence
+  ! compute pressure variation due to mass convergence
   DO ij = ijb, ije
     dp(ij) = convm(ij, 1) / airesurg(ij)
@@ -112 +112 @@
   !$OMP BARRIER
 
-  ! ! compute vertical velocity w()
+  ! compute vertical velocity w()
   CALL vitvert_loc (convm, w)
-  ! ! compute potential vorticity vorpot()
+  ! compute potential vorticity vorpot()
   CALL tourpot_loc (vcov, ucov, massebxy, vorpot)
-  ! ! compute rotation induced du() and dv()
+  ! compute rotation induced du() and dv()
   CALL dudv1_loc   (vorpot, pbaru, pbarv, du, dv)
 
@@ -126 +126 @@
   END IF
 
-  ! ! compute kinetic energy ecin()
+  ! compute kinetic energy ecin()
   CALL enercin_loc (vcov, ucov, vcont, ucont, ecin)
-  ! ! compute Bernouilli function bern()
+  ! compute Bernouilli function bern()
   CALL bernoui_loc (ip1jmp1, llm, phi, ecin, bern)
-  ! ! compute and add du() and dv() contributions from Bernouilli and pressure
+  ! compute and add du() and dv() contributions from Bernouilli and pressure
   CALL dudv2_loc   (teta, pkf, bern, du, dv)
 
@@ -155 +155 @@
   !$OMP END DO
 
-  ! ! compute vertical advection contributions to du(), dv() and dteta()
+  ! compute vertical advection contributions to du(), dv() and dteta()
   CALL advect_new_loc(ang, vcov, teta, w, massebx, masseby, du, dv, dteta)