Ignore:
Timestamp:
Apr 12, 2001, 11:30:31 AM (23 years ago)
Author:
lmdz
Message:

Mise en oeuvre du nouveau zoom. Modifs P. Levan
LF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ.3.3/trunk/libf/dyn3d/fxyhyper.F

    r2 r203  
    1        SUBROUTINE fxyhyper ( yzoom, grossy, dzoomy,tauy,deltay  ,
    2      ,                      xzoom, grossx, dzoomx,taux,
     1c
     2c $Header
     3c
     4       SUBROUTINE fxyhyper ( yzoom, grossy, dzoomy,tauy  ,   
     5     ,                       xzoom, grossx, dzoomx,taux  ,
    36     , rlatu,yprimu,rlatv,yprimv,rlatu1,  yprimu1,  rlatu2,  yprimu2  ,
    47     , rlonu,xprimu,rlonv,xprimv,rlonm025,xprimm025,rlonp025,xprimp025)
     
    1922c     a) le grossissement du zoom  :  grossy  ( en y ) et grossx ( en x )
    2023c     b) l' extension     du zoom  :  dzoomy  ( en y ) et dzoomx ( en x )
    21 c     c) la raideur       du zoom  :   tau   , ici  =  1.
     24c     c) la raideur de la transition du zoom  :   taux et tauy   
    2225c
    23 c     N.B :   le produit ( grossissement * extension )  doit etre <  1.
    24 c    *******
    25 c     En plus , il y a un autre parametre , moins important mais quand
    26 c     meme utile , c'est  deltay , deplacement de la zone du zoom en
    27 c     latitude   :  Si  on deplace  de 10. degres vers le nord  ( deltay
    28 c      = 10.  dans inigeom ) .
    29 c
     26c  N.B : Il vaut mieux avoir   :   grossx * dzoomx <  pi    ( radians )
     27c ******
     28c                  et              grossy * dzoomy <  pi/2  ( radians )
    3029c
    3130#include "dimensions.h"
     
    4039       REAL rlonu(iip1),xprimu(iip1),rlonv(iip1),xprimv(iip1),
    4140     , rlonm025(iip1),xprimm025(iip1), rlonp025(iip1),xprimp025(iip1)
    42        REAL deltay
    4341
    4442c   ....   var. locales   .....
     
    4745c
    4846
    49        CALL fyhyp ( yzoom, grossy, dzoomy,tauy, deltay ,
     47       CALL fyhyp ( yzoom, grossy, dzoomy,tauy  ,
    5048     ,  rlatu, yprimu,rlatv,yprimv,rlatu2,yprimu2,rlatu1,yprimu1 )
    51 
    5249
    5350       CALL fxhyp(xzoom,grossx,dzoomx,taux,rlonm025,xprimm025,rlonv,
     
    5552
    5653
    57         DO i = 1, iim
     54        DO i = 1, iip1
    5855          IF(rlonp025(i).LT.rlonv(i))  THEN
    59            PRINT *,' Attention !  rlonp025 < rlonv',i
     56           WRITE(6,*) ' Attention !  rlonp025 < rlonv',i
    6057            STOP
    6158          ENDIF
    6259
    6360          IF(rlonv(i).LT.rlonm025(i))  THEN
    64            PRINT *,' Attention !  rlonm025 > rlonv',i
     61           WRITE(6,*) ' Attention !  rlonm025 > rlonv',i
    6562            STOP
    6663          ENDIF
    6764
    6865          IF(rlonp025(i).GT.rlonu(i))  THEN
    69            PRINT *,' Attention !  rlonp025 > rlonu',i
     66           WRITE(6,*) ' Attention !  rlonp025 > rlonu',i
    7067            STOP
    7168          ENDIF
    7269        ENDDO
    7370
    74         PRINT *,'  *** TEST DE COHERENCE  OK    POUR   FX **** '
     71        WRITE(6,*) '  *** TEST DE COHERENCE  OK    POUR   FX **** '
    7572
    7673c
     
    7875c
    7976       IF(rlatu1(j).LE.rlatu2(j))   THEN
    80          PRINT *,' Attention ! rlatu1 < rlatu2 ',rlatu1(j), rlatu2(j),j
     77         WRITE(6,*)'Attention ! rlatu1 < rlatu2 ',rlatu1(j), rlatu2(j),j
    8178         STOP 13
    8279       ENDIF
    8380c
    8481       IF(rlatu2(j).LE.rlatu(j+1))  THEN
    85         PRINT *,' Attention ! rlatu2 < rlatup1 ',rlatu2(j),rlatu(j+1),j
     82        WRITE(6,*)'Attention ! rlatu2 < rlatup1 ',rlatu2(j),rlatu(j+1),j
    8683        STOP 14
    8784       ENDIF
    8885c
    8986       IF(rlatu(j).LE.rlatu1(j))    THEN
    90         PRINT *,' Attention ! rlatu < rlatu1 ',rlatu(j),rlatu1(j),j
     87        WRITE(6,*)' Attention ! rlatu < rlatu1 ',rlatu(j),rlatu1(j),j
    9188        STOP 15
    9289       ENDIF
    9390c
    9491       IF(rlatv(j).LE.rlatu2(j))    THEN
    95         PRINT *,' Attention ! rlatv < rlatu2 ',rlatv(j),rlatu2(j),j
     92        WRITE(6,*)' Attention ! rlatv < rlatu2 ',rlatv(j),rlatu2(j),j
    9693        STOP 16
    9794       ENDIF
    9895c
    9996       IF(rlatv(j).ge.rlatu1(j))    THEN
    100         PRINT *,' Attention ! rlatv > rlatu1 ',rlatv(j),rlatu1(j),j
     97        WRITE(6,*)' Attention ! rlatv > rlatu1 ',rlatv(j),rlatu1(j),j
    10198        STOP 17
    10299       ENDIF
    103100c
    104101       IF(rlatv(j).ge.rlatu(j))     THEN
    105         PRINT *,'Attention ! rlatv > rlatu ',rlatv(j),rlatu(j),j
     102        WRITE(6,*) ' Attention ! rlatv > rlatu ',rlatv(j),rlatu(j),j
    106103        STOP 18
    107104       ENDIF
     
    109106       ENDDO
    110107c
    111        PRINT *,'  *** TEST DE COHERENCE  OK    POUR   FY **** '
    112 
     108       WRITE(6,*) '  *** TEST DE COHERENCE  OK    POUR   FY **** '
     109       WRITE(6,25)
     11025     FORMAT(//)
    113111c
    114112
Note: See TracChangeset for help on using the changeset viewer.