source: LMDZ6/trunk/libf/grid/fxy_reg.h @ 5248

Last change on this file since 5248 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.5 KB
RevLine 
[524]1!
2! $Header$
3!
[5246]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 ) =    pi/REAL(jjm) * ( 0.5 * REAL(jjm) +  1. - rj  )
16fyprim( rj ) =    pi/REAL(jjm)
[524]17
[5246]18! fy(rj)=ASIN(1.+2.*((1.-rj)/REAL(jjm)))
19! fyprim(rj)=1./SQRT((rj-1.)*(jjm+1.-rj))
[524]20
[5246]21fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5*  REAL(iim) - 1. )
22! fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5* ( REAL(iim) + 1.) )
23fxprim( ri ) = 2.*pi/REAL(iim)
24!
25!
26!    La valeur de pi est passee par le common/const/ou /const2/ .
27!    Sinon, il faut la calculer avant d'appeler ces fonctions .
28!
29!   ----------------------------------------------------------------
30! Fonctions a changer eventuellement, selon x(x) et y(y) choisis .
31!   -----------------------------------------------------------------
32!
33!    .....  ici, on a l'application particuliere suivante   ........
34!
35!            **************************************
36!            **     x = 2. * pi/iim *  X         **
37!            **     y =      pi/jjm *  Y         **
38!            **************************************
39!
40!   ..................................................................
41!   ..................................................................
42!
43!
44!
45!-----------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.