Changeset 2382

Timestamp:

Nov 2, 2015, 2:52:09 PM (9 years ago)

Author:

jbmadeleine

Message:

Updated phymar routines and made sure it compiles

• Changed X_AT in Mod_PHY_AT_kkl.f90 from REAL*16 (obsolete) to double precision
• The number of tracers nqtot is now given to MAR via the infotrac_phy routine
• Changed iim and jjm into ii and jj in dynlonlat_phylonlat/phymar/iniphysiq_mod.F90
• Added the declaration of the rotational in phymar/physiq.F90

Location:

LMDZ5/trunk/libf/phymar

Files:

• Mod_PHY_AT_kkl.f90 (modified) (1 diff)
• infotrac_phy.F90 (added)
• physiq.F90 (modified) (1 diff)

LMDZ5/trunk/libf/phymar/Mod_PHY_AT_kkl.f90

@@ -88 +88 @@
       real(kind=real8), SAVE,ALLOCATABLE ,dimension(:)    ::  P___AT  ! Auxiliary   Term               P                           [-]
       real(kind=real8), SAVE,ALLOCATABLE ,dimension(:)    ::  Q___AT  ! Auxiliary   Term               Q                           [-]
-      real*16         ,ALLOCATABLE ,dimension(:)    ::  X___AT  ! Auxiliary   Unknown            X                           [x]
+!      real*16         ,ALLOCATABLE ,dimension(:)    ::  X___AT  ! Auxiliary   Unknown            X                           [x]
+      double precision,ALLOCATABLE ,dimension(:)    ::  X___AT  ! Auxiliary   Unknown            X                           [x]
 
       real(kind=real8), SAVE,ALLOCATABLE ,dimension(:)    ::  LMO_AT  ! Monin-Obukhov     Length         (Grid Cell Average)       [m]

LMDZ5/trunk/libf/phymar/physiq.F90

@@ -120 +120 @@
       real,intent(in) :: u(klon,klev) ! eastward zonal wind (m/s)
       real,intent(in) :: v(klon,klev) ! northward meridional wind (m/s)
+      REAL, intent(in):: rot(klon, klev) ! relative vorticity,
+                                         ! in s-1, needed for frontal waves
       real,intent(in) :: t(klon,klev) ! temperature (K)
       real,intent(in) :: qx(klon,klev,nqtot) ! tracers (.../kg_air)