/[lmdze]/trunk/dyn3d/comgeom.f
ViewVC logotype

Diff of /trunk/dyn3d/comgeom.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 112 by guez, Thu Sep 18 13:36:51 2014 UTC revision 113 by guez, Thu Sep 18 19:56:46 2014 UTC
# Line 163  contains Line 163  contains
163    
164      USE comconst, ONLY : g, omeg, rad      USE comconst, ONLY : g, omeg, rad
165      USE comdissnew, ONLY : coefdis, nitergdiv, nitergrot, niterh      USE comdissnew, ONLY : coefdis, nitergdiv, nitergrot, niterh
     use conf_gcm_m, ONLY : fxyhypb, ysinus  
     use fxy_m, only: fxy  
166      use fxyhyper_m, only: fxyhyper      use fxyhyper_m, only: fxyhyper
     use fxysinus_m, only: fxysinus  
167      use jumble, only: new_unit      use jumble, only: new_unit
168      use nr_util, only: pi      use nr_util, only: pi
169      USE paramet_m, ONLY : iip1, jjp1      USE paramet_m, ONLY : iip1, jjp1
     USE serre, ONLY : alphax, alphay, clat, clon, pxo, pyo, transx, transy  
     ! Modifiés pxo, pyo, transx, transy  
170    
171      ! Local:      ! Local:
172      INTEGER i, j, itmax, itmay, iter, unit      INTEGER i, j, unit
173      REAL cvu(iip1, jjp1), cuv(iip1, jjm)      REAL cvu(iip1, jjp1), cuv(iip1, jjm)
174      REAL ai14, ai23, airez, un4rad2      REAL ai14, ai23, airez, un4rad2
     REAL eps, x1, xo1, f, df, xdm, y1, yo1, ydm  
175      REAL coslatm, coslatp, radclatm, radclatp      REAL coslatm, coslatp, radclatm, radclatp
176      REAL, dimension(iip1, jjp1):: cuij1, cuij2, cuij3, cuij4 ! in m      REAL, dimension(iip1, jjp1):: cuij1, cuij2, cuij3, cuij4 ! in m
177      REAL, dimension(iip1, jjp1):: cvij1, cvij2, cvij3, cvij4 ! in m      REAL, dimension(iip1, jjp1):: cvij1, cvij2, cvij3, cvij4 ! in m
# Line 215  contains Line 209  contains
209      print *, "gamdi_grot = ", gamdi_grot      print *, "gamdi_grot = ", gamdi_grot
210      print *, "gamdi_h = ", gamdi_h      print *, "gamdi_h = ", gamdi_h
211    
212      IF (fxyhypb) THEN      print *, 'inigeom: Y = latitude, dérivée tangente hyperbolique'
213         print *, 'inigeom: Y = latitude, dérivée tangente hyperbolique'      CALL fxyhyper(rlatu, yprimu, rlatv, yprimv, rlatu1, yprimu1, rlatu2, &
214         CALL fxyhyper(rlatu, yprimu, rlatv, yprimv, rlatu1, yprimu1, rlatu2, &           yprimu2, rlonu, xprimu, rlonv, xprimv, rlonm025, xprimm025, &
215              yprimu2, rlonu, xprimu, rlonv, xprimv, rlonm025, xprimm025, &           rlonp025, xprimp025)
             rlonp025, xprimp025)  
     ELSE  
        IF (ysinus) THEN  
           print *, 'inigeom: Y = sin(latitude)'  
           ! Utilisation de f(x, y) avec y = sinus de la latitude  
           CALL fxysinus(rlatu, yprimu, rlatv, yprimv, rlatu1, yprimu1, &  
                rlatu2, yprimu2, rlonu, xprimu, rlonv, xprimv, rlonm025, &  
                xprimm025, rlonp025, xprimp025)  
        ELSE  
           print *, 'Inigeom, Y = Latitude, der. sinusoid .'  
           ! utilisation de f(x, y) a tangente sinusoidale, y etant la latit  
   
           pxo = clon * pi / 180.  
           pyo = 2. * clat * pi / 180.  
   
           ! determination de transx (pour le zoom) par Newton-Raphson  
   
           itmax = 10  
           eps = .1E-7  
   
           xo1 = 0.  
           DO iter = 1, itmax  
              x1 = xo1  
              f = x1 + alphax * sin(x1-pxo)  
              df = 1. + alphax * cos(x1-pxo)  
              x1 = x1 - f / df  
              xdm = abs(x1-xo1)  
              IF (xdm<=eps) EXIT  
              xo1 = x1  
           END DO  
   
           transx = xo1  
   
           itmay = 10  
           eps = .1E-7  
   
           yo1 = 0.  
           DO iter = 1, itmay  
              y1 = yo1  
              f = y1 + alphay * sin(y1-pyo)  
              df = 1. + alphay * cos(y1-pyo)  
              y1 = y1 - f / df  
              ydm = abs(y1-yo1)  
              IF (ydm<=eps) EXIT  
              yo1 = y1  
           END DO  
   
           transy = yo1  
   
           CALL fxy(rlatu, yprimu, rlatv, yprimv, rlatu1, yprimu1, rlatu2, &  
                yprimu2, rlonu, xprimu, rlonv, xprimv, rlonm025, xprimm025, &  
                rlonp025, xprimp025)  
        END IF  
     END IF  
216    
217      rlatu(1) = pi / 2.      rlatu(1) = pi / 2.
218      rlatu(jjp1) = -rlatu(1)      rlatu(jjp1) = -rlatu(1)
Legend:
Removed from v.112  
changed lines
  Added in v.113
  ViewVC Help
Powered by ViewVC 1.1.21