Changeset 4200

Timestamp:

Dec 20, 2019, 1:34:11 PM (5 years ago)

Author:

dubos

Message:

simple_physics : cleanup solar angles

Location:

dynamico_lmdz/simple_physics/phyparam

Files:

• param/mufract.F (deleted)
• param/phyparam.F (modified) (1 diff)
• physics/solar.F90 (moved) (moved from dynamico_lmdz/simple_physics/phyparam/param/solang.F) (3 diffs)

dynamico_lmdz/simple_physics/phyparam/param/phyparam.F

@@ -14 +14 @@
       USE astronomy
       USE vdif_mod, ONLY : vdif
+      USE solar, ONLY : solang
       USE radiative, ONLY : mucorr, sw
       USE radiative_lw, ONLY : lw

dynamico_lmdz/simple_physics/phyparam/physics/solar.F90

@@ -1 @@
-      subroutine solang ( kgrid,psilon,pcolon,psilat,pcolat,
-     &                    ptim1,ptim2,ptim3,pmu0,pfract )
-      IMPLICIT NONE
+MODULE solar
+  IMPLICIT NONE
 
-C
-C**** *LW*   - ORGANIZES THE LONGWAVE CALCULATIONS
-C
-C     PURPOSE.
-C     --------
-C          CALCULATES THE SOLAR ANGLE FOR ALL THE POINTS OF THE GRID
-C
-C**   INTERFACE.
-C     ----------
-C      SUBROUTINE SOLANG ( KGRID )
-C
-C        EXPLICIT ARGUMENTS :
-C        --------------------
-C     ==== INPUTS  ===
-C
-C PSILON(KGRID)   : SINUS OF THE LONGITUDE
-C PCOLON(KGRID)   : COSINUS OF THE LONGITUDE
-C PSILAT(KGRID)   : SINUS OF THE LATITUDE
-C PCOLAT(KGRID)   : COSINUS OF THE LATITUDE
-C PTIM1           : SIN(DECLI)
-C PTIM2           : COS(DECLI)*COS(TIME)
-C PTIM3           : SIN(DECLI)*SIN(TIME)
-C
-C     ==== OUTPUTS ===
-C
-C PMU0 (KGRID)    : SOLAR ANGLE
-C PFRACT(KGRID)   : DAY FRACTION OF THE TIME INTERVAL
-C
-C        IMPLICIT ARGUMENTS :   NONE
-C        --------------------
-C
-C     METHOD.
-C     -------
-C
-C     EXTERNALS.
-C     ----------
-C
-C         NONE
-C
-C     REFERENCE.
-C     ----------
-C
-C         RADIATIVE PROCESSES IN METEOROLOGIE AND CLIMATOLOGIE
-C         PALTRIDGE AND PLATT
-C
-C     AUTHOR.
-C     -------
-C        FREDERIC HOURDIN
-C
-C     MODIFICATIONS.
-C     --------------
-C        ORIGINAL :90-01-14
-C                  92-02-14 CALCULATIONS DONE THE ENTIER GRID (J.Polcher)
-C-----------------------------------------------------------------------
-C
-C     ------------------------------------------------------------------
+CONTAINS
 
-C-----------------------------------------------------------------------
-C
-C*      0.1   ARGUMENTS
-C             ---------
-C
+  subroutine solang ( kgrid,psilon,pcolon,psilat,pcolat, &
+       ptim1,ptim2,ptim3,pmu0,pfract )
+
+!
+!**** *LW*   - ORGANIZES THE LONGWAVE CALCULATIONS
+!
+!     PURPOSE.
+!     --------
+!          CALCULATES THE SOLAR ANGLE FOR ALL THE POINTS OF THE GRID
+!
+!**   INTERFACE.
+!     ----------
+!      SUBROUTINE SOLANG ( KGRID )
+!
+!        EXPLICIT ARGUMENTS :
+!        --------------------
+!     ==== INPUTS  ===
+!
+! PSILON(KGRID)   : SINUS OF THE LONGITUDE
+! PCOLON(KGRID)   : COSINUS OF THE LONGITUDE
+! PSILAT(KGRID)   : SINUS OF THE LATITUDE
+! PCOLAT(KGRID)   : COSINUS OF THE LATITUDE
+! PTIM1           : SIN(DECLI)
+! PTIM2           : COS(DECLI)*COS(TIME)
+! PTIM3           : SIN(DECLI)*SIN(TIME)
+!
+!     ==== OUTPUTS ===
+!
+! PMU0 (KGRID)    : SOLAR ANGLE
+! PFRACT(KGRID)   : DAY FRACTION OF THE TIME INTERVAL
+!
+!        IMPLICIT ARGUMENTS :   NONE
+!        --------------------
+!
+!     METHOD.
+!     -------
+!
+!     EXTERNALS.
+!     ----------
+!
+!         NONE
+!
+!     REFERENCE.
+!     ----------
+!
+!         RADIATIVE PROCESSES IN METEOROLOGIE AND CLIMATOLOGIE
+!         PALTRIDGE AND PLATT
+!
+!     AUTHOR.
+!     -------
+!        FREDERIC HOURDIN
+!
+!     MODIFICATIONS.
+!     --------------
+!        ORIGINAL :90-01-14
+!                  92-02-14 CALCULATIONS DONE THE ENTIER GRID (J.Polcher)
+!-----------------------------------------------------------------------
+!
+!     ------------------------------------------------------------------
+!-----------------------------------------------------------------------
+!
+!*      0.1   ARGUMENTS
+!             ---------
+!
       INTEGER kgrid
       REAL ptim1,ptim2,ptim3
       REAL psilon(kgrid),pcolon(kgrid),pmu0(kgrid),pfract(kgrid)
       REAL psilat(kgrid), pcolat(kgrid)
-C
+
       INTEGER jl
       REAL ztim1,ztim2,ztim3
-C------------------------------------------------------------------------
-C------------------------------------------------------------------------
-C------------------------------------------------------------------------
-C
-C------------------------------------------------------------------------
-C
-C*     1.     INITIALISATION
-C             --------------
-C
+!------------------------------------------------------------------------
+!------------------------------------------------------------------------
+!------------------------------------------------------------------------
+!
+!------------------------------------------------------------------------
+!
+!*     1.     INITIALISATION
+!             --------------
+!
  100  CONTINUE
-C
+!
       DO jl=1,kgrid
         pmu0(jl)=0.
         pfract(jl)=0.
       ENDDO
-C
-C*     1.1     COMPUTATION OF THE SOLAR ANGLE
-C              ------------------------------
-C
+!
+!*     1.1     COMPUTATION OF THE SOLAR ANGLE
+!              ------------------------------
+!
       DO jl=1,kgrid
         ztim1=psilat(jl)*ptim1
@@ -97 +100 @@
         pmu0(jl)=ztim1+ztim2*pcolon(jl)+ztim3*psilon(jl)
       ENDDO
-C
-C*     1.2      DISTINCTION BETWEEN DAY AND NIGHT
-C               ---------------------------------
-C
+!
+!*     1.2      DISTINCTION BETWEEN DAY AND NIGHT
+!               ---------------------------------
+!
       DO jl=1,kgrid
         IF (pmu0(jl).gt.0.) THEN
@@ -109 +112 @@
         ENDIF
       ENDDO
-C
-      RETURN
-      END
+!
+
+    END subroutine solang
+  END MODULE solar