Changeset 1523 for trunk/LMDZ.COMMON/libf/dyn3d_common/convmas.F90

Timestamp:

Mar 28, 2016, 5:27:51 PM (9 years ago)

Author:

emillour

Message:

All models: More updates to make planetary codes (+Earth) setups converge.

• in dyn3d_common:
• convmas.F => convmas.F90
• enercin.F => enercin.F90
• flumass.F => flumass.F90
• massbar.F => massbar.F90
• tourpot.F => tourpot.F90
• vitvert.F => vitvert.F90

• in misc:
• move "q_sat" from "dyn3d_common" to "misc" (in Earth model, it is also called by the physics)
• move "write_field" from "dyn3d_common" to "misc"(may be called from physics or dynamics and depends on neither).

• in phy_common:
• move "write_field_phy" here since it may be called from any physics package)
• add module "regular_lonlat_mod" to store global information on lon-lat grid

• in dynlonlat_phylonlat/phy*:
• turn "iniphysiq.F90" into module "iniphysiq_mod.F90" (and of course adapt gcm.F[90] and 1D models accordingly)

EM

File:

• trunk/LMDZ.COMMON/libf/dyn3d_common/convmas.F90 (moved) (moved from trunk/LMDZ.COMMON/libf/dyn3d_common/convmas.F) (1 diff)

trunk/LMDZ.COMMON/libf/dyn3d_common/convmas.F90

@@ +1 @@
+SUBROUTINE convmas (pbaru, pbarv, convm)
 !
-! $Header$
-!
-      SUBROUTINE convmas (pbaru, pbarv, convm )
-c
-      IMPLICIT NONE
+!-------------------------------------------------------------------------------
+! Authors: P. Le Van , Fr. Hourdin.
+!-------------------------------------------------------------------------------
+! Purpose: Compute mass flux convergence at p levels.
+  IMPLICIT NONE
+  include "dimensions.h"
+  include "paramet.h"
+  include "comgeom.h"
+!===============================================================================
+! Arguments:
+  REAL, INTENT(IN)  :: pbaru(ip1jmp1,llm)
+  REAL, INTENT(IN)  :: pbarv(ip1jm  ,llm)
+  REAL, INTENT(OUT) :: convm(ip1jmp1,llm)
+!===============================================================================
+! Method used:   Computation from top to bottom.
+!   Mass convergence at level llm is equal to zero and is not stored in convm.
+!===============================================================================
+! Local variables:
+  INTEGER :: l
+!===============================================================================
 
-c=======================================================================
-c
-c   Auteurs:  P. Le Van , F. Hourdin  .
-c   -------
-c
-c   Objet:
-c   ------
-c
-c   ********************************************************************
-c   .... calcul de la convergence du flux de masse aux niveaux p ...
-c   ********************************************************************
-c
-c
-c     pbaru  et  pbarv  sont des arguments d'entree pour le s-pg  ....
-c      .....  convm      est  un argument de sortie pour le s-pg  ....
-c
-c    le calcul se fait de haut en bas, 
-c    la convergence de masse au niveau p(llm+1) est egale a 0. et
-c    n'est pas stockee dans le tableau convm .
-c
-c
-c=======================================================================
-c
-c   Declarations:
-c   -------------
+!--- Computation of - (d(pbaru)/dx + d(pbarv)/dy )
+  CALL convflu( pbaru, pbarv, llm, convm )
 
-#include "dimensions.h"
-#include "paramet.h"
+!--- Filter
+  CALL filtreg( convm, jjp1, llm, 2, 2, .TRUE., 1 )
 
-      REAL pbaru( ip1jmp1,llm ),pbarv( ip1jm,llm ),convm(  ip1jmp1,llm )
-      INTEGER   l,ij
+!--- Mass convergence is integrated from top to bottom
+  DO l=llmm1,1,-1
+    convm(:,l) = convm(:,l) + convm(:,l+1)
+  END DO
 
-
-c-----------------------------------------------------------------------
-c    ....  calcul de - (d(pbaru)/dx + d(pbarv)/dy ) ......
-
-      CALL  convflu( pbaru, pbarv, llm, convm )
-
-c-----------------------------------------------------------------------
-c   filtrage:
-c   ---------
-
-       CALL filtreg( convm, jjp1, llm, 2, 2, .true., 1 )
-
-c    integration de la convergence de masse de haut  en bas ......
-
-      DO      l      = llmm1, 1, -1
-        DO    ij     = 1, ip1jmp1
-         convm(ij,l) = convm(ij,l) + convm(ij,l+1)
-        ENDDO
-      ENDDO
-c
-      RETURN
-      END
+END SUBROUTINE convmas