Changeset 4233 for dynamico_lmdz/simple_physics/phyparam

Timestamp:

Jan 18, 2020, 5:59:32 PM (5 years ago)

Author:

dubos

Message:

simple_physics : enforce F2003 strictly

Location:

dynamico_lmdz/simple_physics/phyparam

Files:

• Makefile (modified) (1 diff)
• check_F2003.sh (added)
• physics/astronomy.F90 (modified) (8 diffs)
• physics/logging.F90 (modified) (1 diff)
• physics/radiative_lw.F90 (modified) (2 diffs)
• physics/radiative_sw.F90 (modified) (4 diffs)
• physics/read_param_mod.F90 (modified) (1 diff)
• physics/solar.F90 (modified) (5 diffs)
• physics/surface.F90 (modified) (2 diffs)
• physics/turbulence.F90 (modified) (1 diff)

dynamico_lmdz/simple_physics/phyparam/Makefile

@@ -22 +22 @@
 clean :
         rm -f obj/* include/* lib/* xml/*
+
+nice:
         ../bash/beautify.sh physics/*.F90
+
 %.so : $(OBJECTS)
         $(F90) -shared $^ -o $@ 

dynamico_lmdz/simple_physics/phyparam/physics/astronomy.F90

@@ -14 +14 @@
 
   SUBROUTINE solarlong(pday,psollong)
-    REAL, INTENT(IN) :: pday               ! jour de l'annee (le jour 0 correspondant a l'equinoxe)
+    REAL, INTENT(IN) :: pday               ! jour de l annee (le jour 0 correspondant a l equinoxe)
     REAL, INTENT(OUT) :: psollong          ! solar longitude
     LOGICAL, PARAMETER ::  lwrite=.TRUE.
@@ -24 +24 @@
 
     !--------------------------------------------------------
-    ! calcul de l'angle polaire et de la distance au soleil :
+    ! calcul de l angle polaire et de la distance au soleil :
     ! -------------------------------------------------------
 
-    !  calcul de l'zanomalie moyenne
+    !  calcul de l zanomalie moyenne
 
     zz=(pday-peri_day)/year_day
@@ -33 +33 @@
     xref=abs(zanom)
 
-    !  resolution de l'equation horaire  zx0 - e * sin (zx0) = xref
+    !  resolution de l equation horaire  zx0 - e * sin (zx0) = xref
     !  methode de Newton
 
@@ -74 +74 @@
     !   ------
     !    Initialisation du sous programme orbite qui calcule
-    !    a une date donnee de l'annee de duree year_day commencant
-    !    a l'equinoxe de printemps et dont le perihelie se situe
+    !    a une date donnee de l annee de duree year_day commencant
+    !    a l equinoxe de printemps et dont le perihelie se situe
     !    a la date peri_day, la distance au soleil et la declinaison.
     !
@@ -81 +81 @@
     !   ----------
     !   - initialise certaines variables de ce module
-    !   - Doit etre appele avant d'utiliser orbite.
+    !   - Doit etre appele avant d utiliser orbite.
     !
     !   Arguments:
@@ -88 +88 @@
     !   Input:
     !   ------
-    !   aphelie       \   aphelie et perihelie de l'orbite
+    !   aphelie       \   aphelie et perihelie de l orbite
     !   periheli      /   en millions de kilometres.
     !
@@ -113 +113 @@
 
     !-----------------------------------------------------------------------
-    ! calcul de l'angle polaire et de la distance au soleil :
+    ! calcul de l angle polaire et de la distance au soleil :
     ! -------------------------------------------------------
 
-    !  calcul de l'zanomalie moyenne
+    !  calcul de l zanomalie moyenne
 
     zz=(year_day-peri_day)/year_day
@@ -123 +123 @@
     WRITELOG(*,*) 'zanom  ',zanom
 
-    !  resolution de l'equation horaire  zx0 - e * sin (zx0) = zxref
+    !  resolution de l equation horaire  zx0 - e * sin (zx0) = zxref
     !  methode de Newton
 

dynamico_lmdz/simple_physics/phyparam/physics/logging.F90

@@ -98 +98 @@
     INTEGER, INTENT(IN) :: index ! index of gridpoint
     CHARACTER(*), INTENT(OUT) :: line
-    line=''
+    line=' '
   END SUBROUTINE default_log_gridpoint
 

dynamico_lmdz/simple_physics/phyparam/physics/radiative_lw.F90

@@ -22 +22 @@
     !=======================================================================
     !
-    !   calcul de l'evolution de la temperature sous l'effet du rayonnement
+    !   calcul de l evolution de la temperature sous l effet du rayonnement
     !   infra-rouge.
     !   Pour simplifier, les transmissions sont precalculees et ne
-    !   dependent que de l'altitude.
+    !   dependent que de l altitude.
     !
     !   arguments:
@@ -76 +76 @@
 
     !-----------------------------------------------------------------------
-    !   2. calcul des quantites d'absorbants:
+    !   2. calcul des quantites d absorbants:
     !   -------------------------------------
 

dynamico_lmdz/simple_physics/phyparam/physics/radiative_sw.F90

@@ -42 +42 @@
     !
     !   Rayonnement solaire en atmosphere non diffusante avec un
-    !   coefficient d'absoprption gris.
+    !   coefficient d absoprption gris.
     !
     !=======================================================================
@@ -113 +113 @@
     tau0=-.5*log(coefvis)
 
-    ! calcul de la partie homogene de l'opacite
+    ! calcul de la partie homogene de l opacite
     tau0=tau0/ps_rad
     DO l=1,nlayer+1
@@ -122 +122 @@
 
     !-----------------------------------------------------------------------
-    !   2. calcul de la transmission depuis le sommet de l'atmosphere:
+    !   2. calcul de la transmission depuis le sommet de l atmosphere:
     !   -----------------------------------------------------------
 
@@ -202 +202 @@
 
     !-----------------------------------------------------------------------
-    !   6.ajout a l'echauffement de la contribution du ray. sol. reflechit:
+    !   6.ajout a l echauffement de la contribution du ray. sol. reflechit:
     !   -------------------------------------------------------------------
 

dynamico_lmdz/simple_physics/phyparam/physics/read_param_mod.F90

@@ -96 +96 @@
     LOG_INFO('read_param')
   END SUBROUTINE read_paramb
-  
+
 END MODULE read_param_mod

dynamico_lmdz/simple_physics/phyparam/physics/solar.F90

@@ -89 +89 @@
           pfract(jl)=1.
        ELSE
-          pmu0(jl)=0.
+          pmu0(jl)=0.
           pfract(jl)=0.
        ENDIF
@@ -102 +102 @@
     !=============================================================
     ! Auteur : O. Boucher (LMD/CNRS)
-    !          d'apres les routines zenith et angle de Z.X. Li
-    ! Objet  : calculer les valeurs moyennes du cos de l'angle zenithal
-    !          et l'ensoleillement moyen entre gmtime1 et gmtime2
+    !          d apres les routines zenith et angle de Z.X. Li
+    ! Objet  : calculer les valeurs moyennes du cos de l angle zenithal
+    !          et l ensoleillement moyen entre gmtime1 et gmtime2
     !          connaissant la declinaison, la latitude et la longitude.
     ! Rque   : Different de la routine angle en ce sens que zenang
@@ -114 +114 @@
     !===============================================================
     ! longi----INPUT : la longitude vraie de la terre dans son plan
-    !                  solaire a partir de l'equinoxe de printemps (degre)
+    !                  solaire a partir de l equinoxe de printemps (degre)
     ! gmtime---INPUT : temps universel en fraction de jour
     ! pdtrad---INPUT : pas de temps du rayonnement (secondes)
@@ -133 +133 @@
     ! omega1, omega2 : temps 1 et 2 exprime en radian avec 0 a midi.
     ! omega : heure en radian du coucher de soleil
-    ! -omega est donc l'heure en radian de lever du soleil
+    ! -omega est donc l heure en radian de lever du soleil
     real omegadeb, omegafin
     real zfrac1, zfrac2, z1_mu, z2_mu
@@ -275 +275 @@
     !         pdeclin          solar declinaison
     !         plat(npts)        latitude
-    !         phaut            hauteur typique de l'atmosphere
+    !         phaut            hauteur typique de l atmosphere
     !         prad             rayon planetaire
     !

dynamico_lmdz/simple_physics/phyparam/physics/surface.F90

@@ -28 +28 @@
   SUBROUTINE init_soil(ngrid,nsoil)
     INTEGER, INTENT(IN) :: ngrid, nsoil
-    REAL min_period,dalph_soil
-    REAL fz,rk,fz1,rk1,rk2
+    REAL :: min_period,dalph_soil, rk,fz1,rk1,rk2
     INTEGER :: jk
-
-    ! this is a function definition
-    fz(rk)=fz1*(dalph_soil**rk-1.)/(dalph_soil-1.)
 
     !-----------------------------------------------------------------------
@@ -71 +67 @@
     ENDDO
     LOG_INFO('init_soil')
+
+  CONTAINS
+
+    FUNCTION fz(rk) RESULT(val)
+      REAL :: val, rk
+      val = fz1*(dalph_soil**rk-1.)/(dalph_soil-1.)
+    END FUNCTION fz
 
   END SUBROUTINE init_soil

dynamico_lmdz/simple_physics/phyparam/physics/turbulence.F90

@@ -388 +388 @@
 
     !-----------------------------------------------------------------------
-    !   calcul le l'evolution de la temperature du sol:
+    !   calcul le l evolution de la temperature du sol:
     !   -----------------------------------------------