Changeset 941

Timestamp:

Apr 9, 2008, 12:36:53 PM (17 years ago)

Author:

lmdzadmin

Message:

Ajout coordonnees geographiques pour les plantages
IM

Location:

LMDZ4/trunk/libf/phylmd

Files:

• add_phys_tend.F90 (modified) (3 diffs)
• hgardfou.F (modified) (6 diffs)
• print_debug_phys.F90 (modified) (2 diffs)

LMDZ4/trunk/libf/phylmd/add_phys_tend.F90

@@ -16 +16 @@
 use dimphy
 use phys_local_var_mod
+use phys_state_var_mod
 IMPLICIT none
 
@@ -78 +79 @@
       DO j = 1, jbad
          i=jadrs(j)
-         print*,'PLANTAGE POUR LE POINT i=',i,text
+         print*,'PLANTAGE POUR LE POINT i rlon rlat =',i,rlon(i),rlat(i),text
          print*,'l    T     dT       Q     dQ    '
          DO k = 1, klev
@@ -93 +94 @@
       DO j = 1, jqbad
          i=jqadrs(j)
-         print*,'WARNING  : EAU POUR LE POINT i=',i,text
+         print*,'WARNING  : EAU POUR LE POINT i rlon rlat =',i,rlon(i),rlat(i),text
          print*,'l    T     dT       Q     dQ    '
          DO k = 1, klev

LMDZ4/trunk/libf/phylmd/hgardfou.F

@@ -4 +4 @@
       SUBROUTINE hgardfou (t,tsol,text)
       use dimphy
+      use phys_state_var_mod
       IMPLICIT none
 c======================================================================
 c Verifier la temperature
 c======================================================================
-cym#include "dimensions.h"
-cym#include "dimphy.h"
+#include "dimensions.h"
 #include "YOMCST.h"
 #include "indicesol.h"
@@ -26 +26 @@
          PRINT*, 'hgardfou garantit la temperature dans [100,370] K'
          firstcall = .FALSE.
+c        DO i = 1, klon
+c         print*,'i=',i,'rlon=',rlon(i),'rlat=',rlat(i)
+c        ENDDO
+c
       ENDIF
 c
@@ -47 +51 @@
            ok = .FALSE.
            DO i = 1, jbad
-             PRINT *,'i,k,temperature =',jadrs(i),k,zt(jadrs(i))
+             PRINT *,'i,k,temperature =',jadrs(i),k,zt(jadrs(i)),
+     $       rlon(jadrs(i)),rlat(jadrs(i))
            ENDDO
          ENDIF
@@ -65 +70 @@
            ok = .FALSE.
            DO i = 1, jbad
-             PRINT *,'i,k,temperature =',jadrs(i),k,zt(jadrs(i))
+             PRINT *,'i,k,temperature =',jadrs(i),k,zt(jadrs(i)),
+     $       rlon(jadrs(i)),rlat(jadrs(i))
            ENDDO
          ENDIF
@@ -88 +94 @@
            ok = .FALSE.
            DO i = 1, jbad
-             PRINT *,'i,nsrf,temperature =',jadrs(i),nsrf,zt(jadrs(i))
+             PRINT *,'i,nsrf,temperature =',jadrs(i),nsrf,zt(jadrs(i)),
+     $       rlon(jadrs(i)),rlat(jadrs(i))
            ENDDO
          ENDIF
@@ -106 +113 @@
            ok = .FALSE.
            DO i = 1, jbad
-             PRINT *,'i,nsrf,temperature =',jadrs(i),nsrf,zt(jadrs(i))
+             PRINT *,'i,nsrf,temperature =',jadrs(i),nsrf,zt(jadrs(i)),
+     $       rlon(jadrs(i)),rlat(jadrs(i))
            ENDDO
          ENDIF

LMDZ4/trunk/libf/phylmd/print_debug_phys.F90

@@ -3 +3 @@
 use dimphy
 use phys_local_var_mod
+use phys_state_var_mod
 IMPLICIT NONE
 integer i,debug_lev
@@ -13 +14 @@
 print*,'l    u, v, T, q, ql'
 DO k = 1, klev
-   write(*,'(i3,3f14.4,2e14.2)') k,u_seri(i,k),v_seri(i,k),t_seri(i,k),q_seri(i,k),ql_seri(i,k)
+   write(*,'(i3,2f8.4,3f14.4,2e14.2)') k,rlon(i),rlat(i),u_seri(i,k),v_seri(i,k),t_seri(i,k),q_seri(i,k),ql_seri(i,k)
 ENDDO