Changeset 5287 for LMDZ6/trunk/libf/grid/fxy_reg_mod_h.f90

Timestamp:

Oct 28, 2024, 2:53:53 PM (7 hours ago)

Author:

abarral

Message:

Turn fxy*.h into modules
Fix clesphys_rrtm.h name

File:

• LMDZ6/trunk/libf/grid/fxy_reg_mod_h.f90 (moved) (moved from LMDZ6/trunk/libf/grid/fxy_reg.h) (1 diff)

LMDZ6/trunk/libf/grid/fxy_reg_mod_h.f90

@@ -1 @@
-!
-! $Header$
-!
-!-----------------------------------------------------------------------
-! INCLUDE 'fxyprim.h'
-!
-!    ................................................................
-!    ................  Fonctions in line  ...........................
-!    ................................................................
-!
-REAL :: fy, fx, fxprim, fyprim
-REAL :: ri, rj
-!
-!
-fy    ( rj ) =    pi/REAL(jjm) * ( 0.5 * REAL(jjm) +  1. - rj  )
-fyprim( rj ) =    pi/REAL(jjm)
+! This module replaces grid/fxy_reg.h
 
-! fy(rj)=ASIN(1.+2.*((1.-rj)/REAL(jjm)))
-! fyprim(rj)=1./SQRT((rj-1.)*(jjm+1.-rj))
+MODULE fxy_reg_mod_h
+  USE comconst_mod, ONLY: pi
+  USE dimensions_mod, ONLY: iim, jjm, llm, ndm
+  IMPLICIT NONE; PRIVATE
+  PUBLIC fx, fxprim, fy, fyprim
 
-fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5*  REAL(iim) - 1. )
-! fx    ( ri ) = 2.*pi/REAL(iim) * ( ri - 0.5* ( REAL(iim) + 1.) )
-fxprim( ri ) = 2.*pi/REAL(iim)
-!
-!
-!    La valeur de pi est passee par le common/const/ou /const2/ .
-!    Sinon, il faut la calculer avant d'appeler ces fonctions .
-!
-!   ----------------------------------------------------------------
-! Fonctions a changer eventuellement, selon x(x) et y(y) choisis .
-!   -----------------------------------------------------------------
-!
-!    .....  ici, on a l'application particuliere suivante   ........
-!
-!            **************************************
-!            **     x = 2. * pi/iim *  X         **
-!            **     y =      pi/jjm *  Y         **
-!            **************************************
-!
-!   ..................................................................
-!   ..................................................................
-!
-!
-!
-!-----------------------------------------------------------------------
+CONTAINS
+
+  REAL FUNCTION fx(ri)
+    REAL, INTENT(IN) :: ri
+    fx = 2. * pi / REAL(iim) * (ri - 0.5 * REAL(iim) - 1.)
+  END FUNCTION fx
+
+  REAL FUNCTION fxprim(ri)
+    REAL, INTENT(IN) :: ri
+    fxprim = 2. * pi / REAL(iim)
+  END FUNCTION fxprim
+
+  REAL FUNCTION fy(rj)
+    REAL, INTENT(IN) :: rj
+    fy = pi / REAL(jjm) * (0.5 * REAL(jjm) + 1. - rj)
+  END FUNCTION fy
+
+  REAL FUNCTION fyprim(rj)
+    REAL, INTENT(IN) :: rj
+    fyprim = pi / REAL(jjm)
+  END FUNCTION fyprim
+END MODULE fxy_reg_mod_h