Changeset 5117 for LMDZ6/branches/Amaury_dev/libf/dyn3dmem/iniacademic_loc.F90

Timestamp:

Jul 24, 2024, 4:23:34 PM (4 months ago)

Author:

abarral

Message:

rename modules properly lmdz_*
move some unused files to obsolete/
(lint) uppercase fortran keywords

File:

• LMDZ6/branches/Amaury_dev/libf/dyn3dmem/iniacademic_loc.F90 (modified) (15 diffs)

LMDZ6/branches/Amaury_dev/libf/dyn3dmem/iniacademic_loc.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 parallel_lmdz, ONLY: ijb_u, ije_u, ijb_v, ije_v
   USE IOIPSL, ONLY: getin
-  USE Write_Field
+  USE lmdz_write_field
   USE comconst_mod, ONLY: cpp, kappa, g, daysec, dtvr, pi, im, jm
   USE logic_mod, ONLY: iflag_phys, read_start
@@ -17 +17 @@
   USE temps_mod, ONLY: annee_ref, day_ini, day_ref
   USE ener_mod, ONLY: etot0,ptot0,ztot0,stot0,ang0
-  USE readTracFiles_mod, ONLY: addPhase
-  use netcdf, ONLY: nf90_nowrite,nf90_open,nf90_noerr,nf90_inq_varid,nf90_close, nf90_get_var
+  USE lmdz_readTracFiles, ONLY: addPhase
+  USE netcdf, ONLY: nf90_nowrite,nf90_open,nf90_noerr,nf90_inq_varid,nf90_close, nf90_get_var
+  USE lmdz_ran1, ONLY: ran1
 
   !   Author:    Frederic Hourdin      original: 15/01/93
@@ -60 +61 @@
   REAL phi(ip1jmp1,llm)                  ! geopotentiel
   REAL ddsin,zsig,tetapv,w_pv  ! variables auxiliaires
-  real tetastrat ! potential temperature in the stratosphere, in K
-  real tetajl(jjp1,llm)
+  REAL tetastrat ! potential temperature in the stratosphere, in K
+  REAL tetajl(jjp1,llm)
   INTEGER i,j,l,lsup,ij, iq, iName, iPhase, iqParent
 
   INTEGER :: nid_relief,varid,ierr
-  real, dimension(iip1,jjp1) :: relief
-
+  REAL, DIMENSION(iip1,jjp1) :: relief
 
   REAL teta0,ttp,delt_y,delt_z,eps ! Constantes pour profil de T
@@ -74 +74 @@
   REAL phi_pv,dphi_pv,gam_pv,tetanoise   ! Constantes pour polar vortex 
 
-  real zz,ran1
-  integer idum
+  REAL zz
+  INTEGER idum
 
   REAL zdtvr, tnat, alpha_ideal
@@ -84 +84 @@
 
   ! Sanity check: verify that options selected by user are not incompatible
-  if ((iflag_phys==1).and. .not. read_start) THEN
+  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"
@@ -90 +90 @@
     " (iflag_phys >= 100)"
     CALL abort_gcm(modname,"incompatible iflag_phys==1 and read_start==.FALSE.",1)
-  endif
+  ENDIF
   
   !-----------------------------------------------------------------------
@@ -114 +114 @@
   ang0       = 0.
 
-  if (llm == 1) THEN
+  IF (llm == 1) THEN
      ! specific initializations for the shallow water case
      kappa=1
-  endif
+  ENDIF
 
   CALL iniconst
@@ -148 +148 @@
      relief=0.
      ierr = nf90_open ('relief_in.nc', nf90_nowrite,nid_relief)
-     if (ierr==nf90_noerr) THEN
+     IF (ierr==nf90_noerr) THEN
          ierr=nf90_inq_varid(nid_relief,'RELIEF',varid)
-         if (ierr==nf90_noerr) THEN
+         IF (ierr==nf90_noerr) THEN
               ierr=nf90_get_var(nid_relief,varid,relief(1:iim,1:jjp1))
               relief(iip1,:)=relief(1,:)
@@ -173 +173 @@
 
      CALL pression ( ip1jmp1, ap, bp, ps_glo, p       )
-     if (pressure_exner) THEN
+     IF (pressure_exner) THEN
        CALL exner_hyb( ip1jmp1, ps_glo, p, pks, pk )
      else
@@ -181 +181 @@
   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(ijb_u:ije_u)=0.
         q(ijb_u:ije_u,1:llm,1:nqtot)=0
         CALL sw_case_williamson91_6_loc(vcov,ucov,teta,masse,ps)
      endif
-  endif
+  ENDIF
 
   academic_case: if (iflag_phys >= 2) THEN
@@ -258 +258 @@
            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)   &
@@ -295 +295 @@
 
         ! winds
-        if (ok_geost) THEN
+        IF (ok_geost) THEN
            CALL ugeostr(phi,ucov_glo)
         else
@@ -303 +303 @@
 
         ! 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
@@ -313 +313 @@
               ! distill de Rayleigh très simplifiée
               iName    = tracers(iq)%iso_iName
-              if (niso <= 0 .OR. iName <= 0) CYCLE
+              IF (niso <= 0 .OR. iName <= 0) CYCLE
               iPhase   = tracers(iq)%iso_iPhase
               iqParent = tracers(iq)%iqParent
               IF(tracers(iq)%iso_iZone == 0) THEN
-                 if (tnat1) THEN
+                 IF (tnat1) THEN
                          tnat=1.0
                          alpha_ideal=1.0
@@ -374 +374 @@
         deallocate(phis_glo)
      ENDIF ! of IF (.NOT. read_start)
-  endif academic_case
+  ENDIF academic_case
 
 END SUBROUTINE iniacademic_loc