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

Timestamp:

Jul 23, 2024, 3:29:36 PM (3 months 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/calltherm.F90 (modified) (14 diffs)

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

@@ -2 +2 @@
 ! $Id$
 
-      subroutine calltherm(dtime  &
+      SUBROUTINE calltherm(dtime  &
         ,pplay,paprs,pphi,weak_inversion  &
         ,u_seri_,v_seri_,t_seri_,q_seri_,t_env,q_env,zqsat,debut  &
@@ -148 +148 @@
 
       integer i,k,isplit
-      logical, save :: first=.true.
+      logical, save :: first=.TRUE.
       logical :: new_thermcell
 
@@ -168 +168 @@
       if (first) then
         itap=0
-        first=.false.
+        first=.FALSE.
       endif
 
@@ -181 +181 @@
 !  Modele du thermique
 !  ===================
-!         print*,'thermiques: WARNING on passe t au lieu de t_seri'
+!         PRINT*,'thermiques: WARNING on passe t au lieu de t_seri'
 
 
@@ -199 +199 @@
          alp_bl(:)=0.
          if (prt_level>=10) then
-          print*,'thermV4 nsplit: ',nsplit_thermals,' weak_inversion'
+          PRINT*,'thermV4 nsplit: ',nsplit_thermals,' weak_inversion'
          endif
 
@@ -224 +224 @@
                endif
 !               if (logexpr0) &
-!    &             print*,'WARN eau<0 avant therm i=',i,'  k=',k  &
+!    &             PRINT*,'WARN eau<0 avant therm i=',i,'  k=',k  &
 !    &         ,' dq,q',d_q_the(i,k),q_seri(i,k)
             enddo
          enddo
-         if(nbptspb>0) print*,'Number of points with q_seri(i,k)<=0 ',nbptspb
+         if(nbptspb>0) PRINT*,'Number of points with q_seri(i,k)<=0 ',nbptspb
 
 
@@ -330 +330 @@
 
 ! Attention : les noms sont contre intuitif.
-! flag_bidouille_stratocu est .true. si on ne fait pas de bidouille.
+! flag_bidouille_stratocu est .TRUE. si on ne fait pas de bidouille.
 ! Il aurait mieux valu avoir un nobidouille_stratocu
 ! Et pour simplifier :
@@ -422 +422 @@
             alp_bl(i)=alp_bl(i)+alp(i)/REAL(nsplit_thermals)
 !            write(23,*)'ALP CALLTHERM',alp_bl(i),alp(i)
-        if(prt_level>=10) print*,'calltherm i alp_bl alp ale_bl ale',i,alp_bl(i),alp(i),ale_bl(i),ale(i)
+        if(prt_level>=10) PRINT*,'calltherm i alp_bl alp ale_bl ale',i,alp_bl(i),alp(i),ale_bl(i),ale(i)
        ENDDO
 
@@ -441 +441 @@
 #endif
 !                if (prt_level.ge.10) then
-!                  print*,'WARN eau<0 apres therm i=',i,'  k=',k  &
+!                  PRINT*,'WARN eau<0 apres therm i=',i,'  k=',k  &
 !    &         ,' dq,q',d_q_the(i,k),q_seri(i,k),  &
 !    &         'fm=',zfm_therm(i,k),'entr=',entr_therm(i,k)
@@ -457 +457 @@
 #endif
 
-        IF(nbptspb>0) print*,'Number of points with q_seri(i,k)<=0 ',nbptspb
+        IF(nbptspb>0) PRINT*,'Number of points with q_seri(i,k)<=0 ',nbptspb
 ! tests sur les valeurs de la temperature
         nbptspb=0
@@ -466 +466 @@
 !              if ((t_seri(i,k).lt.50.) .or.  &
 !    &              (t_seri(i,k).gt.370.)) then
-!                 print*,'WARN temp apres therm i=',i,'  k=',k  &
+!                 PRINT*,'WARN temp apres therm i=',i,'  k=',k  &
 !    &         ,' t_seri',t_seri(i,k)
 !              CALL abort
@@ -472 +472 @@
             ENDDO
             ENDDO
-        IF(nbptspb>0) print*,'Number of points with q_seri(i,k)<=0 ',nbptspb
+        IF(nbptspb>0) PRINT*,'Number of points with q_seri(i,k)<=0 ',nbptspb
          enddo ! isplit
 
 !***************************************************************
 !     calcul du flux ascencant conservatif
-!            print*,'<<<<calcul flux ascendant conservatif'
+!            PRINT*,'<<<<calcul flux ascendant conservatif'
 
       fmc_therm=0.
@@ -495 +495 @@
 !****************************************************************
 !     calcul de l'humidite dans l'ascendance
-!      print*,'<<<<calcul de lhumidite dans thermique'
+!      PRINT*,'<<<<calcul de lhumidite dans thermique'
 !CR:on ne le calcule que pour le cas sec
       if (iflag_thermals<=11) then
@@ -514 +514 @@
 
 !     calcul de l'eau condensee dans l'ascendance
-!             print*,'<<<<calcul de leau condensee dans thermique'
+!             PRINT*,'<<<<calcul de leau condensee dans thermique'
              do i=1,klon
                 do k=1,klev