Changeset 5128 for LMDZ6/branches/Amaury_dev/libf/phylmd/FCTTRE.h

Timestamp:

Jul 25, 2024, 5:47:25 PM (4 months ago)

Author:

abarral

Message:

Correct bug in vlspltqs_loc.f90 from r2270 where we call SSUM with incorrect arguments.
Merge the three different versions of abort_gcm into one
Fix seq, para 3D compilation broken from r5107 onwards
(lint) usual + Remove uneeded fixed-form continuations

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/FCTTRE.h (modified) (1 diff)

LMDZ6/branches/Amaury_dev/libf/phylmd/FCTTRE.h

@@ -22 +22 @@
 
       FOEEW ( PTARG,PDELARG ) = EXP (                                   &
-     &          (R3LES*(1.-PDELARG)+R3IES*PDELARG) * (PTARG-RTT)        &
-     & / (PTARG-(R4LES*(1.-PDELARG)+R4IES*PDELARG)) )
+                (R3LES*(1.-PDELARG)+R3IES*PDELARG) * (PTARG-RTT)        &
+       / (PTARG-(R4LES*(1.-PDELARG)+R4IES*PDELARG)) )
 
       FOEDE ( PTARG,PDELARG,P5ARG,PQSARG,PCOARG ) = PQSARG*PCOARG*P5ARG &
-     & / (PTARG-(R4LES*(1.-PDELARG)+R4IES*PDELARG))**2
+       / (PTARG-(R4LES*(1.-PDELARG)+R4IES*PDELARG))**2
 
       qsats(ptarg) = 100.0 * 0.622 * 10.0                               &
-     &           ** (2.07023 - 0.00320991 * ptarg                       &
-     &           - 2484.896 / ptarg + 3.56654 * LOG10(ptarg))
+                 ** (2.07023 - 0.00320991 * ptarg                       &
+                 - 2484.896 / ptarg + 3.56654 * LOG10(ptarg))
       qsatl(ptarg) = 100.0 * 0.622 * 10.0                               &
-     &           ** (23.8319 - 2948.964 / ptarg                         &
-     &           - 5.028 * LOG10(ptarg)                                 &
-     &           - 29810.16 * EXP( - 0.0699382 * ptarg)                 &
-     &           + 25.21935 * EXP( - 2999.924 / ptarg))
+                 ** (23.8319 - 2948.964 / ptarg                         &
+                 - 5.028 * LOG10(ptarg)                                 &
+                 - 29810.16 * EXP( - 0.0699382 * ptarg)                 &
+                 + 25.21935 * EXP( - 2999.924 / ptarg))
 
       dqsats(ptarg,pqsarg) = RLVTT/RCPD*pqsarg * (3.56654/ptarg         &
-     &                     +2484.896*LOG(10.)/ptarg**2                  &
-     &                     -0.00320991*LOG(10.))
+                           +2484.896*LOG(10.)/ptarg**2                  &
+                           -0.00320991*LOG(10.))
       dqsatl(ptarg,pqsarg) = RLVTT/RCPD*pqsarg*LOG(10.)*                &
-     &                (2948.964/ptarg**2-5.028/LOG(10.)/ptarg           &
-     &                +25.21935*2999.924/ptarg**2*EXP(-2999.924/ptarg)  &
-     &                +29810.16*0.0699382*EXP(-0.0699382*ptarg))
+                      (2948.964/ptarg**2-5.028/LOG(10.)/ptarg           &
+                      +25.21935*2999.924/ptarg**2*EXP(-2999.924/ptarg)  &
+                      +29810.16*0.0699382*EXP(-0.0699382*ptarg))