source: LMDZ6/trunk/libf/grid/fxy_sin.h @ 5246

Last change on this file since 5246 was 5246, checked in by abarral, 23 hours ago

Convert fixed-form to free-form sources .F -> .{f,F}90
(WIP: some .F remain, will be handled in subsequent commits)
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.4 KB
Line 
1!
2! $Header$
3!
4!-----------------------------------------------------------------------
5! INCLUDE 'fxyprim.h'
6!
7!    ................................................................
8!    ................  Fonctions in line  ...........................
9!    ................................................................
10!
11REAL :: fy, fx, fxprim, fyprim
12REAL :: ri, rj
13!
14!
15fy(rj)=ASIN(1.+2.*((1.-rj)/REAL(jjm)))
16fyprim(rj)=1./SQRT((rj-1.)*(jjm+1.-rj))
17
18fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5*  REAL(iim) - 1. )
19! fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5* ( REAL(iim) + 1.) )
20fxprim( ri ) = 2.*pi/REAL(iim)
21!
22!
23!    La valeur de pi est passee par le common/const/ou /const2/ .
24!    Sinon, il faut la calculer avant d'appeler ces fonctions .
25!
26!   ----------------------------------------------------------------
27! Fonctions a changer eventuellement, selon x(x) et y(y) choisis .
28!   -----------------------------------------------------------------
29!
30!    .....  ici, on a l'application particuliere suivante   ........
31!
32!            **************************************
33!            **     x = 2. * pi/iim *  X         **
34!            **     y =      pi/jjm *  Y         **
35!            **************************************
36!
37!   ..................................................................
38!   ..................................................................
39!
40!
41!
42!-----------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.