Changeset 1508 for trunk/LMDZ.COMMON/libf/dyn3d_common/covcont.F90

Timestamp:

Jan 15, 2016, 8:27:16 AM (9 years ago)

Author:

emillour

Message:

Common dynamics:
Updates in the dynamics (seq and ) to keep up with updates
in LMDZ5 (up to LMDZ5 trunk, rev 2325):
IMPORTANT: Modifications for isotopes are only done in dyn3d, not in dyn3dpar

as in LMDZ5 these modifications were done in dyn3dmem.
Related LMDZ5 revisions are r2270 and r2281

• in dynlonlat_phylonlat:
• add module "grid_atob_m.F90" (a regridding utility so far only used by phylmd/ce0l.F90, used to be dyn3d_common/grid_atob.F)

• in misc:
• follow up updates on wxios.F (add missing_val module variable)

• in dyn3d_common:
• pression.F => pression.F90
• misc_mod.F90: moved from misc to dyn3d_common
• added new iso_verif_dyn.F
• covcont.F => covcont.F90
• infotrac.F90 : add handling of isotopes (reading of corresponding traceur.def for planets not implemented)
• dynetat0.F => dynetat0.F90 with some code factorization
• dynredem.F => dynredem.F90 with some code factorization
• added dynredem_mod.F90: routines used by dynredem
• iniacademic.F90 : added isotopes-related initialization for Earth case

• in dyn3d:
• added check_isotopes.F
• modified (isotopes) advtrac.F90, caladvtrac.F
• guide_mod.F90: ported updates
• leapfrog.F : (isotopes) updates (NB: call integrd with nqtot tracers)
• qminimium.F : adaptations for isotopes (copied over, except that #include comvert.h is not needed).
• vlsplt.F: adaptations for isotopes (copied over, except than #include logic.h, comvert.h not needed, and replace "include comconst.h" with use comconst_mod, ONLY: pi)
• vlspltqs.F : same as vlsplt.F, but also keeping added modification for CP(T)

• in dyn3dpar:
• leapfrog_p.F: remove unecessary #ifdef CPP_EARTH cpp flag. and call integrd_p with nqtot tracers (only important for Earth)
• dynredem_p.F => dynredem_p.F90 and some code factorization
• and no isotopes-relates changes in dyn3dpar (since these changes have been made in LMDZ5 dyn3dmem).

EM

File:

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

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

@@ +1 @@
+SUBROUTINE covcont (klevel,ucov, vcov, ucont, vcont )
 !
-! $Header$
-!
-      SUBROUTINE covcont (klevel,ucov, vcov, ucont, vcont )
-      IMPLICIT NONE
+!-------------------------------------------------------------------------------
+! Author: P. Le Van
+!-------------------------------------------------------------------------------
+! Purpose: Compute contravariant components from covariant components.
+!-------------------------------------------------------------------------------
+  IMPLICIT NONE
+  include "dimensions.h"
+  include "paramet.h"
+  include "comgeom.h"
+!===============================================================================
+! Arguments:
+  INTEGER, INTENT(IN)  :: klevel                    !--- VERTICAL LEVELS NUMBER
+  REAL,    INTENT(IN)  :: ucov ( ip1jmp1,klevel )   !--- U COVARIANT WIND
+  REAL,    INTENT(IN)  :: vcov ( ip1jm  ,klevel )   !--- V COVARIANT WIND
+  REAL,    INTENT(OUT) :: ucont( ip1jmp1,klevel )   !--- U CONTRAVAR WIND
+  REAL,    INTENT(OUT) :: vcont( ip1jm  ,klevel )   !--- V CONTRAVAR WIND
+!===============================================================================
+!   Local variables:
+  INTEGER :: l
+!===============================================================================
+  DO l=1,klevel
+    ucont(iip2:ip1jm,l)=ucov(iip2:ip1jm,l) * unscu2(iip2:ip1jm)
+    vcont(   1:ip1jm,l)=vcov(   1:ip1jm,l) * unscv2(   1:ip1jm)
+  END DO
 
-c=======================================================================
-c
-c   Auteur:  P. Le Van
-c   -------
-c
-c   Objet:
-c   ------
-c
-c  *********************************************************************
-c    calcul des compos. contravariantes a partir des comp.covariantes
-c  ********************************************************************
-c
-c=======================================================================
+END SUBROUTINE covcont
 
-#include "dimensions.h"
-#include "paramet.h"
-#include "comgeom.h"
-
-      INTEGER klevel
-      REAL ucov( ip1jmp1,klevel ),  vcov( ip1jm,klevel )
-      REAL ucont( ip1jmp1,klevel ), vcont( ip1jm,klevel )
-      INTEGER   l,ij
-
-
-      DO 10 l = 1,klevel
-
-      DO 2  ij = iip2, ip1jm
-      ucont( ij,l ) = ucov( ij,l ) * unscu2( ij )
-   2  CONTINUE
-
-      DO 4 ij = 1,ip1jm
-      vcont( ij,l ) = vcov( ij,l ) * unscv2( ij )
-   4  CONTINUE
-
-  10  CONTINUE
-      RETURN
-      END