Changeset 212 for LMDZ.3.3/trunk/libf/dyn3d/fxhyp.F

Timestamp:

Apr 24, 2001, 3:12:45 PM (23 years ago)

Author:

lmdz

Message:

Modifications des sorties de controle Levan
LF

File:

• LMDZ.3.3/trunk/libf/dyn3d/fxhyp.F (modified) (13 diffs)

LMDZ.3.3/trunk/libf/dyn3d/fxhyp.F

@@ -3 +3 @@
 c
        SUBROUTINE fxhyp ( xzoomdeg,grossism,dzoom,tau ,
-     , rlonm025,xprimm025,rlonv,xprimv,rlonu,xprimu,rlonp025,xprimp025)
+     , rlonm025,xprimm025,rlonv,xprimv,rlonu,xprimu,rlonp025,xprimp025,
+     , champmin,champmax                                               )
 
 c      Auteur :  P. Le Van 
@@ -49 +50 @@
        INTEGER is2
        SAVE is2
+
        REAL*8 heavyside
        EXTERNAL coefpoly,heavyside
@@ -57 +59 @@
        xzoom    = xzoomdeg * pi/180. 
 
-       WRITE(6,24) xzoomdeg,grossism,tau,dzoom
        IF( dzoom.LT.1.)  THEN
          dzoom = dzoom * depi
@@ -68 +69 @@
        ENDIF
 
+       WRITE(6,*) ' xzoom( rad.),grossism,tau,dzoom (radians)'
+       WRITE(6,24) xzoom,grossism,tau,dzoom
+
        DO i = 0, nmax2 
         xtild(i) = - pi + FLOAT(i) * depi /nmax2
@@ -117 +121 @@
 
         beta  = ( grossism * ffdx - pi ) / ( ffdx - pi )
-ccc        WRITE(6,*) ' ** X  beta **',beta,grossism
 
        IF( 2.*beta - grossism.LE. 0.)  THEN
@@ -178 +181 @@
 c     .....  xuv = 0.5  si  calcul  aux pts      U        ........
 c
+      WRITE(6,18)
 c
       DO 5000  ik = 1, 4
@@ -271 +275 @@
                 GO TO 1600
       ELSE
-       WRITE(6,18) 
        WRITE(6,*) 'Reorganisation des longitudes pour avoir entre - pi',
      ,  ' et pi '
@@ -331 +334 @@
 
          IF( ik.EQ.1 )  THEN
-          WRITE(6,*)  ' XLON aux pts. V-0.25   apres ( en  deg. ) '
-          WRITE(6,18) 
-          WRITE(6,68) xvrai
+c          WRITE(6,*)  ' XLON aux pts. V-0.25   apres ( en  deg. ) '
+c          WRITE(6,18) 
+c          WRITE(6,68) xvrai
 ccc          WRITE(6,*) ' XPRIM k ',ik
 ccc          WRITE(6,566)  xprimm
@@ -343 +346 @@
 
          ELSE IF( ik.EQ.2 )  THEN
-          WRITE(6,18) 
-          WRITE(6,*)  ' XLON aux pts. V   apres ( en  deg. ) '
-          WRITE(6,68) xvrai
+c          WRITE(6,18) 
+c          WRITE(6,*)  ' XLON aux pts. V   apres ( en  deg. ) '
+c          WRITE(6,68) xvrai
 ccc          WRITE(6,*) ' XPRIM k ',ik
 ccc          WRITE(6,566)  xprimm
@@ -355 +358 @@
 
          ELSE IF( ik.EQ.3)   THEN
-          WRITE(6,18) 
-          WRITE(6,*)  ' XLON aux pts. U   apres ( en  deg. ) '
-          WRITE(6,68) xvrai
+c          WRITE(6,18) 
+c          WRITE(6,*)  ' XLON aux pts. U   apres ( en  deg. ) '
+c          WRITE(6,68) xvrai
 ccc          WRITE(6,*) ' XPRIM ik ',ik
 ccc          WRITE(6,566)  xprimm
@@ -367 +370 @@
 
          ELSE IF( ik.EQ.4 )  THEN
-          WRITE(6,18) 
-          WRITE(6,*)  ' XLON aux pts. V+0.25   apres ( en  deg. ) '
-          WRITE(6,68) xvrai
+c          WRITE(6,18) 
+c          WRITE(6,*)  ' XLON aux pts. V+0.25   apres ( en  deg. ) '
+c          WRITE(6,68) xvrai
 ccc          WRITE(6,*) ' XPRIM ik ',ik
 ccc          WRITE(6,566)  xprimm
@@ -381 +384 @@
 
 5000    CONTINUE
+c
+       WRITE(6,18)
 c
 c    ...........  fin  de la boucle  do 5000      ............
@@ -396 +401 @@
          champmax = champmax * 180./pi
 
-        WRITE(6,18)
-        WRITE(6,*) '  Longitudes '
-        WRITE(6,18)
-        WRITE(6,3)  champmin, champmax
-        WRITE(6,*) ' Si cette derniere est trop lache , modifiez les par
-     ,ametres  grossism , tau , dzoom pour X et repasser ! '
-        WRITE(6,18)
-
-3      Format(1x, ' Au centre du zoom , la longueur de la maille est',
-     ,  ' d environ ',f8.2 ,' degres  ',
-     , ' alors que la maille en dehors de la zone du zoom est d environ
-     , ', f8.2,' degres ' )
 18     FORMAT(/)
-24     FORMAT(2x,' Parametres yzoom,gross,tau ,dzoom pour fyhyp ',4f8.3)
+24     FORMAT(2x,'Parametres xzoom,gross,tau ,dzoom pour fxhyp ',4f8.3)
 68     FORMAT(1x,7f9.2)