Ignore:
Timestamp:
Dec 20, 2019, 1:34:11 PM (5 years ago)
Author:
dubos
Message:

simple_physics : cleanup solar angles

File:
1 moved

Legend:

Unmodified
Added
Removed

  • dynamico_lmdz/simple_physics/phyparam/physics/solar.F90

    r4196 r4200  
    1       subroutine solang ( kgrid,psilon,pcolon,psilat,pcolat,
    2      &                    ptim1,ptim2,ptim3,pmu0,pfract )
    3       IMPLICIT NONE
     1MODULE solar
     2  IMPLICIT NONE
    43
    5 C
    6 C**** *LW*   - ORGANIZES THE LONGWAVE CALCULATIONS
    7 C
    8 C     PURPOSE.
    9 C     --------
    10 C          CALCULATES THE SOLAR ANGLE FOR ALL THE POINTS OF THE GRID
    11 C
    12 C**   INTERFACE.
    13 C     ----------
    14 C      SUBROUTINE SOLANG ( KGRID )
    15 C
    16 C        EXPLICIT ARGUMENTS :
    17 C        --------------------
    18 C     ==== INPUTS  ===
    19 C
    20 C PSILON(KGRID)   : SINUS OF THE LONGITUDE
    21 C PCOLON(KGRID)   : COSINUS OF THE LONGITUDE
    22 C PSILAT(KGRID)   : SINUS OF THE LATITUDE
    23 C PCOLAT(KGRID)   : COSINUS OF THE LATITUDE
    24 C PTIM1           : SIN(DECLI)
    25 C PTIM2           : COS(DECLI)*COS(TIME)
    26 C PTIM3           : SIN(DECLI)*SIN(TIME)
    27 C
    28 C     ==== OUTPUTS ===
    29 C
    30 C PMU0 (KGRID)    : SOLAR ANGLE
    31 C PFRACT(KGRID)   : DAY FRACTION OF THE TIME INTERVAL
    32 C
    33 C        IMPLICIT ARGUMENTS :   NONE
    34 C        --------------------
    35 C
    36 C     METHOD.
    37 C     -------
    38 C
    39 C     EXTERNALS.
    40 C     ----------
    41 C
    42 C         NONE
    43 C
    44 C     REFERENCE.
    45 C     ----------
    46 C
    47 C         RADIATIVE PROCESSES IN METEOROLOGIE AND CLIMATOLOGIE
    48 C         PALTRIDGE AND PLATT
    49 C
    50 C     AUTHOR.
    51 C     -------
    52 C        FREDERIC HOURDIN
    53 C
    54 C     MODIFICATIONS.
    55 C     --------------
    56 C        ORIGINAL :90-01-14
    57 C                  92-02-14 CALCULATIONS DONE THE ENTIER GRID (J.Polcher)
    58 C-----------------------------------------------------------------------
    59 C
    60 C     ------------------------------------------------------------------
     4CONTAINS
    615
    62 C-----------------------------------------------------------------------
    63 C
    64 C*      0.1   ARGUMENTS
    65 C             ---------
    66 C
     6  subroutine solang ( kgrid,psilon,pcolon,psilat,pcolat, &
     7       ptim1,ptim2,ptim3,pmu0,pfract )
     8
     9!
     10!**** *LW*   - ORGANIZES THE LONGWAVE CALCULATIONS
     11!
     12!     PURPOSE.
     13!     --------
     14!          CALCULATES THE SOLAR ANGLE FOR ALL THE POINTS OF THE GRID
     15!
     16!**   INTERFACE.
     17!     ----------
     18!      SUBROUTINE SOLANG ( KGRID )
     19!
     20!        EXPLICIT ARGUMENTS :
     21!        --------------------
     22!     ==== INPUTS  ===
     23!
     24! PSILON(KGRID)   : SINUS OF THE LONGITUDE
     25! PCOLON(KGRID)   : COSINUS OF THE LONGITUDE
     26! PSILAT(KGRID)   : SINUS OF THE LATITUDE
     27! PCOLAT(KGRID)   : COSINUS OF THE LATITUDE
     28! PTIM1           : SIN(DECLI)
     29! PTIM2           : COS(DECLI)*COS(TIME)
     30! PTIM3           : SIN(DECLI)*SIN(TIME)
     31!
     32!     ==== OUTPUTS ===
     33!
     34! PMU0 (KGRID)    : SOLAR ANGLE
     35! PFRACT(KGRID)   : DAY FRACTION OF THE TIME INTERVAL
     36!
     37!        IMPLICIT ARGUMENTS :   NONE
     38!        --------------------
     39!
     40!     METHOD.
     41!     -------
     42!
     43!     EXTERNALS.
     44!     ----------
     45!
     46!         NONE
     47!
     48!     REFERENCE.
     49!     ----------
     50!
     51!         RADIATIVE PROCESSES IN METEOROLOGIE AND CLIMATOLOGIE
     52!         PALTRIDGE AND PLATT
     53!
     54!     AUTHOR.
     55!     -------
     56!        FREDERIC HOURDIN
     57!
     58!     MODIFICATIONS.
     59!     --------------
     60!        ORIGINAL :90-01-14
     61!                  92-02-14 CALCULATIONS DONE THE ENTIER GRID (J.Polcher)
     62!-----------------------------------------------------------------------
     63!
     64!     ------------------------------------------------------------------
     65!-----------------------------------------------------------------------
     66!
     67!*      0.1   ARGUMENTS
     68!             ---------
     69!
    6770      INTEGER kgrid
    6871      REAL ptim1,ptim2,ptim3
    6972      REAL psilon(kgrid),pcolon(kgrid),pmu0(kgrid),pfract(kgrid)
    7073      REAL psilat(kgrid), pcolat(kgrid)
    71 C
     74
    7275      INTEGER jl
    7376      REAL ztim1,ztim2,ztim3
    74 C------------------------------------------------------------------------
    75 C------------------------------------------------------------------------
    76 C------------------------------------------------------------------------
    77 C
    78 C------------------------------------------------------------------------
    79 C
    80 C*     1.     INITIALISATION
    81 C             --------------
    82 C
     77!------------------------------------------------------------------------
     78!------------------------------------------------------------------------
     79!------------------------------------------------------------------------
     80!
     81!------------------------------------------------------------------------
     82!
     83!*     1.     INITIALISATION
     84!             --------------
     85!
    8386 100  CONTINUE
    84 C
     87!
    8588      DO jl=1,kgrid
    8689        pmu0(jl)=0.
    8790        pfract(jl)=0.
    8891      ENDDO
    89 C
    90 C*     1.1     COMPUTATION OF THE SOLAR ANGLE
    91 C              ------------------------------
    92 C
     92!
     93!*     1.1     COMPUTATION OF THE SOLAR ANGLE
     94!              ------------------------------
     95!
    9396      DO jl=1,kgrid
    9497        ztim1=psilat(jl)*ptim1
    … …  
    97100        pmu0(jl)=ztim1+ztim2*pcolon(jl)+ztim3*psilon(jl)
    98101      ENDDO
    99 C
    100 C*     1.2      DISTINCTION BETWEEN DAY AND NIGHT
    101 C               ---------------------------------
    102 C
     102!
     103!*     1.2      DISTINCTION BETWEEN DAY AND NIGHT
     104!               ---------------------------------
     105!
    103106      DO jl=1,kgrid
    104107        IF (pmu0(jl).gt.0.) THEN
    … …  
    109112        ENDIF
    110113      ENDDO
    111 C
    112       RETURN
    113       END
     114!
     115
     116    END subroutine solang
     117  END MODULE solar
Note: See TracChangeset for help on using the changeset viewer.