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

Diff of /trunk/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 52 by guez, Fri Sep 23 12:28:01 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
34    c    ................................................................
35    c    ................  Fonctions in line  ...........................
36    c    ................................................................
37    c
38          REAL  fy, fx, fxprim, fyprim
39          REAL  ri, rj
40    c
41    c
42          fy(rj)=ASIN(1.+2.*((1.-rj)/FLOAT(jjm)))
43          fyprim(rj)=1./SQRT((rj-1.)*(jjm+1.-rj))
44    
45          fx    ( ri ) = 2.*pi/FLOAT(iim) * ( ri - 0.5*  FLOAT(iim) - 1. )
46    c     fx    ( ri ) = 2.*pi/FLOAT(iim) * ( ri - 0.5* ( FLOAT(iim) + 1.) )
47          fxprim( ri ) = 2.*pi/FLOAT(iim)
48    c
49    c
50    c    La valeur de pi est passee par le common/const/ou /const2/ .
51    c    Sinon, il faut la calculer avant d'appeler ces fonctions .
52    c
53    c   ----------------------------------------------------------------
54    c     Fonctions a changer eventuellement, selon x(x) et y(y) choisis .
55    c   -----------------------------------------------------------------
56    c
57    c    .....  ici, on a l'application particuliere suivante   ........
58    c
59    c                **************************************
60    c                **     x = 2. * pi/iim *  X         **
61    c                **     y =      pi/jjm *  Y         **
62    c                **************************************
63    c
64    c   ..................................................................
65    c   ..................................................................
66    c
67    c
68    c
69    c-----------------------------------------------------------------------
70    
71  c    ......  calcul  des  latitudes  et de y'   .....  c    ......  calcul  des  latitudes  et de y'   .....
72  c  c
Legend:
Removed from v.3  
changed lines
  Added in v.52
  ViewVC Help
Powered by ViewVC 1.1.21