Changeset 5116 for LMDZ6/branches/Amaury_dev/libf/dyn3d/iniacademic.F90

Timestamp:

Jul 24, 2024, 2:54:37 PM (2 months ago)

Author:

abarral

Message:

rename modules properly lmdz_*
move ismin, ismax, minmax into new lmdz_libmath.f90
(lint) uppercase fortran keywords

File:

• LMDZ6/branches/Amaury_dev/libf/dyn3d/iniacademic.F90 (modified) (11 diffs)

LMDZ6/branches/Amaury_dev/libf/dyn3d/iniacademic.F90

@@ -7 +7 @@
   USE infotrac,    ONLY: nqtot, niso, iqIsoPha, tracers, getKey, isoName
   USE control_mod, ONLY: day_step,planet_type
-  use exner_hyb_m, only: exner_hyb
-  use exner_milieu_m, only: exner_milieu
+  use exner_hyb_m, ONLY: exner_hyb
+  use exner_milieu_m, ONLY: exner_milieu
   USE IOIPSL, ONLY: getin
   USE Write_Field
@@ -60 +60 @@
   INTEGER i,j,l,lsup,ij, iq, iName, iPhase, iqParent
 
-  integer :: nid_relief,varid,ierr
+  INTEGER :: nid_relief,varid,ierr
   real, dimension(iip1,jjp1) :: relief
 
@@ -75 +75 @@
   LOGICAL,PARAMETER :: tnat1=.TRUE.
   
-  character(len=*),parameter :: modname="iniacademic"
-  character(len=80) :: abort_message
+  CHARACTER(LEN=*),parameter :: modname="iniacademic"
+  CHARACTER(LEN=80) :: abort_message
 
   ! Sanity check: verify that options selected by user are not incompatible
-  if ((iflag_phys==1).and. .not. read_start) then
-    write(lunout,*) trim(modname)," error: if read_start is set to ", &
+  if ((iflag_phys==1).and. .not. read_start) THEN
+    WRITE(lunout,*) trim(modname)," error: if read_start is set to ", &
     " false then iflag_phys should not be 1"
-    write(lunout,*) "You most likely want an aquaplanet initialisation", &
+    WRITE(lunout,*) "You most likely want an aquaplanet initialisation", &
     " (iflag_phys >= 100)"
     CALL abort_gcm(modname,"incompatible iflag_phys==1 and read_start==.FALSE.",1)
@@ -109 +109 @@
   ang0       = 0.
 
-  if (llm == 1) then
+  if (llm == 1) THEN
      ! specific initializations for the shallow water case
      kappa=1
@@ -164 +164 @@
      CALL pression ( ip1jmp1, ap, bp, ps, p       )
 
-     if (pressure_exner) then
+     if (pressure_exner) THEN
        CALL exner_hyb( ip1jmp1, ps, p, pks, pk)
      else
@@ -172 +172 @@
   ENDIF
 
-  if (llm == 1) then
+  if (llm == 1) THEN
      ! initialize fields for the shallow water case, if required
-     if (.not.read_start) then
+     if (.not.read_start) THEN
         phis(:)=0.
         q(:,:,:)=0
@@ -181 +181 @@
   endif
 
-  academic_case: if (iflag_phys >= 2) then
+  academic_case: if (iflag_phys >= 2) THEN
      ! initializations
 
@@ -249 +249 @@
            tetajl(j,l)=teta0-delt_y*ddsin*ddsin+eps*ddsin &
                 -delt_z*(1.-ddsin*ddsin)*log(zsig)
-           if (planet_type=="giant") then
+           if (planet_type=="giant") THEN
              tetajl(j,l)=teta0+(delt_y*                   &
                 ((sin(rlatu(j)*3.14159*eps+0.0001))**2)   &
@@ -293 +293 @@
 
         ! winds
-        if (ok_geost) then
+        if (ok_geost) THEN
            CALL ugeostr(phi,ucov)
         else
@@ -301 +301 @@
 
         ! bulk initialization of tracers
-        if (planet_type=="earth") then
+        if (planet_type=="earth") THEN
            ! Earth: first two tracers will be water
            do iq=1,nqtot
@@ -315 +315 @@
               iqParent = tracers(iq)%iqParent
               IF(tracers(iq)%iso_iZone == 0) THEN
-                 if (tnat1) then
+                 if (tnat1) THEN
                          tnat=1.0
                          alpha_ideal=1.0
-                         write(*,*) 'Attention dans iniacademic: alpha_ideal=1'
+                         WRITE(*,*) 'Attention dans iniacademic: alpha_ideal=1'
                  else
                     IF(getKey('tnat', tnat, isoName(iName)) .OR. getKey('alpha', alpha_ideal, isoName(iName))) &