Changeset 5105 for LMDZ6/branches/Amaury_dev/libf/grid/fxy_new.h

Timestamp:

Jul 23, 2024, 7:14:34 PM (8 weeks ago)

Author:

abarral

Message:

Replace 1DUTILS.h by module lmdz_1dutils.f90
Replace 1DConv.h by module lmdz_old_1dconv.f90 (it's only used by old_* files)
Convert *.F to *.f90
Fix gradsdef.h formatting
Remove unnecessary "RETURN" at the end of functions/subroutines

File:

• LMDZ6/branches/Amaury_dev/libf/grid/fxy_new.h (modified) (1 diff)

LMDZ6/branches/Amaury_dev/libf/grid/fxy_new.h

@@ -2 +2 @@
 ! $Header$
 
-c--------------------------------------------------------------
-         REAL ripx
-         REAL fx,fxprim,fy,fyprim,ri,rj,bigy
-c
-c....stretching in x...
-c
-        ripx(  ri )= (ri-1.0) *2.*pi/REAL(iim) 
-        fx  (  ri )= ripx(ri) + transx  +
-     *         alphax * SIN( ripx(ri)+transx-pxo ) - pi
-        fxprim(ri) = 2.*pi/REAL(iim)  *
-     *        ( 1.+ alphax * COS( ripx(ri)+transx-pxo ) )
+!--------------------------------------------------------------
+   REAL :: ripx
+   REAL :: fx,fxprim,fy,fyprim,ri,rj,bigy
+!
+!....stretching in x...
+!
+  ripx(  ri )= (ri-1.0) *2.*pi/REAL(iim)
+  fx  (  ri )= ripx(ri) + transx  + &
+        alphax * SIN( ripx(ri)+transx-pxo ) - pi
+  fxprim(ri) = 2.*pi/REAL(iim)  * &
+        ( 1.+ alphax * COS( ripx(ri)+transx-pxo ) )
 
-c....stretching in y...
-c
-        bigy(rj)   = 2.* (REAL(jjp1)-rj ) *pi/jjm
-        fy(rj)     =  ( bigy(rj) + transy  +
-     *        alphay * SIN( bigy(rj)+transy-pyo ) ) /2.  - pi/2.
-        fyprim(rj) = ( pi/jjm ) * ( 1.+
-     *           alphay * COS( bigy(rj)+transy-pyo ) )
+!....stretching in y...
+!
+  bigy(rj)   = 2.* (REAL(jjp1)-rj ) *pi/jjm
+  fy(rj)     =  ( bigy(rj) + transy  + &
+        alphay * SIN( bigy(rj)+transy-pyo ) ) /2.  - pi/2.
+  fyprim(rj) = ( pi/jjm ) * ( 1.+ &
+        alphay * COS( bigy(rj)+transy-pyo ) )
 
-c       fy(rj)= pyo-pisjjm*(rj-transy)+coefalpha*SIN(depisjm*(rj-
-c     *  transy ))
-c       fyprim(rj)= pisjjm-pisjjm*coefy2* COS(depisjm*(rj-transy)) 
-c--------------------------------------------------------------
+  ! fy(rj)= pyo-pisjjm*(rj-transy)+coefalpha*SIN(depisjm*(rj-
+! *  transy ))
+!   fyprim(rj)= pisjjm-pisjjm*coefy2* COS(depisjm*(rj-transy))
+!--------------------------------------------------------------