Changeset 4197 for dynamico_lmdz/simple_physics/phyparam/physics/radiative.F90

Timestamp:

Dec 20, 2019, 10:45:29 AM (5 years ago)

Author:

dubos

Message:

simple_physics : cleanup radiative transfer

File:

• dynamico_lmdz/simple_physics/phyparam/physics/radiative.F90 (moved) (moved from dynamico_lmdz/simple_physics/phyparam/param/mucorr.F) (4 diffs)

dynamico_lmdz/simple_physics/phyparam/physics/radiative.F90

@@ -1 @@
-      SUBROUTINE mucorr(npts,pdeclin, plat, pmu, pfract,phaut,prad)
-      IMPLICIT NONE
+MODULE radiative
+  IMPLICIT NONE
+  SAVE
 
-c=======================================================================
-c
-c   Calcul of equivalent solar angle and and fraction of day whithout 
-c   diurnal cycle.
-c
-c   parmeters :
-c   -----------
-c
-c      Input :
-c      -------
-c         npts             number of points
-c         pdeclin          solar declinaison
-c         plat(npts)        latitude
-c         phaut            hauteur typique de l'atmosphere
-c         prad             rayon planetaire
-c
-c      Output :
-c      --------
-c         pmu(npts)          equivalent cosinus of the solar angle
-c         pfract(npts)       fractionnal day
-c
-c=======================================================================
+CONTAINS
 
-c-----------------------------------------------------------------------
-c
-c    0. Declarations :
-c    -----------------
+  SUBROUTINE mucorr(npts,pdeclin, plat, pmu, pfract,phaut,prad)
 
-c     Arguments :
-c     -----------
+!=======================================================================
+!
+!   Calcul of equivalent solar angle and and fraction of day whithout 
+!   diurnal cycle.
+!
+!   parmeters :
+!   -----------
+!
+!      Input :
+!      -------
+!         npts             number of points
+!         pdeclin          solar declinaison
+!         plat(npts)        latitude
+!         phaut            hauteur typique de l'atmosphere
+!         prad             rayon planetaire
+!
+!      Output :
+!      --------
+!         pmu(npts)          equivalent cosinus of the solar angle
+!         pfract(npts)       fractionnal day
+!
+!=======================================================================
+
+!-----------------------------------------------------------------------
+!
+!    0. Declarations :
+!    -----------------
+
+!     Arguments :
+!     -----------
       INTEGER npts
       REAL plat(npts),pmu(npts), pfract(npts)
       REAL phaut,prad,pdeclin
-c
-c     Local variables :
-c     -----------------
+!
+!     Local variables :
+!     -----------------
       INTEGER j
       REAL pi
@@ -44 +49 @@
       REAL alph
 
-c-----------------------------------------------------------------------
+!-----------------------------------------------------------------------
 
       print*,'npts,pdeclin'
@@ -57 +62 @@
        print*,'cz,sz',cz,sz
 
-      DO 20 j = 1, npts
+      DO j = 1, npts
 
          phi = plat(j)
@@ -91 +96 @@
          IF (pmu(j).eq.0.) pfract(j) = 0.
 
-   20 CONTINUE
+      END DO
 
-c-----------------------------------------------------------------------
-c   correction de rotondite:
-c   ------------------------
+!-----------------------------------------------------------------------
+!   correction de rotondite:
+!   ------------------------
 
       alph=phaut/prad
-      DO 30 j=1,npts
-c !!!!!!
+      DO j=1,npts
+! !!!!!!
          pmu(j)=sqrt(1224.*pmu(j)*pmu(j)+1.)/35.
-c    $          (sqrt(alph*alph*pmu(j)*pmu(j)+2.*alph+1.)-alph*pmu(j))
-30    CONTINUE
+!    $          (sqrt(alph*alph*pmu(j)*pmu(j)+2.*alph+1.)-alph*pmu(j))
+      END DO
 
-      RETURN
-      END
+    END SUBROUTINE mucorr
+
+  END MODULE radiative