Changeset 3825 for dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/add_phys_tend.F90

Timestamp:

May 6, 2015, 12:14:12 PM (10 years ago)

Author:

ymipsl

Message:

Reorganize geometry and grid modules. Prepare physics for unstructutured grid support. Simplify initialization of physics from dynamic.
Compiled only with dynd3dmem, but not tested for moment.

YM

File:

• dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/add_phys_tend.F90 (modified) (5 diffs)

dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/add_phys_tend.F90

@@ -18 +18 @@
 use phys_state_var_mod
 USE inifis_mod, ONLY: prt_level
+USE geometry_mod, ONLY : lon_degrees, lat_degrees
 IMPLICIT none
   include "YOMCST.h"
@@ -106 +107 @@
          i=jadrs(j)
          if(prt_level.ge.debug_level) THEN
-          print*,'PLANTAGE POUR LE POINT i rlon rlat =',i,rlon(i),rlat(i),text
+          print*,'PLANTAGE POUR LE POINT i lon lat =',i,lon_degrees(i),lat_degrees(i),text
           print*,'l    T     dT       Q     dQ    '
           DO k = 1, klev
@@ -124 +125 @@
         i=jqadrs(j)
           if(prt_level.ge.debug_level) THEN
-           print*,'WARNING  : EAU POUR LE POINT i rlon rlat =',i,rlon(i),rlat(i),text
+           print*,'WARNING  : EAU POUR LE POINT i lon lat =',i,lon_degrees(i),lat_degrees(i),text
            print*,'l    T     dT       Q     dQ    '
            DO k = 1, klev
@@ -205 +206 @@
          k=kadrs(j)
          if(prt_level.ge.debug_level) THEN
-          print*,'PLANTAGE2 POUR LE POINT i itap rlon rlat txt jbad zdt t',i,itap,rlon(i),rlat(i),text,jbad, &
+          print*,'PLANTAGE2 POUR LE POINT i itap lon lat txt jbad zdt t',i,itap,lon_degrees(i),lat_degrees(i),text,jbad, &
        &        zdt(i,k),t_seri(i,k)-zdt(i,k)
 !!!       if(prt_level.ge.10.and.itap.GE.229.and.i.EQ.3027) THEN 
@@ -222 +223 @@
          k=kqadrs(j)
          if(prt_level.ge.debug_level) THEN
-          print*,'WARNING  : EAU2 POUR LE POINT i itap rlon rlat txt jqbad zdq q zdql ql',i,itap,rlon(i),rlat(i),text,jqbad,&
+          print*,'WARNING  : EAU2 POUR LE POINT i itap lon lat txt jqbad zdq q zdql ql',i,itap,lon_degrees(i), lat_degrees(i), text,jqbad,&
        &        zdq(i,k), q_seri(i,k)-zdq(i,k), zdql(i,k), ql_seri(i,k)-zdql(i,k)
 !!!       if(prt_level.ge.10.and.itap.GE.229.and.i.EQ.3027) THEN