Changeset 1146 for LMDZ4/trunk/libf/dyn3dpar/dynredem.F

Timestamp:

Apr 9, 2009, 12:11:35 PM (15 years ago)

Author:

Laurent Fairhead

Message:

Réintegration dans le tronc des modifications issues de la branche LMDZ-dev
comprises entre la révision 1074 et 1145
Validation: une simulation de 1 jour en séquentiel sur PC donne les mêmes
résultats entre la trunk et la dev
LF

Location:

LMDZ4/trunk

Files:

• . (modified) (1 prop)
• libf/dyn3dpar/dynredem.F (modified) (9 diffs)

LMDZ4/trunk/libf/dyn3dpar/dynredem.F

@@ -3 +3 @@
 !
 c
-      SUBROUTINE dynredem0(fichnom,iday_end,phis,nq)
+      SUBROUTINE dynredem0(fichnom,iday_end,phis)
       USE IOIPSL
+      USE infotrac
       IMPLICIT NONE
 c=======================================================================
@@ -22 +23 @@
 #include "description.h"
 #include "serre.h"
-#include "advtrac.h"
 
 c   Arguments:
@@ -29 +29 @@
       REAL phis(ip1jmp1)
       CHARACTER*(*) fichnom
-      INTEGER nq
 
 c   Local:
@@ -458 +457 @@
       dims4(3) = idim_s
       dims4(4) = idim_tim
-      IF(nq.GE.1) THEN
-      DO iq=1,nq
+
+      DO iq=1,nqtot
 cIM 220306 BEG
 #ifdef NC_DOUBLE
@@ -469 +468 @@
       ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 12,ttext(iq))
       ENDDO
-      ENDIF
 c
       dims4(1) = idim_rlonv
@@ -508 +506 @@
       END
       SUBROUTINE dynredem1(fichnom,time,
-     .                     vcov,ucov,teta,q,nq,masse,ps)
+     .                     vcov,ucov,teta,q,masse,ps)
+      USE infotrac
       IMPLICIT NONE
 c=================================================================
@@ -519 +518 @@
 #include "comvert.h"
 #include "comgeom.h"
-#include "advtrac.h"
 #include "temps.h"
 #include "control.h"
 
-      INTEGER nq, l
+      INTEGER l
       REAL vcov(ip1jm,llm),ucov(ip1jmp1,llm) 
       REAL teta(ip1jmp1,llm)                   
       REAL ps(ip1jmp1),masse(ip1jmp1,llm)                   
-      REAL q(ip1jmp1,llm,nq)
+      REAL q(ip1jmp1,llm,nqtot)
       CHARACTER*(*) fichnom
      
@@ -633 +631 @@
       END IF
 
-      IF(nq.GE.1) THEN
-      do iq=1,nq 
+      do iq=1,nqtot 
 
          IF (config_inca == 'none') THEN
@@ -704 +701 @@
       
       ENDDO
-      ENDIF
 c
       ierr = NF_INQ_VARID(nid, "masse", nvarid)