Changeset 5103 for LMDZ6/branches/Amaury_dev/libf/phylmd/stratosphere_mask.F90

Timestamp:

Jul 23, 2024, 3:29:36 PM (8 weeks ago)

Author:

abarral

Message:

Handle CPP_INLANDSIS in lmdz_cppkeys_wrapper.F90
Remove obsolete key wrgrads_thermcell, _ADV_HALO, _ADV_HALLO, isminmax
Remove redundant uses of CPPKEY_INCA (thanks acozic)
Remove obsolete misc/write_field.F90
Remove unused ioipsl_* wrappers
Remove calls to WriteField_u with wrong signature
Convert .F -> .[fF]90
(lint) uppercase fortran operators
[note: 1d and iso still broken - working on it]

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/stratosphere_mask.F90 (modified) (8 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/stratosphere_mask.F90

@@ -50 +50 @@
 REAL, PARAMETER                        :: pliml=7500.
 REAL, PARAMETER                        :: gamma=-0.002
-LOGICAL, PARAMETER                     :: dofill=.true.
+LOGICAL, PARAMETER                     :: dofill=.TRUE.
 REAL,DIMENSION(klon)                   :: tp
 REAL,DIMENSION(klev)                   :: t, p
@@ -66 +66 @@
   psrf=pplay(i,1)
   zsrf=pphis(i)/RG           !--altitude de la surface
-  call twmo(missing_val, klev, t, p, psrf, zsrf, plimu, pliml, gamma, ptrp, ttrp, ztrp)
+  CALL twmo(missing_val, klev, t, p, psrf, zsrf, plimu, pliml, gamma, ptrp, ttrp, ztrp)
   tp(i)=ptrp
   p_tropopause(i)=ptrp
@@ -107 +107 @@
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
-subroutine twmo(missing_val, level, t, p, ps, zs, plimu, pliml, gamma, ptrp, ttrp, ztrp)
+SUBROUTINE twmo(missing_val, level, t, p, ps, zs, plimu, pliml, gamma, ptrp, ttrp, ztrp)
 
 ! reference: Reichler, T., M. Dameris, and R. Sausen (GRL, 10.1029/2003GL018240, 2003)
@@ -137 +137 @@
 faktor = -RG/RD
 
-do j=level,2,-1
+DO j=level,2,-1
 
    ! dt/dz
@@ -206 +206 @@
 
 jj = LEVEL                                       ! bottom
-do while ((P(jj)>PS) .or. (T(jj)<100))     ! T must be valid too
+DO while ((P(jj)>PS) .or. (T(jj)<100))     ! T must be valid too
   jj=jj-1
-enddo
+END DO
 
 DLNP = log(PS/P(jj))                             ! from surface pressure
@@ -214 +214 @@
 TDLNP = TM*DLNP
 
-do while ( (JJ>=2) .and. (PTRP<P(jj-1)) )
+DO while ( (JJ>=2) .and. (PTRP<P(jj-1)) )
   DLNP = log(P(jj)/P(jj-1))
   TM = 0.5 * (T(jj) + T(jj-1))
   TDLNP = TDLNP + TM*DLNP
   JJ=JJ-1
-enddo
+END DO
 
 DLNP = log(P(jj)/PTRP)                           ! up to tropopause pressure
@@ -228 +228 @@
 
 !!if (ZTRP .lt. 0) then
-!!  print*,'ZTRP=',ZTRP
-!!  print*,'PS=',PS
-!!  print*,'P=',P
-!!  print*,'T=',T
-!!  print*,'ZS=',ZS
+!!  PRINT*,'ZTRP=',ZTRP
+!!  PRINT*,'PS=',PS
+!!  PRINT*,'P=',P
+!!  PRINT*,'T=',T
+!!  PRINT*,'ZS=',ZS
 !!  stop
 !!endif
@@ -246 +246 @@
     b0 = b
 
-enddo 
+END DO
 
 ! no tropopouse found
 return
-end subroutine twmo
+END SUBROUTINE  twmo