/[lmdze]/trunk/libf/dyn3d/fxysinus.f
ViewVC logotype

Diff of /trunk/libf/dyn3d/fxysinus.f

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

revision 3 by guez, Wed Feb 27 13:16:39 2008 UTC revision 39 by guez, Tue Jan 25 15:11:05 2011 UTC
# Line 9  Line 9 
9        use dimens_m        use dimens_m
10        use paramet_m        use paramet_m
11        use comconst        use comconst
12          use nr_util, only: pi
13        IMPLICIT NONE        IMPLICIT NONE
14  c  c
15  c     Calcul  des longitudes et des latitudes  pour une fonction f(x,y)  c     Calcul  des longitudes et des latitudes  pour une fonction f(x,y)
# Line 25  c Line 26  c
26         REAL rlonu(iip1),xprimu(iip1),rlonv(iip1),xprimv(iip1),         REAL rlonu(iip1),xprimu(iip1),rlonv(iip1),xprimv(iip1),
27       , rlonm025(iip1),xprimm025(iip1), rlonp025(iip1),xprimp025(iip1)       , rlonm025(iip1),xprimm025(iip1), rlonp025(iip1),xprimp025(iip1)
28    
29         include "fxy_sin.h"  !
30    ! $Header: /home/cvsroot/LMDZ4/libf/grid/fxy_sin.h,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $
31    !
32    c-----------------------------------------------------------------------
33    c INCLUDE 'fxyprim.h'
34    c
35    c    ................................................................
36    c    ................  Fonctions in line  ...........................
37    c    ................................................................
38    c
39          REAL  fy, fx, fxprim, fyprim
40          REAL  ri, rj
41    c
42    c
43          fy(rj)=ASIN(1.+2.*((1.-rj)/FLOAT(jjm)))
44          fyprim(rj)=1./SQRT((rj-1.)*(jjm+1.-rj))
45    
46          fx    ( ri ) = 2.*pi/FLOAT(iim) * ( ri - 0.5*  FLOAT(iim) - 1. )
47    c     fx    ( ri ) = 2.*pi/FLOAT(iim) * ( ri - 0.5* ( FLOAT(iim) + 1.) )
48          fxprim( ri ) = 2.*pi/FLOAT(iim)
49    c
50    c
51    c    La valeur de pi est passee par le common/const/ou /const2/ .
52    c    Sinon, il faut la calculer avant d'appeler ces fonctions .
53    c
54    c   ----------------------------------------------------------------
55    c     Fonctions a changer eventuellement, selon x(x) et y(y) choisis .
56    c   -----------------------------------------------------------------
57    c
58    c    .....  ici, on a l'application particuliere suivante   ........
59    c
60    c                **************************************
61    c                **     x = 2. * pi/iim *  X         **
62    c                **     y =      pi/jjm *  Y         **
63    c                **************************************
64    c
65    c   ..................................................................
66    c   ..................................................................
67    c
68    c
69    c
70    c-----------------------------------------------------------------------
71    
72  c    ......  calcul  des  latitudes  et de y'   .....  c    ......  calcul  des  latitudes  et de y'   .....
73  c  c
Legend:
Removed from v.3  
changed lines
  Added in v.39
  ViewVC Help
Powered by ViewVC 1.1.21